listloadmodule Module Reference

This module contains the ListLoadModule. More...

Data Types
type	listloadtype
	list input loader for dynamic packages. 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	reset (this)
subroutine	rp (this, parser)
subroutine	destroy (this)
subroutine	ts_link_bnd (this, structvector, ts_strloc)
subroutine	ts_link_aux (this, structvector, ts_strloc)
subroutine	ts_update (this, structarray)
subroutine	ts_link (this, structvector, ts_strloc)
subroutine	create_structarray (this)

Detailed Description

This module contains the routines for reading period block list based input.

Function/Subroutine Documentation

ad()

subroutine listloadmodule::ad ( class(listloadtype), intent(inout)  this )

private

Definition at line 117 of file Mf6FileList.f90.

    class(ListLoadType), intent(inout) :: this
    ! advance timeseries
    call this%tsmanager%ad()

ainit()

subroutine listloadmodule::ainit ( class(listloadtype), 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 52 of file Mf6FileList.f90.

    use inputoutputmodule, only: getunit
    use memorymanagermodule, only: get_isize
    use blockparsermodule, only: blockparsertype
    use loadmf6filemodule, only: loadmf6filetype
    class(ListLoadType), 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 :: ts_fnames
    character(len=LINELENGTH) :: fname
    integer(I4B) :: ts6_size, n
 
    ! init loader
    call this%DynamicPkgLoadType%init(mf6_input, component_name, &
                                      component_input_name, input_name, &
                                      iperblock, iout)
    ! initialize scalars
    this%iboundname = 0
    this%ts_active = 0
 
    ! load static input
    call loader%load(parser, mf6_input, this%nc_vars, this%input_name, iout)
 
    ! create tsmanager
    allocate (this%tsmanager)
    call tsmanager_cr(this%tsmanager, iout)
 
    ! determine if TS6 files were provided in OPTIONS block
    call get_isize('TS6_FILENAME', this%mf6_input%mempath, ts6_size)
    if (ts6_size > 0) then
      this%ts_active = 1
      call mem_setptr(ts_fnames, 'TS6_FILENAME', this%mf6_input%mempath)
      do n = 1, size(ts_fnames)
        fname = ts_fnames(n)
        call this%tsmanager%add_tsfile(fname, getunit())
      end do
    end if
 
    ! initialize package input context
    call this%ctx%init(mf6_input)
 
    ! store in scope SA cols for list input
    call this%ctx%tags(this%param_names, this%nparam, this%input_name)
 
    ! construct and set up the struct array object
    call this%create_structarray()
 
    ! finalize input context setup
    call this%ctx%allocate_arrays()

Here is the call graph for this function:

create_structarray()

subroutine listloadmodule::create_structarray

(

class(listloadtype), intent(inout)

this

)

Definition at line 292 of file Mf6FileList.f90.

    use inputdefinitionmodule, only: inputparamdefinitiontype
    use definitionselectmodule, only: get_param_definition_type
    class(ListLoadType), intent(inout) :: this
    type(InputParamDefinitionType), pointer :: idt
    integer(I4B) :: icol
 
    ! construct and set up the struct array object
    this%structarray => constructstructarray(this%mf6_input, this%nparam, &
                                             this%ctx%maxbound, 0, &
                                             this%mf6_input%mempath, &
                                             this%mf6_input%component_mempath)
    ! set up struct array
    do icol = 1, this%nparam
      idt => get_param_definition_type(this%mf6_input%param_dfns, &
                                       this%mf6_input%component_type, &
                                       this%mf6_input%subcomponent_type, &
                                       'PERIOD', &
                                       this%param_names(icol), this%input_name)
      ! allocate variable in memory manager
      call this%structarray%mem_create_vector(icol, idt)
      ! store boundname index when found
      if (idt%mf6varname == 'BOUNDNAME') this%iboundname = icol
    end do

Here is the call graph for this function:

destroy()

subroutine listloadmodule::destroy

(

class(listloadtype), intent(inout)

this

)

Definition at line 168 of file Mf6FileList.f90.

    class(ListLoadType), intent(inout) :: this
    !
    ! deallocate tsmanager
    call this%tsmanager%da()
    deallocate (this%tsmanager)
    nullify (this%tsmanager)
    !
    ! deallocate StructArray
    call destructstructarray(this%structarray)
    call this%ctx%destroy()

Here is the call graph for this function:

df()

subroutine listloadmodule::df

(

class(listloadtype), intent(inout)

this

)

Definition at line 111 of file Mf6FileList.f90.

    class(ListLoadType), intent(inout) :: this
    ! define tsmanager
    call this%tsmanager%tsmanager_df()

reset()

subroutine listloadmodule::reset ( class(listloadtype), intent(inout)  this )

private

Definition at line 123 of file Mf6FileList.f90.

    class(ListLoadType), intent(inout) :: this
    ! reset tsmanager
    call this%tsmanager%reset(this%mf6_input%subcomponent_name)

rp()

subroutine listloadmodule::rp ( class(listloadtype), intent(inout)  this, type(blockparsertype), intent(inout), pointer  parser  )

private

Definition at line 129 of file Mf6FileList.f90.

    use blockparsermodule, only: blockparsertype
    use loadmf6filemodule, only: read_control_record
    use structvectormodule, only: structvectortype
    use idmloggermodule, only: idm_log_header, idm_log_close
    class(ListLoadType), intent(inout) :: this
    type(BlockParserType), pointer, intent(inout) :: parser
    integer(I4B) :: ibinary
    integer(I4B) :: oc_inunit
    logical(LGP) :: ts_active
 
    call this%reset()
    ibinary = read_control_record(parser, oc_inunit, this%iout)
 
    ! log lst file header
    call idm_log_header(this%mf6_input%component_name, &
                        this%mf6_input%subcomponent_name, this%iout)
 
    if (ibinary == 1) then
      this%ctx%nbound = &
        this%structarray%read_from_binary(oc_inunit, this%iout)
      call parser%terminateblock()
      close (oc_inunit)
    else
      ts_active = (this%ts_active /= 0)
      this%ctx%nbound = &
        this%structarray%read_from_parser(parser, ts_active, this%iout)
    end if
 
    ! update ts links
    if (this%ts_active /= 0) then
      call this%ts_update(this%structarray)
    end if
 
    ! close logging statement
    call idm_log_close(this%mf6_input%component_name, &
                       this%mf6_input%subcomponent_name, this%iout)

Here is the call graph for this function:

ts_link()

subroutine listloadmodule::ts_link	(	class(listloadtype), intent(inout)	this,
		type(structvectortype), intent(in), pointer	structvector,
		type(tsstringloctype), intent(in), pointer	ts_strloc
	)

Definition at line 278 of file Mf6FileList.f90.

    use structvectormodule, only: structvectortype, tsstringloctype
    class(ListLoadType), intent(inout) :: this
    type(StructVectorType), pointer, intent(in) :: structvector
    type(TSStringLocType), pointer, intent(in) :: ts_strloc
    select case (structvector%memtype)
    case (2) ! dbl1d
      call this%ts_link_bnd(structvector, ts_strloc)
    case (6) ! dbl2d
      call this%ts_link_aux(structvector, ts_strloc)
    case default
    end select

ts_link_aux()

subroutine listloadmodule::ts_link_aux	(	class(listloadtype), intent(inout)	this,
		type(structvectortype), intent(in), pointer	structvector,
		type(tsstringloctype), intent(in), pointer	ts_strloc
	)

Definition at line 216 of file Mf6FileList.f90.

    use timeserieslinkmodule, only: timeserieslinktype
    use timeseriesmanagermodule, only: read_value_or_time_series
    use structvectormodule, only: structvectortype, tsstringloctype
    class(ListLoadType), intent(inout) :: this
    type(StructVectorType), pointer, intent(in) :: structvector
    type(TSStringLocType), pointer, intent(in) :: ts_strloc
    real(DP), pointer :: bndElem
    type(TimeSeriesLinkType), pointer :: tsLinkAux
    type(StructVectorType), pointer :: sv_bound
    character(len=LENBOUNDNAME) :: boundname
 
    nullify (tslinkaux)
 
    ! set bound element
    bndelem => structvector%dbl2d(ts_strloc%col, ts_strloc%row)
 
    ! set link
    call read_value_or_time_series(ts_strloc%token, ts_strloc%row, &
                                   ts_strloc%structarray_col, bndelem, &
                                   this%mf6_input%subcomponent_name, &
                                   'AUX', this%tsmanager, &
                                   this%ctx%iprpak, tslinkaux)
    if (associated(tslinkaux)) then
      ! set variable name
      tslinkaux%Text = this%ctx%auxname_cst(ts_strloc%col)
      ! set boundname if provided
      if (this%ctx%boundnames > 0) then
        sv_bound => this%structarray%get(this%iboundname)
        boundname = sv_bound%charstr1d(ts_strloc%row)
        tslinkaux%BndName = boundname
      end if
    end if

Here is the call graph for this function:

ts_link_bnd()

subroutine listloadmodule::ts_link_bnd ( class(listloadtype), intent(inout)  this, type(structvectortype), intent(in), pointer  structvector, type(tsstringloctype), intent(in), pointer  ts_strloc  )

private

Definition at line 181 of file Mf6FileList.f90.

    use timeserieslinkmodule, only: timeserieslinktype
    use timeseriesmanagermodule, only: read_value_or_time_series
    use structvectormodule, only: structvectortype, tsstringloctype
    class(ListLoadType), intent(inout) :: this
    type(StructVectorType), pointer, intent(in) :: structvector
    type(TSStringLocType), pointer, intent(in) :: ts_strloc
    real(DP), pointer :: bndElem
    type(TimeSeriesLinkType), pointer :: tsLinkBnd
    type(StructVectorType), pointer :: sv_bound
    character(len=LENBOUNDNAME) :: boundname
 
    nullify (tslinkbnd)
 
    ! set bound element
    bndelem => structvector%dbl1d(ts_strloc%row)
 
    ! set link
    call read_value_or_time_series(ts_strloc%token, ts_strloc%row, &
                                   ts_strloc%structarray_col, bndelem, &
                                   this%mf6_input%subcomponent_name, &
                                   'BND', this%tsmanager, &
                                   this%ctx%iprpak, tslinkbnd)
    if (associated(tslinkbnd)) then
      ! set variable name
      tslinkbnd%Text = structvector%idt%mf6varname
      ! set boundname if provided
      if (this%ctx%boundnames > 0) then
        sv_bound => this%structarray%get(this%iboundname)
        boundname = sv_bound%charstr1d(ts_strloc%row)
        tslinkbnd%BndName = boundname
      end if
    end if

Here is the call graph for this function:

ts_update()

subroutine listloadmodule::ts_update	(	class(listloadtype), intent(inout)	this,
		type(structarraytype), intent(inout), pointer	structarray
	)

Definition at line 251 of file Mf6FileList.f90.

    use simmodule, only: count_errors, store_error_filename
    use structvectormodule, only: tsstringloctype
    use structvectormodule, only: structvectortype
    class(ListLoadType), intent(inout) :: this
    type(StructArrayType), pointer, intent(inout) :: structarray
    integer(I4B) :: n, m
    type(TSStringLocType), pointer :: ts_strloc
    type(StructVectorType), pointer :: sv
 
    do m = 1, structarray%count()
      sv => structarray%get(m)
      if (sv%idt%timeseries) then
        do n = 1, sv%ts_strlocs%count()
          ts_strloc => sv%get_ts_strloc(n)
          call this%ts_link(sv, ts_strloc)
        end do
        call sv%clear()
      end if
    end do
 
    ! terminate if errors were detected
    if (count_errors() > 0) then
      call store_error_filename(this%input_name)
    end if

Here is the call graph for this function: