loadmf6filemodule Module Reference

This module contains the LoadMf6FileModule. More...

Data Types
type	loadmf6filetype
	Static parser based input loader. More...

Functions/Subroutines
subroutine	load (this, parser, mf6_input, nc_vars, filename, iout)
	load all static input blocks More...
subroutine	init (this, parser, mf6_input, filename, iout)
	init More...
subroutine	load_block (this, iblk)
	load a single block More...
subroutine	finalize (this)
	finalize More...
subroutine	block_post_process (this, iblk)
	Post parse block handling. More...
recursive subroutine	parse_block (this, iblk, recursive_call)
	parse block More...
subroutine	parse_io_tag (this, iblk, pkgtype, which, tag)
subroutine	parse_keyword_tag (this, iblk, tag, idt)
recursive subroutine	parse_tag (this, iblk, recursive_call)
	load an individual input record into memory More...
type(inputparamdefinitiontype) function	block_index_dfn (this, iblk)
subroutine	parse_structarray_block (this, iblk)
	parse a structured array record into memory manager More...
subroutine	load_keyword_type (parser, idt, memoryPath, iout)
	load type keyword More...
subroutine	load_string_type (parser, idt, memoryPath, iout)
	load type string More...
subroutine	load_io_tag (parser, idt, memoryPath, which, iout)
	load io tag More...
subroutine	load_auxvar_names (parser, idt, memoryPath, iout)
	load aux variable names More...
subroutine	load_integer_type (parser, idt, memoryPath, iout)
	load type integer More...
subroutine	load_integer1d_type (parser, idt, mf6_input, mshape, export, nc_vars, input_fname, iout)
	load type 1d integer More...
subroutine	load_integer2d_type (parser, idt, mf6_input, mshape, export, nc_vars, input_fname, iout)
	load type 2d integer More...
subroutine	load_integer3d_type (parser, idt, mf6_input, mshape, export, nc_vars, input_fname, iout)
	load type 3d integer More...
subroutine	load_double_type (parser, idt, memoryPath, iout)
	load type double More...
subroutine	load_double1d_type (parser, idt, mf6_input, mshape, export, nc_vars, input_fname, iout)
	load type 1d double More...
subroutine	load_double2d_type (parser, idt, mf6_input, mshape, export, nc_vars, input_fname, iout)
	load type 2d double More...
subroutine	load_double3d_type (parser, idt, mf6_input, mshape, export, nc_vars, input_fname, iout)
	load type 3d double More...
integer(i4b) function, public	read_control_record (parser, oc_inunit, iout)

Detailed Description

This module contains the input data model routines for loading static data from a MODFLOW 6 input file using the block parser.

Function/Subroutine Documentation

block_index_dfn()

type(inputparamdefinitiontype) function loadmf6filemodule::block_index_dfn	(	class(loadmf6filetype)	this,
		integer(i4b), intent(in)	iblk
	)

Returns

input data type object describing this record

Definition at line 519 of file LoadMf6File.f90.

    ! -- modules
    ! -- dummy
    class(LoadMf6FileType) :: this
    integer(I4B), intent(in) :: iblk
    ! -- local
    type(InputParamDefinitionType) :: idt !< input data type object describing this record
    character(len=LENVARNAME) :: varname
    integer(I4B) :: ilen
    character(len=3) :: block_suffix = 'NUM'
    !
    ! -- assign first column as the block number
    ilen = len_trim(this%mf6_input%block_dfns(iblk)%blockname)
    !
    if (ilen > (lenvarname - len(block_suffix))) then
      varname = &
        this%mf6_input%block_dfns(iblk)% &
        blockname(1:(lenvarname - len(block_suffix)))//block_suffix
    else
      varname = trim(this%mf6_input%block_dfns(iblk)%blockname)//block_suffix
    end if
    !
    idt%component_type = trim(this%mf6_input%component_type)
    idt%subcomponent_type = trim(this%mf6_input%subcomponent_type)
    idt%blockname = trim(this%mf6_input%block_dfns(iblk)%blockname)
    idt%tagname = varname
    idt%mf6varname = varname
    idt%datatype = 'INTEGER'

block_post_process()

subroutine loadmf6filemodule::block_post_process	(	class(loadmf6filetype)	this,
		integer(i4b), intent(in)	iblk
	)

Definition at line 214 of file LoadMf6File.f90.

    ! -- modules
    use constantsmodule, only: lenboundname
    use characterstringmodule, only: characterstringtype
    use sourcecommonmodule, only: set_model_shape
    ! -- dummy
    class(LoadMf6FileType) :: this
    integer(I4B), intent(in) :: iblk
    ! -- local
    type(InputParamDefinitionType), pointer :: idt
    integer(I4B) :: iparam
    integer(I4B), pointer :: intptr
    !
    ! -- update state based on read tags
    do iparam = 1, size(this%block_tags)
      select case (this%mf6_input%block_dfns(iblk)%blockname)
      case ('OPTIONS')
        if (this%block_tags(iparam) == 'AUXILIARY') then
          this%iauxiliary = 1
        else if (this%block_tags(iparam) == 'BOUNDNAMES') then
          this%inamedbound = 1
        else if (this%block_tags(iparam) == 'READASARRAYS') then
          this%readasarrays = .true.
        else if (this%block_tags(iparam) == 'TS6') then
          this%ts_active = .true.
        else if (this%block_tags(iparam) == 'EXPORT_ARRAY_ASCII') then
          this%export = .true.
        end if
      case default
      end select
    end do
    !
    ! -- update input context allocations based on dfn set and input
    select case (this%mf6_input%block_dfns(iblk)%blockname)
    case ('OPTIONS')
      ! -- allocate naux and set to 0 if not allocated
      do iparam = 1, size(this%mf6_input%param_dfns)
        idt => this%mf6_input%param_dfns(iparam)
        !
        if (idt%blockname == 'OPTIONS' .and. &
            idt%tagname == 'AUXILIARY') then
          if (this%iauxiliary == 0) then
            call mem_allocate(intptr, 'NAUX', this%mf6_input%mempath)
            intptr = 0
          end if
          exit
        end if
      end do
    case ('DIMENSIONS')
      ! -- set model shape if discretization dimensions have been read
      if (this%mf6_input%pkgtype(1:3) == 'DIS') then
        call set_model_shape(this%mf6_input%pkgtype, this%filename, &
                             this%mf6_input%component_mempath, &
                             this%mf6_input%mempath, this%mshape)
      end if
    case default
    end select

Here is the call graph for this function:

finalize()

subroutine loadmf6filemodule::finalize

(

class(loadmf6filetype)

this

)

init / finalize are only used when load_block() will be called

Definition at line 193 of file LoadMf6File.f90.

    ! -- modules
    use structarraymodule, only: destructstructarray
    ! -- dummy
    class(LoadMf6FileType) :: this
    ! -- local
    !
    ! -- cleanup
    if (associated(this%structarray)) then
      ! -- destroy the structured array reader
      call destructstructarray(this%structarray)
    end if
    !
    ! -- close logging block
    call idm_log_close(this%mf6_input%component_name, &
                       this%mf6_input%subcomponent_name, this%iout)

Here is the call graph for this function:

init()

subroutine loadmf6filemodule::init	(	class(loadmf6filetype)	this,
		type(blockparsertype), intent(inout), target	parser,
		type(modflowinputtype), intent(in)	mf6_input,
		character(len=*), intent(in)	filename,
		integer(i4b), intent(in)	iout
	)

init / finalize are only used when load_block() will be called

Definition at line 122 of file LoadMf6File.f90.

    ! -- modules
    use memorymanagermodule, only: get_isize
    ! -- dummy
    class(LoadMf6FileType) :: this
    type(BlockParserType), target, intent(inout) :: parser
    type(ModflowInputType), intent(in) :: mf6_input
    character(len=*), intent(in) :: filename
    integer(I4B), intent(in) :: iout
    ! -- local
    integer(I4B) :: isize
    !
    this%parser => parser
    this%mf6_input = mf6_input
    this%filename = filename
    this%ts_active = .false.
    this%export = .false.
    this%readasarrays = .false.
    this%inamedbound = 0
    this%iauxiliary = 0
    this%iout = iout
    !
    call get_isize('MODEL_SHAPE', mf6_input%component_mempath, isize)
    !
    if (isize > 0) then
      call mem_setptr(this%mshape, 'MODEL_SHAPE', mf6_input%component_mempath)
    end if
    !
    ! -- log lst file header
    call idm_log_header(this%mf6_input%component_name, &
                        this%mf6_input%subcomponent_name, this%iout)

Here is the call graph for this function:

load()

subroutine loadmf6filemodule::load ( class(loadmf6filetype)  this, type(blockparsertype), intent(inout), target  parser, type(modflowinputtype), intent(in)  mf6_input, type(ncpackagevarstype), intent(in), pointer  nc_vars, character(len=*), intent(in)  filename, integer(i4b), intent(in)  iout  )

private

Invoke this routine to load all static input blocks in single call.

Definition at line 84 of file LoadMf6File.f90.

    ! -- modules
    use memorymanagermodule, only: get_isize
    ! -- dummy
    class(LoadMf6FileType) :: this
    type(BlockParserType), target, intent(inout) :: parser
    type(ModflowInputType), intent(in) :: mf6_input
    type(NCPackageVarsType), pointer, intent(in) :: nc_vars
    character(len=*), intent(in) :: filename
    integer(I4B), intent(in) :: iout
    ! -- local
    integer(I4B) :: iblk
    !
    ! -- initialize static load
    call this%init(parser, mf6_input, filename, iout)
    !
    this%nc_vars => nc_vars
    !
    ! -- process blocks
    do iblk = 1, size(this%mf6_input%block_dfns)
      !
      ! -- don't load dynamic input data
      if (this%mf6_input%block_dfns(iblk)%blockname == 'PERIOD') exit
      !
      ! -- load the block
      call this%load_block(iblk)
      !
    end do
    !
    ! -- finalize static load
    call this%finalize()

Here is the call graph for this function:

load_auxvar_names()

subroutine loadmf6filemodule::load_auxvar_names	(	type(blockparsertype), intent(inout)	parser,
		type(inputparamdefinitiontype), intent(in)	idt,
		character(len=*), intent(in)	memoryPath,
		integer(i4b), intent(in)	iout
	)

Parameters

[in,out]	parser	block parser
[in]	idt	input data type object describing this record
[in]	memorypath	memorypath to put loaded information
[in]	iout	unit number for output

Definition at line 749 of file LoadMf6File.f90.

    use constantsmodule, only: lenauxname, linelength, lenpackagename
    use inputoutputmodule, only: urdaux
    use characterstringmodule, only: characterstringtype
    type(BlockParserType), intent(inout) :: parser !< block parser
    type(InputParamDefinitionType), intent(in) :: idt !< input data type object describing this record
    character(len=*), intent(in) :: memoryPath !< memorypath to put loaded information
    integer(I4B), intent(in) :: iout !< unit number for output
    character(len=:), allocatable :: line
    character(len=LENAUXNAME), dimension(:), allocatable :: caux
    integer(I4B) :: lloc
    integer(I4B) :: istart
    integer(I4B) :: istop
    integer(I4B) :: i
    character(len=LENPACKAGENAME) :: text = ''
    integer(I4B), pointer :: intvar
    type(CharacterStringType), dimension(:), &
      pointer, contiguous :: acharstr1d !< variable for allocation
    call mem_allocate(intvar, idt%shape, memorypath)
    intvar = 0
    call parser%GetRemainingLine(line)
    lloc = 1
    call urdaux(intvar, parser%iuactive, iout, lloc, &
                istart, istop, caux, line, text)
    call mem_allocate(acharstr1d, lenauxname, intvar, idt%mf6varname, memorypath)
    do i = 1, intvar
      acharstr1d(i) = caux(i)
    end do
    deallocate (line)
    deallocate (caux)
    return

Here is the call graph for this function:

Here is the caller graph for this function:

load_block()

subroutine loadmf6filemodule::load_block	(	class(loadmf6filetype)	this,
		integer(i4b), intent(in)	iblk
	)

Assumed in order load of single (next) block. If a StructArray object is allocated to load this block it persists until this routine (or finalize) is called again.

Definition at line 163 of file LoadMf6File.f90.

    ! -- modules
    use structarraymodule, only: destructstructarray
    ! -- dummy
    class(LoadMf6FileType) :: this
    integer(I4B), intent(in) :: iblk
    ! -- local
    !
    ! -- reset structarray if it was created for previous block
    if (associated(this%structarray)) then
      ! -- destroy the structured array reader
      call destructstructarray(this%structarray)
    end if
    !
    allocate (this%block_tags(0))
    !
    ! -- load the block
    call this%parse_block(iblk, .false.)
    !
    ! -- post process block
    call this%block_post_process(iblk)
    !
    deallocate (this%block_tags)

Here is the call graph for this function:

load_double1d_type()

subroutine loadmf6filemodule::load_double1d_type ( type(blockparsertype), intent(inout)  parser, type(inputparamdefinitiontype), intent(in)  idt, type(modflowinputtype), intent(in)  mf6_input, integer(i4b), dimension(:), intent(in), pointer, contiguous  mshape, logical(lgp), intent(in)  export, type(ncpackagevarstype), intent(in), pointer  nc_vars, character(len=*), intent(in)  input_fname, integer(i4b), intent(in)  iout  )

private

Parameters

[in,out]	parser	block parser
[in]	idt	input data type object describing this record
[in]	mf6_input	description of input
[in]	mshape	model shape
[in]	export	export to ascii layer files
[in]	input_fname	ascii input file name
[in]	iout	unit number for output

Definition at line 994 of file LoadMf6File.f90.

    use sourcecommonmodule, only: get_shape_from_string, get_layered_shape
    use loadncinputmodule, only: netcdf_read_array
    type(BlockParserType), intent(inout) :: parser !< block parser
    type(InputParamDefinitionType), intent(in) :: idt !< input data type object describing this record
    type(ModflowInputType), intent(in) :: mf6_input !< description of input
    integer(I4B), dimension(:), contiguous, pointer, intent(in) :: mshape !< model shape
    logical(LGP), intent(in) :: export !< export to ascii layer files
    type(NCPackageVarsType), pointer, intent(in) :: nc_vars
    character(len=*), intent(in) :: input_fname !< ascii input file name
    integer(I4B), intent(in) :: iout !< unit number for output
    real(DP), dimension(:), pointer, contiguous :: dbl1d
    integer(I4B) :: nlay
    integer(I4B) :: nvals
    integer(I4B), dimension(:), allocatable :: array_shape
    integer(I4B), dimension(:), allocatable :: layer_shape
    character(len=LINELENGTH) :: keyword
 
    ! Check if it is a full grid sized array (NODES)
    if (idt%shape == 'NODES') then
      nvals = product(mshape)
    else
      call get_shape_from_string(idt%shape, array_shape, mf6_input%mempath)
      nvals = array_shape(1)
    end if
 
    ! allocate memory for the array
    call mem_allocate(dbl1d, nvals, idt%mf6varname, mf6_input%mempath)
 
    ! read keyword
    keyword = ''
    call parser%GetStringCaps(keyword)
 
    ! check for "NETCDF" and "LAYERED"
    if (keyword == 'NETCDF') then
      call netcdf_read_array(dbl1d, mshape, idt, mf6_input, nc_vars, &
                             input_fname, iout)
    else if (keyword == 'LAYERED' .and. idt%layered) then
      call get_layered_shape(mshape, nlay, layer_shape)
      call read_dbl1d_layered(parser, dbl1d, idt%mf6varname, nlay, layer_shape)
    else
      call read_dbl1d(parser, dbl1d, idt%mf6varname)
    end if
 
    ! log information on the loaded array to the list file
    call idm_log_var(dbl1d, idt%tagname, mf6_input%mempath, iout)
 
    ! create export file for griddata parameters if optioned
    if (export) then
      if (idt%blockname == 'GRIDDATA') then
        call idm_export(dbl1d, idt%tagname, mf6_input%mempath, idt%shape, iout)
      end if
    end if
 
    return

Here is the call graph for this function:

Here is the caller graph for this function:

load_double2d_type()

subroutine loadmf6filemodule::load_double2d_type	(	type(blockparsertype), intent(inout)	parser,
		type(inputparamdefinitiontype), intent(in)	idt,
		type(modflowinputtype), intent(in)	mf6_input,
		integer(i4b), dimension(:), intent(in), pointer, contiguous	mshape,
		logical(lgp), intent(in)	export,
		type(ncpackagevarstype), intent(in), pointer	nc_vars,
		character(len=*), intent(in)	input_fname,
		integer(i4b), intent(in)	iout
	)

Parameters

[in,out]	parser	block parser
[in]	idt	input data type object describing this record
[in]	mf6_input	description of input
[in]	mshape	model shape
[in]	export	export to ascii layer files
[in]	input_fname	ascii input file name
[in]	iout	unit number for output

Definition at line 1054 of file LoadMf6File.f90.

    use sourcecommonmodule, only: get_shape_from_string, get_layered_shape
    use loadncinputmodule, only: netcdf_read_array
    type(BlockParserType), intent(inout) :: parser !< block parser
    type(InputParamDefinitionType), intent(in) :: idt !< input data type object describing this record
    type(ModflowInputType), intent(in) :: mf6_input !< description of input
    integer(I4B), dimension(:), contiguous, pointer, intent(in) :: mshape !< model shape
    logical(LGP), intent(in) :: export !< export to ascii layer files
    type(NCPackageVarsType), pointer, intent(in) :: nc_vars
    character(len=*), intent(in) :: input_fname !< ascii input file name
    integer(I4B), intent(in) :: iout !< unit number for output
    real(DP), dimension(:, :), pointer, contiguous :: dbl2d
    integer(I4B) :: nlay
    integer(I4B) :: nsize1, nsize2
    integer(I4B), dimension(:), allocatable :: array_shape
    integer(I4B), dimension(:), allocatable :: layer_shape
    character(len=LINELENGTH) :: keyword
 
    ! determine the array shape from the input data definition (idt%shape),
    ! which looks like "NCOL, NROW, NLAY"
    call get_shape_from_string(idt%shape, array_shape, mf6_input%mempath)
    nsize1 = array_shape(1)
    nsize2 = array_shape(2)
 
    ! create a new 3d memory managed variable
    call mem_allocate(dbl2d, nsize1, nsize2, idt%mf6varname, mf6_input%mempath)
 
    ! read keyword
    keyword = ''
    call parser%GetStringCaps(keyword)
 
    ! check for "NETCDF" and "LAYERED"
    if (keyword == 'NETCDF') then
      call netcdf_read_array(dbl2d, mshape, idt, mf6_input, nc_vars, &
                             input_fname, iout)
    else if (keyword == 'LAYERED' .and. idt%layered) then
      call get_layered_shape(mshape, nlay, layer_shape)
      call read_dbl2d_layered(parser, dbl2d, idt%mf6varname, nlay, layer_shape)
    else
      call read_dbl2d(parser, dbl2d, idt%mf6varname)
    end if
 
    ! log information on the loaded array to the list file
    call idm_log_var(dbl2d, idt%tagname, mf6_input%mempath, iout)
 
    ! create export file for griddata parameters if optioned
    if (export) then
      if (idt%blockname == 'GRIDDATA') then
        call idm_export(dbl2d, idt%tagname, mf6_input%mempath, idt%shape, iout)
      end if
    end if
 
    return

Here is the call graph for this function:

Here is the caller graph for this function:

load_double3d_type()

subroutine loadmf6filemodule::load_double3d_type	(	type(blockparsertype), intent(inout)	parser,
		type(inputparamdefinitiontype), intent(in)	idt,
		type(modflowinputtype), intent(in)	mf6_input,
		integer(i4b), dimension(:), intent(in), pointer, contiguous	mshape,
		logical(lgp), intent(in)	export,
		type(ncpackagevarstype), intent(in), pointer	nc_vars,
		character(len=*), intent(in)	input_fname,
		integer(i4b), intent(in)	iout
	)

Parameters

[in,out]	parser	block parser
[in]	idt	input data type object describing this record
[in]	mf6_input	description of input
[in]	mshape	model shape
[in]	export	export to ascii layer files
[in]	input_fname	ascii input file name
[in]	iout	unit number for output

Definition at line 1112 of file LoadMf6File.f90.

    use sourcecommonmodule, only: get_shape_from_string, get_layered_shape
    use loadncinputmodule, only: netcdf_read_array
    type(BlockParserType), intent(inout) :: parser !< block parser
    type(InputParamDefinitionType), intent(in) :: idt !< input data type object describing this record
    type(ModflowInputType), intent(in) :: mf6_input !< description of input
    integer(I4B), dimension(:), contiguous, pointer, intent(in) :: mshape !< model shape
    logical(LGP), intent(in) :: export !< export to ascii layer files
    type(NCPackageVarsType), pointer, intent(in) :: nc_vars
    character(len=*), intent(in) :: input_fname !< ascii input file name
    integer(I4B), intent(in) :: iout !< unit number for output
    real(DP), dimension(:, :, :), pointer, contiguous :: dbl3d
    integer(I4B) :: nlay
    integer(I4B) :: nsize1, nsize2, nsize3
    integer(I4B), dimension(:), allocatable :: array_shape
    integer(I4B), dimension(:), allocatable :: layer_shape
    real(DP), dimension(:), pointer, contiguous :: dbl1d_ptr
    character(len=LINELENGTH) :: keyword
 
    ! determine the array shape from the input data definition (idt%shape),
    ! which looks like "NCOL, NROW, NLAY"
    call get_shape_from_string(idt%shape, array_shape, mf6_input%mempath)
    nsize1 = array_shape(1)
    nsize2 = array_shape(2)
    nsize3 = array_shape(3)
 
    ! create a new 3d memory managed variable
    call mem_allocate(dbl3d, nsize1, nsize2, nsize3, idt%mf6varname, &
                      mf6_input%mempath)
 
    ! read keyword
    keyword = ''
    call parser%GetStringCaps(keyword)
 
    ! check for "NETCDF" and "LAYERED"
    if (keyword == 'NETCDF') then
      call netcdf_read_array(dbl3d, mshape, idt, mf6_input, nc_vars, &
                             input_fname, iout)
    else if (keyword == 'LAYERED' .and. idt%layered) then
      call get_layered_shape(mshape, nlay, layer_shape)
      call read_dbl3d_layered(parser, dbl3d, idt%mf6varname, nlay, &
                              layer_shape)
    else
      dbl1d_ptr(1:nsize1 * nsize2 * nsize3) => dbl3d(:, :, :)
      call read_dbl1d(parser, dbl1d_ptr, idt%mf6varname)
    end if
 
    ! log information on the loaded array to the list file
    call idm_log_var(dbl3d, idt%tagname, mf6_input%mempath, iout)
 
    ! create export file for griddata parameters if optioned
    if (export) then
      if (idt%blockname == 'GRIDDATA') then
        call idm_export(dbl3d, idt%tagname, mf6_input%mempath, idt%shape, iout)
      end if
    end if
 
    return

Here is the call graph for this function:

Here is the caller graph for this function:

load_double_type()

subroutine loadmf6filemodule::load_double_type	(	type(blockparsertype), intent(inout)	parser,
		type(inputparamdefinitiontype), intent(in)	idt,
		character(len=*), intent(in)	memoryPath,
		integer(i4b), intent(in)	iout
	)

Parameters

[in,out]	parser	block parser
[in]	idt	input data type object describing this record
[in]	memorypath	memorypath to put loaded information
[in]	iout	unit number for output

Definition at line 980 of file LoadMf6File.f90.

    type(BlockParserType), intent(inout) :: parser !< block parser
    type(InputParamDefinitionType), intent(in) :: idt !< input data type object describing this record
    character(len=*), intent(in) :: memoryPath !< memorypath to put loaded information
    integer(I4B), intent(in) :: iout !< unit number for output
    real(DP), pointer :: dblvar
    call mem_allocate(dblvar, idt%mf6varname, memorypath)
    dblvar = parser%GetDouble()
    call idm_log_var(dblvar, idt%tagname, memorypath, iout)
    return

Here is the caller graph for this function:

load_integer1d_type()

subroutine loadmf6filemodule::load_integer1d_type ( type(blockparsertype), intent(inout)  parser, type(inputparamdefinitiontype), intent(in)  idt, type(modflowinputtype), intent(in)  mf6_input, integer(i4b), dimension(:), intent(in), pointer, contiguous  mshape, logical(lgp), intent(in)  export, type(ncpackagevarstype), intent(in), pointer  nc_vars, character(len=*), intent(in)  input_fname, integer(i4b), intent(in)  iout  )

private

Parameters

[in,out]	parser	block parser
[in]	idt	input data type object describing this record
[in]	mf6_input	description of input
[in]	mshape	model shape
[in]	export	export to ascii layer files
[in]	input_fname	ascii input file name
[in]	iout	unit number for output

Definition at line 798 of file LoadMf6File.f90.

    use sourcecommonmodule, only: get_shape_from_string, get_layered_shape
    use loadncinputmodule, only: netcdf_read_array
    type(BlockParserType), intent(inout) :: parser !< block parser
    type(InputParamDefinitionType), intent(in) :: idt !< input data type object describing this record
    type(ModflowInputType), intent(in) :: mf6_input !< description of input
    integer(I4B), dimension(:), contiguous, pointer, intent(in) :: mshape !< model shape
    logical(LGP), intent(in) :: export !< export to ascii layer files
    type(NCPackageVarsType), pointer, intent(in) :: nc_vars
    character(len=*), intent(in) :: input_fname !< ascii input file name
    integer(I4B), intent(in) :: iout !< unit number for output
    integer(I4B), dimension(:), pointer, contiguous :: int1d
    integer(I4B) :: nlay
    integer(I4B) :: nvals
    integer(I4B), dimension(:), allocatable :: array_shape
    integer(I4B), dimension(:), allocatable :: layer_shape
    character(len=LINELENGTH) :: keyword
 
    ! Check if it is a full grid sized array (NODES), otherwise use
    ! idt%shape to construct shape from variables in memoryPath
    if (idt%shape == 'NODES') then
      nvals = product(mshape)
    else
      call get_shape_from_string(idt%shape, array_shape, mf6_input%mempath)
      nvals = array_shape(1)
    end if
 
    ! allocate memory for the array
    call mem_allocate(int1d, nvals, idt%mf6varname, mf6_input%mempath)
 
    ! read keyword
    keyword = ''
    call parser%GetStringCaps(keyword)
 
    ! check for "NETCDF" and "LAYERED"
    if (keyword == 'NETCDF') then
      call netcdf_read_array(int1d, mshape, idt, mf6_input, nc_vars, &
                             input_fname, iout)
    else if (keyword == 'LAYERED' .and. idt%layered) then
      call get_layered_shape(mshape, nlay, layer_shape)
      call read_int1d_layered(parser, int1d, idt%mf6varname, nlay, layer_shape)
    else
      call read_int1d(parser, int1d, idt%mf6varname)
    end if
 
    ! log information on the loaded array to the list file
    call idm_log_var(int1d, idt%tagname, mf6_input%mempath, iout)
 
    ! create export file for griddata parameters if optioned
    if (export) then
      if (idt%blockname == 'GRIDDATA') then
        call idm_export(int1d, idt%tagname, mf6_input%mempath, idt%shape, iout)
      end if
    end if
 
    return

Here is the call graph for this function:

Here is the caller graph for this function:

load_integer2d_type()

subroutine loadmf6filemodule::load_integer2d_type	(	type(blockparsertype), intent(inout)	parser,
		type(inputparamdefinitiontype), intent(in)	idt,
		type(modflowinputtype), intent(in)	mf6_input,
		integer(i4b), dimension(:), intent(in), pointer, contiguous	mshape,
		logical(lgp), intent(in)	export,
		type(ncpackagevarstype), intent(in), pointer	nc_vars,
		character(len=*), intent(in)	input_fname,
		integer(i4b), intent(in)	iout
	)

Parameters

[in,out]	parser	block parser
[in]	idt	input data type object describing this record
[in]	mf6_input	description of input
[in]	mshape	model shape
[in]	export	export to ascii layer files
[in]	input_fname	ascii input file name
[in]	iout	unit number for output

Definition at line 859 of file LoadMf6File.f90.

    use sourcecommonmodule, only: get_shape_from_string, get_layered_shape
    use loadncinputmodule, only: netcdf_read_array
    type(BlockParserType), intent(inout) :: parser !< block parser
    type(InputParamDefinitionType), intent(in) :: idt !< input data type object describing this record
    type(ModflowInputType), intent(in) :: mf6_input !< description of input
    integer(I4B), dimension(:), contiguous, pointer, intent(in) :: mshape !< model shape
    logical(LGP), intent(in) :: export !< export to ascii layer files
    type(NCPackageVarsType), pointer, intent(in) :: nc_vars
    character(len=*), intent(in) :: input_fname !< ascii input file name
    integer(I4B), intent(in) :: iout !< unit number for output
    integer(I4B), dimension(:, :), pointer, contiguous :: int2d
    integer(I4B) :: nlay
    integer(I4B) :: nsize1, nsize2
    integer(I4B), dimension(:), allocatable :: array_shape
    integer(I4B), dimension(:), allocatable :: layer_shape
    character(len=LINELENGTH) :: keyword
 
    ! determine the array shape from the input data definition (idt%shape),
    ! which looks like "NCOL, NROW, NLAY"
    call get_shape_from_string(idt%shape, array_shape, mf6_input%mempath)
    nsize1 = array_shape(1)
    nsize2 = array_shape(2)
 
    ! create a new 3d memory managed variable
    call mem_allocate(int2d, nsize1, nsize2, idt%mf6varname, mf6_input%mempath)
 
    ! read keyword
    keyword = ''
    call parser%GetStringCaps(keyword)
 
    ! check for "NETCDF" and "LAYERED"
    if (keyword == 'NETCDF') then
      call netcdf_read_array(int2d, mshape, idt, mf6_input, nc_vars, &
                             input_fname, iout)
    else if (keyword == 'LAYERED' .and. idt%layered) then
      call get_layered_shape(mshape, nlay, layer_shape)
      call read_int2d_layered(parser, int2d, idt%mf6varname, nlay, layer_shape)
    else
      call read_int2d(parser, int2d, idt%mf6varname)
    end if
 
    ! log information on the loaded array to the list file
    call idm_log_var(int2d, idt%tagname, mf6_input%mempath, iout)
 
    ! create export file for griddata parameters if optioned
    if (export) then
      if (idt%blockname == 'GRIDDATA') then
        call idm_export(int2d, idt%tagname, mf6_input%mempath, idt%shape, iout)
      end if
    end if
 
    return

Here is the call graph for this function:

Here is the caller graph for this function:

load_integer3d_type()

subroutine loadmf6filemodule::load_integer3d_type	(	type(blockparsertype), intent(inout)	parser,
		type(inputparamdefinitiontype), intent(in)	idt,
		type(modflowinputtype), intent(in)	mf6_input,
		integer(i4b), dimension(:), intent(in), pointer, contiguous	mshape,
		logical(lgp), intent(in)	export,
		type(ncpackagevarstype), intent(in), pointer	nc_vars,
		character(len=*), intent(in)	input_fname,
		integer(i4b), intent(in)	iout
	)

Parameters

[in,out]	parser	block parser
[in]	idt	input data type object describing this record
[in]	mf6_input	description of input
[in]	mshape	model shape
[in]	export	export to ascii layer files
[in]	input_fname	ascii input file name
[in]	iout	unit number for output

Definition at line 917 of file LoadMf6File.f90.

    use sourcecommonmodule, only: get_shape_from_string, get_layered_shape
    use loadncinputmodule, only: netcdf_read_array
    type(BlockParserType), intent(inout) :: parser !< block parser
    type(InputParamDefinitionType), intent(in) :: idt !< input data type object describing this record
    type(ModflowInputType), intent(in) :: mf6_input !< description of input
    integer(I4B), dimension(:), contiguous, pointer, intent(in) :: mshape !< model shape
    logical(LGP), intent(in) :: export !< export to ascii layer files
    type(NCPackageVarsType), pointer, intent(in) :: nc_vars
    character(len=*), intent(in) :: input_fname !< ascii input file name
    integer(I4B), intent(in) :: iout !< unit number for output
    integer(I4B), dimension(:, :, :), pointer, contiguous :: int3d
    integer(I4B) :: nlay
    integer(I4B) :: nsize1, nsize2, nsize3
    integer(I4B), dimension(:), allocatable :: array_shape
    integer(I4B), dimension(:), allocatable :: layer_shape
    integer(I4B), dimension(:), pointer, contiguous :: int1d_ptr
    character(len=LINELENGTH) :: keyword
 
    ! determine the array shape from the input data definition (idt%shape),
    ! which looks like "NCOL, NROW, NLAY"
    call get_shape_from_string(idt%shape, array_shape, mf6_input%mempath)
    nsize1 = array_shape(1)
    nsize2 = array_shape(2)
    nsize3 = array_shape(3)
 
    ! create a new 3d memory managed variable
    call mem_allocate(int3d, nsize1, nsize2, nsize3, idt%mf6varname, &
                      mf6_input%mempath)
 
    ! read keyword
    keyword = ''
    call parser%GetStringCaps(keyword)
 
    ! check for "NETCDF" and "LAYERED"
    if (keyword == 'NETCDF') then
      call netcdf_read_array(int3d, mshape, idt, mf6_input, nc_vars, &
                             input_fname, iout)
    else if (keyword == 'LAYERED' .and. idt%layered) then
      call get_layered_shape(mshape, nlay, layer_shape)
      call read_int3d_layered(parser, int3d, idt%mf6varname, nlay, &
                              layer_shape)
    else
      int1d_ptr(1:nsize1 * nsize2 * nsize3) => int3d(:, :, :)
      call read_int1d(parser, int1d_ptr, idt%mf6varname)
    end if
 
    ! log information on the loaded array to the list file
    call idm_log_var(int3d, idt%tagname, mf6_input%mempath, iout)
 
    ! create export file for griddata parameters if optioned
    if (export) then
      if (idt%blockname == 'GRIDDATA') then
        call idm_export(int3d, idt%tagname, mf6_input%mempath, idt%shape, iout)
      end if
    end if
 
    return

Here is the call graph for this function:

Here is the caller graph for this function:

load_integer_type()

subroutine loadmf6filemodule::load_integer_type	(	type(blockparsertype), intent(inout)	parser,
		type(inputparamdefinitiontype), intent(in)	idt,
		character(len=*), intent(in)	memoryPath,
		integer(i4b), intent(in)	iout
	)

Parameters

[in,out]	parser	block parser
[in]	idt	input data type object describing this record
[in]	memorypath	memorypath to put loaded information
[in]	iout	unit number for output

Definition at line 784 of file LoadMf6File.f90.

    type(BlockParserType), intent(inout) :: parser !< block parser
    type(InputParamDefinitionType), intent(in) :: idt !< input data type object describing this record
    character(len=*), intent(in) :: memoryPath !< memorypath to put loaded information
    integer(I4B), intent(in) :: iout !< unit number for output
    integer(I4B), pointer :: intvar
    call mem_allocate(intvar, idt%mf6varname, memorypath)
    intvar = parser%GetInteger()
    call idm_log_var(intvar, idt%tagname, memorypath, idt%datatype, iout)
    return

Here is the caller graph for this function:

load_io_tag()

subroutine loadmf6filemodule::load_io_tag	(	type(blockparsertype), intent(inout)	parser,
		type(inputparamdefinitiontype), intent(in)	idt,
		character(len=*), intent(in)	memoryPath,
		character(len=*), intent(in)	which,
		integer(i4b), intent(in)	iout
	)

Parameters

[in,out]	parser	block parser
[in]	idt	input data type object describing this record
[in]	memorypath	memorypath to put loaded information
[in]	iout	unit number for output

Definition at line 714 of file LoadMf6File.f90.

    use memorymanagermodule, only: mem_allocate, mem_reallocate, &
                                   mem_setptr, get_isize
    use characterstringmodule, only: characterstringtype
    type(BlockParserType), intent(inout) :: parser !< block parser
    type(InputParamDefinitionType), intent(in) :: idt !< input data type object describing this record
    character(len=*), intent(in) :: memoryPath !< memorypath to put loaded information
    character(len=*), intent(in) :: which
    integer(I4B), intent(in) :: iout !< unit number for output
    character(len=LINELENGTH) :: cstr
    type(CharacterStringType), dimension(:), pointer, contiguous :: charstr1d
    integer(I4B) :: ilen, isize, idx
    ilen = linelength
    if (which == 'FILEIN') then
      call get_isize(idt%mf6varname, memorypath, isize)
      if (isize < 0) then
        call mem_allocate(charstr1d, ilen, 1, idt%mf6varname, memorypath)
        idx = 1
      else
        call mem_setptr(charstr1d, idt%mf6varname, memorypath)
        call mem_reallocate(charstr1d, ilen, isize + 1, idt%mf6varname, &
                            memorypath)
        idx = isize + 1
      end if
      call parser%GetString(cstr, (.not. idt%preserve_case))
      charstr1d(idx) = cstr
    else if (which == 'FILEOUT') then
      call load_string_type(parser, idt, memorypath, iout)
    end if
    return

Here is the call graph for this function:

Here is the caller graph for this function:

load_keyword_type()

subroutine loadmf6filemodule::load_keyword_type	(	type(blockparsertype), intent(inout)	parser,
		type(inputparamdefinitiontype), intent(in)	idt,
		character(len=*), intent(in)	memoryPath,
		integer(i4b), intent(in)	iout
	)

Parameters

[in,out]	parser	block parser
[in]	idt	input data type object describing this record
[in]	memorypath	memorypath to put loaded information
[in]	iout	unit number for output

Definition at line 674 of file LoadMf6File.f90.

    type(BlockParserType), intent(inout) :: parser !< block parser
    type(InputParamDefinitionType), intent(in) :: idt !< input data type object describing this record
    character(len=*), intent(in) :: memoryPath !< memorypath to put loaded information
    integer(I4B), intent(in) :: iout !< unit number for output
    integer(I4B), pointer :: intvar
    call mem_allocate(intvar, idt%mf6varname, memorypath)
    intvar = 1
    call idm_log_var(intvar, idt%tagname, memorypath, idt%datatype, iout)
    return

Here is the caller graph for this function:

load_string_type()

subroutine loadmf6filemodule::load_string_type ( type(blockparsertype), intent(inout)  parser, type(inputparamdefinitiontype), intent(in)  idt, character(len=*), intent(in)  memoryPath, integer(i4b), intent(in)  iout  )

private

Parameters

[in,out]	parser	block parser
[in]	idt	input data type object describing this record
[in]	memorypath	memorypath to put loaded information
[in]	iout	unit number for output

Definition at line 688 of file LoadMf6File.f90.

    use constantsmodule, only: lenbigline
    type(BlockParserType), intent(inout) :: parser !< block parser
    type(InputParamDefinitionType), intent(in) :: idt !< input data type object describing this record
    character(len=*), intent(in) :: memoryPath !< memorypath to put loaded information
    integer(I4B), intent(in) :: iout !< unit number for output
    character(len=LINELENGTH), pointer :: cstr
    character(len=LENBIGLINE), pointer :: bigcstr
    integer(I4B) :: ilen
    select case (idt%shape)
    case ('LENBIGLINE')
      ilen = lenbigline
      call mem_allocate(bigcstr, ilen, idt%mf6varname, memorypath)
      call parser%GetString(bigcstr, (.not. idt%preserve_case))
      call idm_log_var(bigcstr, idt%tagname, memorypath, iout)
    case default
      ilen = linelength
      call mem_allocate(cstr, ilen, idt%mf6varname, memorypath)
      call parser%GetString(cstr, (.not. idt%preserve_case))
      call idm_log_var(cstr, idt%tagname, memorypath, iout)
    end select
    return

Here is the caller graph for this function:

parse_block()

recursive subroutine loadmf6filemodule::parse_block	(	class(loadmf6filetype)	this,
		integer(i4b), intent(in)	iblk,
		logical(lgp), intent(in)	recursive_call
	)

Parameters

[in]

recursive_call

true if recursive call

Definition at line 276 of file LoadMf6File.f90.

    ! -- modules
    use memorytypemodule, only: memorytype
    use memorymanagermodule, only: get_from_memorystore
    ! -- dummy
    class(LoadMf6FileType) :: this
    integer(I4B), intent(in) :: iblk
    logical(LGP), intent(in) :: recursive_call !< true if recursive call
    ! -- local
    logical(LGP) :: isblockfound
    logical(LGP) :: endOfBlock
    logical(LGP) :: supportOpenClose
    integer(I4B) :: ierr
    logical(LGP) :: found, required
    type(MemoryType), pointer :: mt
    !
    ! -- disu vertices/cell2d blocks are contingent on NVERT dimension
    if (this%mf6_input%pkgtype == 'DISU6' .or. &
        this%mf6_input%pkgtype == 'DISV1D6' .or. &
        this%mf6_input%pkgtype == 'DISV2D6') then
      if (this%mf6_input%block_dfns(iblk)%blockname == 'VERTICES' .or. &
          this%mf6_input%block_dfns(iblk)%blockname == 'CELL2D') then
        call get_from_memorystore('NVERT', this%mf6_input%mempath, mt, found, &
                                  .false.)
        if (.not. found) return
        if (mt%intsclr == 0) return
      end if
    end if
    !
    ! -- block open/close support
    supportopenclose = (this%mf6_input%block_dfns(iblk)%blockname /= 'GRIDDATA')
    !
    ! -- parser search for block
    required = this%mf6_input%block_dfns(iblk)%required .and. .not. recursive_call
    call this%parser%GetBlock(this%mf6_input%block_dfns(iblk)%blockname, &
                              isblockfound, ierr, &
                              supportopenclose=supportopenclose, &
                              blockrequired=required)
    !
    ! -- process block
    if (isblockfound) then
      if (this%mf6_input%block_dfns(iblk)%aggregate) then
        !
        ! -- process block recarray type, set of variable 1d/2d types
        call this%parse_structarray_block(iblk)
        !
      else
        do
          ! process each line in block
          call this%parser%GetNextLine(endofblock)
          if (endofblock) exit
          !
          ! -- process line as tag(s)
          call this%parse_tag(iblk, .false.)
          !
        end do
      end if
    end if
    !
    ! -- recurse if block is reloadable and was just read
    if (this%mf6_input%block_dfns(iblk)%block_variable) then
      if (isblockfound) then
        call this%parse_block(iblk, .true.)
      end if
    end if

Here is the call graph for this function:

parse_io_tag()

subroutine loadmf6filemodule::parse_io_tag	(	class(loadmf6filetype)	this,
		integer(i4b), intent(in)	iblk,
		character(len=*), intent(in)	pkgtype,
		character(len=*), intent(in)	which,
		character(len=*), intent(in)	tag
	)

Definition at line 343 of file LoadMf6File.f90.

    ! -- modules
    ! -- dummy
    class(LoadMf6FileType) :: this
    integer(I4B), intent(in) :: iblk
    character(len=*), intent(in) :: pkgtype
    character(len=*), intent(in) :: which
    character(len=*), intent(in) :: tag
    ! -- local
    type(InputParamDefinitionType), pointer :: idt !< input data type object describing this record
    !
    ! -- matches, read and load file name
    idt => &
      get_param_definition_type(this%mf6_input%param_dfns, &
                                this%mf6_input%component_type, &
                                this%mf6_input%subcomponent_type, &
                                this%mf6_input%block_dfns(iblk)%blockname, &
                                tag, this%filename)
    !
    ! -- load io tag
    call load_io_tag(this%parser, idt, this%mf6_input%mempath, which, this%iout)

Here is the call graph for this function:

parse_keyword_tag()

subroutine loadmf6filemodule::parse_keyword_tag ( class(loadmf6filetype)  this, integer(i4b), intent(in)  iblk, character(len=linelength), intent(in)  tag, type(inputparamdefinitiontype), intent(in), pointer  idt  )

private

Definition at line 366 of file LoadMf6File.f90.

    ! -- modules
    use definitionselectmodule, only: split_record_definition
    ! -- dummy
    class(LoadMf6FileType) :: this
    integer(I4B), intent(in) :: iblk
    character(len=LINELENGTH), intent(in) :: tag
    type(InputParamDefinitionType), pointer, intent(in) :: idt
    ! -- local
    character(len=40), dimension(:), allocatable :: words
    integer(I4B) :: nwords
    character(len=LINELENGTH) :: io_tag
    logical(LGP) :: found
    !
    ! -- initialization
    found = .false.
    !
    ! -- if in record tag check and load if input/output file
    if (idt%in_record) then
      !
      ! -- get tokens in matching definition
      call split_record_definition(this%mf6_input%param_dfns, &
                                   this%mf6_input%component_type, &
                                   this%mf6_input%subcomponent_type, &
                                   tag, nwords, words)
      !
      ! -- a filein/fileout record tag definition has 4 tokens
      if (nwords == 4) then
        !
        ! -- verify third definition token is FILEIN/FILEOUT
        if (words(3) == 'FILEIN' .or. words(3) == 'FILEOUT') then
          !
          ! -- read 3rd token
          call this%parser%GetStringCaps(io_tag)
          !
          ! -- check if 3rd token matches definition
          if (io_tag == words(3)) then
            call this%parse_io_tag(iblk, words(2), words(3), words(4))
            found = .true.
          else
            errmsg = 'Expected "'//trim(words(3))//'" following keyword "'// &
                     trim(tag)//'" but instead found "'//trim(io_tag)//'"'
            call store_error(errmsg)
            call this%parser%StoreErrorUnit()
          end if
          !
        end if
      end if
      !
      ! -- deallocate words
      if (allocated(words)) deallocate (words)
    end if
    !
    if (.not. found) then
      ! -- load standard keyword tag
      call load_keyword_type(this%parser, idt, this%mf6_input%mempath, this%iout)
      !
      ! -- check/set as dev option
      if (idt%tagname(1:4) == 'DEV_' .and. &
          this%mf6_input%block_dfns(iblk)%blockname == 'OPTIONS') then
        call this%parser%DevOpt()
      end if
    end if

Here is the call graph for this function:

parse_structarray_block()

subroutine loadmf6filemodule::parse_structarray_block ( class(loadmf6filetype)  this, integer(i4b), intent(in)  iblk  )

private

A structarray is similar to a numpy recarray. It it used to load a list of data in which each column in the list may be a different type. Each column in the list is stored as a 1d vector.

Definition at line 557 of file LoadMf6File.f90.

    ! -- modules
    use structarraymodule, only: structarraytype, constructstructarray
    use dynamicpackageparamsmodule, only: dynamicpackageparamstype
    ! -- dummy
    class(LoadMf6FileType) :: this
    integer(I4B), intent(in) :: iblk
    ! -- local
    type(DynamicPackageParamsType) :: block_params
    type(InputParamDefinitionType), pointer :: idt !< input data type object describing this record
    type(InputParamDefinitionType), target :: blockvar_idt
    integer(I4B) :: blocknum
    integer(I4B), pointer :: nrow
    integer(I4B) :: nrows, nrowsread
    integer(I4B) :: ibinary, oc_inunit
    integer(I4B) :: icol, iparam
    integer(I4B) :: ncol
    !
    ! -- initialize package params object
    call block_params%init(this%mf6_input, &
                           this%mf6_input%block_dfns(iblk)%blockname, &
                           this%readasarrays, this%iauxiliary, this%inamedbound)
    !
    ! -- set input definition for this block
    idt => &
      get_aggregate_definition_type(this%mf6_input%aggregate_dfns, &
                                    this%mf6_input%component_type, &
                                    this%mf6_input%subcomponent_type, &
                                    this%mf6_input%block_dfns(iblk)%blockname)
    !
    ! -- if block is reloadable read the block number
    if (this%mf6_input%block_dfns(iblk)%block_variable) then
      blocknum = this%parser%GetInteger()
    else
      blocknum = 0
    end if
    !
    ! -- set ncol
    ncol = block_params%nparam
    !
    ! -- add col if block is reloadable
    if (blocknum > 0) ncol = ncol + 1
    !
    ! -- use shape to set the max num of rows
    if (idt%shape /= '') then
      call mem_setptr(nrow, idt%shape, this%mf6_input%mempath)
      nrows = nrow
    else
      nrows = -1
    end if
    !
    ! -- create a structured array
    this%structarray => constructstructarray(this%mf6_input, ncol, nrows, &
                                             blocknum, this%mf6_input%mempath, &
                                             this%mf6_input%component_mempath)
    !
    ! -- create structarray vectors for each column
    do icol = 1, ncol
      !
      ! -- if block is reloadable, block number is first column
      if (blocknum > 0) then
        if (icol == 1) then
          !
          blockvar_idt = this%block_index_dfn(iblk)
          idt => blockvar_idt
          !
          call this%structarray%mem_create_vector(icol, idt)
          !
          ! -- continue as this column managed by internally SA object
          cycle
        end if
        !
        ! -- set indexes (where first column is blocknum)
        iparam = icol - 1
      else
        !
        ! -- set indexes (no blocknum column)
        iparam = icol
      end if
      !
      ! -- set pointer to input definition for this 1d vector
      idt => &
        get_param_definition_type(this%mf6_input%param_dfns, &
                                  this%mf6_input%component_type, &
                                  this%mf6_input%subcomponent_type, &
                                  this%mf6_input%block_dfns(iblk)%blockname, &
                                  block_params%params(iparam), this%filename)
      !
      ! -- allocate variable in memory manager
      call this%structarray%mem_create_vector(icol, idt)
    end do
    !
    ! -- read the block control record
    ibinary = read_control_record(this%parser, oc_inunit, this%iout)
    !
    if (ibinary == 1) then
      !
      ! -- read from binary
      nrowsread = this%structarray%read_from_binary(oc_inunit, this%iout)
      !
      call this%parser%terminateblock()
      !
      close (oc_inunit)
      !
    else
      !
      ! -- read from ascii
      nrowsread = this%structarray%read_from_parser(this%parser, this%ts_active, &
                                                    this%iout)
    end if
    !
    ! -- clean up
    call block_params%destroy()

Here is the call graph for this function:

parse_tag()

recursive subroutine loadmf6filemodule::parse_tag	(	class(loadmf6filetype)	this,
		integer(i4b), intent(in)	iblk,
		logical(lgp), intent(in)	recursive_call
	)

Load an individual input record into the memory manager. Allow for recursive calls in the case that multiple tags are on a single line.

Parameters

[in]

recursive_call

true if recursive call

Definition at line 438 of file LoadMf6File.f90.

    ! -- modules
    use arrayhandlersmodule, only: expandarray
    ! -- dummy
    class(LoadMf6FileType) :: this
    integer(I4B), intent(in) :: iblk
    logical(LGP), intent(in) :: recursive_call !< true if recursive call
    ! -- local
    character(len=LINELENGTH) :: tag
    type(InputParamDefinitionType), pointer :: idt !< input data type object describing this record
    !
    ! -- read tag name
    call this%parser%GetStringCaps(tag)
    if (recursive_call) then
      if (tag == '') then
        ! no data on line so return
        return
      end if
    end if
    !
    ! -- find keyword in input definition
    idt => get_param_definition_type(this%mf6_input%param_dfns, &
                                     this%mf6_input%component_type, &
                                     this%mf6_input%subcomponent_type, &
                                     this%mf6_input%block_dfns(iblk)%blockname, &
                                     tag, this%filename)
    !
    ! -- allocate and load data type
    select case (idt%datatype)
    case ('KEYWORD')
      call this%parse_keyword_tag(iblk, tag, idt)
    case ('STRING')
      if (idt%shape == 'NAUX') then
        call load_auxvar_names(this%parser, idt, this%mf6_input%mempath, &
                               this%iout)
      else
        call load_string_type(this%parser, idt, this%mf6_input%mempath, this%iout)
      end if
    case ('INTEGER')
      call load_integer_type(this%parser, idt, this%mf6_input%mempath, this%iout)
    case ('INTEGER1D')
      call load_integer1d_type(this%parser, idt, this%mf6_input, this%mshape, &
                               this%export, this%nc_vars, this%filename, &
                               this%iout)
    case ('INTEGER2D')
      call load_integer2d_type(this%parser, idt, this%mf6_input, this%mshape, &
                               this%export, this%nc_vars, this%filename, &
                               this%iout)
    case ('INTEGER3D')
      call load_integer3d_type(this%parser, idt, this%mf6_input, this%mshape, &
                               this%export, this%nc_vars, this%filename, &
                               this%iout)
    case ('DOUBLE')
      call load_double_type(this%parser, idt, this%mf6_input%mempath, this%iout)
    case ('DOUBLE1D')
      call load_double1d_type(this%parser, idt, this%mf6_input, this%mshape, &
                              this%export, this%nc_vars, this%filename, this%iout)
    case ('DOUBLE2D')
      call load_double2d_type(this%parser, idt, this%mf6_input, this%mshape, &
                              this%export, this%nc_vars, this%filename, this%iout)
    case ('DOUBLE3D')
      call load_double3d_type(this%parser, idt, this%mf6_input, this%mshape, &
                              this%export, this%nc_vars, this%filename, this%iout)
    case default
      write (errmsg, '(a,a)') 'Failure reading data for tag: ', trim(tag)
      call store_error(errmsg)
      call this%parser%StoreErrorUnit()
    end select
    !
    ! -- continue line if in same record
    if (idt%in_record) then
      !
      ! recursively call parse tag again to read rest of line
      call this%parse_tag(iblk, .true.)
    end if
    !
    !
    call expandarray(this%block_tags)
    this%block_tags(size(this%block_tags)) = trim(idt%tagname)

Here is the call graph for this function:

read_control_record()

integer(i4b) function, public loadmf6filemodule::read_control_record	(	type(blockparsertype), intent(inout)	parser,
		integer(i4b), intent(inout)	oc_inunit,
		integer(i4b), intent(in)	iout
	)

Definition at line 1173 of file LoadMf6File.f90.

    ! -- modules
    use simmodule, only: store_error_unit
    use inputoutputmodule, only: urword
    use inputoutputmodule, only: openfile
    use openspecmodule, only: form, access
    use constantsmodule, only: linelength
    use blockparsermodule, only: blockparsertype
    ! -- dummy
    type(BlockParserType), intent(inout) :: parser
    integer(I4B), intent(inout) :: oc_inunit
    integer(I4B), intent(in) :: iout
    ! -- return
    integer(I4B) :: ibinary
    ! -- local
    integer(I4B) :: lloc, istart, istop, idum, inunit, itmp, ierr
    integer(I4B) :: nunopn = 99
    character(len=:), allocatable :: line
    character(len=LINELENGTH) :: fname
    logical :: exists
    real(DP) :: r
    ! -- formats
    character(len=*), parameter :: fmtocne = &
      &"('Specified OPEN/CLOSE file ',(A),' does not exist')"
    character(len=*), parameter :: fmtobf = &
      &"(1X,/1X,'OPENING BINARY FILE ON UNIT ',I0,':',/1X,A)"
    !
    ! -- initialize oc_inunit and ibinary
    oc_inunit = 0
    ibinary = 0
    !
    inunit = parser%getunit()
    !
    ! -- Read to the first non-commented line
    lloc = 1
    call parser%line_reader%rdcom(inunit, iout, line, ierr)
    call urword(line, lloc, istart, istop, 1, idum, r, iout, inunit)
    !
    if (line(istart:istop) == 'OPEN/CLOSE') then
      !
      ! -- get filename
      call urword(line, lloc, istart, istop, 0, idum, r, &
                  iout, inunit)
      !
      fname = line(istart:istop)
      !
      ! -- check to see if file OPEN/CLOSE file exists
      inquire (file=fname, exist=exists)
      !
      if (.not. exists) then
        write (errmsg, fmtocne) line(istart:istop)
        call store_error(errmsg)
        call store_error('Specified OPEN/CLOSE file does not exist')
        call store_error_unit(inunit)
      end if
      !
      ! -- Check for (BINARY) keyword
      call urword(line, lloc, istart, istop, 1, idum, r, &
                  iout, inunit)
      !
      if (line(istart:istop) == '(BINARY)') ibinary = 1
      !
      ! -- Open the file depending on ibinary flag
      if (ibinary == 1) then
        oc_inunit = nunopn
        itmp = iout
        !
        if (iout > 0) then
          itmp = 0
          write (iout, fmtobf) oc_inunit, trim(adjustl(fname))
        end if
        !
        call openfile(oc_inunit, itmp, fname, 'OPEN/CLOSE', &
                      fmtarg_opt=form, accarg_opt=access)
      end if
    end if
    !
    if (ibinary == 0) then
      call parser%line_reader%bkspc(parser%getunit())
    end if

Here is the call graph for this function:

Here is the caller graph for this function: