tspmvtmodule Module Reference

Data Types
type	mvrtermtype
	Type to store mover-transport flow and concentration for one budget term. More...
type	tspmvttype

Functions/Subroutines
subroutine, public	mvt_cr (mvt, name_model, inunit, iout, fmi1, eqnsclfac, depvartype, gwfmodelname1, gwfmodelname2, fmi2)
subroutine	mvt_init (this, name_model, inunit, iout, fmi1, eqnsclfac, depvartype, gwfmodelname1, gwfmodelname2, fmi2)
	Create a new mover transport object. More...
subroutine	mvt_df (this, dis)
	Define mover transport object. More...
subroutine	set_pointer_mvrbudobj (this, mvrbudobj)
	@ brief Set pointer to mvrbudobj More...
subroutine	mvt_ar (this)
	Allocate and read mover-for-transport information. More...
subroutine	mvt_rp (this)
	Read and prepare mover transport object. More...
subroutine	mvt_fc (this, cnew1, cnew2)
	Calculate coefficients and fill amat and rhs. More...
subroutine	mvt_fill_mvrterm (this, cnew1, cnew2)
	Fill mvrterm with provider concentration and flow for each budget term. More...
subroutine	mvt_update_qmfrommvr (this)
	Assign mover transport mass into each receiver's qmfrommvr. More...
subroutine	set_fmi (this, ibudterm, fmi, set_provider)
subroutine	mvt_cc (this, kiter, iend, icnvgmod, cpak, dpak)
	Extra convergence check for mover. More...
subroutine	mvt_bd (this)
	Write mover terms to listing file. More...
subroutine	mvt_ot_saveflow (this, icbcfl, ibudfl)
	Write mover budget terms. More...
subroutine	mvt_ot_printflow (this, icbcfl, ibudfl)
	Print mover flow table. More...
subroutine	mvt_ot_bdsummary (this, ibudfl)
	Write mover budget to listing file. More...
subroutine	mvt_da (this)
	@ brief Deallocate memory More...
subroutine	allocate_scalars (this)
	@ brief Allocate scalar variables for package More...
subroutine	read_options (this)
	Read mover-for-transport options block. More...
subroutine	mvt_setup_budobj (this)
	Set up the budget object that stores all the mvr flows. More...
subroutine	mvt_fill_budobj (this)
	Copy mover-for-transport flow terms into thisbudobj. More...
subroutine	mvt_scan_mvrbudobj (this)
	Determine max number of packages in use. More...
subroutine	mvt_setup_outputtab (this)
	Set up the mover-for-transport output table. More...
subroutine	mvt_print_outputtab (this)
	Set up mover-for-transport output table. More...

Function/Subroutine Documentation

allocate_scalars()

subroutine tspmvtmodule::allocate_scalars

(

class(tspmvttype)

this

)

Method to allocate scalar variables for the MVT package.

Definition at line 569 of file tsp-mvt.f90.

    ! -- modules
    use memorymanagermodule, only: mem_allocate, mem_setptr
    ! -- dummy
    class(TspMvtType) :: this
    !
    ! -- Allocate scalars in NumericalPackageType
    call this%NumericalPackageType%allocate_scalars()
    !
    ! -- Allocate
    call mem_allocate(this%maxpackages, 'MAXPACKAGES', this%memoryPath)
    call mem_allocate(this%ibudgetout, 'IBUDGETOUT', this%memoryPath)
    call mem_allocate(this%ibudcsv, 'IBUDCSV', this%memoryPath)
    !
    ! -- Initialize
    this%maxpackages = 0
    this%ibudgetout = 0
    this%ibudcsv = 0

mvt_ar()

subroutine tspmvtmodule::mvt_ar ( class(tspmvttype)  this )

private

Definition at line 197 of file tsp-mvt.f90.

    ! -- dummy
    class(TspMvtType) :: this
    !
    ! -- Setup the output table
    call this%mvt_setup_outputtab()

mvt_bd()

subroutine tspmvtmodule::mvt_bd ( class(tspmvttype)  this )

private

Definition at line 402 of file tsp-mvt.f90.

    ! -- dummy
    class(TspMvtType) :: this
    !
    ! -- Fill the budget object
    call this%mvt_fill_budobj()

mvt_cc()

subroutine tspmvtmodule::mvt_cc ( class(tspmvttype)  this, integer(i4b), intent(in)  kiter, integer(i4b), intent(in)  iend, integer(i4b), intent(in)  icnvgmod, character(len=lenpakloc), intent(inout)  cpak, real(dp), intent(inout)  dpak  )

private

Definition at line 377 of file tsp-mvt.f90.

    ! -- dummy
    class(TspMvtType) :: this
    integer(I4B), intent(in) :: kiter
    integer(I4B), intent(in) :: iend
    integer(I4B), intent(in) :: icnvgmod
    character(len=LENPAKLOC), intent(inout) :: cpak
    real(DP), intent(inout) :: dpak
    ! -- formats
    character(len=*), parameter :: fmtmvrcnvg = &
      "(/,1x,'MOVER PACKAGE REQUIRES AT LEAST TWO OUTER ITERATIONS. CONVERGE &
      &FLAG HAS BEEN RESET TO FALSE.')"
    !
    ! -- If there are active movers, then at least 2 outers required
    if (associated(this%mvrbudobj)) then
      if (icnvgmod == 1 .and. kiter == 1) then
        dpak = dnodata
        cpak = trim(this%packName)
        write (this%iout, fmtmvrcnvg)
      end if
    end if

mvt_cr()

subroutine, public tspmvtmodule::mvt_cr	(	type(tspmvttype), pointer	mvt,
		character(len=*), intent(in)	name_model,
		integer(i4b), intent(in)	inunit,
		integer(i4b), intent(in)	iout,
		type(tspfmitype), intent(in), pointer	fmi1,
		real(dp), intent(in), pointer	eqnsclfac,
		character(len=lenvarname), intent(in)	depvartype,
		character(len=*), intent(in), optional	gwfmodelname1,
		character(len=*), intent(in), optional	gwfmodelname2,
		type(tspfmitype), intent(in), optional, pointer	fmi2
	)

Parameters

[in]	eqnsclfac	governing equation scale factor
[in]	depvartype	dependent variable type ('concentration' or 'temperature')

Definition at line 83 of file tsp-mvt.f90.

    type(TspMvtType), pointer :: mvt
    character(len=*), intent(in) :: name_model
    integer(I4B), intent(in) :: inunit
    integer(I4B), intent(in) :: iout
    type(TspFmiType), intent(in), pointer :: fmi1
    real(DP), intent(in), pointer :: eqnsclfac !< governing equation scale factor
    character(len=LENVARNAME), intent(in) :: depvartype !< dependent variable type ('concentration' or 'temperature')
    character(len=*), intent(in), optional :: gwfmodelname1
    character(len=*), intent(in), optional :: gwfmodelname2
    type(TspFmiType), intent(in), pointer, optional :: fmi2
 
    allocate (mvt)
    call mvt%mvt_init(name_model, inunit, iout, fmi1, eqnsclfac, &
                      depvartype, gwfmodelname1, gwfmodelname2, fmi2)
 

Here is the caller graph for this function:

mvt_da()

subroutine tspmvtmodule::mvt_da

(

class(tspmvttype)

this

)

Method to deallocate memory for the package.

Definition at line 514 of file tsp-mvt.f90.

    ! -- modules
    use memorymanagermodule, only: mem_deallocate
    ! -- dummy
    class(TspMvtType) :: this
    ! -- local
    integer(I4B) :: i
    !
    ! -- Deallocate arrays if package was active
    if (this%inunit > 0) then
      !
      ! -- mvrterm
      if (associated(this%mvrterm)) then
        do i = 1, size(this%mvrterm)
          deallocate (this%mvrterm(i)%qty)
        end do
        deallocate (this%mvrterm)
      end if
      !
      ! -- Character array
      deallocate (this%paknames)
      !
      ! -- Budget object
      call this%budget%budget_da()
      deallocate (this%budget)
      !
      ! -- Budobj
      call this%budobj%budgetobject_da()
      deallocate (this%budobj)
      nullify (this%budobj)
      !
      ! -- Output table object
      if (associated(this%outputtab)) then
        call this%outputtab%table_da()
        deallocate (this%outputtab)
        nullify (this%outputtab)
      end if
    end if
    !
    ! -- Scalars
    this%fmi1 => null()
    this%fmi1 => null()
    this%mvrbudobj => null()
    call mem_deallocate(this%maxpackages)
    call mem_deallocate(this%ibudgetout)
    call mem_deallocate(this%ibudcsv)
    !
    ! -- Deallocate scalars in NumericalPackageType
    call this%NumericalPackageType%da()

mvt_df()

subroutine tspmvtmodule::mvt_df ( class(tspmvttype)  this, class(disbasetype), intent(in), pointer  dis  )

private

Definition at line 158 of file tsp-mvt.f90.

    ! -- dummy
    class(TspMvtType) :: this
    class(DisBaseType), pointer, intent(in) :: dis
    ! -- formats
    character(len=*), parameter :: fmtmvt = &
      "(1x,/1x,'MVT -- MOVER TRANSPORT PACKAGE, VERSION 1, 4/15/2020', &
      &' INPUT READ FROM UNIT ', i0, //)"
    !
    ! -- Set pointer to dis
    this%dis => dis
    !
    ! -- Print a message identifying the MVT package.
    write (this%iout, fmtmvt) this%inunit
    !
    ! -- Initialize block parser
    call this%parser%Initialize(this%inunit, this%iout)
    !
    ! -- Initialize the budget table writer
    call budget_cr(this%budget, this%memoryPath)
    !
    ! -- Read mvt options
    call this%read_options()

Here is the call graph for this function:

mvt_fc()

subroutine tspmvtmodule::mvt_fc	(	class(tspmvttype)	this,
		real(dp), dimension(:), intent(in), target, contiguous	cnew1,
		real(dp), dimension(:), intent(in), target, contiguous	cnew2
	)

The mvt package adds the mass flow rate to the provider qmfrommvr array. The advanced packages know enough to subtract any mass that is leaving, so the mvt just adds mass coming in from elsewhere. Because the movers change by stress period, their solute effects must be added to the right- hand side of the transport matrix equations.

Definition at line 240 of file tsp-mvt.f90.

    class(TspMvtType) :: this
    real(DP), intent(in), dimension(:), contiguous, target :: cnew1
    real(DP), intent(in), dimension(:), contiguous, target :: cnew2
 
    call this%mvt_fill_mvrterm(cnew1, cnew2)
    call this%mvt_update_qmfrommvr()
 

mvt_fill_budobj()

subroutine tspmvtmodule::mvt_fill_budobj

(

class(tspmvttype)

this

)

Definition at line 724 of file tsp-mvt.f90.

    ! -- dummy
    class(TspMvtType) :: this
    ! -- local
    type(TspFmiType), pointer :: fmi_rc
    integer(I4B) :: nbudterm
    integer(I4B) :: nlist
    integer(I4B) :: irc
    integer(I4B) :: i
    integer(I4B) :: j
    integer(I4B) :: n1, n2
    real(DP) :: cp
    real(DP) :: q
    real(DP) :: rate
    !
    ! -- Go through the water mover budget terms and set up the transport
    !    mover budget terms
    nbudterm = this%mvrbudobj%nbudterm
    do i = 1, nbudterm
      nlist = this%mvrbudobj%budterm(i)%nlist
      call this%set_fmi(i, fmi_rc, .false.)
      if (.not. associated(fmi_rc)) cycle
 
      call fmi_rc%get_package_index(this%mvrbudobj%budterm(i)%text2id2, irc)
      call this%budobj%budterm(i)%reset(nlist)
      do j = 1, nlist
        n1 = this%mvrbudobj%budterm(i)%id1(j)
        n2 = this%mvrbudobj%budterm(i)%id2(j)
        q = this%mvrbudobj%budterm(i)%flow(j)
        cp = this%mvrterm(i)%qty(j)
 
        ! calculate solute mover rate
        rate = dzero
        if (fmi_rc%iatp(irc) /= 0) then
          rate = -q * cp * this%eqnsclfac
        end if
        !
        ! -- Add the rate to the budterm
        call this%budobj%budterm(i)%update_term(n1, n2, rate)
      end do
    end do
    !
    ! --Terms are filled, now accumulate them for this time step
    call this%budobj%accumulate_terms()

mvt_fill_mvrterm()

subroutine tspmvtmodule::mvt_fill_mvrterm ( class(tspmvttype)  this, real(dp), dimension(:), intent(in), target, contiguous  cnew1, real(dp), dimension(:), intent(in), target, contiguous  cnew2  )

private

Definition at line 252 of file tsp-mvt.f90.

    class(TspMvtType) :: this
    real(DP), intent(in), dimension(:), contiguous, target :: cnew1
    real(DP), intent(in), dimension(:), contiguous, target :: cnew2
    ! local
    integer(I4B) :: i, n
    integer(I4B) :: id1, nlist
    integer(I4B) :: ipr
    integer(I4B) :: igwtnode
    integer(I4B) :: nbudterm
    real(DP) :: cp
    real(DP), dimension(:), pointer :: concpak
    real(DP), dimension(:), contiguous, pointer :: cnew_pr !, pointer to provider dep. var.
    type(TspFmiType), pointer :: fmi_pr !< pointer to provider model fmi package
 
    nbudterm = this%mvrbudobj%nbudterm
    do i = 1, nbudterm
      nlist = this%mvrbudobj%budterm(i)%nlist
      if (nlist == 0) cycle
 
      ! fetch fmi package
      call this%set_fmi(i, fmi_pr, .true.)
      if (.not. associated(fmi_pr)) cycle
 
      cnew_pr => cnew1
      if (associated(fmi_pr, this%fmi2)) then
        cnew_pr => cnew2
      end if
 
      call fmi_pr%get_package_index(this%mvrbudobj%budterm(i)%text2id1, ipr)
 
      ! use package concentration for APT
      if (fmi_pr%iatp(ipr) /= 0) then
        concpak => fmi_pr%datp(ipr)%concpack
      end if
 
      do n = 1, nlist
 
        ! lak/sfr/maw/uzf id1 (provider)
        id1 = this%mvrbudobj%budterm(i)%id1(n)
        cp = dzero
        if (fmi_pr%iatp(ipr) /= 0) then
          ! concentration from APT
          cp = concpak(id1)
        else
          igwtnode = fmi_pr%gwfpackages(ipr)%nodelist(id1)
          cp = cnew_pr(igwtnode)
        end if
 
        this%mvrterm(i)%qty(n) = cp
      end do
    end do

mvt_init()

subroutine tspmvtmodule::mvt_init ( class(tspmvttype)  this, character(len=*), intent(in)  name_model, integer(i4b), intent(in)  inunit, integer(i4b), intent(in)  iout, type(tspfmitype), intent(in), pointer  fmi1, real(dp), intent(in), pointer  eqnsclfac, character(len=lenvarname), intent(in)  depvartype, character(len=*), intent(in), optional  gwfmodelname1, character(len=*), intent(in), optional  gwfmodelname2, type(tspfmitype), intent(in), optional, pointer  fmi2  )

private

Parameters

[in]	eqnsclfac	governing equation scale factor
[in]	depvartype	dependent variable type ('concentration' or 'temperature')

Definition at line 104 of file tsp-mvt.f90.

    ! -- dummy
    class(TspMvtType) :: this
    character(len=*), intent(in) :: name_model
    integer(I4B), intent(in) :: inunit
    integer(I4B), intent(in) :: iout
    type(TspFmiType), intent(in), pointer :: fmi1
    real(DP), intent(in), pointer :: eqnsclfac !< governing equation scale factor
    character(len=LENVARNAME), intent(in) :: depvartype !< dependent variable type ('concentration' or 'temperature')
    character(len=*), intent(in), optional :: gwfmodelname1
    character(len=*), intent(in), optional :: gwfmodelname2
    type(TspFmiType), intent(in), pointer, optional :: fmi2
    !
    ! -- Create name and memory path
    call this%set_names(1, name_model, 'MVT', 'MVT')
    !
    ! -- Allocate scalars
    call this%allocate_scalars()
    !
    this%inunit = inunit
    this%iout = iout
    !
    ! -- Assume that this MVT is owned by a GWT Model
    this%fmi1 => fmi1
    this%fmi2 => fmi1
    !
    ! -- Set pointers
    if (present(fmi2)) then
      this%fmi2 => fmi2
    end if
    !
    ! -- Set model names
    if (present(gwfmodelname1)) then
      this%gwfmodelname1 = gwfmodelname1
    end if
    if (present(gwfmodelname2)) then
      this%gwfmodelname2 = gwfmodelname2
    end if
    !
    ! -- Create the budget object
    call budgetobject_cr(this%budobj, 'TRANSPORT MOVER')
    !
    ! -- Store pointer to governing equation scale factor
    this%eqnsclfac => eqnsclfac
    !
    ! -- Store pointer to labels associated with the current model so that the
    !    package has access to the corresponding dependent variable type
    this%depvartype = depvartype
 

Here is the call graph for this function:

mvt_ot_bdsummary()

subroutine tspmvtmodule::mvt_ot_bdsummary ( class(tspmvttype)  this, integer(i4b), intent(in)  ibudfl  )

private

Definition at line 450 of file tsp-mvt.f90.

    ! -- modules
    use tdismodule, only: kstp, kper, delt, totim
    ! -- dummy
    class(TspMvtType) :: this
    integer(I4B), intent(in) :: ibudfl
    ! -- locals
    integer(I4B) :: i, j, n
    real(DP), allocatable, dimension(:) :: ratin, ratout
    !
    ! -- Allocate and initialize ratin/ratout
    allocate (ratin(this%maxpackages), ratout(this%maxpackages))
    do j = 1, this%maxpackages
      ratin(j) = dzero
      ratout(j) = dzero
    end do
    !
    ! -- Accumulate the rates
    do i = 1, this%maxpackages
      do j = 1, this%budobj%nbudterm
        do n = 1, this%budobj%budterm(j)%nlist
          !
          ! -- Provider is inflow to mover
          if (this%paknames(i) == this%budobj%budterm(j)%text2id1) then
            ratin(i) = ratin(i) + this%budobj%budterm(j)%flow(n)
          end if
          !
          ! -- Receiver is outflow from mover
          if (this%paknames(i) == this%budobj%budterm(j)%text2id2) then
            ratout(i) = ratout(i) + this%budobj%budterm(j)%flow(n)
          end if
        end do
      end do
    end do
    !
    ! -- Send rates to budget object
    call this%budget%reset()
    do j = 1, this%maxpackages
      call this%budget%addentry(ratin(j), ratout(j), delt, this%paknames(j))
    end do
    !
    ! -- Write the budget
    call this%budget%finalize_step(delt)
    if (ibudfl /= 0) then
      call this%budget%budget_ot(kstp, kper, this%iout)
    end if
    !
    ! -- Write budget csv
    call this%budget%writecsv(totim)
    !
    ! -- Deallocate
    deallocate (ratin, ratout)
    !
    ! -- Output mvr budget
    !    Not using budobj write_table here because it would result
    !    in a table that has one entry.  A custom table looks
    !    better here with a row for each package.
    !call this%budobj%write_budtable(kstp, kper, this%iout)

mvt_ot_printflow()

subroutine tspmvtmodule::mvt_ot_printflow	(	class(tspmvttype)	this,
		integer(i4b), intent(in)	icbcfl,
		integer(i4b), intent(in)	ibudfl
	)

Definition at line 436 of file tsp-mvt.f90.

    ! -- dummy
    class(TspMvtType) :: this
    integer(I4B), intent(in) :: icbcfl
    integer(I4B), intent(in) :: ibudfl
    !
    ! -- Print the mover flow table
    if (ibudfl /= 0 .and. this%iprflow /= 0) then
      call this%mvt_print_outputtab()
    end if

mvt_ot_saveflow()

subroutine tspmvtmodule::mvt_ot_saveflow ( class(tspmvttype)  this, integer(i4b), intent(in)  icbcfl, integer(i4b), intent(in)  ibudfl  )

private

Definition at line 412 of file tsp-mvt.f90.

    ! -- modules
    use tdismodule, only: kstp, kper, delt, pertim, totim
    ! -- dummy
    class(TspMvttype) :: this
    integer(I4B), intent(in) :: icbcfl
    integer(I4B), intent(in) :: ibudfl
    ! -- locals
    integer(I4B) :: ibinun
    !
    ! -- Save the mover flows from the budobj to a mover binary file
    ibinun = 0
    if (this%ibudgetout /= 0) then
      ibinun = this%ibudgetout
    end if
    if (icbcfl == 0) ibinun = 0
    if (ibinun > 0) then
      call this%budobj%save_flows(this%dis, ibinun, kstp, kper, delt, &
                                  pertim, totim, this%iout)
    end if

mvt_print_outputtab()

subroutine tspmvtmodule::mvt_print_outputtab ( class(tspmvttype), intent(inout)  this )

private

Definition at line 864 of file tsp-mvt.f90.

    ! -- module
    use tdismodule, only: kstp, kper
    ! -- dummy
    class(TspMvttype), intent(inout) :: this
    ! -- local
    character(len=LINELENGTH) :: title
    character(len=LENMODELNAME + LENPACKAGENAME + 1) :: cloc1, cloc2
    integer(I4B) :: i
    integer(I4B) :: n
    integer(I4B) :: inum
    integer(I4B) :: ntabrows
    integer(I4B) :: nlist
    !
    ! -- Determine number of table rows
    ntabrows = 0
    do i = 1, this%budobj%nbudterm
      nlist = this%budobj%budterm(i)%nlist
      ntabrows = ntabrows + nlist
    end do
    !
    ! -- Set table kstp and kper
    call this%outputtab%set_kstpkper(kstp, kper)
    !
    ! -- Add terms and print the table
    title = 'TRANSPORT MOVER PACKAGE ('//trim(this%packName)// &
            ') FLOW RATES'
    call this%outputtab%set_title(title)
    call this%outputtab%set_maxbound(ntabrows)
    !
    ! -- Process each table row
    inum = 1
    do i = 1, this%budobj%nbudterm
      nlist = this%budobj%budterm(i)%nlist
      do n = 1, nlist
        cloc1 = trim(adjustl(this%budobj%budterm(i)%text1id1))//' '// &
                trim(adjustl(this%budobj%budterm(i)%text2id1))
        cloc2 = trim(adjustl(this%budobj%budterm(i)%text1id2))//' '// &
                trim(adjustl(this%budobj%budterm(i)%text2id2))
        call this%outputtab%add_term(inum)
        call this%outputtab%add_term(cloc1)
        call this%outputtab%add_term(this%budobj%budterm(i)%id1(n))
        call this%outputtab%add_term(-this%mvrbudobj%budterm(i)%flow(n))
        call this%outputtab%add_term(this%budobj%budterm(i)%flow(n))
        call this%outputtab%add_term(cloc2)
        call this%outputtab%add_term(this%budobj%budterm(i)%id2(n))
        inum = inum + 1
      end do
    end do

mvt_rp()

subroutine tspmvtmodule::mvt_rp ( class(tspmvttype)  this )

private

Definition at line 207 of file tsp-mvt.f90.

    ! -- modules
    use tdismodule, only: kper, kstp
    ! -- dummy
    class(TspMvtType) :: this
    !
    ! -- At this point, the mvrbudobj is available to set up the mvt budobj
    if (kper * kstp == 1) then
      !
      ! -- If mvt is for a single model then point to fmi1
      !cdl todo: this needs to be called from GwtGwtExg somehow for the 2 model case
      if (associated(this%fmi1, this%fmi2)) then
        call this%set_pointer_mvrbudobj(this%fmi1%mvrbudobj)
      end if
      !
      ! -- Set up the mvt budobject
      call this%mvt_scan_mvrbudobj()
      call this%mvt_setup_budobj()
      !
      ! -- Define the budget object to be the size of maxpackages
      call this%budget%budget_df(this%maxpackages, 'TRANSPORT MOVER', bddim='M')
      call this%budget%set_ibudcsv(this%ibudcsv)
    end if

mvt_scan_mvrbudobj()

subroutine tspmvtmodule::mvt_scan_mvrbudobj ( class(tspmvttype)  this )

private

Scan through the gwf water mover budget object and determine the maximum number of packages and unique package names

Definition at line 775 of file tsp-mvt.f90.

    class(TspMvtType) :: this
    integer(I4B) :: nbudterm
    integer(I4B) :: maxpackages
    integer(I4B) :: i, j
    integer(I4B) :: ipos
    logical :: found
    !
    ! -- Calculate maxpackages, which is the the square of nbudterm
    nbudterm = this%mvrbudobj%nbudterm
    do i = 1, nbudterm
      if (i * i == nbudterm) then
        maxpackages = i
        exit
      end if
    end do
    this%maxpackages = maxpackages
    !
    ! -- Allocate paknames
    if (associated(this%paknames)) deallocate (this%paknames)
    allocate (this%paknames(this%maxpackages))
    do i = 1, this%maxpackages
      this%paknames(i) = ''
    end do
    !
    ! -- Scan through mvrbudobj and create unique paknames
    ipos = 1
    do i = 1, nbudterm
      found = .false.
      do j = 1, ipos
        if (this%mvrbudobj%budterm(i)%text2id1 == this%paknames(j)) then
          found = .true.
          exit
        end if
      end do
      if (.not. found) then
        this%paknames(ipos) = this%mvrbudobj%budterm(i)%text2id1
        ipos = ipos + 1
      end if
    end do

mvt_setup_budobj()

subroutine tspmvtmodule::mvt_setup_budobj

(

class(tspmvttype)

this

)

Definition at line 670 of file tsp-mvt.f90.

    ! -- modules
    use constantsmodule, only: lenbudtxt
    ! -- dummy
    class(TspMvtType) :: this
    ! -- local
    integer(I4B) :: nbudterm
    integer(I4B) :: ncv
    integer(I4B) :: maxlist
    integer(I4B) :: i
    integer(I4B) :: naux
    character(len=LENMODELNAME) :: modelname1, modelname2
    character(len=LENPACKAGENAME) :: packagename1, packagename2
    character(len=LENBUDTXT) :: text
    !
    ! -- Assign terms to set up the mover budget object
    nbudterm = this%mvrbudobj%nbudterm
    ncv = 0
    naux = 0
    if (this%depvartype == 'CONCENTRATION') then
      text = '        MVT-FLOW'
    else
      text = '        MVE-FLOW'
    end if
    !
    ! -- Set up budobj
    call this%budobj%budgetobject_df(ncv, nbudterm, 0, 0, bddim_opt='M')
    !
    ! -- Allocate mvrterm
    allocate (this%mvrterm(nbudterm))
    !
    ! -- Go through the water mover budget terms and set up the transport
    !    mover budget terms
    do i = 1, nbudterm
      modelname1 = this%mvrbudobj%budterm(i)%text1id1
      packagename1 = this%mvrbudobj%budterm(i)%text2id1
      modelname2 = this%mvrbudobj%budterm(i)%text1id2
      packagename2 = this%mvrbudobj%budterm(i)%text2id2
      maxlist = this%mvrbudobj%budterm(i)%maxlist
      call this%budobj%budterm(i)%initialize(text, &
                                             modelname1, &
                                             packagename1, &
                                             modelname2, &
                                             packagename2, &
                                             maxlist, .false., .false., &
                                             naux)
      this%mvrterm(i)%prov_is_m1 = (this%gwfmodelname1 == modelname1)
      allocate (this%mvrterm(i)%qty(maxlist))
      this%mvrterm(i)%qty = dnodata
    end do

mvt_setup_outputtab()

subroutine tspmvtmodule::mvt_setup_outputtab ( class(tspmvttype), intent(inout)  this )

private

Definition at line 819 of file tsp-mvt.f90.

    ! -- dummy
    class(TspMvtType), intent(inout) :: this
    ! -- local
    character(len=LINELENGTH) :: title
    character(len=LINELENGTH) :: text
    integer(I4B) :: ntabcol
    integer(I4B) :: maxrow
    integer(I4B) :: ilen
    !
    ! -- Allocate and initialize the output table
    if (this%iprflow /= 0) then
      !
      ! -- Dimension table
      ntabcol = 7
      maxrow = 0
      !
      ! -- Initialize the output table object
      title = 'TRANSPORT MOVER PACKAGE ('//trim(this%packName)// &
              ') FLOW RATES'
      call table_cr(this%outputtab, this%packName, title)
      call this%outputtab%table_df(maxrow, ntabcol, this%iout, &
                                   transient=.true.)
      text = 'NUMBER'
      call this%outputtab%initialize_column(text, 10, alignment=tabcenter)
      text = 'PROVIDER LOCATION'
      ilen = lenmodelname + lenpackagename + 1
      call this%outputtab%initialize_column(text, ilen)
      text = 'PROVIDER ID'
      call this%outputtab%initialize_column(text, 10)
      text = 'PROVIDER FLOW RATE'
      call this%outputtab%initialize_column(text, 10)
      text = 'PROVIDER TRANSPORT RATE'
      call this%outputtab%initialize_column(text, 10)
      text = 'RECEIVER LOCATION'
      ilen = lenmodelname + lenpackagename + 1
      call this%outputtab%initialize_column(text, ilen)
      text = 'RECEIVER ID'
      call this%outputtab%initialize_column(text, 10)
      !
    end if

Here is the call graph for this function:

mvt_update_qmfrommvr()

subroutine tspmvtmodule::mvt_update_qmfrommvr ( class(tspmvttype)  this )

private

Definition at line 308 of file tsp-mvt.f90.

    class(TspMvtType) :: this
    ! local
    integer(I4B) :: i, n
    integer(I4B) :: id2, nlist
    integer(I4B) :: irc
    integer(I4B) :: nbudterm
    type(TspFmiType), pointer :: fmi_rc !< pointer to receiver model fmi package
 
    nbudterm = this%mvrbudobj%nbudterm
    do i = 1, nbudterm
      nlist = this%mvrbudobj%budterm(i)%nlist
      if (nlist == 0) cycle
 
      call this%set_fmi(i, fmi_rc, .false.)
      if (.not. associated(fmi_rc)) cycle
 
      call fmi_rc%get_package_index(this%mvrbudobj%budterm(i)%text2id2, irc)
 
      do n = 1, nlist
 
        ! lak/sfr/maw/uzf id2 (receiver)
        id2 = this%mvrbudobj%budterm(i)%id2(n)
 
        ! add the mover rate times the provider concentration into the receiver;
        ! accumulated since multiple providers can move water into the same receiver
        if (fmi_rc%iatp(irc) /= 0) then
          fmi_rc%datp(irc)%qmfrommvr(id2) = fmi_rc%datp(irc)%qmfrommvr(id2) - &
                                            this%mvrbudobj%budterm(i)%flow(n) * &
                                            this%mvrterm(i)%qty(n) * &
                                            this%eqnsclfac
        end if
      end do
    end do

read_options()

subroutine tspmvtmodule::read_options

(

class(tspmvttype)

this

)

Definition at line 591 of file tsp-mvt.f90.

    ! -- modules
    use openspecmodule, only: access, form
    use inputoutputmodule, only: assign_iounit, openfile
    ! -- dummy
    class(TspMvtType) :: this
    ! -- local
    character(len=LINELENGTH) :: errmsg, keyword
    character(len=MAXCHARLEN) :: fname
    integer(I4B) :: ierr
    logical :: isfound, endOfBlock
    ! -- formats
    character(len=*), parameter :: fmtflow = &
      "(4x, a, 1x, a, 1x, ' WILL BE SAVED TO FILE: ', a, &
      &/4x, 'OPENED ON UNIT: ', I0)"
    character(len=*), parameter :: fmtflow2 = &
      &"(4x, 'FLOWS WILL BE SAVED TO BUDGET FILE')"
    !
    ! -- Get options block
    call this%parser%GetBlock('OPTIONS', isfound, ierr, blockrequired=.false., &
                              supportopenclose=.true.)
    !
    ! -- Parse options block if detected
    if (isfound) then
      write (this%iout, '(1x,a)') 'PROCESSING MVT OPTIONS'
      do
        call this%parser%GetNextLine(endofblock)
        if (endofblock) exit
        call this%parser%GetStringCaps(keyword)
        select case (keyword)
        case ('SAVE_FLOWS')
          this%ipakcb = -1
          write (this%iout, fmtflow2)
        case ('PRINT_INPUT')
          this%iprpak = 1
          write (this%iout, '(4x,a)') 'MVT INPUT WILL BE PRINTED.'
        case ('PRINT_FLOWS')
          this%iprflow = 1
          write (this%iout, '(4x,a)') &
            'MVT FLOWS WILL BE PRINTED TO LISTING FILE.'
        case ('BUDGET')
          call this%parser%GetStringCaps(keyword)
          if (keyword == 'FILEOUT') then
            call this%parser%GetString(fname)
            call assign_iounit(this%ibudgetout, this%inunit, "BUDGET fileout")
            call openfile(this%ibudgetout, this%iout, fname, 'DATA(BINARY)', &
                          form, access, 'REPLACE')
            write (this%iout, fmtflow) 'MVT', 'BUDGET', trim(adjustl(fname)), &
              this%ibudgetout
          else
            call store_error('Optional BUDGET keyword must &
                             &be followed by FILEOUT')
          end if
        case ('BUDGETCSV')
          call this%parser%GetStringCaps(keyword)
          if (keyword == 'FILEOUT') then
            call this%parser%GetString(fname)
            call assign_iounit(this%ibudcsv, this%inunit, "BUDGETCSV fileout")
            call openfile(this%ibudcsv, this%iout, fname, 'CSV', &
                          filstat_opt='REPLACE')
            write (this%iout, fmtflow) 'MVT', 'BUDGET CSV', &
              trim(adjustl(fname)), this%ibudcsv
          else
            call store_error('Optional BUDGETCSV keyword must be followed by &
              &FILEOUT')
          end if
        case default
          write (errmsg, '(a,a)') 'Unknown MVT option: ', &
            trim(keyword)
          call store_error(errmsg)
          call this%parser%StoreErrorUnit()
        end select
      end do
      write (this%iout, '(1x,a)') 'END OF MVT OPTIONS'
    end if

Here is the call graph for this function:

set_fmi()

subroutine tspmvtmodule::set_fmi ( class(tspmvttype)  this, integer(i4b), intent(in)  ibudterm, type(tspfmitype), pointer  fmi, logical(lgp)  set_provider  )

private

Parameters

set_provider

set provider fmi when true, receiver fmi otherwise

Definition at line 344 of file tsp-mvt.f90.

    class(TspMvtType) :: this
    integer(I4B), intent(in) :: ibudterm
    type(TspFmiType), pointer :: fmi
    logical(LGP) :: set_provider !< set provider fmi when true, receiver fmi otherwise
    ! local
    character(len=LENMODELNAME) :: model_name
 
    fmi => null()
    if (this%gwfmodelname1 == '' .and. this%gwfmodelname2 == '') then
      fmi => this%fmi1 ! fmi2 equals fmi1 here
      return
    end if
 
    ! model name for provider is text1id1, receiver is text1id2
    if (set_provider) then
      model_name = this%mvrbudobj%budterm(ibudterm)%text1id1
    else
      model_name = this%mvrbudobj%budterm(ibudterm)%text1id2
    end if
 
    if (model_name == this%gwfmodelname1) then
      ! model 1 is the provider/receiver
      fmi => this%fmi1
    else if (model_name == this%gwfmodelname2) then
      ! model 2 is the provider
      fmi => this%fmi2
    end if
 

set_pointer_mvrbudobj()

subroutine tspmvtmodule::set_pointer_mvrbudobj ( class(tspmvttype)  this, type(budgetobjecttype), intent(in), target  mvrbudobj  )

private

Store a pointer to mvrbudobj, which contains the simulated water mover flows from either a gwf model MVR package or from a gwf-gwf exchange MVR package.

Definition at line 189 of file tsp-mvt.f90.

    class(TspMvtType) :: this
    type(BudgetObjectType), intent(in), target :: mvrbudobj
    this%mvrbudobj => mvrbudobj