Data Types
type	gwtlkttype

Functions/Subroutines
subroutine, public	lkt_create (packobj, id, ibcnum, inunit, iout, namemodel, pakname, fmi, eqnsclfac, dvt, dvu, dvua)
	Create a new lkt package. More...

subroutine	find_lkt_package (this)
	Find corresponding lkt package. More...

subroutine	lkt_fc_expanded (this, rhs, ia, idxglo, matrix_sln)
	Add matrix terms related to LKT. More...

subroutine	lkt_solve (this)
	Add terms specific to lakes to the explicit lake solve. More...

integer(i4b) function	lkt_get_nbudterms (this)
	Function to return the number of budget terms just for this package. More...

subroutine	lkt_setup_budobj (this, idx)
	Set up the budget object that stores all the lake flows. More...

subroutine	lkt_fill_budobj (this, idx, x, flowja, ccratin, ccratout)
	Copy flow terms into thisbudobj. More...

subroutine	allocate_scalars (this)
	Allocate scalars specific to the lake mass transport (LKT) package. More...

subroutine	lkt_allocate_arrays (this)
	Allocate arrays specific to the lake mass transport (LKT) package. More...

subroutine	lkt_da (this)
	Deallocate memory. More...

subroutine	lkt_rain_term (this, ientry, n1, n2, rrate, rhsval, hcofval)
	Rain term. More...

subroutine	lkt_evap_term (this, ientry, n1, n2, rrate, rhsval, hcofval)
	Evaporative term. More...

subroutine	lkt_roff_term (this, ientry, n1, n2, rrate, rhsval, hcofval)
	Runoff term. More...

subroutine	lkt_iflw_term (this, ientry, n1, n2, rrate, rhsval, hcofval)
	Inflow Term. More...

subroutine	lkt_wdrl_term (this, ientry, n1, n2, rrate, rhsval, hcofval)
	Specified withdrawal term. More...

subroutine	lkt_outf_term (this, ientry, n1, n2, rrate, rhsval, hcofval)
	Outflow term. More...

subroutine	lkt_df_obs (this)
	Defined observation types. More...

subroutine	lkt_rp_obs (this, obsrv, found)
	Process package specific obs. More...

subroutine	lkt_bd_obs (this, obstypeid, jj, v, found)
	Calculate observation value and pass it back to APT. More...

subroutine	lkt_set_stressperiod (this, itemno, keyword, found)
	Sets the stress period attributes for keyword use. More...

Variables
character(len= *), parameter	ftype = 'LKT'

character(len= *), parameter	flowtype = 'LAK'

character(len=16)	text = ' LKT'

Function/Subroutine Documentation

◆ allocate_scalars()

subroutine gwtlktmodule::allocate_scalars ( class(gwtlkttype) this )

Definition at line 606 of file gwt-lkt.f90.

     ! -- modules
     use memorymanagermodule, only: mem_allocate
     ! -- dummy
     class(GwtLktType) :: this
     ! -- local
     !
     ! -- allocate scalars in TspAptType
     call this%TspAptType%allocate_scalars()
     !
     ! -- Allocate
     call mem_allocate(this%idxbudrain, 'IDXBUDRAIN', this%memoryPath)
     call mem_allocate(this%idxbudevap, 'IDXBUDEVAP', this%memoryPath)
     call mem_allocate(this%idxbudroff, 'IDXBUDROFF', this%memoryPath)
     call mem_allocate(this%idxbudiflw, 'IDXBUDIFLW', this%memoryPath)
     call mem_allocate(this%idxbudwdrl, 'IDXBUDWDRL', this%memoryPath)
     call mem_allocate(this%idxbudoutf, 'IDXBUDOUTF', this%memoryPath)
     !
     ! -- Initialize
     this%idxbudrain = 0
     this%idxbudevap = 0
     this%idxbudroff = 0
     this%idxbudiflw = 0
     this%idxbudwdrl = 0
     this%idxbudoutf = 0

◆ find_lkt_package()

subroutine gwtlktmodule::find_lkt_package ( class(gwtlkttype) this )

Definition at line 152 of file gwt-lkt.f90.

     ! -- modules
     use memorymanagermodule, only: mem_allocate
     ! -- dummy
     class(GwtLktType) :: this
     ! -- local
     character(len=LINELENGTH) :: errmsg
     class(BndType), pointer :: packobj
     integer(I4B) :: ip, icount
     integer(I4B) :: nbudterm
     logical :: found
     !
     ! -- Initialize found to false, and error later if flow package cannot
     !    be found
     found = .false.
     !
     ! -- If user is specifying flows in a binary budget file, then set up
     !    the budget file reader, otherwise set a pointer to the flow package
     !    budobj
     if (this%fmi%flows_from_file) then
       call this%fmi%set_aptbudobj_pointer(this%flowpackagename, this%flowbudptr)
       if (associated(this%flowbudptr)) found = .true.
       !
     else
       if (associated(this%fmi%gwfbndlist)) then
         ! -- Look through gwfbndlist for a flow package with the same name as
         !    this transport package name
         do ip = 1, this%fmi%gwfbndlist%Count()
           packobj => getbndfromlist(this%fmi%gwfbndlist, ip)
           if (packobj%packName == this%flowpackagename) then
             found = .true.
             !
             ! -- store BndType pointer to packobj, and then
             !    use the select type to point to the budobj in flow package
             this%flowpackagebnd => packobj
             select type (packobj)
             type is (laktype)
               this%flowbudptr => packobj%budobj
             end select
           end if
           if (found) exit
         end do
       end if
     end if
     !
     ! -- error if flow package not found
     if (.not. found) then
       write (errmsg, '(a)') 'Could not find flow package with name '&
                             &//trim(adjustl(this%flowpackagename))//'.'
       call store_error(errmsg)
       call this%parser%StoreErrorUnit()
     end if
     !
     ! -- allocate space for idxbudssm, which indicates whether this is a
     !    special budget term or one that is a general source and sink
     nbudterm = this%flowbudptr%nbudterm
     call mem_allocate(this%idxbudssm, nbudterm, 'IDXBUDSSM', this%memoryPath)
     !
     ! -- Process budget terms and identify special budget terms
     write (this%iout, '(/, a, a)') &
       'PROCESSING '//ftype//' INFORMATION FOR ', this%packName
     write (this%iout, '(a)') '  IDENTIFYING FLOW TERMS IN '//flowtype//' PACKAGE'
     write (this%iout, '(a, i0)') &
       '  NUMBER OF '//flowtype//' = ', this%flowbudptr%ncv
     icount = 1
     do ip = 1, this%flowbudptr%nbudterm
       select case (trim(adjustl(this%flowbudptr%budterm(ip)%flowtype)))
       case ('FLOW-JA-FACE')
         this%idxbudfjf = ip
         this%idxbudssm(ip) = 0
       case ('GWF')
         this%idxbudgwf = ip
         this%idxbudssm(ip) = 0
       case ('STORAGE')
         this%idxbudsto = ip
         this%idxbudssm(ip) = 0
       case ('RAINFALL')
         this%idxbudrain = ip
         this%idxbudssm(ip) = 0
       case ('EVAPORATION')
         this%idxbudevap = ip
         this%idxbudssm(ip) = 0
       case ('RUNOFF')
         this%idxbudroff = ip
         this%idxbudssm(ip) = 0
       case ('EXT-INFLOW')
         this%idxbudiflw = ip
         this%idxbudssm(ip) = 0
       case ('WITHDRAWAL')
         this%idxbudwdrl = ip
         this%idxbudssm(ip) = 0
       case ('EXT-OUTFLOW')
         this%idxbudoutf = ip
         this%idxbudssm(ip) = 0
       case ('TO-MVR')
         this%idxbudtmvr = ip
         this%idxbudssm(ip) = 0
       case ('FROM-MVR')
         this%idxbudfmvr = ip
         this%idxbudssm(ip) = 0
       case ('AUXILIARY')
         this%idxbudaux = ip
         this%idxbudssm(ip) = 0
       case default
         !
         ! -- set idxbudssm equal to a column index for where the concentrations
         !    are stored in the concbud(nbudssm, ncv) array
         this%idxbudssm(ip) = icount
         icount = icount + 1
       end select
       write (this%iout, '(a, i0, " = ", a,/, a, i0)') &
         '  TERM ', ip, trim(adjustl(this%flowbudptr%budterm(ip)%flowtype)), &
         '   MAX NO. OF ENTRIES = ', this%flowbudptr%budterm(ip)%maxlist
     end do
     write (this%iout, '(a, //)') 'DONE PROCESSING '//ftype//' INFORMATION'

Here is the call graph for this function:

◆ lkt_allocate_arrays()

subroutine gwtlktmodule::lkt_allocate_arrays ( class(gwtlkttype), intent(inout) this )

Definition at line 636 of file gwt-lkt.f90.

     ! -- modules
     use memorymanagermodule, only: mem_allocate
     ! -- dummy
     class(GwtLktType), intent(inout) :: this
     ! -- local
     integer(I4B) :: n
     !
     ! -- time series
     call mem_allocate(this%concrain, this%ncv, 'CONCRAIN', this%memoryPath)
     call mem_allocate(this%concevap, this%ncv, 'CONCEVAP', this%memoryPath)
     call mem_allocate(this%concroff, this%ncv, 'CONCROFF', this%memoryPath)
     call mem_allocate(this%conciflw, this%ncv, 'CONCIFLW', this%memoryPath)
     !
     ! -- call standard TspAptType allocate arrays
     call this%TspAptType%apt_allocate_arrays()
     !
     ! -- Initialize
     do n = 1, this%ncv
       this%concrain(n) = dzero
       this%concevap(n) = dzero
       this%concroff(n) = dzero
       this%conciflw(n) = dzero
     end do
     !

◆ lkt_bd_obs()

subroutine gwtlktmodule::lkt_bd_obs	(	class(gwtlkttype), intent(inout)	this,
		character(len=*), intent(in)	obstypeid,
		integer(i4b), intent(in)	jj,
		real(dp), intent(inout)	v,
		logical, intent(inout)	found
	)

Definition at line 971 of file gwt-lkt.f90.

     ! -- dummy
     class(GwtLktType), intent(inout) :: this
     character(len=*), intent(in) :: obstypeid
     real(DP), intent(inout) :: v
     integer(I4B), intent(in) :: jj
     logical, intent(inout) :: found
     ! -- local
     integer(I4B) :: n1, n2
     !
     found = .true.
     select case (obstypeid)
     case ('RAINFALL')
       if (this%iboundpak(jj) /= 0) then
         call this%lkt_rain_term(jj, n1, n2, v)
       end if
     case ('EVAPORATION')
       if (this%iboundpak(jj) /= 0) then
         call this%lkt_evap_term(jj, n1, n2, v)
       end if
     case ('RUNOFF')
       if (this%iboundpak(jj) /= 0) then
         call this%lkt_roff_term(jj, n1, n2, v)
       end if
     case ('EXT-INFLOW')
       if (this%iboundpak(jj) /= 0) then
         call this%lkt_iflw_term(jj, n1, n2, v)
       end if
     case ('WITHDRAWAL')
       if (this%iboundpak(jj) /= 0) then
         call this%lkt_wdrl_term(jj, n1, n2, v)
       end if
     case ('EXT-OUTFLOW')
       if (this%iboundpak(jj) /= 0) then
         call this%lkt_outf_term(jj, n1, n2, v)
       end if
     case default
       found = .false.
     end select

◆ lkt_create()

subroutine, public gwtlktmodule::lkt_create	(	class(bndtype), pointer	packobj,
		integer(i4b), intent(in)	id,
		integer(i4b), intent(in)	ibcnum,
		integer(i4b), intent(in)	inunit,
		integer(i4b), intent(in)	iout,
		character(len=*), intent(in)	namemodel,
		character(len=*), intent(in)	pakname,
		type(tspfmitype), pointer	fmi,
		real(dp), intent(in), pointer	eqnsclfac,
		character(len=*), intent(in)	dvt,
		character(len=*), intent(in)	dvu,
		character(len=*), intent(in)	dvua
	)

Parameters

[in]	eqnsclfac	governing equation scale factor
[in]	dvt	For GWT, set to "CONCENTRATION" in TspAptType
[in]	dvu	For GWT, set to "mass" in TspAptType
[in]	dvua	For GWT, set to "M" in TspAptType

Definition at line 97 of file gwt-lkt.f90.

     ! -- dummy
     class(BndType), pointer :: packobj
     integer(I4B), intent(in) :: id
     integer(I4B), intent(in) :: ibcnum
     integer(I4B), intent(in) :: inunit
     integer(I4B), intent(in) :: iout
     character(len=*), intent(in) :: namemodel
     character(len=*), intent(in) :: pakname
     type(TspFmiType), pointer :: fmi
     real(DP), intent(in), pointer :: eqnsclfac !< governing equation scale factor
     character(len=*), intent(in) :: dvt !< For GWT, set to "CONCENTRATION" in TspAptType
     character(len=*), intent(in) :: dvu !< For GWT, set to "mass" in TspAptType
     character(len=*), intent(in) :: dvua !< For GWT, set to "M" in TspAptType
     ! -- local
     type(GwtLktType), pointer :: lktobj
     !
     ! -- allocate the object and assign values to object variables
     allocate (lktobj)
     packobj => lktobj
     !
     ! -- create name and memory path
     call packobj%set_names(ibcnum, namemodel, pakname, ftype)
     packobj%text = text
     !
     ! -- allocate scalars
     call lktobj%allocate_scalars()
     !
     ! -- initialize package
     call packobj%pack_initialize()
  
     packobj%inunit = inunit
     packobj%iout = iout
     packobj%id = id
     packobj%ibcnum = ibcnum
     packobj%ncolbnd = 1
     packobj%iscloc = 1
  
     ! -- Store pointer to flow model interface.  When the GwfGwt exchange is
     !    created, it sets fmi%bndlist so that the GWT model has access to all
     !    the flow packages
     lktobj%fmi => fmi
     !
     ! -- Store pointer to governing equation scale factor
     lktobj%eqnsclfac => eqnsclfac
     !
     ! -- Set labels that will be used in generalized APT class
     lktobj%depvartype = dvt
     lktobj%depvarunit = dvu
     lktobj%depvarunitabbrev = dvua

Here is the caller graph for this function:

◆ lkt_da()

subroutine gwtlktmodule::lkt_da ( class(gwtlkttype) this )

Definition at line 665 of file gwt-lkt.f90.

     ! -- modules
     use memorymanagermodule, only: mem_deallocate
     ! -- dummy
     class(GwtLktType) :: this
     ! -- local
     !
     ! -- deallocate scalars
     call mem_deallocate(this%idxbudrain)
     call mem_deallocate(this%idxbudevap)
     call mem_deallocate(this%idxbudroff)
     call mem_deallocate(this%idxbudiflw)
     call mem_deallocate(this%idxbudwdrl)
     call mem_deallocate(this%idxbudoutf)
     !
     ! -- deallocate time series
     call mem_deallocate(this%concrain)
     call mem_deallocate(this%concevap)
     call mem_deallocate(this%concroff)
     call mem_deallocate(this%conciflw)
     !
     ! -- deallocate scalars in TspAptType
     call this%TspAptType%bnd_da()

◆ lkt_df_obs()

subroutine gwtlktmodule::lkt_df_obs ( class(gwtlkttype) this )

Store the observation type supported by the APT package and override BndTypebnd_df_obs

Definition at line 863 of file gwt-lkt.f90.

     ! -- modules
     ! -- dummy
     class(GwtLktType) :: this
     ! -- local
     integer(I4B) :: indx
     !
     ! -- Store obs type and assign procedure pointer
     !    for concentration observation type.
     call this%obs%StoreObsType('concentration', .false., indx)
     this%obs%obsData(indx)%ProcessIdPtr => apt_process_obsid
     !
     ! -- Store obs type and assign procedure pointer
     !    for flow between features, such as lake to lake.
     call this%obs%StoreObsType('flow-ja-face', .true., indx)
     this%obs%obsData(indx)%ProcessIdPtr => apt_process_obsid12
     !
     ! -- Store obs type and assign procedure pointer
     !    for from-mvr observation type.
     call this%obs%StoreObsType('from-mvr', .true., indx)
     this%obs%obsData(indx)%ProcessIdPtr => apt_process_obsid
     !
     ! -- Store obs type and assign procedure pointer
     !    for to-mvr observation type.
     call this%obs%StoreObsType('to-mvr', .true., indx)
     this%obs%obsData(indx)%ProcessIdPtr => apt_process_obsid
     !
     ! -- Store obs type and assign procedure pointer
     !    for storage observation type.
     call this%obs%StoreObsType('storage', .true., indx)
     this%obs%obsData(indx)%ProcessIdPtr => apt_process_obsid
     !
     ! -- Store obs type and assign procedure pointer
     !    for constant observation type.
     call this%obs%StoreObsType('constant', .true., indx)
     this%obs%obsData(indx)%ProcessIdPtr => apt_process_obsid
     !
     ! -- Store obs type and assign procedure pointer
     !    for observation type: lkt
     call this%obs%StoreObsType('lkt', .true., indx)
     this%obs%obsData(indx)%ProcessIdPtr => apt_process_obsid12
     !
     ! -- Store obs type and assign procedure pointer
     !    for rainfall observation type.
     call this%obs%StoreObsType('rainfall', .true., indx)
     this%obs%obsData(indx)%ProcessIdPtr => apt_process_obsid
     !
     ! -- Store obs type and assign procedure pointer
     !    for evaporation observation type.
     call this%obs%StoreObsType('evaporation', .true., indx)
     this%obs%obsData(indx)%ProcessIdPtr => apt_process_obsid
     !
     ! -- Store obs type and assign procedure pointer
     !    for runoff observation type.
     call this%obs%StoreObsType('runoff', .true., indx)
     this%obs%obsData(indx)%ProcessIdPtr => apt_process_obsid
     !
     ! -- Store obs type and assign procedure pointer
     !    for inflow observation type.
     call this%obs%StoreObsType('ext-inflow', .true., indx)
     this%obs%obsData(indx)%ProcessIdPtr => apt_process_obsid
     !
     ! -- Store obs type and assign procedure pointer
     !    for withdrawal observation type.
     call this%obs%StoreObsType('withdrawal', .true., indx)
     this%obs%obsData(indx)%ProcessIdPtr => apt_process_obsid
     !
     ! -- Store obs type and assign procedure pointer
     !    for ext-outflow observation type.
     call this%obs%StoreObsType('ext-outflow', .true., indx)
     this%obs%obsData(indx)%ProcessIdPtr => apt_process_obsid

Here is the call graph for this function:

◆ lkt_evap_term()

subroutine gwtlktmodule::lkt_evap_term	(	class(gwtlkttype)	this,
		integer(i4b), intent(in)	ientry,
		integer(i4b), intent(inout)	n1,
		integer(i4b), intent(inout)	n2,
		real(dp), intent(inout), optional	rrate,
		real(dp), intent(inout), optional	rhsval,
		real(dp), intent(inout), optional	hcofval
	)

Definition at line 717 of file gwt-lkt.f90.

     ! -- dummy
     class(GwtLktType) :: this
     integer(I4B), intent(in) :: ientry
     integer(I4B), intent(inout) :: n1
     integer(I4B), intent(inout) :: n2
     real(DP), intent(inout), optional :: rrate
     real(DP), intent(inout), optional :: rhsval
     real(DP), intent(inout), optional :: hcofval
     ! -- local
     real(DP) :: qbnd
     real(DP) :: ctmp
     real(DP) :: omega
     !
     n1 = this%flowbudptr%budterm(this%idxbudevap)%id1(ientry)
     n2 = this%flowbudptr%budterm(this%idxbudevap)%id2(ientry)
     ! -- note that qbnd is negative for evap
     qbnd = this%flowbudptr%budterm(this%idxbudevap)%flow(ientry)
     ctmp = this%concevap(n1)
     if (this%xnewpak(n1) < ctmp) then
       omega = done
     else
       omega = dzero
     end if
     if (present(rrate)) &
       rrate = omega * qbnd * this%xnewpak(n1) + &
               (done - omega) * qbnd * ctmp
     if (present(rhsval)) rhsval = -(done - omega) * qbnd * ctmp
     if (present(hcofval)) hcofval = omega * qbnd

◆ lkt_fc_expanded()

subroutine gwtlktmodule::lkt_fc_expanded	(	class(gwtlkttype)	this,
		real(dp), dimension(:), intent(inout)	rhs,
		integer(i4b), dimension(:), intent(in)	ia,
		integer(i4b), dimension(:), intent(in)	idxglo,
		class(matrixbasetype), pointer	matrix_sln
	)

This will be called from TspAptTypeapt_fc_expanded() in order to add matrix terms specifically for LKT

Definition at line 274 of file gwt-lkt.f90.

     ! -- modules
     ! -- dummy
     class(GwtLktType) :: this
     real(DP), dimension(:), intent(inout) :: rhs
     integer(I4B), dimension(:), intent(in) :: ia
     integer(I4B), dimension(:), intent(in) :: idxglo
     class(MatrixBaseType), pointer :: matrix_sln
     ! -- local
     integer(I4B) :: j, n1, n2
     integer(I4B) :: iloc
     integer(I4B) :: iposd
     real(DP) :: rrate
     real(DP) :: rhsval
     real(DP) :: hcofval
     !
     ! -- add rainfall contribution
     if (this%idxbudrain /= 0) then
       do j = 1, this%flowbudptr%budterm(this%idxbudrain)%nlist
         call this%lkt_rain_term(j, n1, n2, rrate, rhsval, hcofval)
         iloc = this%idxlocnode(n1)
         iposd = this%idxpakdiag(n1)
         call matrix_sln%add_value_pos(iposd, hcofval)
         rhs(iloc) = rhs(iloc) + rhsval
       end do
     end if
     !
     ! -- add evaporation contribution
     if (this%idxbudevap /= 0) then
       do j = 1, this%flowbudptr%budterm(this%idxbudevap)%nlist
         call this%lkt_evap_term(j, n1, n2, rrate, rhsval, hcofval)
         iloc = this%idxlocnode(n1)
         iposd = this%idxpakdiag(n1)
         call matrix_sln%add_value_pos(iposd, hcofval)
         rhs(iloc) = rhs(iloc) + rhsval
       end do
     end if
     !
     ! -- add runoff contribution
     if (this%idxbudroff /= 0) then
       do j = 1, this%flowbudptr%budterm(this%idxbudroff)%nlist
         call this%lkt_roff_term(j, n1, n2, rrate, rhsval, hcofval)
         iloc = this%idxlocnode(n1)
         iposd = this%idxpakdiag(n1)
         call matrix_sln%add_value_pos(iposd, hcofval)
         rhs(iloc) = rhs(iloc) + rhsval
       end do
     end if
     !
     ! -- add inflow contribution
     if (this%idxbudiflw /= 0) then
       do j = 1, this%flowbudptr%budterm(this%idxbudiflw)%nlist
         call this%lkt_iflw_term(j, n1, n2, rrate, rhsval, hcofval)
         iloc = this%idxlocnode(n1)
         iposd = this%idxpakdiag(n1)
         call matrix_sln%add_value_pos(iposd, hcofval)
         rhs(iloc) = rhs(iloc) + rhsval
       end do
     end if
     !
     ! -- add withdrawal contribution
     if (this%idxbudwdrl /= 0) then
       do j = 1, this%flowbudptr%budterm(this%idxbudwdrl)%nlist
         call this%lkt_wdrl_term(j, n1, n2, rrate, rhsval, hcofval)
         iloc = this%idxlocnode(n1)
         iposd = this%idxpakdiag(n1)
         call matrix_sln%add_value_pos(iposd, hcofval)
         rhs(iloc) = rhs(iloc) + rhsval
       end do
     end if
     !
     ! -- add outflow contribution
     if (this%idxbudoutf /= 0) then
       do j = 1, this%flowbudptr%budterm(this%idxbudoutf)%nlist
         call this%lkt_outf_term(j, n1, n2, rrate, rhsval, hcofval)
         iloc = this%idxlocnode(n1)
         iposd = this%idxpakdiag(n1)
         call matrix_sln%add_value_pos(iposd, hcofval)
         rhs(iloc) = rhs(iloc) + rhsval
       end do
     end if

◆ lkt_fill_budobj()

subroutine gwtlktmodule::lkt_fill_budobj	(	class(gwtlkttype)	this,
		integer(i4b), intent(inout)	idx,
		real(dp), dimension(:), intent(in)	x,
		real(dp), dimension(:), intent(inout), contiguous	flowja,
		real(dp), intent(inout)	ccratin,
		real(dp), intent(inout)	ccratout
	)

Definition at line 527 of file gwt-lkt.f90.

     ! -- modules
     ! -- dummy
     class(GwtLktType) :: this
     integer(I4B), intent(inout) :: idx
     real(DP), dimension(:), intent(in) :: x
     real(DP), dimension(:), contiguous, intent(inout) :: flowja
     real(DP), intent(inout) :: ccratin
     real(DP), intent(inout) :: ccratout
     ! -- local
     integer(I4B) :: j, n1, n2
     integer(I4B) :: nlist
     real(DP) :: q
     ! -- formats
     !
     ! -- RAIN
     idx = idx + 1
     nlist = this%flowbudptr%budterm(this%idxbudrain)%nlist
     call this%budobj%budterm(idx)%reset(nlist)
     do j = 1, nlist
       call this%lkt_rain_term(j, n1, n2, q)
       call this%budobj%budterm(idx)%update_term(n1, n2, q)
       call this%apt_accumulate_ccterm(n1, q, ccratin, ccratout)
     end do
     !
     ! -- EVAPORATION
     idx = idx + 1
     nlist = this%flowbudptr%budterm(this%idxbudevap)%nlist
     call this%budobj%budterm(idx)%reset(nlist)
     do j = 1, nlist
       call this%lkt_evap_term(j, n1, n2, q)
       call this%budobj%budterm(idx)%update_term(n1, n2, q)
       call this%apt_accumulate_ccterm(n1, q, ccratin, ccratout)
     end do
     !
     ! -- RUNOFF
     idx = idx + 1
     nlist = this%flowbudptr%budterm(this%idxbudroff)%nlist
     call this%budobj%budterm(idx)%reset(nlist)
     do j = 1, nlist
       call this%lkt_roff_term(j, n1, n2, q)
       call this%budobj%budterm(idx)%update_term(n1, n2, q)
       call this%apt_accumulate_ccterm(n1, q, ccratin, ccratout)
     end do
     !
     ! -- EXT-INFLOW
     idx = idx + 1
     nlist = this%flowbudptr%budterm(this%idxbudiflw)%nlist
     call this%budobj%budterm(idx)%reset(nlist)
     do j = 1, nlist
       call this%lkt_iflw_term(j, n1, n2, q)
       call this%budobj%budterm(idx)%update_term(n1, n2, q)
       call this%apt_accumulate_ccterm(n1, q, ccratin, ccratout)
     end do
     !
     ! -- WITHDRAWAL
     idx = idx + 1
     nlist = this%flowbudptr%budterm(this%idxbudwdrl)%nlist
     call this%budobj%budterm(idx)%reset(nlist)
     do j = 1, nlist
       call this%lkt_wdrl_term(j, n1, n2, q)
       call this%budobj%budterm(idx)%update_term(n1, n2, q)
       call this%apt_accumulate_ccterm(n1, q, ccratin, ccratout)
     end do
     !
     ! -- EXT-OUTFLOW
     idx = idx + 1
     nlist = this%flowbudptr%budterm(this%idxbudoutf)%nlist
     call this%budobj%budterm(idx)%reset(nlist)
     do j = 1, nlist
       call this%lkt_outf_term(j, n1, n2, q)
       call this%budobj%budterm(idx)%update_term(n1, n2, q)
       call this%apt_accumulate_ccterm(n1, q, ccratin, ccratout)
     end do

◆ lkt_get_nbudterms()

integer(i4b) function gwtlktmodule::lkt_get_nbudterms ( class(gwtlkttype) this )

This overrides a function in the parent class.

Definition at line 420 of file gwt-lkt.f90.

     ! -- modules
     ! -- dummy
     class(GwtLktType) :: this
     ! -- return
     integer(I4B) :: nbudterms
     ! -- local
     !
     ! -- Number of budget terms is 6
     nbudterms = 6

◆ lkt_iflw_term()

subroutine gwtlktmodule::lkt_iflw_term	(	class(gwtlkttype)	this,
		integer(i4b), intent(in)	ientry,
		integer(i4b), intent(inout)	n1,
		integer(i4b), intent(inout)	n2,
		real(dp), intent(inout), optional	rrate,
		real(dp), intent(inout), optional	rhsval,
		real(dp), intent(inout), optional	hcofval
	)

Accounts for mass flowing into a lake from a connected stream, for example.

Definition at line 779 of file gwt-lkt.f90.

     ! -- dummy
     class(GwtLktType) :: this
     integer(I4B), intent(in) :: ientry
     integer(I4B), intent(inout) :: n1
     integer(I4B), intent(inout) :: n2
     real(DP), intent(inout), optional :: rrate
     real(DP), intent(inout), optional :: rhsval
     real(DP), intent(inout), optional :: hcofval
     ! -- local
     real(DP) :: qbnd
     real(DP) :: ctmp
     !
     n1 = this%flowbudptr%budterm(this%idxbudiflw)%id1(ientry)
     n2 = this%flowbudptr%budterm(this%idxbudiflw)%id2(ientry)
     qbnd = this%flowbudptr%budterm(this%idxbudiflw)%flow(ientry)
     ctmp = this%conciflw(n1)
     if (present(rrate)) rrate = ctmp * qbnd
     if (present(rhsval)) rhsval = -rrate
     if (present(hcofval)) hcofval = dzero

◆ lkt_outf_term()

subroutine gwtlktmodule::lkt_outf_term	(	class(gwtlkttype)	this,
		integer(i4b), intent(in)	ientry,
		integer(i4b), intent(inout)	n1,
		integer(i4b), intent(inout)	n2,
		real(dp), intent(inout), optional	rrate,
		real(dp), intent(inout), optional	rhsval,
		real(dp), intent(inout), optional	hcofval
	)

Accounts for the mass leaving a lake, for example, mass exiting a lake via a flow into a draining stream channel.

Definition at line 835 of file gwt-lkt.f90.

     ! -- dummy
     class(GwtLktType) :: this
     integer(I4B), intent(in) :: ientry
     integer(I4B), intent(inout) :: n1
     integer(I4B), intent(inout) :: n2
     real(DP), intent(inout), optional :: rrate
     real(DP), intent(inout), optional :: rhsval
     real(DP), intent(inout), optional :: hcofval
     ! -- local
     real(DP) :: qbnd
     real(DP) :: ctmp
     !
     n1 = this%flowbudptr%budterm(this%idxbudoutf)%id1(ientry)
     n2 = this%flowbudptr%budterm(this%idxbudoutf)%id2(ientry)
     qbnd = this%flowbudptr%budterm(this%idxbudoutf)%flow(ientry)
     ctmp = this%xnewpak(n1)
     if (present(rrate)) rrate = ctmp * qbnd
     if (present(rhsval)) rhsval = dzero
     if (present(hcofval)) hcofval = qbnd

◆ lkt_rain_term()

subroutine gwtlktmodule::lkt_rain_term	(	class(gwtlkttype)	this,
		integer(i4b), intent(in)	ientry,
		integer(i4b), intent(inout)	n1,
		integer(i4b), intent(inout)	n2,
		real(dp), intent(inout), optional	rrate,
		real(dp), intent(inout), optional	rhsval,
		real(dp), intent(inout), optional	hcofval
	)

Definition at line 692 of file gwt-lkt.f90.

     ! -- dummy
     class(GwtLktType) :: this
     integer(I4B), intent(in) :: ientry
     integer(I4B), intent(inout) :: n1
     integer(I4B), intent(inout) :: n2
     real(DP), intent(inout), optional :: rrate
     real(DP), intent(inout), optional :: rhsval
     real(DP), intent(inout), optional :: hcofval
     ! -- local
     real(DP) :: qbnd
     real(DP) :: ctmp
     !
     n1 = this%flowbudptr%budterm(this%idxbudrain)%id1(ientry)
     n2 = this%flowbudptr%budterm(this%idxbudrain)%id2(ientry)
     qbnd = this%flowbudptr%budterm(this%idxbudrain)%flow(ientry)
     ctmp = this%concrain(n1)
     if (present(rrate)) rrate = ctmp * qbnd
     if (present(rhsval)) rhsval = -rrate
     if (present(hcofval)) hcofval = dzero

◆ lkt_roff_term()

subroutine gwtlktmodule::lkt_roff_term	(	class(gwtlkttype)	this,
		integer(i4b), intent(in)	ientry,
		integer(i4b), intent(inout)	n1,
		integer(i4b), intent(inout)	n2,
		real(dp), intent(inout), optional	rrate,
		real(dp), intent(inout), optional	rhsval,
		real(dp), intent(inout), optional	hcofval
	)

Definition at line 751 of file gwt-lkt.f90.

     ! -- dummy
     class(GwtLktType) :: this
     integer(I4B), intent(in) :: ientry
     integer(I4B), intent(inout) :: n1
     integer(I4B), intent(inout) :: n2
     real(DP), intent(inout), optional :: rrate
     real(DP), intent(inout), optional :: rhsval
     real(DP), intent(inout), optional :: hcofval
     ! -- local
     real(DP) :: qbnd
     real(DP) :: ctmp
     !
     n1 = this%flowbudptr%budterm(this%idxbudroff)%id1(ientry)
     n2 = this%flowbudptr%budterm(this%idxbudroff)%id2(ientry)
     qbnd = this%flowbudptr%budterm(this%idxbudroff)%flow(ientry)
     ctmp = this%concroff(n1)
     if (present(rrate)) rrate = ctmp * qbnd
     if (present(rhsval)) rhsval = -rrate
     if (present(hcofval)) hcofval = dzero

◆ lkt_rp_obs()

subroutine gwtlktmodule::lkt_rp_obs	(	class(gwtlkttype), intent(inout)	this,
		type(observetype), intent(inout)	obsrv,
		logical, intent(inout)	found
	)

Method to process specific observations for this package.

Parameters

[in,out]	this	package class
[in,out]	obsrv	observation object
[in,out]	found	indicate whether observation was found

Definition at line 940 of file gwt-lkt.f90.

     ! -- dummy
     class(GwtLktType), intent(inout) :: this !< package class
     type(ObserveType), intent(inout) :: obsrv !< observation object
     logical, intent(inout) :: found !< indicate whether observation was found
     ! -- local
     !
     found = .true.
     select case (obsrv%ObsTypeId)
     case ('RAINFALL')
       call this%rp_obs_byfeature(obsrv)
     case ('EVAPORATION')
       call this%rp_obs_byfeature(obsrv)
     case ('RUNOFF')
       call this%rp_obs_byfeature(obsrv)
     case ('EXT-INFLOW')
       call this%rp_obs_byfeature(obsrv)
     case ('WITHDRAWAL')
       call this%rp_obs_byfeature(obsrv)
     case ('EXT-OUTFLOW')
       call this%rp_obs_byfeature(obsrv)
     case ('TO-MVR')
       call this%rp_obs_budterm(obsrv, &
                                this%flowbudptr%budterm(this%idxbudtmvr))
     case default
       found = .false.
     end select

◆ lkt_set_stressperiod()

subroutine gwtlktmodule::lkt_set_stressperiod	(	class(gwtlkttype), intent(inout)	this,
		integer(i4b), intent(in)	itemno,
		character(len=*), intent(in)	keyword,
		logical, intent(inout)	found
	)

Definition at line 1014 of file gwt-lkt.f90.

     use timeseriesmanagermodule, only: read_value_or_time_series_adv
     ! -- dummy
     class(GwtLktType), intent(inout) :: this
     integer(I4B), intent(in) :: itemno
     character(len=*), intent(in) :: keyword
     logical, intent(inout) :: found
     ! -- local
     character(len=LINELENGTH) :: text
     integer(I4B) :: ierr
     integer(I4B) :: jj
     real(DP), pointer :: bndElem => null()
     ! -- formats
     !
     ! RAINFALL <rainfall>
     ! EVAPORATION <evaporation>
     ! RUNOFF <runoff>
     ! EXT-INFLOW <inflow>
     ! WITHDRAWAL <withdrawal>
     !
     found = .true.
     select case (keyword)
     case ('RAINFALL')
       ierr = this%apt_check_valid(itemno)
       if (ierr /= 0) then
         goto 999
       end if
       call this%parser%GetString(text)
       jj = 1
       bndelem => this%concrain(itemno)
       call read_value_or_time_series_adv(text, itemno, jj, bndelem, &
                                          this%packName, 'BND', this%tsManager, &
                                          this%iprpak, 'RAINFALL')
     case ('EVAPORATION')
       ierr = this%apt_check_valid(itemno)
       if (ierr /= 0) then
         goto 999
       end if
       call this%parser%GetString(text)
       jj = 1
       bndelem => this%concevap(itemno)
       call read_value_or_time_series_adv(text, itemno, jj, bndelem, &
                                          this%packName, 'BND', this%tsManager, &
                                          this%iprpak, 'EVAPORATION')
     case ('RUNOFF')
       ierr = this%apt_check_valid(itemno)
       if (ierr /= 0) then
         goto 999
       end if
       call this%parser%GetString(text)
       jj = 1
       bndelem => this%concroff(itemno)
       call read_value_or_time_series_adv(text, itemno, jj, bndelem, &
                                          this%packName, 'BND', this%tsManager, &
                                          this%iprpak, 'RUNOFF')
     case ('EXT-INFLOW')
       ierr = this%apt_check_valid(itemno)
       if (ierr /= 0) then
         goto 999
       end if
       call this%parser%GetString(text)
       jj = 1
       bndelem => this%conciflw(itemno)
       call read_value_or_time_series_adv(text, itemno, jj, bndelem, &
                                          this%packName, 'BND', this%tsManager, &
                                          this%iprpak, 'EXT-INFLOW')
     case default
       !
       ! -- keyword not recognized so return to caller with found = .false.
       found = .false.
     end select
     !
 999 continue

Here is the call graph for this function:

◆ lkt_setup_budobj()

subroutine gwtlktmodule::lkt_setup_budobj	(	class(gwtlkttype)	this,
		integer(i4b), intent(inout)	idx
	)

Definition at line 434 of file gwt-lkt.f90.

     ! -- modules
     use constantsmodule, only: lenbudtxt
     ! -- dummy
     class(GwtLktType) :: this
     integer(I4B), intent(inout) :: idx
     ! -- local
     integer(I4B) :: maxlist, naux
     character(len=LENBUDTXT) :: text
     !
     ! -- Addition of mass associated with rainfall directly on lake surface
     text = '        RAINFALL'
     idx = idx + 1
     maxlist = this%flowbudptr%budterm(this%idxbudrain)%maxlist
     naux = 0
     call this%budobj%budterm(idx)%initialize(text, &
                                              this%name_model, &
                                              this%packName, &
                                              this%name_model, &
                                              this%packName, &
                                              maxlist, .false., .false., &
                                              naux)
     !
     ! -- Loss of dissolved mass associated with evaporation when a non-zero
     !    evaporative concentration is specified
     text = '     EVAPORATION'
     idx = idx + 1
     maxlist = this%flowbudptr%budterm(this%idxbudevap)%maxlist
     naux = 0
     call this%budobj%budterm(idx)%initialize(text, &
                                              this%name_model, &
                                              this%packName, &
                                              this%name_model, &
                                              this%packName, &
                                              maxlist, .false., .false., &
                                              naux)
     !
     ! -- Addition of mass associated with runoff that flows to the lake
     text = '          RUNOFF'
     idx = idx + 1
     maxlist = this%flowbudptr%budterm(this%idxbudroff)%maxlist
     naux = 0
     call this%budobj%budterm(idx)%initialize(text, &
                                              this%name_model, &
                                              this%packName, &
                                              this%name_model, &
                                              this%packName, &
                                              maxlist, .false., .false., &
                                              naux)
     !
     ! -- Addition of mass associated with user-specified inflow to the lake
     text = '      EXT-INFLOW'
     idx = idx + 1
     maxlist = this%flowbudptr%budterm(this%idxbudiflw)%maxlist
     naux = 0
     call this%budobj%budterm(idx)%initialize(text, &
                                              this%name_model, &
                                              this%packName, &
                                              this%name_model, &
                                              this%packName, &
                                              maxlist, .false., .false., &
                                              naux)
     !
     ! -- Removal of mass associated with user-specified withdrawal from lake
     text = '      WITHDRAWAL'
     idx = idx + 1
     maxlist = this%flowbudptr%budterm(this%idxbudwdrl)%maxlist
     naux = 0
     call this%budobj%budterm(idx)%initialize(text, &
                                              this%name_model, &
                                              this%packName, &
                                              this%name_model, &
                                              this%packName, &
                                              maxlist, .false., .false., &
                                              naux)
     !
     ! -- Removal of heat associated with outflow from lake that leaves
     !    model domain
     text = '     EXT-OUTFLOW'
     idx = idx + 1
     maxlist = this%flowbudptr%budterm(this%idxbudoutf)%maxlist
     naux = 0
     call this%budobj%budterm(idx)%initialize(text, &
                                              this%name_model, &
                                              this%packName, &
                                              this%name_model, &
                                              this%packName, &
                                              maxlist, .false., .false., &
                                              naux)

◆ lkt_solve()

subroutine gwtlktmodule::lkt_solve ( class(gwtlkttype) this )

Definition at line 359 of file gwt-lkt.f90.

     ! -- dummy
     class(GwtLktType) :: this
     ! -- local
     integer(I4B) :: j
     integer(I4B) :: n1, n2
     real(DP) :: rrate
     !
     ! -- add rainfall contribution
     if (this%idxbudrain /= 0) then
       do j = 1, this%flowbudptr%budterm(this%idxbudrain)%nlist
         call this%lkt_rain_term(j, n1, n2, rrate)
         this%dbuff(n1) = this%dbuff(n1) + rrate
       end do
     end if
     !
     ! -- add evaporation contribution
     if (this%idxbudevap /= 0) then
       do j = 1, this%flowbudptr%budterm(this%idxbudevap)%nlist
         call this%lkt_evap_term(j, n1, n2, rrate)
         this%dbuff(n1) = this%dbuff(n1) + rrate
       end do
     end if
     !
     ! -- add runoff contribution
     if (this%idxbudroff /= 0) then
       do j = 1, this%flowbudptr%budterm(this%idxbudroff)%nlist
         call this%lkt_roff_term(j, n1, n2, rrate)
         this%dbuff(n1) = this%dbuff(n1) + rrate
       end do
     end if
     !
     ! -- add inflow contribution
     if (this%idxbudiflw /= 0) then
       do j = 1, this%flowbudptr%budterm(this%idxbudiflw)%nlist
         call this%lkt_iflw_term(j, n1, n2, rrate)
         this%dbuff(n1) = this%dbuff(n1) + rrate
       end do
     end if
     !
     ! -- add withdrawal contribution
     if (this%idxbudwdrl /= 0) then
       do j = 1, this%flowbudptr%budterm(this%idxbudwdrl)%nlist
         call this%lkt_wdrl_term(j, n1, n2, rrate)
         this%dbuff(n1) = this%dbuff(n1) + rrate
       end do
     end if
     !
     ! -- add outflow contribution
     if (this%idxbudoutf /= 0) then
       do j = 1, this%flowbudptr%budterm(this%idxbudoutf)%nlist
         call this%lkt_outf_term(j, n1, n2, rrate)
         this%dbuff(n1) = this%dbuff(n1) + rrate
       end do
     end if

◆ lkt_wdrl_term()

subroutine gwtlktmodule::lkt_wdrl_term	(	class(gwtlkttype)	this,
		integer(i4b), intent(in)	ientry,
		integer(i4b), intent(inout)	n1,
		integer(i4b), intent(inout)	n2,
		real(dp), intent(inout), optional	rrate,
		real(dp), intent(inout), optional	rhsval,
		real(dp), intent(inout), optional	hcofval
	)

Accounts for mass associated with a withdrawal of water from a lake or group of lakes.

Definition at line 807 of file gwt-lkt.f90.

     ! -- dummy
     class(GwtLktType) :: this
     integer(I4B), intent(in) :: ientry
     integer(I4B), intent(inout) :: n1
     integer(I4B), intent(inout) :: n2
     real(DP), intent(inout), optional :: rrate
     real(DP), intent(inout), optional :: rhsval
     real(DP), intent(inout), optional :: hcofval
     ! -- local
     real(DP) :: qbnd
     real(DP) :: ctmp
     !
     n1 = this%flowbudptr%budterm(this%idxbudwdrl)%id1(ientry)
     n2 = this%flowbudptr%budterm(this%idxbudwdrl)%id2(ientry)
     qbnd = this%flowbudptr%budterm(this%idxbudwdrl)%flow(ientry)
     ctmp = this%xnewpak(n1)
     if (present(rrate)) rrate = ctmp * qbnd
     if (present(rhsval)) rhsval = dzero
     if (present(hcofval)) hcofval = qbnd

Variable Documentation

◆ flowtype

character(len=*), parameter gwtlktmodule::flowtype = 'LAK'

Definition at line 52 of file gwt-lkt.f90.

52 character(len=*), parameter :: flowtype = 'LAK'

◆ ftype

character(len=*), parameter gwtlktmodule::ftype = 'LKT'

Definition at line 51 of file gwt-lkt.f90.

51 character(len=*), parameter :: ftype = 'LKT'

◆ text

character(len=16) gwtlktmodule::text = ' LKT'

Definition at line 53 of file gwt-lkt.f90.

53 character(len=16) :: text = ' LKT'

Data Types

Functions/Subroutines

Variables

Function/Subroutine Documentation

◆ allocate_scalars()

◆ find_lkt_package()

◆ lkt_allocate_arrays()

◆ lkt_bd_obs()

◆ lkt_create()

◆ lkt_da()

◆ lkt_df_obs()

◆ lkt_evap_term()

◆ lkt_fc_expanded()

◆ lkt_fill_budobj()

◆ lkt_get_nbudterms()

◆ lkt_iflw_term()

◆ lkt_outf_term()

◆ lkt_rain_term()

◆ lkt_roff_term()

◆ lkt_rp_obs()

◆ lkt_set_stressperiod()

◆ lkt_setup_budobj()

◆ lkt_solve()

◆ lkt_wdrl_term()

Variable Documentation

◆ flowtype

◆ ftype

◆ text