layerarrayloadmodule Module Reference

This module contains the LayerArrayLoadModule. More...

Data Types
type	layerarrayloadtype
	Ascii array layer dynamic loader type. More...

Functions/Subroutines
subroutine	ainit (this, mf6_input, component_name, component_input_name, input_name, iperblock, parser, iout)
subroutine	df (this)
subroutine	ts_advance (this)
subroutine	rp (this, parser)
subroutine	destroy (this)
subroutine	reset (this)
subroutine	init_charstr1d (this, varname, input_name)
subroutine	params_alloc (this)
subroutine	param_load (this, parser, idt, mempath, netcdf, iaux)
subroutine	tas_arrays_alloc (this)
subroutine	tas_links_create (this, inunit)

Detailed Description

This module contains the routines for reading period block array based input that is associated with a layer and an layer index array, such as with the EVTA and RCHA packages.

Function/Subroutine Documentation

ainit()

subroutine layerarrayloadmodule::ainit ( class(layerarrayloadtype), 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), intent(inout), pointer  parser, integer(i4b), intent(in)  iout  )

private

Definition at line 56 of file Mf6FileLayerArray.f90.

    use memorymanagermodule, only: get_isize
    use blockparsermodule, only: blockparsertype
    use loadmf6filemodule, only: loadmf6filetype
    class(LayerArrayLoadType), 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
    type(CharacterStringType), dimension(:), pointer, &
      contiguous :: tas_fnames
    character(len=LINELENGTH) :: fname
    integer(I4B) :: tas6_size, n
 
    ! initialize base type
    call this%DynamicPkgLoadType%init(mf6_input, component_name, &
                                      component_input_name, &
                                      input_name, iperblock, iout)
    ! initialize
    nullify (this%aux_tasnames)
    nullify (this%param_tasnames)
    this%tas_active = 0
    this%iout = iout
 
    ! load static input
    call loader%load(parser, mf6_input, this%nc_vars, this%input_name, iout)
 
    ! create tasmanager
    allocate (this%tasmanager)
    call tasmanager_cr(this%tasmanager, modelname=this%mf6_input%component_name, &
                       iout=this%iout)
 
    ! determine if TAS6 files were provided in OPTIONS block
    call get_isize('TAS6_FILENAME', this%mf6_input%mempath, tas6_size)
    if (tas6_size > 0) then
      this%tas_active = 1
      call mem_setptr(tas_fnames, 'TAS6_FILENAME', this%mf6_input%mempath)
      ! add files to tasmanager
      do n = 1, size(tas_fnames)
        fname = tas_fnames(n)
        call this%tasmanager%add_tasfile(fname)
      end do
    end if
 
    ! initialize input context memory
    call this%ctx%init(mf6_input)
 
    ! allocate dfn params
    call this%params_alloc()
 
    ! allocate memory for storing TAS strings
    call this%tas_arrays_alloc()

Here is the call graph for this function:

destroy()

subroutine layerarrayloadmodule::destroy

(

class(layerarrayloadtype), intent(inout)

this

)

Definition at line 226 of file Mf6FileLayerArray.f90.

    class(LayerArrayLoadType), intent(inout) :: this
    !
    ! nullify ctx pointers (including mshape) before deallocate
    call this%ctx%destroy()
    !
    ! deallocate tasmanager
    call this%tasmanager%da()
    deallocate (this%tasmanager)
    nullify (this%tasmanager)

df()

subroutine layerarrayloadmodule::df

(

class(layerarrayloadtype), intent(inout)

this

)

Definition at line 116 of file Mf6FileLayerArray.f90.

    class(LayerArrayLoadType), intent(inout) :: this
    call this%tasmanager%tasmanager_df()

init_charstr1d()

subroutine layerarrayloadmodule::init_charstr1d ( class(layerarrayloadtype)  this, character(len=*), intent(in)  varname, character(len=*), intent(in)  input_name  )

private

Definition at line 257 of file Mf6FileLayerArray.f90.

    use memorymanagermodule, only: mem_setptr
    class(LayerArrayLoadType) :: this
    character(len=*), intent(in) :: varname
    character(len=*), intent(in) :: input_name
    type(CharacterStringType), dimension(:), pointer, &
      contiguous :: charstr1d
    integer(I4B) :: n
    call mem_setptr(charstr1d, varname, this%mf6_input%mempath)
    do n = 1, size(charstr1d)
      charstr1d(n) = ''
    end do

param_load()

subroutine layerarrayloadmodule::param_load ( class(layerarrayloadtype), intent(inout)  this, type(blockparsertype), intent(in)  parser, type(inputparamdefinitiontype), intent(in)  idt, character(len=*), intent(in)  mempath, logical(lgp), intent(in)  netcdf, integer(i4b), intent(in)  iaux  )

private

Definition at line 295 of file Mf6FileLayerArray.f90.

    use tdismodule, only: kper
    use memorymanagermodule, only: mem_setptr
    use arrayhandlersmodule, only: ifind
    use inputdefinitionmodule, only: inputparamdefinitiontype
    use definitionselectmodule, only: get_param_definition_type
    use double1dreadermodule, only: read_dbl1d
    use double2dreadermodule, only: read_dbl2d
    use integer1dreadermodule, only: read_int1d
    use loadncinputmodule, only: netcdf_read_array
    use idmloggermodule, only: idm_log_var
    class(LayerArrayLoadType), intent(inout) :: this
    type(BlockParserType), intent(in) :: parser
    type(InputParamDefinitionType), intent(in) :: idt
    character(len=*), intent(in) :: mempath
    logical(LGP), intent(in) :: netcdf
    integer(I4B), intent(in) :: iaux
    integer(I4B), dimension(:), pointer, contiguous :: int1d
    real(DP), dimension(:), pointer, contiguous :: dbl1d
    real(DP), dimension(:, :), pointer, contiguous :: dbl2d
    integer(I4B) :: iparam, n
 
    select case (idt%datatype)
    case ('INTEGER1D')
      call mem_setptr(int1d, idt%mf6varname, mempath)
      if (netcdf) then
        call netcdf_read_array(int1d, this%ctx%mshape, idt, &
                               this%mf6_input, this%nc_vars, this%input_name, &
                               this%iout, kper)
      else
        call read_int1d(parser, int1d, idt%mf6varname)
      end if
      call idm_log_var(int1d, idt%tagname, mempath, this%iout)
    case ('DOUBLE1D')
      call mem_setptr(dbl1d, idt%mf6varname, mempath)
      if (netcdf) then
        call netcdf_read_array(dbl1d, this%ctx%mshape, idt, &
                               this%mf6_input, this%nc_vars, this%input_name, &
                               this%iout, kper)
      else
        call read_dbl1d(parser, dbl1d, idt%mf6varname)
      end if
      call idm_log_var(dbl1d, idt%tagname, mempath, this%iout)
    case ('DOUBLE2D')
      call mem_setptr(dbl2d, idt%mf6varname, mempath)
      allocate (dbl1d(this%ctx%ncpl))
      if (netcdf) then
        call netcdf_read_array(dbl1d, this%ctx%mshape, idt, &
                               this%mf6_input, this%nc_vars, this%input_name, &
                               this%iout, kper, iaux)
      else
        call read_dbl1d(parser, dbl1d, idt%mf6varname)
      end if
      do n = 1, this%ctx%ncpl
        dbl2d(iaux, n) = dbl1d(n)
      end do
      call idm_log_var(dbl1d, idt%tagname, mempath, this%iout)
      deallocate (dbl1d)
    case default
      errmsg = 'IDM unimplemented. LayerArrayLoad::param_load &
               &datatype='//trim(idt%datatype)
      call store_error(errmsg)
      call store_error_filename(this%input_name)
    end select
 
    ! if param is tracked set read state
    iparam = ifind(this%param_names, idt%tagname)
    if (iparam > 0) then
      this%param_reads(iparam)%invar = 1
    end if

Here is the call graph for this function:

params_alloc()

subroutine layerarrayloadmodule::params_alloc

(

class(layerarrayloadtype), intent(inout)

this

)

Definition at line 271 of file Mf6FileLayerArray.f90.

    class(LayerArrayLoadType), intent(inout) :: this
    character(len=LENVARNAME) :: rs_varname
    integer(I4B), pointer :: intvar
    integer(I4B) :: iparam
 
    ! set in scope param names
    call this%ctx%tags(this%param_names, this%nparam, this%input_name, &
                       create=.true.)
    call this%ctx%allocate_arrays()
 
    ! allocate and set param_reads pointer array
    allocate (this%param_reads(this%nparam))
 
    ! store read state variable pointers
    do iparam = 1, this%nparam
      ! allocate and store name of read state variable
      rs_varname = this%ctx%rsv_alloc(this%param_names(iparam))
      call mem_setptr(intvar, rs_varname, this%mf6_input%mempath)
      this%param_reads(iparam)%invar => intvar
      this%param_reads(iparam)%invar = 0
    end do

reset()

subroutine layerarrayloadmodule::reset ( class(layerarrayloadtype), intent(inout)  this )

private

Definition at line 238 of file Mf6FileLayerArray.f90.

    class(LayerArrayLoadType), intent(inout) :: this
    integer(I4B) :: n
 
    if (this%tas_active /= 0) then
      ! reset tasmanager
      call this%tasmanager%reset(this%mf6_input%subcomponent_name)
      ! reinitialize tas name arrays
      call this%init_charstr1d('AUXTASNAME', this%input_name)
      call this%init_charstr1d('PARAMTASNAME', this%input_name)
    end if
 
    do n = 1, this%nparam
      ! reset read state
      this%param_reads(n)%invar = 0
    end do
 

rp()

subroutine layerarrayloadmodule::rp ( class(layerarrayloadtype), intent(inout)  this, type(blockparsertype), intent(inout), pointer  parser  )

private

Definition at line 126 of file Mf6FileLayerArray.f90.

    use memorymanagermodule, only: mem_setptr
    use blockparsermodule, only: blockparsertype
    use inputdefinitionmodule, only: inputparamdefinitiontype
    use definitionselectmodule, only: get_param_definition_type
    use arrayhandlersmodule, only: ifind
    use sourcecommonmodule, only: ifind_charstr
    use idmloggermodule, only: idm_log_header, idm_log_close, idm_log_var
    class(LayerArrayLoadType), intent(inout) :: this
    type(BlockParserType), pointer, intent(inout) :: parser
    logical(LGP) :: endOfBlock, netcdf
    character(len=LINELENGTH) :: keyword, param_tag
    type(InputParamDefinitionType), pointer :: idt
    integer(I4B) :: iaux, iparam
    character(len=LENTIMESERIESNAME) :: tas_name
    integer(I4B), dimension(:), pointer, contiguous :: int1d
 
    ! reset for this period
    call this%reset()
 
    ! log lst file header
    call idm_log_header(this%mf6_input%component_name, &
                        this%mf6_input%subcomponent_name, this%iout)
 
    ! read array block
    do
      ! initialize
      iaux = 0
      netcdf = .false.
 
      ! read next line
      call parser%GetNextLine(endofblock)
      if (endofblock) exit
      ! read param_tag
      call parser%GetStringCaps(param_tag)
 
      ! is param tag an auxvar?
      iaux = ifind_charstr(this%ctx%auxname_cst, param_tag)
      ! any auvxar corresponds to the definition tag 'AUX'
      if (iaux > 0) param_tag = 'AUX'
 
      ! set input definition
      idt => get_param_definition_type(this%mf6_input%param_dfns, &
                                       this%mf6_input%component_type, &
                                       this%mf6_input%subcomponent_type, &
                                       'PERIOD', param_tag, this%input_name)
      ! look for TAS and NetCDF keywords
      call parser%GetStringCaps(keyword)
      if (keyword == 'TIMEARRAYSERIES') then
        if (this%tas_active /= 0) then
          call parser%GetStringCaps(tas_name)
          if (param_tag == 'AUX') then
            this%aux_tasnames(iaux) = tas_name
          else
            iparam = ifind(this%param_names, param_tag)
            this%param_tasnames(iparam) = tas_name
            this%param_reads(iparam)%invar = 2
          end if
          ! log variable
          call idm_log_var(param_tag, this%mf6_input%mempath, this%iout, .true., &
                           trim(tas_name))
          ! cycle to next input param
          cycle
        else
          ! TODO: throw error
        end if
      else if (keyword == 'NETCDF') then
        netcdf = .true.
      end if
 
      ! read and load the parameter
      call this%param_load(parser, idt, this%mf6_input%mempath, netcdf, iaux)
    end do
 
    ! check if layer index variable was read; default to 1 if not provided.
    ! the layer index variable follows the I<TYPE3> naming convention (e.g.
    ! IEVT, IRCH), consistent with the in_scope check in LoadContext.
    if (this%param_reads(1)%invar == 0) then
      if (this%param_names(1) == &
          'I'//trim(this%mf6_input%subcomponent_type(1:3))) then
        idt => get_param_definition_type(this%mf6_input%param_dfns, &
                                         this%mf6_input%component_type, &
                                         this%mf6_input%subcomponent_type, &
                                         'PERIOD', this%param_names(1), &
                                         this%input_name)
        ! set to default of 1 without updating invar
        call mem_setptr(int1d, idt%mf6varname, this%mf6_input%mempath)
        int1d = 1
      end if
    end if
 
    if (this%tas_active /= 0) then
      call this%tas_links_create(parser%iuactive)
    end if
 
    ! log lst file header
    call idm_log_close(this%mf6_input%component_name, &
                       this%mf6_input%subcomponent_name, this%iout)

Here is the call graph for this function:

tas_arrays_alloc()

subroutine layerarrayloadmodule::tas_arrays_alloc

(

class(layerarrayloadtype), intent(inout)

this

)

Definition at line 367 of file Mf6FileLayerArray.f90.

    use memorymanagermodule, only: mem_allocate
    class(LayerArrayLoadType), intent(inout) :: this
 
    ! count params other than AUX
    if (this%tas_active /= 0) then
      call mem_allocate(this%aux_tasnames, lentimeseriesname, &
                        this%ctx%naux, 'AUXTASNAME', &
                        this%mf6_input%mempath)
      call mem_allocate(this%param_tasnames, lentimeseriesname, this%nparam, &
                        'PARAMTASNAME', this%mf6_input%mempath)
      call this%init_charstr1d('AUXTASNAME', this%input_name)
      call this%init_charstr1d('PARAMTASNAME', this%input_name)
    else
      call mem_allocate(this%aux_tasnames, lentimeseriesname, 0, &
                        'AUXTASNAME', this%mf6_input%mempath)
      call mem_allocate(this%param_tasnames, lentimeseriesname, 0, &
                        'PARAMTASNAME', this%mf6_input%mempath)
    end if

tas_links_create()

subroutine layerarrayloadmodule::tas_links_create	(	class(layerarrayloadtype), intent(inout)	this,
		integer(i4b), intent(in)	inunit
	)

Definition at line 389 of file Mf6FileLayerArray.f90.

    use inputdefinitionmodule, only: inputparamdefinitiontype
    use definitionselectmodule, only: get_param_definition_type
    class(LayerArrayLoadType), intent(inout) :: this
    integer(I4B), intent(in) :: inunit
    type(InputParamDefinitionType), pointer :: idt
    ! non-contiguous because a slice of bound is passed
    real(DP), dimension(:), pointer :: auxArrayPtr, bndArrayPtr
    real(DP), dimension(:), pointer, contiguous :: bound
    integer(I4B), dimension(:), pointer, contiguous :: nodelist
    character(len=LENTIMESERIESNAME) :: tas_name
    character(len=LENAUXNAME) :: aux_name
    logical :: convertFlux
    integer(I4B) :: n
 
    ! initialize
    nullify (auxarrayptr)
    nullify (bndarrayptr)
    nullify (nodelist)
    convertflux = .false.
 
    ! Create AUX Time Array Series links
    do n = 1, this%ctx%naux
      tas_name = this%aux_tasnames(n)
      if (tas_name /= '') then
        ! set auxvar pointer
        auxarrayptr => this%ctx%auxvar(n, :)
        aux_name = this%ctx%auxname_cst(n)
        call this%tasmanager%MakeTasLink(this%mf6_input%subcomponent_name, &
                                         auxarrayptr, this%ctx%iprpak, &
                                         tas_name, aux_name, convertflux, &
                                         nodelist, inunit)
      end if
    end do
 
    ! Create BND Time Array Series links
    do n = 1, this%nparam
      ! assign param definition pointer
      idt => get_param_definition_type(this%mf6_input%param_dfns, &
                                       this%mf6_input%component_type, &
                                       this%mf6_input%subcomponent_type, &
                                       'PERIOD', this%param_names(n), &
                                       this%input_name)
      if (idt%timeseries) then
        if (this%param_reads(n)%invar == 2) then
          tas_name = this%param_tasnames(n)
          call mem_setptr(bound, idt%mf6varname, this%mf6_input%mempath)
          ! set bound pointer
          bndarrayptr => bound(:)
          call this%tasmanager%MakeTasLink(this%mf6_input%subcomponent_name, &
                                           bndarrayptr, &
                                           this%ctx%iprpak, &
                                           tas_name, idt%mf6varname, &
                                           convertflux, nodelist, inunit)
        end if
      end if
    end do

Here is the call graph for this function:

ts_advance()

subroutine layerarrayloadmodule::ts_advance ( class(layerarrayloadtype), intent(inout)  this )

private

Definition at line 121 of file Mf6FileLayerArray.f90.

    class(LayerArrayLoadType), intent(inout) :: this
    call this%tasmanager%ad()