timeseriesmanagermodule Module Reference

Data Types
type	timeseriesmanagertype

Functions/Subroutines
subroutine, public	tsmanager_cr (this, iout, removeTsLinksOnCompletion, extendTsToEndOfSimulation)
	Create the tsmanager. More...
subroutine	tsmanager_df (this)
	Define time series manager object. More...
subroutine	add_tsfile (this, fname, inunit)
	Add a time series file to this manager. More...
subroutine	tsmgr_ad (this)
	Time step (or subtime step) advance. Call this each time step or subtime step. More...
subroutine	tsmgr_da (this)
	Deallocate memory. More...
subroutine	reset (this, pkgName)
	Call this when a new BEGIN PERIOD block is read for a new stress period. More...
subroutine	make_link (this, timeSeries, pkgName, auxOrBnd, bndElem, irow, jcol, iprpak, tsLink, text, bndName)
	Make link. More...
type(timeserieslinktype) function, pointer	getlink (this, auxOrBnd, indx)
	Get link. More...
integer(i4b) function	countlinks (this, auxOrBnd)
	Count links. More...
type(timeseriestype) function, pointer	get_time_series (this, name)
	Get time series. More...
subroutine	hashbndtimeseries (this)
	Store all boundary (stress) time series links in TsContainers and construct hash table BndTsHashTable. More...
subroutine, public	read_value_or_time_series (textInput, ii, jj, bndElem, pkgName, auxOrBnd, tsManager, iprpak, tsLink)
	Call this subroutine if the time-series link is available or needed. More...
subroutine, public	read_value_or_time_series_adv (textInput, ii, jj, bndElem, pkgName, auxOrBnd, tsManager, iprpak, varName)
	Call this subroutine from advanced packages to define timeseries link for a variable (varName). More...
logical function	remove_existing_link (tsManager, ii, jj, pkgName, auxOrBnd, varName)
	Remove an existing timeseries link if it is defined. More...
logical function, public	var_timeseries (tsManager, pkgName, varName, auxOrBnd)
	Determine if a timeseries link with varName is defined. More...

Function/Subroutine Documentation

add_tsfile()

subroutine timeseriesmanagermodule::add_tsfile ( class(timeseriesmanagertype)  this, character(len=*), intent(in)  fname, integer(i4b), intent(in)  inunit  )

private

Definition at line 102 of file TimeSeriesManager.f90.

    ! -- modules
    use simmodule, only: store_error, store_error_unit
    use arrayhandlersmodule, only: expandarray
    ! -- dummy
    class(TimeSeriesManagerType) :: this
    character(len=*), intent(in) :: fname
    integer(I4B), intent(in) :: inunit
    ! -- local
    integer(I4B) :: isize
    integer(I4B) :: i
    class(TimeSeriesFileType), pointer :: tsfile => null()
    !
    ! -- Check for fname duplicates
    if (this%numtsfiles > 0) then
      do i = 1, this%numtsfiles
        if (this%tsfiles(i) == fname) then
          call store_error('Found duplicate time-series file name: '//trim(fname))
          call store_error_unit(inunit)
        end if
      end do
    end if
    !
    ! -- Save fname
    this%numtsfiles = this%numtsfiles + 1
    isize = size(this%tsfiles)
    if (this%numtsfiles > isize) then
      call expandarray(this%tsfiles, 1000)
    end if
    this%tsfiles(this%numtsfiles) = fname
    !
    ! --
    call this%tsfileList%Add(fname, this%iout, tsfile)
    !
    ! -- Return
    return

Here is the call graph for this function:

countlinks()

integer(i4b) function timeseriesmanagermodule::countlinks ( class(timeseriesmanagertype)  this, character(len=3), intent(in)  auxOrBnd  )

private

Definition at line 461 of file TimeSeriesManager.f90.

    ! -- return
    integer(I4B) :: CountLinks
    ! -- dummy
    class(TimeSeriesManagerType) :: this
    character(len=3), intent(in) :: auxOrBnd
    !
    countlinks = 0
    if (auxorbnd == 'BND') then
      countlinks = this%boundTsLinks%Count()
    elseif (auxorbnd == 'AUX') then
      countlinks = this%auxvarTsLinks%count()
    end if
    !
    ! -- Return
    return

get_time_series()

type(timeseriestype) function, pointer timeseriesmanagermodule::get_time_series ( class(timeseriesmanagertype)  this, character(len=*), intent(in)  name  )

private

Definition at line 481 of file TimeSeriesManager.f90.

    ! -- dummy
    class(TimeSeriesManagerType) :: this
    character(len=*), intent(in) :: name
    ! -- result
    type(TimeSeriesType), pointer :: res
    ! -- local
    integer(I4B) :: indx
    !
    ! Get index from hash table, get time series from TsContainers,
    !     and assign result to time series contained in link.
    res => null()
    indx = this%BndTsHashTable%get(name)
    if (indx > 0) then
      res => this%TsContainers(indx)%timeSeries
    end if
    !
    ! -- Return
    return

getlink()

type(timeserieslinktype) function, pointer timeseriesmanagermodule::getlink ( class(timeseriesmanagertype)  this, character(len=3), intent(in)  auxOrBnd, integer(i4b), intent(in)  indx  )

private

Definition at line 432 of file TimeSeriesManager.f90.

    ! -- dummy
    class(TimeSeriesManagerType) :: this
    character(len=3), intent(in) :: auxOrBnd
    integer(I4B), intent(in) :: indx
    type(TimeSeriesLinkType), pointer :: tsLink
    ! -- local
    type(ListType), pointer :: list
    !
    list => null()
    tslink => null()
    !
    select case (auxorbnd)
    case ('AUX')
      list => this%auxvarTsLinks
    case ('BND')
      list => this%boundTsLinks
    end select
    !
    if (associated(list)) then
      tslink => gettimeserieslinkfromlist(list, indx)
    end if
    !
    ! -- Return
    return

Here is the call graph for this function:

hashbndtimeseries()

subroutine timeseriesmanagermodule::hashbndtimeseries ( class(timeseriesmanagertype), intent(inout)  this )

private

Definition at line 505 of file TimeSeriesManager.f90.

    ! -- dummy
    class(TimeSeriesManagerType), intent(inout) :: this
    ! -- local
    integer(I4B) :: i, j, k, numtsfiles, numts
    character(len=LENTIMESERIESNAME) :: name
    type(TimeSeriesFileType), pointer :: tsfile => null()
    !
    ! Initialize the hash table
    call hash_table_cr(this%BndTsHashTable)
    !
    ! Allocate the TsContainers array to accommodate all time-series links.
    numts = this%tsfileList%CountTimeSeries()
    allocate (this%TsContainers(numts))
    !
    ! Store a pointer to each time series in the TsContainers array
    ! and put its key (time-series name) and index in the hash table.
    numtsfiles = this%tsfileList%Counttsfiles()
    k = 0
    do i = 1, numtsfiles
      tsfile => this%tsfileList%Gettsfile(i)
      numts = tsfile%Count()
      do j = 1, numts
        k = k + 1
        this%TsContainers(k)%timeSeries => tsfile%GetTimeSeries(j)
        if (associated(this%TsContainers(k)%timeSeries)) then
          name = this%TsContainers(k)%timeSeries%Name
          call this%BndTsHashTable%add(name, k)
        end if
      end do
    end do
    !
    ! -- Return
    return

Here is the call graph for this function:

make_link()

subroutine timeseriesmanagermodule::make_link ( class(timeseriesmanagertype), intent(inout)  this, type(timeseriestype), intent(inout), pointer  timeSeries, character(len=*), intent(in)  pkgName, character(len=3), intent(in)  auxOrBnd, real(dp), intent(inout), pointer  bndElem, integer(i4b), intent(in)  irow, integer(i4b), intent(in)  jcol, integer(i4b), intent(in)  iprpak, type(timeserieslinktype), intent(inout), pointer  tsLink, character(len=*), intent(in)  text, character(len=*), intent(in)  bndName  )

private

Definition at line 397 of file TimeSeriesManager.f90.

    ! -- dummy
    class(TimeSeriesManagerType), intent(inout) :: this
    type(TimeSeriesType), pointer, intent(inout) :: timeSeries
    character(len=*), intent(in) :: pkgName
    character(len=3), intent(in) :: auxOrBnd
    real(DP), pointer, intent(inout) :: bndElem
    integer(I4B), intent(in) :: irow, jcol
    integer(I4B), intent(in) :: iprpak
    type(TimeSeriesLinkType), pointer, intent(inout) :: tsLink
    character(len=*), intent(in) :: text
    character(len=*), intent(in) :: bndName
    !
    tslink => null()
    call constructtimeserieslink(tslink, timeseries, pkgname, &
                                 auxorbnd, bndelem, irow, jcol, iprpak)
    if (associated(tslink)) then
      if (auxorbnd == 'BND') then
        call addtimeserieslinktolist(this%boundTsLinks, tslink)
      elseif (auxorbnd == 'AUX') then
        call addtimeserieslinktolist(this%auxvarTsLinks, tslink)
      else
        call store_error('programmer error in make_link', terminate=.true.)
      end if
      tslink%Text = text
      tslink%BndName = bndname
    end if
    !
    ! -- Return
    return

Here is the call graph for this function:

read_value_or_time_series()

subroutine, public timeseriesmanagermodule::read_value_or_time_series	(	character(len=*), intent(in)	textInput,
		integer(i4b), intent(in)	ii,
		integer(i4b), intent(in)	jj,
		real(dp), intent(inout), pointer	bndElem,
		character(len=*), intent(in)	pkgName,
		character(len=3), intent(in)	auxOrBnd,
		type(timeseriesmanagertype), intent(inout)	tsManager,
		integer(i4b), intent(in)	iprpak,
		type(timeserieslinktype), intent(inout), pointer	tsLink
	)

This routine assumes that there is not an existing link for the specified package and array row and column. For the standard boundary package, all links are removed by calling the tsmanagerreset() method. For advanced packages, there is a separate routine called read_value_or_time_series_adv, which should be used instead of this one.

Definition at line 551 of file TimeSeriesManager.f90.

    ! dummy
    character(len=*), intent(in) :: textInput
    integer(I4B), intent(in) :: ii
    integer(I4B), intent(in) :: jj
    real(DP), pointer, intent(inout) :: bndElem
    character(len=*), intent(in) :: pkgName
    character(len=3), intent(in) :: auxOrBnd
    type(TimeSeriesManagerType), intent(inout) :: tsManager
    integer(I4B), intent(in) :: iprpak
    type(TimeSeriesLinkType), pointer, intent(inout) :: tsLink
    ! local
    type(TimeSeriesType), pointer :: timeseries => null()
    integer(I4B) :: istat
    real(DP) :: r
    character(len=LINELENGTH) :: errmsg
    character(len=LENTIMESERIESNAME) :: tsNameTemp
 
    read (textinput, *, iostat=istat) r
    if (istat == 0) then
      bndelem = r
    else
 
      ! Check to see if this is a time series name
      tsnametemp = textinput
      call upcase(tsnametemp)
      timeseries => tsmanager%get_time_series(tsnametemp)
 
      ! If this is a time series, then create a link between
      ! the time series and the row and column position of
      ! bndElem
      if (associated(timeseries)) then
        ! -- Assign value from time series to current
        !    array element
        r = timeseries%GetValue(totimsav, totim, &
                                tsmanager%extendTsToEndOfSimulation)
        bndelem = r
        ! Make a new link between the time series and the row and
        ! column position in the array
        call tsmanager%make_link(timeseries, pkgname, auxorbnd, bndelem, &
                                 ii, jj, iprpak, tslink, '', '')
      else
        errmsg = 'Error in list input. Expected numeric value or '// &
                 "time-series name, but found '"//trim(textinput)//"'."
        call store_error(errmsg)
      end if
    end if
 

Here is the call graph for this function:

Here is the caller graph for this function:

read_value_or_time_series_adv()

subroutine, public timeseriesmanagermodule::read_value_or_time_series_adv	(	character(len=*), intent(in)	textInput,
		integer(i4b), intent(in)	ii,
		integer(i4b), intent(in)	jj,
		real(dp), intent(inout), pointer	bndElem,
		character(len=*), intent(in)	pkgName,
		character(len=3), intent(in)	auxOrBnd,
		type(timeseriesmanagertype), intent(inout)	tsManager,
		integer(i4b), intent(in)	iprpak,
		character(len=*), intent(in)	varName
	)

Arguments are as follows: textInput : string that is either a float or a string name ii : column number jj : row number bndElem : pointer to a position in an array in package pkgName pkgName : package name auxOrBnd : 'AUX' or 'BND' keyword tsManager : timeseries manager object for package iprpak : integer flag indicating if interpolated timeseries values should be printed to package iout during TsManagerad() varName : variable name

Definition at line 617 of file TimeSeriesManager.f90.

    ! -- dummy
    character(len=*), intent(in) :: textInput
    integer(I4B), intent(in) :: ii
    integer(I4B), intent(in) :: jj
    real(DP), pointer, intent(inout) :: bndElem
    character(len=*), intent(in) :: pkgName
    character(len=3), intent(in) :: auxOrBnd
    type(TimeSeriesManagerType), intent(inout) :: tsManager
    integer(I4B), intent(in) :: iprpak
    character(len=*), intent(in) :: varName
    ! -- local
    integer(I4B) :: istat
    real(DP) :: v
    character(len=LINELENGTH) :: errmsg
    character(len=LENTIMESERIESNAME) :: tsNameTemp
    logical :: found
    type(TimeSeriesType), pointer :: timeseries => null()
    type(TimeSeriesLinkType), pointer :: tsLink => null()
    !
    ! -- attempt to read textInput as a real value
    read (textinput, *, iostat=istat) v
    !
    ! -- numeric value
    if (istat == 0) then
      !
      ! -- Numeric value was successfully read.
      bndelem = v
      !
      ! -- remove existing link if it exists for this boundary element
      found = remove_existing_link(tsmanager, ii, jj, pkgname, &
                                   auxorbnd, varname)
      !
      ! -- timeseries
    else
      !
      ! -- attempt to read numeric value from textInput failed.
      !    Text should be a time-series name.
      tsnametemp = textinput
      call upcase(tsnametemp)
      !
      ! -- if textInput is a time-series name, get average value
      !    from time series.
      timeseries => tsmanager%get_time_series(tsnametemp)
      !
      ! -- create a time series link and add it to the package
      !    list of time series links used by the array.
      if (associated(timeseries)) then
        !
        ! -- Assign average value from time series to current array element
        v = timeseries%GetValue(totimsav, totim, &
                                tsmanager%extendTsToEndOfSimulation)
        bndelem = v
        !
        ! -- remove existing link if it exists for this boundary element
        found = remove_existing_link(tsmanager, ii, jj, &
                                     pkgname, auxorbnd, varname)
        !
        ! -- Add link to the list.
        call tsmanager%make_link(timeseries, pkgname, auxorbnd, bndelem, &
                                 ii, jj, iprpak, tslink, varname, '')
        !
        ! -- not a valid timeseries name
      else
        errmsg = 'Error in list input. Expected numeric value or '// &
                 "time-series name, but found '"//trim(textinput)//"'."
        call store_error(errmsg)
      end if
    end if
    !
    ! -- Return
    return

Here is the call graph for this function:

Here is the caller graph for this function:

remove_existing_link()

logical function timeseriesmanagermodule::remove_existing_link ( type(timeseriesmanagertype), intent(inout)  tsManager, integer(i4b), intent(in)  ii, integer(i4b), intent(in)  jj, character(len=*), intent(in)  pkgName, character(len=3), intent(in)  auxOrBnd, character(len=*), intent(in)  varName  )

private

Arguments are as follows: tsManager : timeseries manager object for package ii : column number jj : row number pkgName : package name auxOrBnd : 'AUX' or 'BND' keyword varName : variable name

Definition at line 704 of file TimeSeriesManager.f90.

    ! -- return
    logical :: found
    ! -- dummy
    type(TimeSeriesManagerType), intent(inout) :: tsManager
    integer(I4B), intent(in) :: ii
    integer(I4B), intent(in) :: jj
    character(len=*), intent(in) :: pkgName
    character(len=3), intent(in) :: auxOrBnd
    character(len=*), intent(in) :: varName
    ! -- local
    integer(I4B) :: i
    integer(I4B) :: nlinks
    integer(I4B) :: removeLink
    type(TimeSeriesLinkType), pointer :: tslTemp => null()
    !
    ! -- determine if link exists
    nlinks = tsmanager%CountLinks(auxorbnd)
    found = .false.
    removelink = -1
    csearchlinks: do i = 1, nlinks
      tsltemp => tsmanager%GetLink(auxorbnd, i)
      !
      ! -- Check ii against iRow, jj against jCol, and varName
      !    against Text member of link
      if (tsltemp%PackageName == pkgname) then
        !
        ! -- This array element is already linked to a time series.
        if (tsltemp%IRow == ii .and. tsltemp%JCol == jj .and. &
            same_word(tsltemp%Text, varname)) then
          found = .true.
          removelink = i
          exit csearchlinks
        end if
      end if
    end do csearchlinks
    !
    ! -- remove link if it was found
    if (removelink > 0) then
      if (auxorbnd == 'BND') then
        call tsmanager%boundTsLinks%RemoveNode(removelink, .true.)
      else if (auxorbnd == 'AUX') then
        call tsmanager%auxvarTsLinks%RemoveNode(removelink, .true.)
      end if
    end if
    !
    ! -- Return
    return

Here is the call graph for this function:

Here is the caller graph for this function:

reset()

subroutine timeseriesmanagermodule::reset ( class(timeseriesmanagertype)  this, character(len=*), intent(in)  pkgName  )

private

Definition at line 346 of file TimeSeriesManager.f90.

    ! -- dummy
    class(TimeSeriesManagerType) :: this
    character(len=*), intent(in) :: pkgName
    ! -- local
    integer(I4B) :: i, nlinks
    type(TimeSeriesLinkType), pointer :: tslink
    !
    ! Zero out values for time-series controlled stresses.
    ! Also deallocate all tslinks too.
    ! Then when time series are
    ! specified in this or another stress period,
    ! a new tslink would be set up.
    !
    ! Reassign all linked elements to zero
    nlinks = this%boundTsLinks%Count()
    do i = 1, nlinks
      tslink => gettimeserieslinkfromlist(this%boundTsLinks, i)
      if (associated(tslink)) then
        if (tslink%PackageName == pkgname) then
          tslink%BndElement = dzero
        end if
      end if
    end do
    !
    ! Remove links belonging to calling package
    nlinks = this%boundTsLinks%Count()
    do i = nlinks, 1, -1
      tslink => gettimeserieslinkfromlist(this%boundTsLinks, i)
      if (associated(tslink)) then
        if (tslink%PackageName == pkgname) then
          call this%boundTsLinks%RemoveNode(i, .true.)
        end if
      end if
    end do
    nlinks = this%auxvarTsLinks%Count()
    do i = nlinks, 1, -1
      tslink => gettimeserieslinkfromlist(this%auxvarTsLinks, i)
      if (associated(tslink)) then
        if (tslink%PackageName == pkgname) then
          call this%auxvarTsLinks%RemoveNode(i, .true.)
        end if
      end if
    end do
    !
    ! -- Return
    return

Here is the call graph for this function:

tsmanager_cr()

subroutine, public timeseriesmanagermodule::tsmanager_cr	(	type(timeseriesmanagertype)	this,
		integer(i4b), intent(in)	iout,
		logical, intent(in), optional	removeTsLinksOnCompletion,
		logical, intent(in), optional	extendTsToEndOfSimulation
	)

Definition at line 62 of file TimeSeriesManager.f90.

    ! -- dummy
    type(TimeSeriesManagerType) :: this
    integer(I4B), intent(in) :: iout
    logical, intent(in), optional :: removeTsLinksOnCompletion
    logical, intent(in), optional :: extendTsToEndOfSimulation
    !
    this%iout = iout
    if (present(removetslinksoncompletion)) then
      this%removeTsLinksOnCompletion = removetslinksoncompletion
    end if
    if (present(extendtstoendofsimulation)) then
      this%extendTsToEndOfSimulation = extendtstoendofsimulation
    end if
    allocate (this%boundTsLinks)
    allocate (this%auxvarTsLinks)
    allocate (this%tsfileList)
    allocate (this%tsfiles(1000))
    !
    ! -- Return
    return

Here is the caller graph for this function:

tsmanager_df()

subroutine timeseriesmanagermodule::tsmanager_df ( class(timeseriesmanagertype)  this )

private

Definition at line 88 of file TimeSeriesManager.f90.

    ! -- dummy
    class(TimeSeriesManagerType) :: this
    !
    if (this%numtsfiles > 0) then
      call this%HashBndTimeSeries()
    end if
    !
    ! -- Return
    return

tsmgr_ad()

subroutine timeseriesmanagermodule::tsmgr_ad

(

class(timeseriesmanagertype)

this

)

Definition at line 143 of file TimeSeriesManager.f90.

    ! -- dummy
    class(TimeSeriesManagerType) :: this
    ! -- local
    type(TimeSeriesLinkType), pointer :: tsLink => null()
    type(TimeSeriesType), pointer :: timeseries => null()
    integer(I4B) :: i, nlinks, nauxlinks
    real(DP) :: begintime, endtime, tsendtime
    character(len=LENPACKAGENAME + 2) :: pkgID
    ! -- formats
    character(len=*), parameter :: fmt5 = &
      "(/,'Time-series controlled values in stress period: ', i0, &
      &', time step ', i0, ':')"
10  format(a, ' package: Boundary ', i0, ', entry ', i0, &
           ' value from time series "', a, '" = ', g12.5)
15  format(a, ' package: Boundary ', i0, ', entry ', i0, &
           ' value from time series "', a, '" = ', g12.5, ' (', a, ')')
20  format(a, ' package: Boundary ', i0, ', ', a, &
           ' value from time series "', a, '" = ', g12.5)
25  format(a, ' package: Boundary ', i0, ', ', a, &
           ' value from time series "', a, '" = ', g12.5, ' (', a, ')')
    !
    ! -- Initialize time variables
    begintime = totimc
    endtime = begintime + delt
    !
    ! -- Determine number of ts links
    nlinks = this%boundtslinks%Count()
    nauxlinks = this%auxvartslinks%Count()
    !
    ! -- Iterate through auxvartslinks and replace specified
    !    elements of auxvar with average value obtained from
    !    appropriate time series.  Need to do auxvartslinks
    !    first because they may be a multiplier column
    i = 1
    do while (i <= nauxlinks)
      tslink => gettimeserieslinkfromlist(this%auxvarTsLinks, i)
      timeseries => tslink%timeSeries
      !
      ! -- Remove time series link once its end time has passed, if requested
      if (this%removeTsLinksOnCompletion) then
        tsendtime = timeseries%FindLatestTime(.true.)
        if (tsendtime < begintime) then
          call this%auxvarTsLinks%RemoveNode(i, .true.)
          nauxlinks = this%auxvartslinks%Count()
          cycle
        end if
      end if
      !
      if (i == 1) then
        if (tslink%Iprpak == 1) then
          write (this%iout, fmt5) kper, kstp
        end if
      end if
      tslink%BndElement = timeseries%GetValue(begintime, endtime, &
                                              this%extendTsToEndOfSimulation)
      !
      ! -- Write time series values to output file
      if (tslink%Iprpak == 1) then
        pkgid = '"'//trim(tslink%PackageName)//'"'
        if (tslink%Text == '') then
          if (tslink%BndName == '') then
            write (this%iout, 10) trim(pkgid), tslink%IRow, tslink%JCol, &
              trim(tslink%timeSeries%Name), &
              tslink%BndElement
          else
            write (this%iout, 15) trim(pkgid), tslink%IRow, tslink%JCol, &
              trim(tslink%timeSeries%Name), &
              tslink%BndElement, trim(tslink%BndName)
          end if
        else
          if (tslink%BndName == '') then
            write (this%iout, 20) trim(pkgid), tslink%IRow, trim(tslink%Text), &
              trim(tslink%timeSeries%Name), &
              tslink%BndElement
          else
            write (this%iout, 25) trim(pkgid), tslink%IRow, trim(tslink%Text), &
              trim(tslink%timeSeries%Name), &
              tslink%BndElement, trim(tslink%BndName)
          end if
        end if
      end if
      !
      i = i + 1
    end do
    !
    ! -- Iterate through boundtslinks and replace specified
    !    elements of bound with average value obtained from
    !    appropriate time series. (For list-type packages)
    i = 1
    do while (i <= nlinks)
      tslink => gettimeserieslinkfromlist(this%boundTsLinks, i)
      timeseries => tslink%timeSeries
      !
      ! -- Remove time series link once its end time has passed, if requested
      if (this%removeTsLinksOnCompletion) then
        tsendtime = timeseries%FindLatestTime(.true.)
        if (tsendtime < begintime) then
          call this%boundTsLinks%RemoveNode(i, .true.)
          nlinks = this%boundTsLinks%Count()
          cycle
        end if
      end if
      !
      if (i == 1 .and. nauxlinks == 0) then
        if (tslink%Iprpak == 1) then
          write (this%iout, fmt5) kper, kstp
        end if
      end if
      ! this part needs to be different for MAW because MAW does not use
      ! bound array for well rate (although rate is stored in
      ! this%bound(4,ibnd)), it uses this%mawwells(n)%rate%value
      if (tslink%UseDefaultProc) then
        timeseries => tslink%timeSeries
        tslink%BndElement = timeseries%GetValue(begintime, endtime, &
                                                this%extendTsToEndOfSimulation)
        !
        ! -- If multiplier is active and it applies to this element,
        !    do the multiplication.  This must be done after the auxlinks
        !    have been calculated in case iauxmultcol is being used.
        if (associated(tslink%RMultiplier)) then
          tslink%BndElement = tslink%BndElement * tslink%RMultiplier
        end if
        !
        ! -- Write time series values to output files
        if (tslink%Iprpak == 1) then
          pkgid = '"'//trim(tslink%PackageName)//'"'
          if (tslink%Text == '') then
            if (tslink%BndName == '') then
              write (this%iout, 10) trim(pkgid), tslink%IRow, tslink%JCol, &
                trim(tslink%timeSeries%Name), &
                tslink%BndElement
            else
              write (this%iout, 15) trim(pkgid), tslink%IRow, tslink%JCol, &
                trim(tslink%timeSeries%Name), &
                tslink%BndElement, trim(tslink%BndName)
            end if
          else
            if (tslink%BndName == '') then
              write (this%iout, 20) trim(pkgid), tslink%IRow, trim(tslink%Text), &
                trim(tslink%timeSeries%Name), &
                tslink%BndElement
            else
              write (this%iout, 25) trim(pkgid), tslink%IRow, trim(tslink%Text), &
                trim(tslink%timeSeries%Name), &
                tslink%BndElement, trim(tslink%BndName)
            end if
          end if
        end if
        !
        ! -- If conversion from flux to flow is required, multiply by cell area
        if (tslink%ConvertFlux) then
          tslink%BndElement = tslink%BndElement * tslink%CellArea
        end if
      end if
      !
      i = i + 1
    end do
    !
    ! -- Finish with ending line
    if (nlinks + nauxlinks > 0) then
      if (tslink%Iprpak == 1) then
        write (this%iout, '()')
      end if
    end if
    !
    ! -- Return
    return

Here is the call graph for this function:

tsmgr_da()

subroutine timeseriesmanagermodule::tsmgr_da ( class(timeseriesmanagertype)  this )

private

Definition at line 315 of file TimeSeriesManager.f90.

    ! -- dummy
    class(TimeSeriesManagerType) :: this
    ! -- local
    !
    ! -- Deallocate time-series links in boundTsLinks
    call this%boundTsLinks%Clear(.true.)
    deallocate (this%boundTsLinks)
    !
    ! -- Deallocate time-series links in auxvarTsLinks
    call this%auxvarTsLinks%Clear(.true.)
    deallocate (this%auxvarTsLinks)
    !
    ! -- Deallocate tsfileList
    call this%tsfileList%da()
    deallocate (this%tsfileList)
    !
    ! -- Deallocate the hash table
    if (associated(this%BndTsHashTable)) then
      call hash_table_da(this%BndTsHashTable)
    end if
    !
    deallocate (this%tsfiles)
    !
    ! -- Return
    return

Here is the call graph for this function:

var_timeseries()

logical function, public timeseriesmanagermodule::var_timeseries	(	type(timeseriesmanagertype), intent(inout)	tsManager,
		character(len=*), intent(in)	pkgName,
		character(len=*), intent(in)	varName,
		character(len=3), intent(in), optional	auxOrBnd
	)

Arguments are as follows: tsManager : timeseries manager object for package pkgName : package name varName : variable name auxOrBnd : optional 'AUX' or 'BND' keyword

Definition at line 763 of file TimeSeriesManager.f90.

    ! -- return
    logical :: tsexists
    ! -- dummy
    type(TimeSeriesManagerType), intent(inout) :: tsManager
    character(len=*), intent(in) :: pkgName
    character(len=*), intent(in) :: varName
    character(len=3), intent(in), optional :: auxOrBnd
    ! -- local
    character(len=3) :: ctstype
    integer(I4B) :: i
    integer(I4B) :: nlinks
    type(TimeSeriesLinkType), pointer :: tslTemp => null()
    !
    ! -- process optional variables
    if (present(auxorbnd)) then
      ctstype = auxorbnd
    else
      ctstype = 'BND'
    end if
    !
    ! -- initialize the return variable and the number of timeseries links
    tsexists = .false.
    nlinks = tsmanager%CountLinks(ctstype)
    !
    ! -- determine if link exists
    csearchlinks: do i = 1, nlinks
      tsltemp => tsmanager%GetLink(ctstype, i)
      if (tsltemp%PackageName == pkgname) then
        !
        ! -- Check varName against Text member of link
        if (same_word(tsltemp%Text, varname)) then
          tsexists = .true.
          exit csearchlinks
        end if
      end if
    end do csearchlinks
    !
    ! -- Return
    return

Here is the call graph for this function:

Here is the caller graph for this function: