swfstomodule Module Reference

This module contains the storage package methods. More...

Data Types
type	stooptionsfoundtype
	flags indicating which STO options were found in the input More...
type	swfstotype

Functions/Subroutines
subroutine, public	sto_cr (stoobj, name_model, mempath, inunit, iout, cxs)
	@ brief Create a new package object More...
subroutine	sto_ar (this, dis, ibound)
	@ brief Allocate and read method for package More...
subroutine	sto_rp (this)
	@ brief Read and prepare method for package More...
subroutine	sto_ad (this)
	@ brief Advance the package More...
subroutine	sto_fc (this, kiter, stage_old, stage_new, matrix_sln, idxglo, rhs)
	@ brief Fill A and right-hand side for the package More...
subroutine	sto_fc_dis1d (this, kiter, stage_old, stage_new, matrix_sln, idxglo, rhs)
	@ brief Fill A and right-hand side for the package More...
subroutine	sto_fc_dis2d (this, kiter, stage_old, stage_new, matrix_sln, idxglo, rhs)
	@ brief Fill A and right-hand side for the package More...
subroutine	sto_cq (this, flowja, stage_new, stage_old)
	@ brief Calculate flows for package More...
subroutine	calc_storage_dis1d (this, n, stage_new, stage_old, dx, qsto, derv)
subroutine	calc_storage_dis2d (this, n, stage_new, stage_old, qsto, derv)
subroutine	sto_bd (this, isuppress_output, model_budget)
	@ brief Model budget calculation for package More...
subroutine	sto_save_model_flows (this, icbcfl, icbcun)
	@ brief Save model flows for package More...
subroutine	sto_da (this)
	@ brief Deallocate package memory More...
subroutine	allocate_scalars (this)
	@ brief Allocate scalars More...
subroutine	allocate_arrays (this, nodes)
	@ brief Allocate package arrays More...
subroutine	source_options (this)
	@ brief Source input options for package More...
subroutine	log_options (this, found)
	@ brief Log found options for package More...
subroutine	source_data (this)
	@ brief Source input data for package More...
subroutine	set_dfw_pointers (this)
	Set pointers to channel properties in DFW Package. More...
real(dp) function, dimension(:), pointer	reach_length_pointer (this)

Variables
character(len=lenbudtxt), dimension(1)	budtxt = [' STORAGE']

Detailed Description

This module contains the methods used to add the effects of storage on the surface water flow equation.

Function/Subroutine Documentation

allocate_arrays()

subroutine swfstomodule::allocate_arrays	(	class(swfstotype), target	this,
		integer(i4b), intent(in)	nodes
	)

Allocate and initialize STO package arrays.

Parameters

	this	SwfStoType object
[in]	nodes	active model nodes

Definition at line 548 of file swf-sto.f90.

    ! -- modules
    use memorymanagermodule, only: mem_allocate, mem_setptr
    ! -- dummy variables
    class(SwfStoType), target :: this !< SwfStoType object
    integer(I4B), intent(in) :: nodes !< active model nodes
    ! -- local variables
    integer(I4B) :: n
    !
    ! -- Allocate arrays
    call mem_allocate(this%qsto, nodes, 'STRGSS', this%memoryPath)
    !
    ! -- set input context pointers
    if (this%inunit > 0) then
      call mem_setptr(this%iper, 'IPER', this%input_mempath)
      call mem_setptr(this%storage, 'STORAGE', this%input_mempath)
    end if
    !
    ! -- Initialize arrays
    this%iss = 0
    do n = 1, nodes
      this%qsto(n) = dzero
    end do

allocate_scalars()

subroutine swfstomodule::allocate_scalars

(

class(swfstotype)

this

)

Allocate and initialize scalars for the STO package. The base numerical package allocate scalars method is also called.

Parameters

this	SwfStoType object

Definition at line 527 of file swf-sto.f90.

    ! -- modules
    use memorymanagermodule, only: mem_allocate
    ! -- dummy variables
    class(SwfStoType) :: this !< SwfStoType object
    !
    ! -- allocate scalars in NumericalPackageType
    call this%NumericalPackageType%allocate_scalars()
    !
    ! -- allocate scalars
    !call mem_allocate(this%xxx, 'XXX', this%memoryPath)
    !
    ! -- initialize scalars
    !this%xxx = 0

calc_storage_dis1d()

subroutine swfstomodule::calc_storage_dis1d	(	class(swfstotype)	this,
		integer(i4b), intent(in)	n,
		real(dp), intent(in)	stage_new,
		real(dp), intent(in)	stage_old,
		real(dp), intent(in)	dx,
		real(dp), intent(inout)	qsto,
		real(dp), intent(inout), optional	derv
	)

Definition at line 363 of file swf-sto.f90.

    ! module
    use tdismodule, only: delt
    use mathutilmodule, only: get_perturbation
    ! dummy
    class(SwfStoType) :: this
    integer(I4B), intent(in) :: n
    real(DP), intent(in) :: stage_new
    real(DP), intent(in) :: stage_old
    real(DP), intent(in) :: dx
    real(DP), intent(inout) :: qsto
    real(DP), intent(inout), optional :: derv
    ! local
    real(DP) :: depth_new
    real(DP) :: depth_old
    real(DP) :: width_n
    real(DP) :: width_m
    real(DP) :: cxs_area_new
    real(DP) :: cxs_area_old
    real(DP) :: cxs_area_eps
    real(DP) :: eps
 
    call this%dis%get_flow_width(n, n, 0, width_n, width_m)
    depth_new = stage_new - this%dis%bot(n)
    depth_old = stage_old - this%dis%bot(n)
    cxs_area_new = this%cxs%get_area(this%idcxs(n), width_n, depth_new)
    cxs_area_old = this%cxs%get_area(this%idcxs(n), width_n, depth_old)
    qsto = (cxs_area_new - cxs_area_old) * dx / delt
    if (present(derv)) then
      eps = get_perturbation(depth_new)
      cxs_area_eps = this%cxs%get_area(this%idcxs(n), width_n, depth_new + eps)
      derv = (cxs_area_eps - cxs_area_new) * dx / delt / eps
    end if
 

Here is the call graph for this function:

calc_storage_dis2d()

subroutine swfstomodule::calc_storage_dis2d	(	class(swfstotype)	this,
		integer(i4b), intent(in)	n,
		real(dp), intent(in)	stage_new,
		real(dp), intent(in)	stage_old,
		real(dp), intent(inout)	qsto,
		real(dp), intent(inout), optional	derv
	)

Definition at line 399 of file swf-sto.f90.

    ! module
    use tdismodule, only: delt
    use mathutilmodule, only: get_perturbation
    ! dummy
    class(SwfStoType) :: this
    integer(I4B), intent(in) :: n
    real(DP), intent(in) :: stage_new
    real(DP), intent(in) :: stage_old
    real(DP), intent(inout) :: qsto
    real(DP), intent(inout), optional :: derv
    ! local
    real(DP) :: area
    real(DP) :: depth_new
    real(DP) :: depth_old
    real(DP) :: depth_eps
    real(DP) :: volume_new
    real(DP) :: volume_old
    real(DP) :: eps
 
    area = this%dis%get_area(n)
    depth_new = stage_new - this%dis%bot(n)
    depth_old = stage_old - this%dis%bot(n)
    volume_new = area * depth_new
    volume_old = area * depth_old
    qsto = (volume_new - volume_old) / delt
 
    if (present(derv)) then
      eps = get_perturbation(depth_new)
      depth_eps = depth_new + eps
      derv = (depth_eps - depth_new) * area / delt / eps
    end if
 

Here is the call graph for this function:

log_options()

subroutine swfstomodule::log_options	(	class(swfstotype)	this,
		type(stooptionsfoundtype), intent(in)	found
	)

Log options block for STO package.

Parameters

this	SwfStoType object

Definition at line 603 of file swf-sto.f90.

    ! -- modules
    use memorymanagerextmodule, only: mem_set_value
    ! -- dummy variables
    class(SwfStoType) :: this !< SwfStoType object
    type(StoOptionsFoundType), intent(in) :: found
    ! -- local variables
    ! -- formats
    character(len=*), parameter :: fmtisvflow = &
      "(4x,'CELL-BY-CELL FLOW INFORMATION WILL BE SAVED TO BINARY FILE &
      &WHENEVER ICBCFL IS NOT ZERO.')"
    !
    write (this%iout, '(1x,a)') 'PROCESSING STORAGE OPTIONS'
    !
    if (found%ipakcb) then
      write (this%iout, fmtisvflow)
    end if
    !
    write (this%iout, '(1x,a)') 'END OF STORAGE OPTIONS'

reach_length_pointer()

real(dp) function, dimension(:), pointer swfstomodule::reach_length_pointer

(

class(swfstotype)

this

)

Parameters

this	this instance

Definition at line 665 of file swf-sto.f90.

    ! dummy
    class(SwfStoType) :: this !< this instance
    ! return
    real(DP), dimension(:), pointer :: ptr
    ! local
    class(DisBaseType), pointer :: dis
 
    ptr => null()
    dis => this%dis
    select type (dis)
    type is (disv1dtype)
      ptr => dis%length
    end select
 

set_dfw_pointers()

subroutine swfstomodule::set_dfw_pointers

(

class(swfstotype)

this

)

Parameters

this	this instance

Definition at line 652 of file swf-sto.f90.

    ! -- modules
    use memorymanagermodule, only: mem_setptr, mem_allocate
    ! -- dummy
    class(SwfStoType) :: this !< this instance
    ! -- local
    character(len=LENMEMPATH) :: dfw_mem_path
 
    dfw_mem_path = create_mem_path(this%name_model, 'DFW')
    call mem_setptr(this%idcxs, 'IDCXS', dfw_mem_path)
 

Here is the call graph for this function:

source_data()

subroutine swfstomodule::source_data

(

class(swfstotype)

this

)

Source griddata block parameters for STO package.

Parameters

this	SwfStoType object

Definition at line 629 of file swf-sto.f90.

    ! -- modules
    use memorymanagerextmodule, only: mem_set_value
    ! -- dummy variables
    class(SwfStotype) :: this !< SwfStoType object
    ! -- local variables
    character(len=24), dimension(1) :: aname
    integer(I4B), dimension(:), pointer, contiguous :: map
    !
    ! -- initialize data
    data aname(1)/'                     XXX'/
    !
    ! -- set map to reduce data input arrays
    map => null()
    if (this%dis%nodes < this%dis%nodesuser) map => this%dis%nodeuser
    !
    ! -- log griddata
    write (this%iout, '(1x,a)') 'PROCESSING GRIDDATA'
    write (this%iout, '(1x,a)') 'END PROCESSING GRIDDATA'

source_options()

subroutine swfstomodule::source_options

(

class(swfstotype)

this

)

Source options block parameters for STO package.

Parameters

this	SwfStoType object

Definition at line 578 of file swf-sto.f90.

    ! -- modules
    use memorymanagerextmodule, only: mem_set_value
    use sourcecommonmodule, only: filein_fname
    ! -- dummy variables
    class(SwfStoType) :: this !< SwfStoType object
    ! -- local variables
    type(StoOptionsFoundType) :: found
    !
    ! -- source package input
    call mem_set_value(this%ipakcb, 'IPAKCB', this%input_mempath, found%ipakcb)
    !
    if (found%ipakcb) then
      this%ipakcb = -1
    end if
    !
    ! -- log found options
    call this%log_options(found)

Here is the call graph for this function:

sto_ad()

subroutine swfstomodule::sto_ad

(

class(swfstotype)

this

)

Advance data in the STO package.

Parameters

this	SwfStoType object

Definition at line 194 of file swf-sto.f90.

    ! -- modules
    ! -- dummy variables
    class(SwfStoType) :: this !< SwfStoType object

sto_ar()

subroutine swfstomodule::sto_ar	(	class(swfstotype)	this,
		class(disbasetype), intent(in), pointer	dis,
		integer(i4b), dimension(:), pointer, contiguous	ibound
	)

Method to allocate and read static data for the STO package.

Parameters

	this	SwfStoType object
[in]	dis	model discretization object
	ibound	model ibound array

Definition at line 107 of file swf-sto.f90.

    ! -- modules
    use memorymanagermodule, only: mem_setptr
    use memoryhelpermodule, only: create_mem_path
    ! -- dummy variables
    class(SwfStoType) :: this !< SwfStoType object
    class(DisBaseType), pointer, intent(in) :: dis !< model discretization object
    integer(I4B), dimension(:), pointer, contiguous :: ibound !< model ibound array
    ! -- local variables
    ! -- formats
    character(len=*), parameter :: fmtsto = &
      "(1x,/1x,'STO -- STORAGE PACKAGE, VERSION 1, 10/27/2023', &
      &' INPUT READ FROM UNIT ', i0, //)"
    !
    ! --print a message identifying the storage package.
    write (this%iout, fmtsto) this%inunit
 
    ! -- set pointers to data in dfw package
    call this%set_dfw_pointers()
 
    !
    ! -- store pointers to arguments that were passed in
    this%dis => dis
    this%ibound => ibound
    !
    ! -- set pointer to model iss
    call mem_setptr(this%iss, 'ISS', create_mem_path(this%name_model))
    !
    ! -- Allocate arrays
    call this%allocate_arrays(dis%nodes)
    !
    ! -- Read storage options
    call this%source_options()
    !
    ! -- read the data block
    ! no griddata at the moment for SWF Storage Package
    ! call this%source_data()

Here is the call graph for this function:

sto_bd()

subroutine swfstomodule::sto_bd	(	class(swfstotype)	this,
		integer(i4b), intent(in)	isuppress_output,
		type(budgettype), intent(inout)	model_budget
	)

Budget calculation for the STO package components. Components include specific storage and specific yield storage.

Parameters

	this	SwfStoType object
[in]	isuppress_output	flag to suppress model output
[in,out]	model_budget	model budget object

Definition at line 440 of file swf-sto.f90.

    ! -- modules
    use tdismodule, only: delt
    use budgetmodule, only: budgettype, rate_accumulator
    ! -- dummy variables
    class(SwfStoType) :: this !< SwfStoType object
    integer(I4B), intent(in) :: isuppress_output !< flag to suppress model output
    type(BudgetType), intent(inout) :: model_budget !< model budget object
    ! -- local variables
    real(DP) :: rin
    real(DP) :: rout
    !
    ! -- Add storage rates to model budget
    call rate_accumulator(this%qsto, rin, rout)
    call model_budget%addentry(rin, rout, delt, '             STO', &
                               isuppress_output, '         STORAGE')

Here is the call graph for this function:

sto_cq()

subroutine swfstomodule::sto_cq	(	class(swfstotype)	this,
		real(dp), dimension(:), intent(inout)	flowja,
		real(dp), dimension(:), intent(inout)	stage_new,
		real(dp), dimension(:), intent(inout)	stage_old
	)

Definition at line 323 of file swf-sto.f90.

    ! -- dummy
    class(SwfStoType) :: this
    real(DP), intent(inout), dimension(:) :: flowja
    real(DP), intent(inout), dimension(:) :: stage_new
    real(DP), intent(inout), dimension(:) :: stage_old
    ! -- local
    real(DP), dimension(:), pointer :: reach_length
    integer(I4B) :: n
    integer(I4B) :: idiag
    real(DP) :: dx
    real(DP) :: q
 
    ! -- test if steady-state stress period
    if (this%iss /= 0) return
 
    ! Set pointer to reach_length for 1d
    reach_length => this%reach_length_pointer()
 
    ! -- Calculate storage term
    do n = 1, this%dis%nodes
      !
      ! -- skip if constant stage
      if (this%ibound(n) < 0) cycle
 
      ! Calculate storage for either the DIS1D or DIS2D cases and
      ! add to flowja
      if (associated(reach_length)) then
        dx = reach_length(n)
        call this%calc_storage_dis1d(n, stage_new(n), stage_old(n), dx, q)
      else
        call this%calc_storage_dis2d(n, stage_new(n), stage_old(n), q)
      end if
      this%qsto(n) = -q
      idiag = this%dis%con%ia(n)
      flowja(idiag) = flowja(idiag) + this%qsto(n)
 
    end do

sto_cr()

subroutine, public swfstomodule::sto_cr	(	type(swfstotype), pointer	stoobj,
		character(len=*), intent(in)	name_model,
		character(len=*), intent(in)	mempath,
		integer(i4b), intent(in)	inunit,
		integer(i4b), intent(in)	iout,
		type(swfcxstype), intent(in), pointer	cxs
	)

Create a new storage (STO) object

Parameters

	stoobj	SwfStoType object
[in]	name_model	name of model
[in]	mempath	input context mem path
[in]	inunit	package input file unit
[in]	iout	model listing file unit
[in]	cxs	the pointer to the cxs package

Definition at line 75 of file swf-sto.f90.

    ! -- dummy variables
    type(SwfStoType), pointer :: stoobj !< SwfStoType object
    character(len=*), intent(in) :: name_model !< name of model
    character(len=*), intent(in) :: mempath !< input context mem path
    integer(I4B), intent(in) :: inunit !< package input file unit
    integer(I4B), intent(in) :: iout !< model listing file unit
    type(SwfCxsType), pointer, intent(in) :: cxs !< the pointer to the cxs package
    !
    ! -- Create the object
    allocate (stoobj)
    !
    ! -- create name and memory path
    call stoobj%set_names(1, name_model, 'STO', 'STO', mempath)
    !
    ! -- Allocate scalars
    call stoobj%allocate_scalars()
    !
    ! -- Set variables
    stoobj%inunit = inunit
    stoobj%iout = iout
 
    ! -- store pointers
    stoobj%cxs => cxs
 

Here is the caller graph for this function:

sto_da()

subroutine swfstomodule::sto_da

(

class(swfstotype)

this

)

Deallocate STO package scalars and arrays.

Parameters

this	SwfStoType object

Definition at line 501 of file swf-sto.f90.

    ! -- modules
    use memorymanagermodule, only: mem_deallocate
    ! -- dummy variables
    class(SwfStoType) :: this !< SwfStoType object
    !
    ! -- Deallocate arrays if package is active
    if (this%inunit > 0) then
      call mem_deallocate(this%qsto)
      nullify (this%idcxs)
      nullify (this%iper)
      nullify (this%storage)
    end if
    !
    ! -- Deallocate scalars
    !
    ! -- deallocate parent
    call this%NumericalPackageType%da()

sto_fc()

subroutine swfstomodule::sto_fc	(	class(swfstotype)	this,
		integer(i4b)	kiter,
		real(dp), dimension(:), intent(inout)	stage_old,
		real(dp), dimension(:), intent(inout)	stage_new,
		class(matrixbasetype), pointer	matrix_sln,
		integer(i4b), dimension(:), intent(in)	idxglo,
		real(dp), dimension(:), intent(inout)	rhs
	)

Fill the coefficient matrix and right-hand side with the STO package terms.

Definition at line 205 of file swf-sto.f90.

    ! -- modules
    use tdismodule, only: delt
    ! -- dummy
    class(SwfStoType) :: this
    integer(I4B) :: kiter
    real(DP), intent(inout), dimension(:) :: stage_old
    real(DP), intent(inout), dimension(:) :: stage_new
    class(MatrixBaseType), pointer :: matrix_sln
    integer(I4B), intent(in), dimension(:) :: idxglo
    real(DP), intent(inout), dimension(:) :: rhs
    ! -- local
    ! -- formats
    character(len=*), parameter :: fmtsperror = &
      &"('Detected time step length of zero.  SWF Storage Package cannot be ', &
      &'used unless delt is non-zero.')"
    !
    ! -- test if steady-state stress period
    if (this%iss /= 0) return
    !
    ! -- Ensure time step length is not zero
    if (delt == dzero) then
      write (errmsg, fmtsperror)
      call store_error(errmsg, terminate=.true.)
    end if
 
    if (this%dis%is_1d()) then
      call this%sto_fc_dis1d(kiter, stage_old, stage_new, matrix_sln, idxglo, rhs)
    else
      call this%sto_fc_dis2d(kiter, stage_old, stage_new, matrix_sln, idxglo, rhs)
    end if
 

Here is the call graph for this function:

sto_fc_dis1d()

subroutine swfstomodule::sto_fc_dis1d	(	class(swfstotype)	this,
		integer(i4b)	kiter,
		real(dp), dimension(:), intent(inout)	stage_old,
		real(dp), dimension(:), intent(inout)	stage_new,
		class(matrixbasetype), pointer	matrix_sln,
		integer(i4b), dimension(:), intent(in)	idxglo,
		real(dp), dimension(:), intent(inout)	rhs
	)

Fill the coefficient matrix and right-hand side with the STO package terms.

Definition at line 244 of file swf-sto.f90.

    ! -- modules
    ! -- dummy
    class(SwfStoType) :: this
    integer(I4B) :: kiter
    real(DP), intent(inout), dimension(:) :: stage_old
    real(DP), intent(inout), dimension(:) :: stage_new
    class(MatrixBaseType), pointer :: matrix_sln
    integer(I4B), intent(in), dimension(:) :: idxglo
    real(DP), intent(inout), dimension(:) :: rhs
    ! -- local
    integer(I4B) :: n, idiag
    real(DP) :: derv
    real(DP) :: qsto
    real(DP), dimension(:), pointer :: reach_length
 
    ! Set pointer to reach_length for 1d
    reach_length => this%reach_length_pointer()
 
    ! -- Calculate coefficients and put into amat
    do n = 1, this%dis%nodes
 
      ! -- skip if constant stage
      if (this%ibound(n) < 0) cycle
 
      call this%calc_storage_dis1d(n, stage_new(n), stage_old(n), &
                                   reach_length(n), qsto, derv)
 
      ! -- Fill amat and rhs
      idiag = this%dis%con%ia(n)
      call matrix_sln%add_value_pos(idxglo(idiag), -derv)
      rhs(n) = rhs(n) + qsto - derv * stage_new(n)
 
    end do

sto_fc_dis2d()

subroutine swfstomodule::sto_fc_dis2d	(	class(swfstotype)	this,
		integer(i4b)	kiter,
		real(dp), dimension(:), intent(inout)	stage_old,
		real(dp), dimension(:), intent(inout)	stage_new,
		class(matrixbasetype), pointer	matrix_sln,
		integer(i4b), dimension(:), intent(in)	idxglo,
		real(dp), dimension(:), intent(inout)	rhs
	)

Fill the coefficient matrix and right-hand side with the STO package terms.

Definition at line 286 of file swf-sto.f90.

    ! -- modules
    ! -- dummy
    class(SwfStoType) :: this
    integer(I4B) :: kiter
    real(DP), intent(inout), dimension(:) :: stage_old
    real(DP), intent(inout), dimension(:) :: stage_new
    class(MatrixBaseType), pointer :: matrix_sln
    integer(I4B), intent(in), dimension(:) :: idxglo
    real(DP), intent(inout), dimension(:) :: rhs
    ! -- local
    integer(I4B) :: n, idiag
    real(DP) :: derv
    real(DP) :: qsto
 
    ! -- Calculate coefficients and put into amat
    do n = 1, this%dis%nodes
      !
      ! -- skip if constant stage
      if (this%ibound(n) < 0) cycle
 
      ! Calculate storage and derivative term
      call this%calc_storage_dis2d(n, stage_new(n), stage_old(n), &
                                   qsto, derv)
 
      ! -- Fill amat and rhs
      idiag = this%dis%con%ia(n)
      call matrix_sln%add_value_pos(idxglo(idiag), -derv)
      rhs(n) = rhs(n) + qsto - derv * stage_new(n)
 
    end do
 

sto_rp()

subroutine swfstomodule::sto_rp

(

class(swfstotype)

this

)

Method to read and prepare stress period data for the STO package.

Parameters

this	SwfStoType object

Definition at line 151 of file swf-sto.f90.

    ! -- modules
    use tdismodule, only: kper
    implicit none
    ! -- dummy variables
    class(SwfStoType) :: this !< SwfStoType object
    ! -- local variables
    character(len=16) :: css(0:1)
    ! -- data
    data css(0)/'       TRANSIENT'/
    data css(1)/'    STEADY-STATE'/
    !
    ! -- confirm package is active
    if (this%inunit <= 0) return
    !
    ! -- confirm loaded iper
    if (this%iper /= kper) return
    !
    write (this%iout, '(//,1x,a)') 'PROCESSING STORAGE PERIOD DATA'
    !
    ! -- set period iss
    if (this%storage == 'STEADY-STATE') then
      this%iss = 1
    else if (this%storage == 'TRANSIENT') then
      this%iss = 0
    else
      write (errmsg, '(a,a)') 'Unknown STORAGE data tag: ', &
        trim(this%storage)
      call store_error(errmsg)
      call store_error_filename(this%input_fname)
    end if
    !
    write (this%iout, '(1x,a)') 'END PROCESSING STORAGE PERIOD DATA'
    !
    write (this%iout, '(//1X,A,I0,A,A,/)') &
      'STRESS PERIOD ', kper, ' IS ', trim(adjustl(css(this%iss)))

Here is the call graph for this function:

sto_save_model_flows()

subroutine swfstomodule::sto_save_model_flows	(	class(swfstotype)	this,
		integer(i4b), intent(in)	icbcfl,
		integer(i4b), intent(in)	icbcun
	)

Save cell-by-cell budget terms for the STO package.

Parameters

	this	SwfStoType object
[in]	icbcfl	flag to output budget data
[in]	icbcun	cell-by-cell file unit number

Definition at line 463 of file swf-sto.f90.

    ! -- dummy variables
    class(SwfStoType) :: this !< SwfStoType object
    integer(I4B), intent(in) :: icbcfl !< flag to output budget data
    integer(I4B), intent(in) :: icbcun !< cell-by-cell file unit number
    ! -- local variables
    integer(I4B) :: ibinun
    integer(I4B) :: iprint, nvaluesp, nwidthp
    character(len=1) :: cdatafmp = ' ', editdesc = ' '
    real(DP) :: dinact
    !
    ! -- Set unit number for binary output
    if (this%ipakcb < 0) then
      ibinun = icbcun
    elseif (this%ipakcb == 0) then
      ibinun = 0
    else
      ibinun = this%ipakcb
    end if
    if (icbcfl == 0) ibinun = 0
    !
    ! -- Record the storage rates if requested
    if (ibinun /= 0) then
      iprint = 0
      dinact = dzero
      !
      ! -- qsto
      call this%dis%record_array(this%qsto, this%iout, iprint, -ibinun, &
                                 budtxt(1), cdatafmp, nvaluesp, &
                                 nwidthp, editdesc, dinact)
    end if

Variable Documentation

budtxt

character(len=lenbudtxt), dimension(1) swfstomodule::budtxt = [' STORAGE']

Definition at line 23 of file swf-sto.f90.

  character(len=LENBUDTXT), dimension(1) :: budtxt = & !< text labels for budget terms
    &['         STORAGE']