gridarrayloadmodule Module Reference

This module contains the GridArrayLoadModule. More...

Data Types
type	gridarrayloadtype
	Ascii grid based 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	ad (this)
subroutine	rp (this, parser)
subroutine	destroy (this)
subroutine	reset (this)
subroutine	params_alloc (this)
subroutine	param_load (this, parser, idt, mempath, layered, netcdf, iaux)

Detailed Description

This module contains the routines for reading period block array based input associated with the full grid, such as with the GHBA package.

Function/Subroutine Documentation

ad()

subroutine gridarrayloadmodule::ad ( class(gridarrayloadtype), intent(inout)  this )

private

Definition at line 97 of file Mf6FileGridArray.f90.

    class(GridArrayLoadType), intent(inout) :: this

ainit()

subroutine gridarrayloadmodule::ainit ( class(gridarrayloadtype), 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 46 of file Mf6FileGridArray.f90.

    use memorymanagermodule, only: get_isize
    use blockparsermodule, only: blockparsertype
    use loadmf6filemodule, only: loadmf6filetype
    class(GridArrayLoadType), 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
    integer(I4B) :: isize
    integer(I4B), pointer :: maxbound
 
    ! initialize base type
    call this%DynamicPkgLoadType%init(mf6_input, component_name, &
                                      component_input_name, &
                                      input_name, iperblock, iout)
    ! initialize
    this%iout = iout
 
    ! load static input
    call loader%load(parser, mf6_input, this%nc_vars, this%input_name, iout)
 
    ! maxbound is optional
    call get_isize('MAXBOUND', mf6_input%mempath, isize)
    if (isize < 0) then
      ! set maxbound to grid nodes
      call mem_allocate(maxbound, 'MAXBOUND', mf6_input%mempath)
      maxbound = product(loader%mshape)
    end if
 
    ! initialize input context memory
    call this%ctx%init(mf6_input)
 
    ! allocate user nodelist
    call mem_allocate(this%nodeulist, this%ctx%maxbound, &
                      'NODEULIST', mf6_input%mempath)
 
    ! allocate dfn params
    call this%params_alloc()

Here is the call graph for this function:

destroy()

subroutine gridarrayloadmodule::destroy

(

class(gridarrayloadtype), intent(inout)

this

)

Definition at line 164 of file Mf6FileGridArray.f90.

    use memorymanagermodule, only: mem_deallocate
    class(GridArrayLoadType), intent(inout) :: this
    call mem_deallocate(this%nodeulist)

df()

subroutine gridarrayloadmodule::df

(

class(gridarrayloadtype), intent(inout)

this

)

Definition at line 93 of file Mf6FileGridArray.f90.

    class(GridArrayLoadType), intent(inout) :: this

param_load()

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

private

Definition at line 207 of file Mf6FileGridArray.f90.

    use tdismodule, only: kper
    use constantsmodule, only: dnodata
    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 layeredarrayreadermodule, only: read_dbl1d_layered
    use loadncinputmodule, only: netcdf_read_array
    use sourcecommonmodule, only: get_shape_from_string, get_layered_shape
    use idmloggermodule, only: idm_log_var
    class(GridArrayLoadType), intent(inout) :: this
    type(BlockParserType), intent(in) :: parser
    type(InputParamDefinitionType), intent(in) :: idt
    character(len=*), intent(in) :: mempath
    logical(LGP), intent(in) :: layered
    logical(LGP), intent(in) :: netcdf
    real(DP), dimension(:), pointer, contiguous :: dbl1d, nodes
    real(DP), dimension(:, :), pointer, contiguous :: dbl2d
    integer(I4B), dimension(:), allocatable :: layer_shape
    integer(I4B) :: iaux, iparam, n, nlay, nnode
 
    nnode = 0
 
    select case (idt%datatype)
    case ('DOUBLE1D')
      call mem_setptr(dbl1d, idt%mf6varname, mempath)
      allocate (nodes(this%ctx%nodes))
      if (netcdf) then
        call netcdf_read_array(nodes, this%ctx%mshape, idt, &
                               this%mf6_input, this%nc_vars, this%input_name, &
                               this%iout, kper)
      else if (layered) then
        call get_layered_shape(this%ctx%mshape, nlay, layer_shape)
        call read_dbl1d_layered(parser, nodes, idt%mf6varname, nlay, layer_shape)
      else
        call read_dbl1d(parser, nodes, idt%mf6varname)
      end if
 
      call idm_log_var(nodes, idt%tagname, mempath, this%iout)
 
      do n = 1, this%ctx%nodes
        if (nodes(n) /= dnodata) then
          nnode = nnode + 1
          dbl1d(nnode) = nodes(n)
          if (this%ctx%nbound == 0) then
            this%nodeulist(nnode) = n
          else if (this%nodeulist(nnode) /= n) then
            write (errmsg, '(a,i0)') 'Grid input position mismatch param='// &
              trim(idt%tagname)//', period=', kper
            call store_error(errmsg)
            call store_error_filename(this%input_name)
          end if
        end if
      end do
      deallocate (nodes)
    case ('DOUBLE2D')
      call mem_setptr(dbl2d, idt%mf6varname, mempath)
      allocate (nodes(this%ctx%nodes))
 
      if (netcdf) then
        call netcdf_read_array(nodes, this%ctx%mshape, idt, &
                               this%mf6_input, this%nc_vars, this%input_name, &
                               this%iout, kper, iaux)
      else if (layered) then
        call get_layered_shape(this%ctx%mshape, nlay, layer_shape)
        call read_dbl1d_layered(parser, nodes, idt%mf6varname, nlay, layer_shape)
      else
        call read_dbl1d(parser, nodes, idt%mf6varname)
      end if
 
      call idm_log_var(nodes, idt%tagname, mempath, this%iout)
 
      do n = 1, this%ctx%nodes
        if (nodes(n) /= dnodata) then
          nnode = nnode + 1
          dbl2d(iaux, nnode) = nodes(n)
          if (this%ctx%nbound == 0) then
            this%nodeulist(nnode) = n
          else if (this%nodeulist(nnode) /= n) then
            write (errmsg, '(a,i0)') 'Grid input position mismatch param='// &
              trim(idt%tagname)//', period=', kper
            call store_error(errmsg)
            call store_error_filename(this%input_name)
          end if
        end if
      end do
      deallocate (nodes)
    case default
      errmsg = 'IDM unimplemented. GridArrayLoad::param_load &
               &datatype='//trim(idt%datatype)
      call store_error(errmsg)
      call store_error_filename(this%input_name)
    end select
 
    ! set nbound
    if (this%ctx%nbound == 0) this%ctx%nbound = nnode
 
    ! 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 gridarrayloadmodule::params_alloc

(

class(gridarrayloadtype), intent(inout)

this

)

Definition at line 183 of file Mf6FileGridArray.f90.

    class(GridArrayLoadType), 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 gridarrayloadmodule::reset

(

class(gridarrayloadtype), intent(inout)

this

)

Definition at line 170 of file Mf6FileGridArray.f90.

    use constantsmodule, only: dnodata
    class(GridArrayLoadType), intent(inout) :: this
    integer(I4B) :: n
 
    this%ctx%nbound = 0
 
    do n = 1, this%nparam
      ! reset read state
      this%param_reads(n)%invar = 0
    end do

rp()

subroutine gridarrayloadmodule::rp ( class(gridarrayloadtype), intent(inout)  this, type(blockparsertype), intent(inout), pointer  parser  )

private

Definition at line 101 of file Mf6FileGridArray.f90.

    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(GridArrayLoadType), intent(inout) :: this
    type(BlockParserType), pointer, intent(inout) :: parser
    logical(LGP) :: endOfBlock, netcdf, layered
    character(len=LINELENGTH) :: keyword, param_tag
    type(InputParamDefinitionType), pointer :: idt
    integer(I4B) :: iaux
 
    ! 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.
      layered = .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 Layered and NetCDF keywords
      call parser%GetStringCaps(keyword)
      if (keyword == 'LAYERED' .and. idt%layered) then
        layered = .true.
      else if (keyword == 'NETCDF') then
        netcdf = .true.
      end if
 
      ! read and load the parameter
      call this%param_load(parser, idt, this%mf6_input%mempath, layered, &
                           netcdf, iaux)
    end do
 
    ! 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: