This module contains the TspSsm Module. More...

Data Types
type	tspssmtype
	Derived type for the SSM Package. More...

Functions/Subroutines
subroutine, public	ssm_cr (ssmobj, name_model, input_mempath, inunit, iout, fmi, eqnsclfac, depvartype)
	@ brief Create a new SSM package More...

subroutine	ssm_df (this)
	@ brief Define SSM Package More...

subroutine	ssm_ar (this, dis, ibound, cnew)
	@ brief Allocate and read SSM Package More...

subroutine	ssm_rp (this)
	@ brief Read and prepare this SSM Package More...

subroutine	ssm_ad (this)
	@ brief Advance the SSM Package More...

subroutine	ssm_term (this, ipackage, ientry, rrate, rhsval, hcofval, cssm, qssm)
	@ brief Calculate the SSM mass flow rate and hcof and rhs values More...

subroutine	get_ssm_conc (this, ipackage, ientry, nbound_flow, conc, lauxmixed)
	@ brief Provide bound concentration (or temperature) and mixed flag More...

subroutine	ssm_fc (this, matrix_sln, idxglo, rhs)
	@ brief Fill coefficients More...

subroutine	ssm_cq (this, flowja)
	@ brief Calculate flow More...

subroutine	ssm_bd (this, isuppress_output, model_budget)
	@ brief Calculate the global SSM budget terms More...

subroutine	ssm_ot_flow (this, icbcfl, ibudfl, icbcun)
	@ brief Output flows More...

subroutine	ssm_da (this)
	@ brief Deallocate More...

subroutine	allocate_scalars (this)
	@ brief Allocate scalars More...

subroutine	allocate_arrays (this)
	@ brief Allocate arrays More...

subroutine	source_options (this)
	Source input options. More...

subroutine	source_sources (this)
	Source sources input block. More...

subroutine	source_fileinput (this)
	Source fileinput input block. More...

subroutine	set_iauxpak (this, ip, packname, auxname)
	@ brief Set iauxpak array value for package ip More...

subroutine	set_ssmivec (this, ip, packname, spc_mempath, input_fname)
	Set ssmivec array value for package ip. More...

subroutine	pak_setup_outputtab (this)
	@ brief Setup the output table More...

Variables
character(len=lenftype)	ftype = 'SSM'

character(len=lenpackagename)	text = ' SOURCE-SINK MIX'

Detailed Description

This module contains the code for handling sources and sinks associated with groundwater flow model stress packages.

todo: need observations for SSM terms

Function/Subroutine Documentation

◆ allocate_arrays()

subroutine tspssmmodule::allocate_arrays ( class(tspssmtype) this )

Allocate array variables for this derived type

Parameters

this	TspSsmType object

Definition at line 711 of file tsp-ssm.f90.

     ! -- modules
     use memorymanagermodule, only: mem_allocate, mem_setptr
     ! -- dummy
     class(TspSsmType) :: this !< TspSsmType object
     ! -- local
     integer(I4B) :: nflowpack
     integer(I4B) :: i
     !
     ! -- Allocate
     nflowpack = this%fmi%nflowpack
     call mem_allocate(this%iauxpak, nflowpack, 'IAUXPAK', this%memoryPath)
     call mem_allocate(this%isrctype, nflowpack, 'ISRCTYPE', this%memoryPath)
     !
     ! -- Initialize
     do i = 1, nflowpack
       this%iauxpak(i) = 0
       this%isrctype(i) = 0
     end do
     !
     ! -- Allocate the ssmivec array
     allocate (this%ssmivec(nflowpack))

◆ allocate_scalars()

subroutine tspssmmodule::allocate_scalars ( class(tspssmtype) this )

Allocate scalar variables for this derived type

Parameters

this	TspSsmType object

Definition at line 691 of file tsp-ssm.f90.

     ! -- modules
     use memorymanagermodule, only: mem_allocate, mem_setptr
     ! -- dummy
     class(TspSsmType) :: this !< TspSsmType object
     !
     ! -- allocate scalars in NumericalPackageType
     call this%NumericalPackageType%allocate_scalars()
     !
     ! -- Allocate
     call mem_allocate(this%nbound, 'NBOUND', this%memoryPath)
     !
     ! -- Initialize
     this%nbound = 0

◆ get_ssm_conc()

subroutine tspssmmodule::get_ssm_conc	(	class(tspssmtype)	this,
		integer(i4b), intent(in)	ipackage,
		integer(i4b), intent(in)	ientry,
		integer(i4b), intent(in)	nbound_flow,
		real(dp), intent(out)	conc,
		logical(lgp), intent(out)	lauxmixed
	)

SSM concentrations and temperatures can be provided in auxiliary variables or through separate SPC files. If not provided, the default concentration (or temperature) is zero. This single routine provides the SSM bound concentration (or temperature) based on these different approaches. The mixed flag indicates whether or not the boundary as a mixed type.

Parameters

	this	TspSsmType
[in]	ipackage	package number
[in]	ientry	bound number
[in]	nbound_flow	size of flow package bound list
[out]	conc	user-specified concentration/temperature for this bound
[out]	lauxmixed	user-specified flag for marking this as a mixed boundary

Definition at line 355 of file tsp-ssm.f90.

     ! -- dummy
     class(TspSsmType) :: this !< TspSsmType
     integer(I4B), intent(in) :: ipackage !< package number
     integer(I4B), intent(in) :: ientry !< bound number
     integer(I4B), intent(in) :: nbound_flow !< size of flow package bound list
     real(DP), intent(out) :: conc !< user-specified concentration/temperature for this bound
     logical(LGP), intent(out) :: lauxmixed !< user-specified flag for marking this as a mixed boundary
     ! -- local
     integer(I4B) :: isrctype
     integer(I4B) :: iauxpos
     !
     conc = dzero
     lauxmixed = .false.
     isrctype = this%isrctype(ipackage)
     !
     select case (isrctype)
     case (1, 2)
       iauxpos = this%iauxpak(ipackage)
       conc = this%fmi%gwfpackages(ipackage)%auxvar(iauxpos, ientry)
       if (isrctype == 2) lauxmixed = .true.
     case (3, 4)
       conc = this%ssmivec(ipackage)%get_value(ientry, nbound_flow)
       if (isrctype == 4) lauxmixed = .true.
     end select

◆ pak_setup_outputtab()

subroutine tspssmmodule::pak_setup_outputtab ( class(tspssmtype), intent(inout) this )

Setup the output table by creating the column headers.

Definition at line 1017 of file tsp-ssm.f90.

     ! -- dummy
     class(TspSsmType), intent(inout) :: this
     ! -- local
     character(len=LINELENGTH) :: title
     character(len=LINELENGTH) :: text
     integer(I4B) :: ntabcol
     !
     ! -- allocate and initialize the output table
     if (this%iprflow /= 0) then
       !
       ! -- dimension table
       ntabcol = 6
       !if (this%inamedbound > 0) then
       !  ntabcol = ntabcol + 1
       !end if
       !
       ! -- initialize the output table object
       title = 'SSM PACKAGE ('//trim(this%packName)// &
               ') FLOW RATES'
       call table_cr(this%outputtab, this%packName, title)
       call this%outputtab%table_df(1, ntabcol, this%iout, transient=.true.)
       text = 'NUMBER'
       call this%outputtab%initialize_column(text, 10, alignment=tabcenter)
       text = 'CELLID'
       call this%outputtab%initialize_column(text, 20, alignment=tableft)
       text = 'BOUND Q'
       call this%outputtab%initialize_column(text, 15, alignment=tabcenter)
       text = 'SSM CONC'
       call this%outputtab%initialize_column(text, 15, alignment=tabcenter)
       text = 'RATE'
       call this%outputtab%initialize_column(text, 15, alignment=tabcenter)
       text = 'PACKAGE NAME'
       call this%outputtab%initialize_column(text, 16, alignment=tabcenter)
       !if (this%inamedbound > 0) then
       !  text = 'NAME'
       !  call this%outputtab%initialize_column(text, 20, alignment=TABLEFT)
       !end if
     end if

Here is the call graph for this function:

◆ set_iauxpak()

subroutine tspssmmodule::set_iauxpak	(	class(tspssmtype), intent(inout)	this,
		integer(i4b), intent(in)	ip,
		character(len=*), intent(in)	packname,
		character(len=*), intent(in)	auxname
	)

The routine searches through the auxiliary names in package ip and sets iauxpak to the column number corresponding to the correct auxiliary column.

Parameters

[in,out]	this	TspSsmType
[in]	ip	package number
[in]	packname	name of package
[in]	auxname	name of aux

Definition at line 960 of file tsp-ssm.f90.

     ! -- dummy
     class(TspSsmType), intent(inout) :: this !< TspSsmType
     integer(I4B), intent(in) :: ip !< package number
     character(len=*), intent(in) :: packname !< name of package
     character(len=*), intent(in) :: auxname !< name of aux
     ! -- local
     logical :: auxfound
     integer(I4B) :: iaux
     !
     ! -- match package auxiliary
     auxfound = .false.
     do iaux = 1, this%fmi%gwfpackages(ip)%naux
       if (trim(this%fmi%gwfpackages(ip)%auxname(iaux)) == &
           trim(auxname)) then
         auxfound = .true.
         exit
       end if
     end do
     if (.not. auxfound) then
       write (errmsg, '(a, a)') &
         'Auxiliary name cannot be found: ', trim(auxname)
       call store_error(errmsg)
       call store_error_filename(this%input_fname)
     end if
     !
     ! -- set iauxpak and write message
     this%iauxpak(ip) = iaux
     write (this%iout, '(4x, a, i0, a, a)') 'USING AUX COLUMN ', &
       iaux, ' IN PACKAGE ', trim(packname)

Here is the call graph for this function:

◆ set_ssmivec()

subroutine tspssmmodule::set_ssmivec	(	class(tspssmtype), intent(inout)	this,
		integer(i4b), intent(in)	ip,
		character(len=*), intent(in)	packname,
		character(len=*), intent(in)	spc_mempath,
		character(len=*), intent(in)	input_fname
	)

Initialize TspSpcType from the input mempath.

Parameters

[in,out]	this	TspSsmType
[in]	ip	package number
[in]	packname	name of corresponding flow package
[in]	spc_mempath	input mempath for this SPC instance
[in]	input_fname	SPC input file name (for error messages)

Definition at line 996 of file tsp-ssm.f90.

     ! -- dummy
     class(TspSsmType), intent(inout) :: this !< TspSsmType
     integer(I4B), intent(in) :: ip !< package number
     character(len=*), intent(in) :: packname !< name of corresponding flow package
     character(len=*), intent(in) :: spc_mempath !< input mempath for this SPC instance
     character(len=*), intent(in) :: input_fname !< SPC input file name (for error messages)
     ! -- local
     type(TspSpcType), pointer :: ssmiptr
     !
     ! -- initialize the TspSpcType reader for this package
     ssmiptr => this%ssmivec(ip)
     call ssmiptr%initialize(this%dis, ip, spc_mempath, this%iout, &
                             this%name_model, trim(packname), &
                             this%depvartype, input_fname)

◆ source_fileinput()

subroutine tspssmmodule::source_fileinput ( class(tspssmtype) this )

Initialize an SPC input file reader for each input entry.

Definition at line 851 of file tsp-ssm.f90.

     ! -- modules
     use memorymanagermodule, only: mem_setptr, get_isize
     use characterstringmodule, only: characterstringtype
     ! -- dummy
     class(TspSsmType) :: this
     ! -- locals
     type(CharacterStringType), dimension(:), pointer, &
       contiguous :: pnames, ftypes, iotypes, fnames, conditions, spc6_mempaths
     character(len=LINELENGTH) :: pname, ftype, iotype, fname, condition
     character(len=LENMEMPATH) :: spc_mempath
     integer(I4B) :: n, ip, isize
     logical(LGP) :: found
     ! -- formats
  
     call get_isize('PNAME', this%input_mempath, isize)
     if (isize <= 0) then
       write (this%iout, '(/1x,a)') &
         'OPTIONAL SSM FILEINPUT BLOCK NOT FOUND OR EMPTY.'
       return
     end if
  
     ! set pointers to input context
     call mem_setptr(pnames, 'PNAME', this%input_mempath)
     call mem_setptr(ftypes, 'SPC6', this%input_mempath)
     call mem_setptr(iotypes, 'FILEIN', this%input_mempath)
     call mem_setptr(fnames, 'SPC6_FILENAME', this%input_mempath)
     call mem_setptr(conditions, 'MIXED', this%input_mempath)
     call mem_setptr(spc6_mempaths, 'SPC6_MEMPATH', this%input_mempath)
  
     write (this%iout, '(/1x,a)') 'PROCESSING SSM FILEINPUT'
     do n = 1, size(pnames)
  
       pname = pnames(n)
       ftype = ftypes(n)
       iotype = iotypes(n)
       fname = fnames(n)
       condition = conditions(n)
       spc_mempath = spc6_mempaths(n)
       found = .false.
  
       do ip = 1, this%fmi%nflowpack
         if (trim(adjustl(this%fmi%gwfpackages(ip)%name)) == pname) then
           found = .true.
           exit
         end if
       end do
  
       if (.not. found) then
         write (errmsg, '(a,a)') 'Flow package cannot be found: ', &
           trim(pname)
         call store_error(errmsg)
         cycle
       end if
  
       ! -- Ensure package was not specified more than once in SOURCES block
       if (this%isrctype(ip) /= 0) then
         write (errmsg, '(a, a)') &
           'A package cannot be specified more than once in the SSM SOURCES &
           &block.  The following package was specified more than once: ', &
           trim(pname)
         call store_error(errmsg)
         cycle
       end if
  
       ! verify ftype
       if (ftype == 'SPC6') then
         write (this%iout, '(4x,a)') 'SPC6 SOURCE DETECTED:'
       else
         write (errmsg, '(a,a)') &
           'SRCTYPE must be SPC6. Found: ', trim(ftype)
         call store_error(errmsg)
         cycle
       end if
  
       ! verify iotype
       if (iotype /= 'FILEIN') then
         errmsg = 'SPC6 keyword must be followed by "FILEIN" '// &
                  'then by filename and optionally by <MIXED>.'
         call store_error(errmsg)
         cycle
       end if
  
       ! -- Use set_ssmivec to initialize the ssmi object from the
       !    input mempath
       call this%set_ssmivec(ip, pname, spc_mempath, trim(fname))
  
       if (condition == 'MIXED') then
         this%isrctype(ip) = 4
         write (this%iout, '(4x,a,a)') 'ASSIGNED MIXED SSM TYPE TO PACKAGE ', &
           trim(pname)
       else
         this%isrctype(ip) = 3
       end if
     end do
     write (this%iout, '(1x,a)') 'END PROCESSING SSM FILEINPUT'
  
     ! -- terminate if errors
     if (count_errors() > 0) then
       call store_error_filename(this%input_fname)
     end if

Here is the call graph for this function:

◆ source_options()

subroutine tspssmmodule::source_options ( class(tspssmtype) this )

Definition at line 737 of file tsp-ssm.f90.

     ! -- modules
     !use KindModule, only: LGP
     use memorymanagerextmodule, only: mem_set_value
     use gwtssminputmodule, only: gwtssmparamfoundtype
     ! -- dummy
     class(TspSsmType) :: this
     ! -- locals
     type(GwtSsmParamFoundType) :: found
     ! -- formats
     character(len=*), parameter :: fmtiprflow = &
       "(4x,'SSM FLOW INFORMATION WILL BE PRINTED TO LISTING FILE &
       &WHENEVER ICBCFL IS NOT ZERO.')"
     character(len=*), parameter :: fmtisvflow = &
       "(4x,'CELL-BY-CELL FLOW INFORMATION WILL BE SAVED TO BINARY FILE &
       &WHENEVER ICBCFL IS NOT ZERO.')"
     !
     ! -- update defaults with idm sourced values
     call mem_set_value(this%iprflow, 'PRINT_FLOWS', this%input_mempath, &
                        found%print_flows)
     call mem_set_value(this%ipakcb, 'SAVE_FLOWS', this%input_mempath, &
                        found%save_flows)
  
     if (found%save_flows) this%ipakcb = -1
  
     ! log options
     write (this%iout, '(1x,a)') 'PROCESSING SSM OPTIONS'
     if (found%print_flows) write (this%iout, fmtiprflow)
     if (found%save_flows) write (this%iout, fmtisvflow)
     write (this%iout, '(1x,a)') 'END OF SSM OPTIONS'

◆ source_sources()

subroutine tspssmmodule::source_sources ( class(tspssmtype) this )

Definition at line 771 of file tsp-ssm.f90.

     ! -- modules
     !use KindModule, only: LGP
     use memorymanagermodule, only: mem_setptr
     use characterstringmodule, only: characterstringtype
     ! -- dummy
     class(TspSsmType) :: this
     ! -- locals
     type(CharacterStringType), dimension(:), pointer, &
       contiguous :: pnames, srctypes, auxnames
     character(len=LINELENGTH) :: pname, srctype, auxname
     integer(I4B) :: n, ip
     logical(LGP) :: found
     ! -- formats
  
     ! set pointers to input context
     call mem_setptr(pnames, 'PNAME_SOURCES', this%input_mempath)
     call mem_setptr(srctypes, 'SRCTYPE', this%input_mempath)
     call mem_setptr(auxnames, 'AUXNAME', this%input_mempath)
  
     write (this%iout, '(/1x,a)') 'PROCESSING SSM SOURCES'
     do n = 1, size(pnames)
  
       pname = pnames(n)
       srctype = srctypes(n)
       auxname = auxnames(n)
       found = .false.
  
       do ip = 1, this%fmi%nflowpack
         if (trim(adjustl(this%fmi%gwfpackages(ip)%name)) == pname) then
           found = .true.
           exit
         end if
       end do
  
       if (.not. found) then
         write (errmsg, '(a,a)') 'Flow package cannot be found: ', &
           trim(pname)
         call store_error(errmsg)
         cycle
       end if
  
       ! -- Ensure package was not specified more than once in SOURCES block
       if (this%isrctype(ip) /= 0) then
         write (errmsg, '(a, a)') &
           'A package cannot be specified more than once in the SSM SOURCES &
           &block.  The following package was specified more than once: ', &
           trim(pname)
         call store_error(errmsg)
         cycle
       end if
  
       if (srctype == 'AUX') then
         this%isrctype(ip) = 1
         write (this%iout, '(4x,a)') 'AUX SOURCE DETECTED.'
       else if (srctype == 'AUXMIXED') then
         this%isrctype(ip) = 2
         write (this%iout, '(4x,a)') 'AUXMIXED SOURCE DETECTED.'
       else
         write (errmsg, '(a, a)') &
           'SRCTYPE must be AUX or AUXMIXED.  Found: ', trim(srctype)
         call store_error(errmsg)
         cycle
       end if
  
       ! -- Find and store the auxiliary column
       call this%set_iauxpak(ip, srctype, auxname)
     end do
     write (this%iout, '(1x,a)') 'END PROCESSING SSM SOURCES'
  
     ! -- terminate if errors
     if (count_errors() > 0) then
       call store_error_filename(this%input_fname)
     end if

Here is the call graph for this function:

◆ ssm_ad()

subroutine tspssmmodule::ssm_ad ( class(tspssmtype) this )

This routine is called from gwt_ad(). It is called at the beginning of each time step. The total number of flow boundaries is counted and stored in thisnbound. Also, if any SPC input files are used to provide source and sink concentrations (or temperatures) and time series are referenced in those files, then ssm concenrations must be interpolated for the time step.

Parameters

this	TspSsmType object

Definition at line 213 of file tsp-ssm.f90.

     ! -- dummy
     class(TspSsmType) :: this !< TspSsmType object
     ! -- local
     integer(I4B) :: ip
     type(TspSpcType), pointer :: ssmiptr
     integer(I4B) :: i
     integer(I4B) :: node
     !
     ! -- Calculate total number of existing flow boundaries. It is possible
     !    that a node may equal zero.  In this case, the bound should be
     !    skipped and not written to ssm output.
     this%nbound = 0
     do ip = 1, this%fmi%nflowpack
       if (this%fmi%iatp(ip) /= 0) cycle
       do i = 1, this%fmi%gwfpackages(ip)%nbound
         node = this%fmi%gwfpackages(ip)%nodelist(i)
         if (node > 0) then
           this%nbound = this%nbound + 1
         end if
       end do
     end do
     !
     ! -- Call the ad method on any ssm input files so that values for
     !    time-series are interpolated
     do ip = 1, this%fmi%nflowpack
       if (this%fmi%iatp(ip) /= 0) cycle
       if (this%isrctype(ip) == 3 .or. this%isrctype(ip) == 4) then
         ssmiptr => this%ssmivec(ip)
         call ssmiptr%spc_ad(this%fmi%gwfpackages(ip)%nbound, &
                             this%fmi%gwfpackages(ip)%budtxt)
       end if
     end do

◆ ssm_ar()

subroutine tspssmmodule::ssm_ar	(	class(tspssmtype)	this,
		class(disbasetype), intent(in), pointer	dis,
		integer(i4b), dimension(:), pointer, contiguous	ibound,
		real(dp), dimension(:), pointer, contiguous	cnew
	)

This routine is called from gwt_ar(). It allocates arrays, reads options and data, and sets up the output table.

Parameters

	this	TspSsmType object
[in]	dis	discretization package
	ibound	GWT model ibound
	cnew	GWT model dependent variable

Definition at line 130 of file tsp-ssm.f90.

     ! -- modules
     use memorymanagermodule, only: mem_setptr
     ! -- dummy
     class(TspSsmType) :: this !< TspSsmType object
     class(DisBaseType), pointer, intent(in) :: dis !< discretization package
     integer(I4B), dimension(:), pointer, contiguous :: ibound !< GWT model ibound
     real(DP), dimension(:), pointer, contiguous :: cnew !< GWT model dependent variable
     ! -- formats
     character(len=*), parameter :: fmtssm = &
       "(1x,/1x,'SSM -- SOURCE-SINK MIXING PACKAGE, VERSION 1, 8/25/2017', &
       &' INPUT READ FROM MEMPATH: ', A, //)"
     !
     ! -- print a message identifying the storage package.
     write (this%iout, fmtssm) this%input_mempath
     !
     ! -- store pointers to arguments that were passed in
     this%dis => dis
     this%ibound => ibound
     this%cnew => cnew
     !
     ! -- Check to make sure that there are flow packages
     if (this%fmi%nflowpack == 0) then
       write (errmsg, '(a)') 'SSM package does not detect any boundary flows &
                             &that require SSM terms.  Activate GWF-GWT (or &
                             &GWF-GWE, as appropriate) exchange or activate &
                             &FMI package and provide a budget file that &
                             &contains boundary flows.  If no boundary flows &
                             &are present in corresponding GWF model then this &
                             &SSM package should be removed.'
       call store_error(errmsg)
       call store_error_filename(this%input_fname)
     end if
     !
     ! -- Allocate arrays
     call this%allocate_arrays()
     !
     ! -- read ssm options
     call this%source_options()
     !
     ! -- read data blocks
     call this%source_sources()
     call this%source_fileinput()
     !
     ! -- setup the output table
     call this%pak_setup_outputtab()

Here is the call graph for this function:

◆ ssm_bd()

subroutine tspssmmodule::ssm_bd	(	class(tspssmtype)	this,
		integer(i4b), intent(in)	isuppress_output,
		type(budgettype), intent(inout)	model_budget
	)

Calculate the global SSM budget terms using separate in and out entries for each flow package.

Parameters

	this	TspSsmType object
[in]	isuppress_output	flag to suppress output
[in,out]	model_budget	budget object for the GWT model

Definition at line 465 of file tsp-ssm.f90.

     ! -- modules
     use tdismodule, only: delt
     use budgetmodule, only: budgettype
     ! -- dummy
     class(TspSsmType) :: this !< TspSsmType object
     integer(I4B), intent(in) :: isuppress_output !< flag to suppress output
     type(BudgetType), intent(inout) :: model_budget !< budget object for the GWT model
     ! -- local
     character(len=LENBUDROWLABEL) :: rowlabel
     integer(I4B) :: ip
     integer(I4B) :: i
     integer(I4B) :: n
     real(DP) :: rate
     real(DP) :: rin
     real(DP) :: rout
     !
     ! -- do for each flow package, unless it is being handled by an advanced
     !    transport package
     do ip = 1, this%fmi%nflowpack
       !
       ! -- cycle if package is being managed as an advanced package
       if (this%fmi%iatp(ip) /= 0) cycle
       !
       ! -- Initialize the rate accumulators
       rin = dzero
       rout = dzero
       !
       ! -- do for each boundary
       do i = 1, this%fmi%gwfpackages(ip)%nbound
         n = this%fmi%gwfpackages(ip)%nodelist(i)
         if (n <= 0) cycle
         call this%ssm_term(ip, i, rrate=rate)
         if (rate < dzero) then
           rout = rout - rate
         else
           rin = rin + rate
         end if
         !
       end do
       !
       rowlabel = 'SSM_'//adjustl(trim(this%fmi%flowpacknamearray(ip)))
       call model_budget%addentry(rin, rout, delt, &
                                  this%fmi%gwfpackages(ip)%budtxt, &
                                  isuppress_output, rowlabel=rowlabel)
     end do

◆ ssm_cq()

subroutine tspssmmodule::ssm_cq	(	class(tspssmtype)	this,
		real(dp), dimension(:), intent(inout), contiguous	flowja
	)

Calculate the resulting mass flow between the boundary and the connected GWT/GWE model cell. Update the diagonal position of the flowja array so that it ultimately contains the solute balance residual.

Parameters

	this	TspSsmType object
[in,out]	flowja	flow across each face in the model grid

Definition at line 430 of file tsp-ssm.f90.

     ! -- dummy
     class(TspSsmType) :: this !< TspSsmType object
     real(DP), dimension(:), contiguous, intent(inout) :: flowja !< flow across each face in the model grid
     ! -- local
     integer(I4B) :: ip
     integer(I4B) :: i
     integer(I4B) :: n
     integer(I4B) :: idiag
     real(DP) :: rate
     !
     ! -- do for each flow package
     do ip = 1, this%fmi%nflowpack
       !
       ! -- cycle if package is being managed as an advanced package
       if (this%fmi%iatp(ip) /= 0) cycle
       !
       ! -- do for each boundary
       do i = 1, this%fmi%gwfpackages(ip)%nbound
         n = this%fmi%gwfpackages(ip)%nodelist(i)
         if (n <= 0) cycle
         call this%ssm_term(ip, i, rrate=rate)
         idiag = this%dis%con%ia(n)
         flowja(idiag) = flowja(idiag) + rate
         !
       end do
       !
     end do

◆ ssm_cr()

subroutine, public tspssmmodule::ssm_cr	(	type(tspssmtype), pointer	ssmobj,
		character(len=*), intent(in)	name_model,
		character(len=*), intent(in)	input_mempath,
		integer(i4b), intent(in)	inunit,
		integer(i4b), intent(in)	iout,
		type(tspfmitype), intent(in), target	fmi,
		real(dp), intent(in), pointer	eqnsclfac,
		character(len=lenvarname), intent(in)	depvartype
	)

Create a new SSM package by defining names, allocating scalars

Parameters

	ssmobj	TspSsmType object
[in]	name_model	name of the model
[in]	input_mempath	input mempath
[in]	inunit	fortran unit for input
[in]	iout	fortran unit for output
[in]	fmi	Transport FMI package
[in]	eqnsclfac	governing equation scale factor
[in]	depvartype	dependent variable type ('concentration' or 'temperature')

Definition at line 80 of file tsp-ssm.f90.

     ! -- dummy
     type(TspSsmType), pointer :: ssmobj !< TspSsmType object
     character(len=*), intent(in) :: name_model !< name of the model
     character(len=*), intent(in) :: input_mempath !< input mempath
     integer(I4B), intent(in) :: inunit !< fortran unit for input
     integer(I4B), intent(in) :: iout !< fortran unit for output
     type(TspFmiType), intent(in), target :: fmi !< Transport FMI package
     real(DP), intent(in), pointer :: eqnsclfac !< governing equation scale factor
     character(len=LENVARNAME), intent(in) :: depvartype !< dependent variable type ('concentration' or 'temperature')
     !
     ! -- Create the object
     allocate (ssmobj)
     !
     ! -- create name and memory path
     call ssmobj%set_names(1, name_model, 'SSM', 'SSM', input_mempath)
     !
     ! -- Allocate scalars
     call ssmobj%allocate_scalars()
     !
     ! -- Set variables
     ssmobj%inunit = inunit
     ssmobj%iout = iout
     ssmobj%fmi => fmi
     ssmobj%eqnsclfac => eqnsclfac
     !
     ! -- Store pointer to labels associated with the current model so that the
     !    package has access to the corresponding dependent variable type
     ssmobj%depvartype = depvartype

Here is the caller graph for this function:

◆ ssm_da()

subroutine tspssmmodule::ssm_da ( class(tspssmtype) this )

Deallocate the memory associated with this derived type

Parameters

this	TspSsmType object

Definition at line 645 of file tsp-ssm.f90.

     ! -- modules
     use memorymanagermodule, only: mem_deallocate
     ! -- dummy
     class(TspSsmType) :: this !< TspSsmType object
     ! -- local
     integer(I4B) :: ip
     type(TspSpcType), pointer :: ssmiptr
     !
     ! -- Deallocate the ssmi objects if package was active
     if (this%inunit > 0) then
       do ip = 1, size(this%ssmivec)
         if (this%isrctype(ip) == 3 .or. this%isrctype(ip) == 4) then
           ssmiptr => this%ssmivec(ip)
           call ssmiptr%spc_da()
         end if
       end do
       deallocate (this%ssmivec)
     end if
     !
     ! -- Deallocate arrays if package was active
     if (this%inunit > 0) then
       call mem_deallocate(this%iauxpak)
       call mem_deallocate(this%isrctype)
       this%ibound => null()
       this%fmi => null()
     end if
     !
     ! -- output table object
     if (associated(this%outputtab)) then
       call this%outputtab%table_da()
       deallocate (this%outputtab)
       nullify (this%outputtab)
     end if
     !
     ! -- Scalars
     call mem_deallocate(this%nbound)
     !
     ! -- deallocate parent
     call this%NumericalPackageType%da()

◆ ssm_df()

subroutine tspssmmodule::ssm_df ( class(tspssmtype) this )

This routine is called from gwt_df(), but does not do anything because df is typically used to set up dimensions. For the ssm package, the total number of ssm entries is defined by the flow model.

Parameters

this	TspSsmType object

Definition at line 118 of file tsp-ssm.f90.

     ! -- modules
     use memorymanagermodule, only: mem_setptr
     ! -- dummy
     class(TspSsmType) :: this !< TspSsmType object

◆ ssm_fc()

subroutine tspssmmodule::ssm_fc	(	class(tspssmtype)	this,
		class(matrixbasetype), pointer	matrix_sln,
		integer(i4b), dimension(:), intent(in)	idxglo,
		real(dp), dimension(:), intent(inout)	rhs
	)

This routine adds the effects of the SSM to the matrix equations by updating the a matrix and right-hand side vector.

Definition at line 388 of file tsp-ssm.f90.

     ! -- dummy
     class(TspSsmType) :: this
     class(MatrixBaseType), pointer :: matrix_sln
     integer(I4B), intent(in), dimension(:) :: idxglo
     real(DP), intent(inout), dimension(:) :: rhs
     ! -- local
     integer(I4B) :: ip
     integer(I4B) :: i
     integer(I4B) :: n
     integer(I4B) :: idiag
     integer(I4B) :: nflowpack
     integer(I4B) :: nbound
     real(DP) :: hcofval
     real(DP) :: rhsval
     !
     ! -- do for each flow package
     nflowpack = this%fmi%nflowpack
     do ip = 1, nflowpack
       if (this%fmi%iatp(ip) /= 0) cycle
       !
       ! -- do for each entry in package (ip)
       nbound = this%fmi%gwfpackages(ip)%nbound
       do i = 1, nbound
         n = this%fmi%gwfpackages(ip)%nodelist(i)
         if (n <= 0) cycle
         call this%ssm_term(ip, i, rhsval=rhsval, hcofval=hcofval)
         idiag = idxglo(this%dis%con%ia(n))
         call matrix_sln%add_value_pos(idiag, hcofval)
         rhs(n) = rhs(n) + rhsval
         !
       end do
       !
     end do

◆ ssm_ot_flow()

subroutine tspssmmodule::ssm_ot_flow	(	class(tspssmtype)	this,
		integer(i4b), intent(in)	icbcfl,
		integer(i4b), intent(in)	ibudfl,
		integer(i4b), intent(in)	icbcun
	)

Based on user-specified controls, print SSM mass flow rates to the GWT listing file and/or write the SSM mass flow rates to the GWT binary budget file.

Parameters

	this	TspSsmType object
[in]	icbcfl	flag for writing binary budget terms
[in]	ibudfl	flag for printing budget terms to list file
[in]	icbcun	fortran unit number for binary budget file

Definition at line 519 of file tsp-ssm.f90.

     ! -- modules
     use tdismodule, only: kstp, kper
     use constantsmodule, only: lenpackagename, lenboundname, lenauxname, dzero
     ! -- dummy
     class(TspSsmType) :: this !< TspSsmType object
     integer(I4B), intent(in) :: icbcfl !< flag for writing binary budget terms
     integer(I4B), intent(in) :: ibudfl !< flag for printing budget terms to list file
     integer(I4B), intent(in) :: icbcun !< fortran unit number for binary budget file
     ! -- local
     character(len=LINELENGTH) :: title
     integer(I4B) :: node, nodeu
     character(len=20) :: nodestr
     integer(I4B) :: maxrows
     integer(I4B) :: ip
     integer(I4B) :: i, n2, ibinun
     real(DP) :: rrate
     real(DP) :: qssm
     real(DP) :: cssm
     integer(I4B) :: naux
     real(DP), dimension(0) :: auxrow
     character(len=LENAUXNAME), dimension(0) :: auxname
     ! -- for observations
     character(len=LENBOUNDNAME) :: bname
     ! -- formats
     character(len=*), parameter :: fmttkk = &
       &"(1X,/1X,A,'   PERIOD ',I0,'   STEP ',I0)"
     !
     ! -- initialize
     naux = 0
     maxrows = 0
     if (ibudfl /= 0 .and. this%iprflow /= 0) then
       call this%outputtab%set_kstpkper(kstp, kper)
       do ip = 1, this%fmi%nflowpack
         if (this%fmi%iatp(ip) /= 0) cycle
         !
         ! -- do for each boundary
         do i = 1, this%fmi%gwfpackages(ip)%nbound
           node = this%fmi%gwfpackages(ip)%nodelist(i)
           if (node > 0) then
             maxrows = maxrows + 1
           end if
         end do
       end do
       if (maxrows > 0) then
         call this%outputtab%set_maxbound(maxrows)
       end if
       title = 'SSM PACKAGE ('//trim(this%packName)// &
               ') FLOW RATES'
       call this%outputtab%set_title(title)
     end if
     !
     ! -- Set unit number for binary output
     if (this%ipakcb < 0) then
       ibinun = icbcun
     else if (this%ipakcb == 0) then
       ibinun = 0
     else
       ibinun = this%ipakcb
     end if
     if (icbcfl == 0) ibinun = 0
     !
     ! -- If cell-by-cell flows will be saved as a list, write header.
     if (ibinun /= 0) then
       call this%dis%record_srcdst_list_header(text, this%name_model, &
                                               this%name_model, this%name_model, &
                                               this%packName, naux, auxname, &
                                               ibinun, this%nbound, this%iout)
     end if
     !
     ! -- If no boundaries, skip flow calculations.
     if (this%nbound > 0) then
       !
       ! -- Loop through each boundary calculating flow.
       do ip = 1, this%fmi%nflowpack
         if (this%fmi%iatp(ip) /= 0) cycle
         !
         ! -- do for each boundary
         do i = 1, this%fmi%gwfpackages(ip)%nbound
           !
           ! -- Calculate rate for this entry
           node = this%fmi%gwfpackages(ip)%nodelist(i)
           if (node <= 0) cycle
           call this%ssm_term(ip, i, rrate=rrate, qssm=qssm, cssm=cssm)
           !
           ! -- Print the individual rates if the budget is being printed
           !    and PRINT_FLOWS was specified (this%iprflow<0)
           if (ibudfl /= 0) then
             if (this%iprflow /= 0) then
               !
               ! -- set nodestr and write outputtab table
               nodeu = this%dis%get_nodeuser(node)
               call this%dis%nodeu_to_string(nodeu, nodestr)
               bname = this%fmi%gwfpackages(ip)%name
               call this%outputtab%add_term(i)
               call this%outputtab%add_term(trim(adjustl(nodestr)))
               call this%outputtab%add_term(qssm)
               call this%outputtab%add_term(cssm)
               call this%outputtab%add_term(rrate)
               call this%outputtab%add_term(bname)
             end if
           end if
           !
           ! -- If saving cell-by-cell flows in list, write flow
           if (ibinun /= 0) then
             n2 = i
             call this%dis%record_mf6_list_entry(ibinun, node, n2, rrate, &
                                                 naux, auxrow, &
                                                 olconv2=.false.)
           end if
           !
         end do
         !
       end do
     end if
     if (ibudfl /= 0) then
       if (this%iprflow /= 0) then
         write (this%iout, '(1x)')
       end if
     end if

◆ ssm_rp()

subroutine tspssmmodule::ssm_rp ( class(tspssmtype) this )

This routine is called from gwt_rp(). It is called at the beginning of each stress period. If any SPC input files are used to provide source and sink concentrations (or temperatures), then period blocks for the current stress period are read.

Parameters

this	TspSsmType object

Definition at line 185 of file tsp-ssm.f90.

     ! -- dummy
     class(TspSsmType) :: this !< TspSsmType object
     ! -- local
     integer(I4B) :: ip
     type(TspSpcType), pointer :: ssmiptr
     ! -- formats
     !
     ! -- Call the rp method on any ssm input files
     do ip = 1, this%fmi%nflowpack
       if (this%fmi%iatp(ip) /= 0) cycle
       if (this%isrctype(ip) == 3 .or. this%isrctype(ip) == 4) then
         ssmiptr => this%ssmivec(ip)
         call ssmiptr%spc_rp(this%fmi%gwfpackages(ip)%nbound, &
                             this%fmi%gwfpackages(ip)%budtxt)
       end if
     end do

◆ ssm_term()

subroutine tspssmmodule::ssm_term	(	class(tspssmtype)	this,
		integer(i4b), intent(in)	ipackage,
		integer(i4b), intent(in)	ientry,
		real(dp), intent(out), optional	rrate,
		real(dp), intent(out), optional	rhsval,
		real(dp), intent(out), optional	hcofval,
		real(dp), intent(out), optional	cssm,
		real(dp), intent(out), optional	qssm
	)

This is the primary SSM routine that calculates the matrix coefficient and right-hand-side value for any package and package entry. It returns several different optional variables that are used throughout this package to update matrix terms, budget calculations, and output tables.

Parameters

	this	TspSsmType
[in]	ipackage	package number
[in]	ientry	bound number
[out]	rrate	calculated mass flow rate
[out]	rhsval	calculated rhs value
[out]	hcofval	calculated hcof value
[out]	cssm	calculated source concentration/temperature depending on flow direction
[out]	qssm	water flow rate into model cell from boundary package

Definition at line 255 of file tsp-ssm.f90.

     ! -- dummy
     class(TspSsmType) :: this !< TspSsmType
     integer(I4B), intent(in) :: ipackage !< package number
     integer(I4B), intent(in) :: ientry !< bound number
     real(DP), intent(out), optional :: rrate !< calculated mass flow rate
     real(DP), intent(out), optional :: rhsval !< calculated rhs value
     real(DP), intent(out), optional :: hcofval !< calculated hcof value
     real(DP), intent(out), optional :: cssm !< calculated source concentration/temperature depending on flow direction
     real(DP), intent(out), optional :: qssm !< water flow rate into model cell from boundary package
     ! -- local
     logical(LGP) :: lauxmixed
     integer(I4B) :: n
     integer(I4B) :: nbound_flow
     real(DP) :: qbnd
     real(DP) :: ctmp
     real(DP) :: omega
     real(DP) :: hcoftmp
     real(DP) :: rhstmp
     !
     ! -- retrieve node number, qbnd and iauxpos
     hcoftmp = dzero
     rhstmp = dzero
     ctmp = dzero
     qbnd = dzero
     nbound_flow = this%fmi%gwfpackages(ipackage)%nbound
     n = this%fmi%gwfpackages(ipackage)%nodelist(ientry)
     !
     ! -- If cell is active (ibound > 0) then calculate values
     if (this%ibound(n) > 0) then
       !
       ! -- retrieve qbnd and iauxpos
       qbnd = this%fmi%gwfpackages(ipackage)%get_flow(ientry)
       call this%get_ssm_conc(ipackage, ientry, nbound_flow, ctmp, lauxmixed)
       !
       ! -- assign values for hcoftmp, rhstmp, and ctmp for subsequent assignment
       !    of hcof, rhs, and rate
       if (.not. lauxmixed) then
         !
         ! -- If qbnd is positive, then concentration (or temperature) represents
         !    the inflow concentration (or temperature).  If qbnd is negative,
         !    then the outflow concentration (or temperature) is set equal to the
         !    simulated cell's concentration (or temperature).
         if (qbnd >= dzero) then
           omega = dzero ! rhs
         else
           ctmp = this%cnew(n)
           omega = done ! lhs
           if (ctmp < dzero) then
             omega = dzero ! concentration/temperature is negative, so set mass flux to zero
           end if
         end if
       else
         !
         ! -- lauxmixed value indicates that this is a mixed sink type where
         !    the concentration (or temperature) value represents the injected
         !    concentration (or temperature) if qbnd is positive. If qbnd is
         !    negative, then the withdrawn water is equal to the minimum of the
         !    aux concentration (or temperature) and the cell concentration
         !    (or temperature).
         if (qbnd >= dzero) then
           omega = dzero ! rhs (ctmp is aux value)
         else
           if (ctmp < this%cnew(n)) then
             omega = dzero ! rhs (ctmp is aux value)
           else
             omega = done ! lhs (ctmp is cell concentration)
             ctmp = this%cnew(n)
           end if
         end if
       end if
       !
       ! -- Add terms based on qbnd sign
       if (qbnd <= dzero) then
         hcoftmp = qbnd * omega * this%eqnsclfac
       else
         rhstmp = -qbnd * ctmp * (done - omega) * this%eqnsclfac
       end if
       !
       ! -- end of active ibound
     end if
     !
     ! -- set requested values
     if (present(hcofval)) hcofval = hcoftmp
     if (present(rhsval)) rhsval = rhstmp
     if (present(rrate)) rrate = hcoftmp * ctmp - rhstmp
     if (present(cssm)) cssm = ctmp
     if (present(qssm)) qssm = qbnd

Variable Documentation

◆ ftype

character(len=lenftype) tspssmmodule::ftype = 'SSM'

Definition at line 28 of file tsp-ssm.f90.

28 character(len=LENFTYPE) :: ftype = 'SSM'

◆ text

character(len=lenpackagename) tspssmmodule::text = ' SOURCE-SINK MIX'

Definition at line 29 of file tsp-ssm.f90.

29 character(len=LENPACKAGENAME) :: text = ' SOURCE-SINK MIX'

Data Types

Functions/Subroutines

Variables

Detailed Description

Function/Subroutine Documentation

◆ allocate_arrays()

◆ allocate_scalars()

◆ get_ssm_conc()

◆ pak_setup_outputtab()

◆ set_iauxpak()

◆ set_ssmivec()

◆ source_fileinput()

◆ source_options()

◆ source_sources()

◆ ssm_ad()

◆ ssm_ar()

◆ ssm_bd()

◆ ssm_cq()

◆ ssm_cr()

◆ ssm_da()

◆ ssm_df()

◆ ssm_fc()

◆ ssm_ot_flow()

◆ ssm_rp()

◆ ssm_term()

Variable Documentation

◆ ftype

◆ text