Mf6FileStoInput.f90

Go to the documentation of this file.

!> @brief This module contains the Mf6FileStoInputModule
!!
!! This module contains the routines for reading STO period block input
!!
!<
module mf6filestoinputmodule
 
  use kindmodule, only: i4b, dp, lgp
  use constantsmodule, only: linelength
  use inputdefinitionmodule, only: inputparamdefinitiontype
  use memorymanagermodule, only: mem_setptr, mem_allocate
  use modflowinputmodule, only: modflowinputtype
  use asciiinputloadtypemodule, only: asciidynamicpkgloadbasetype
 
  implicit none
  private
  public :: stoinputtype
 
  !> @brief STO package loader
  !!
  !<
  type, extends(asciidynamicpkgloadbasetype) :: stoinputtype
    character(len=LINELENGTH), pointer :: storage => null()
  contains
    procedure :: ainit => sto_init
    procedure :: rp => sto_rp
    procedure :: destroy => sto_destroy
  end type stoinputtype
 
contains
 
  subroutine sto_init(this, mf6_input, component_name, component_input_name, &
                      input_name, iperblock, parser, iout)
    use blockparsermodule, only: blockparsertype
    use loadmf6filemodule, only: loadmf6filetype
    class(stoinputtype), intent(inout) :: this
    type(modflowinputtype), intent(in) :: mf6_input
    character(len=*), intent(in) :: component_name
    character(len=*), intent(in) :: component_input_name
    character(len=*), intent(in) :: input_name
    integer(I4B), intent(in) :: iperblock
    type(blockparsertype), pointer, intent(inout) :: parser
    integer(I4B), intent(in) :: iout
    type(loadmf6filetype) :: loader
 
    ! init loader
    call this%DynamicPkgLoadType%init(mf6_input, component_name, &
                                      component_input_name, input_name, &
                                      iperblock, iout)
    ! initialize static loader
    call loader%load(parser, mf6_input, this%nc_vars, this%input_name, iout)
 
    ! allocate storage string
    call mem_allocate(this%storage, linelength, 'STORAGE', this%mf6_input%mempath)
 
    ! initialize storage to TRANSIENT (model iss=0)
    this%storage = 'TRANSIENT'
  end subroutine sto_init
 
  subroutine sto_rp(this, parser)
    use blockparsermodule, only: blockparsertype
    use definitionselectmodule, only: get_param_definition_type
    use idmloggermodule, only: idm_log_header, idm_log_close, idm_log_var
    class(stoinputtype), intent(inout) :: this
    type(blockparsertype), pointer, intent(inout) :: parser
    character(len=LINELENGTH) :: tagname
    type(inputparamdefinitiontype), pointer :: idt
    logical(LGP) :: endOfBlock
 
    ! read next line
    call parser%GetNextLine(endofblock)
 
    ! return if no input
    if (endofblock) return
 
    ! read the tag
    call parser%GetStringCaps(tagname)
 
    ! verify tag is supported
    idt => get_param_definition_type(this%mf6_input%param_dfns, &
                                     this%mf6_input%component_type, &
                                     this%mf6_input%subcomponent_type, &
                                     'PERIOD', tagname, this%input_name)
    ! set storage
    this%storage = idt%tagname
 
    ! only one input line is expected, terminate block
    call parser%terminateblock()
 
    ! log lst file header
    call idm_log_header(this%mf6_input%component_name, &
                        this%mf6_input%subcomponent_name, this%iout)
 
    call idm_log_var(this%storage, tagname, this%mf6_input%mempath, this%iout)
 
    ! close logging statement
    call idm_log_close(this%mf6_input%component_name, &
                       this%mf6_input%subcomponent_name, this%iout)
  end subroutine sto_rp
 
  subroutine sto_destroy(this)
    use memorymanagermodule, only: mem_deallocate
    class(stoinputtype), intent(inout) :: this
    call mem_deallocate(this%storage, 'STORAGE', this%mf6_input%mempath)
    call this%DynamicPkgLoadType%destroy()
  end subroutine sto_destroy
 
end module mf6filestoinputmodule