swfstomodule Module Reference

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

Data Types
type	swfstotype

Functions/Subroutines
subroutine, public	sto_cr (stoobj, name_model, 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	read_options (this)
	@ brief Read options for package More...
subroutine	read_data (this)
	@ brief Read 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 605 of file swf-sto.f90.

    ! -- modules
    use memorymanagermodule, only: mem_allocate
    ! -- 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)
    !
    ! -- Initialize arrays
    this%iss = 0
    do n = 1, nodes
      this%qsto(n) = dzero
    end do
    !
    ! -- return
    return

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 581 of file swf-sto.f90.

    ! -- modules
    use memorymanagermodule, only: mem_allocate, mem_setptr
    ! -- 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
    !
    ! -- return
    return

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 410 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 446 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:

reach_length_pointer()

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

(

class(swfstotype)

this

)

Parameters

this	this instance

Definition at line 745 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
 

read_data()

subroutine swfstomodule::read_data

(

class(swfstotype)

this

)

Read griddata block for STO package.

Parameters

this	SwfStoType object

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

    ! -- modules
    ! -- dummy variables
    class(SwfStotype) :: this !< SwfStoType object
    ! -- local variables
    character(len=LINELENGTH) :: keyword
    character(len=:), allocatable :: line
    integer(I4B) :: lloc, ierr
    logical :: isfound, endOfBlock
    character(len=24), dimension(1) :: aname
    ! -- formats
    !data
    data aname(1)/'                     XXX'/
    !
    ! -- initialize
    isfound = .false.
    !
    ! -- get stodata block
    call this%parser%GetBlock('GRIDDATA', isfound, ierr)
    if (isfound) then
      write (this%iout, '(1x,a)') 'PROCESSING GRIDDATA'
      do
        call this%parser%GetNextLine(endofblock)
        if (endofblock) exit
        call this%parser%GetStringCaps(keyword)
        call this%parser%GetRemainingLine(line)
        lloc = 1
        select case (keyword)
        case default
          write (errmsg, '(a,a)') 'Unknown GRIDDATA tag: ', &
            trim(keyword)
          call store_error(errmsg)
          call this%parser%StoreErrorUnit()
        end select
      end do
      write (this%iout, '(1x,a)') 'END PROCESSING GRIDDATA'
    else
      write (errmsg, '(a)') 'Required GRIDDATA block not found.'
      call store_error(errmsg)
      call this%parser%StoreErrorUnit()
    end if
    !
    if (count_errors() > 0) then
      call this%parser%StoreErrorUnit()
    end if
    !
    ! -- return
    return

Here is the call graph for this function:

read_options()

subroutine swfstomodule::read_options

(

class(swfstotype)

this

)

Read options block for STO package.

Parameters

this	SwfStoType object

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

    ! -- modules
    ! -- dummy variables
    class(SwfStoType) :: this !< SwfStoType object
    ! -- local variables
    character(len=LINELENGTH) :: keyword
    integer(I4B) :: ierr
    logical :: isfound, endOfBlock
    ! -- formats
    character(len=*), parameter :: fmtisvflow = &
      "(4x,'CELL-BY-CELL FLOW INFORMATION WILL BE SAVED TO BINARY FILE &
      &WHENEVER ICBCFL IS NOT ZERO.')"
    character(len=*), parameter :: fmtflow = &
      &"(4x, 'FLOWS WILL BE SAVED TO FILE: ', a, /4x, 'OPENED ON UNIT: ', I7)"
    !
    ! -- get options block
    call this%parser%GetBlock('OPTIONS', isfound, ierr, &
                              supportopenclose=.true., blockrequired=.false.)
    !
    ! -- parse options block if detected
    if (isfound) then
      write (this%iout, '(1x,a)') 'PROCESSING STORAGE 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, fmtisvflow)
        case default
          write (errmsg, '(a,a)') 'Unknown STO option: ', &
            trim(keyword)
          call store_error(errmsg, terminate=.true.)
        end select
      end do
      write (this%iout, '(1x,a)') 'END OF STORAGE OPTIONS'
    end if
    !
    ! -- return
    return

Here is the call graph for this function:

set_dfw_pointers()

subroutine swfstomodule::set_dfw_pointers

(

class(swfstotype)

this

)

Parameters

this	this instance

Definition at line 732 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:

sto_ad()

subroutine swfstomodule::sto_ad

(

class(swfstotype)

this

)

Advance data in the STO package.

Parameters

this	SwfStoType object

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

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

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 103 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%read_options()
    !
    ! -- read the data block
    ! no griddata at the moment for SWF Storage Package
    ! call this%read_data()
    !
    ! -- return
    return

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 487 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')
    !
    ! -- return
    return

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 367 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
    !
    ! -- Return
    return

sto_cr()

subroutine, public swfstomodule::sto_cr	(	type(swfstotype), pointer	stoobj,
		character(len=*), intent(in)	name_model,
		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]	inunit	package input file unit
[in]	iout	model listing file unit
[in]	cxs	the pointer to the cxs package

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

    ! -- dummy variables
    type(SwfStoType), pointer :: stoobj !< SwfStoType object
    character(len=*), intent(in) :: name_model !< name of model
    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')
    !
    ! -- Allocate scalars
    call stoobj%allocate_scalars()
    !
    ! -- Set variables
    stoobj%inunit = inunit
    stoobj%iout = iout
 
    ! -- store pointers
    stoobj%cxs => cxs
 
    ! -- Initialize block parser
    call stoobj%parser%Initialize(stoobj%inunit, stoobj%iout)
 

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 554 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)
    end if
    !
    ! -- Deallocate scalars
    !
    ! -- deallocate parent
    call this%NumericalPackageType%da()
    !
    ! -- return
    return

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 244 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
    character(len=LINELENGTH) :: distype = ''
    ! -- 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
 
    call this%dis%get_dis_type(distype)
    if (distype == 'DISV1D') 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 285 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
    !
    ! -- Return
    return

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 330 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 150 of file swf-sto.f90.

    ! -- modules
    use tdismodule, only: kper, nper
    implicit none
    ! -- dummy variables
    class(SwfStoType) :: this !< SwfStoType object
    ! -- local variables
    integer(I4B) :: ierr
    logical :: isfound, endOfBlock
    character(len=16) :: css(0:1)
    character(len=LINELENGTH) :: line, keyword
    ! -- formats
    character(len=*), parameter :: fmtlsp = &
      &"(1X,/1X,'REUSING ',A,' FROM LAST STRESS PERIOD')"
    character(len=*), parameter :: fmtblkerr = &
      &"('Error.  Looking for BEGIN PERIOD iper.  Found ', a, ' instead.')"
    ! -- data
    data css(0)/'       TRANSIENT'/
    data css(1)/'    STEADY-STATE'/
    !
    ! -- get stress period data
    if (this%ionper < kper) then
      !
      ! -- get period block
      call this%parser%GetBlock('PERIOD', isfound, ierr, &
                                supportopenclose=.true., &
                                blockrequired=.false.)
      if (isfound) then
        !
        ! -- read ionper and check for increasing period numbers
        call this%read_check_ionper()
      else
        !
        ! -- PERIOD block not found
        if (ierr < 0) then
          ! -- End of file found; data applies for remainder of simulation.
          this%ionper = nper + 1
        else
          ! -- Found invalid block
          call this%parser%GetCurrentLine(line)
          write (errmsg, fmtblkerr) adjustl(trim(line))
          call store_error(errmsg)
          call this%parser%StoreErrorUnit()
        end if
      end if
    end if
 
    if (this%ionper == kper) then
      write (this%iout, '(//,1x,a)') 'PROCESSING STORAGE PERIOD DATA'
      do
        call this%parser%GetNextLine(endofblock)
        if (endofblock) exit
        call this%parser%GetStringCaps(keyword)
        select case (keyword)
        case ('STEADY-STATE')
          this%iss = 1
        case ('TRANSIENT')
          this%iss = 0
        case default
          write (errmsg, '(a,a)') 'Unknown STORAGE data tag: ', &
            trim(keyword)
          call store_error(errmsg)
          call this%parser%StoreErrorUnit()
        end select
      end do
      write (this%iout, '(1x,a)') 'END PROCESSING STORAGE PERIOD DATA'
    end if
 
    write (this%iout, '(//1X,A,I0,A,A,/)') &
      'STRESS PERIOD ', kper, ' IS ', trim(adjustl(css(this%iss)))
    !
    ! -- return
    return

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 513 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
    !
    ! -- return
    return

Variable Documentation

budtxt

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

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

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