sourceloadmodule Module Reference

This module contains the SourceLoadModule. More...

Functions/Subroutines
class(staticpkgloadbasetype) function, pointer, public	create_input_loader (component_type, subcomponent_type, component_name, subcomponent_name, input_type, input_fname, component_fname, nc_vars)
	factory function to create and setup model package static loader More...
class(staticpkgloadbasetype) function, pointer	package_loader (source_type)
	allocate source model package static loader More...
integer(i4b) function, public	open_source_file (pkgtype, filename, modelfname, iout)
subroutine, public	load_modelnam (mtype, mfname, mname, iout)
subroutine, public	load_simnam ()
subroutine, public	load_simtdis ()
integer(i4b) function, public	remote_model_ndim (mtype, mfname)
subroutine, public	export_cr ()
	create model exports list More...
subroutine, public	export_post_prepare ()
	model exports post prepare step actions More...
subroutine, public	export_post_step ()
	model exports post step actions More...
subroutine, public	export_da ()
	deallocate model export objects and list More...
subroutine, public	nc_close (ncid, nc_fname)
	close an open netcdf file More...
type(ncfilevarstype) function, pointer, public	netcdf_context (modeltype, component_type, modelname, modelfname, iout)
	create model netcdf context More...

Detailed Description

This module contains the routines needed to generate a loader object for an input source and routines that distribute processing to a particular source.

Function/Subroutine Documentation

create_input_loader()

class(staticpkgloadbasetype) function, pointer, public sourceloadmodule::create_input_loader	(	character(len=*), intent(in)	component_type,
		character(len=*), intent(in)	subcomponent_type,
		character(len=*), intent(in)	component_name,
		character(len=*), intent(in)	subcomponent_name,
		character(len=*), intent(in)	input_type,
		character(len=*), intent(in)	input_fname,
		character(len=*), intent(in)	component_fname,
		type(ncfilevarstype), intent(in), optional, pointer	nc_vars
	)

Definition at line 33 of file SourceLoad.F90.

    use sourcecommonmodule, only: package_source_type, idm_subcomponent_name
    use inputloadtypemodule, only: staticpkgloadbasetype
    character(len=*), intent(in) :: component_type
    character(len=*), intent(in) :: subcomponent_type
    character(len=*), intent(in) :: component_name
    character(len=*), intent(in) :: subcomponent_name
    character(len=*), intent(in) :: input_type
    character(len=*), intent(in) :: input_fname
    character(len=*), intent(in) :: component_fname
    type(NCFileVarsType), pointer, optional, intent(in) :: nc_vars
    class(StaticPkgLoadBaseType), pointer :: loader
    type(ModflowInputType) :: mf6_input
    character(len=LENPACKAGENAME) :: source_type
    character(len=LENPACKAGENAME) :: sc_name
 
    ! set subcomponent name
    sc_name = idm_subcomponent_name(component_type, subcomponent_type, &
                                    subcomponent_name)
    ! create description of input
    mf6_input = getmodflowinput(input_type, component_type, subcomponent_type, &
                                component_name, sc_name, input_fname)
    ! set package source
    source_type = package_source_type(input_fname)
 
    ! set source loader for model package
    loader => package_loader(source_type)
 
    ! initialize loader
    call loader%init(mf6_input, component_name, component_fname, input_fname)
 
    ! initialize loader netcdf variables data structure
    if (present(nc_vars)) then
      call nc_vars%create_varlists(component_name, sc_name, loader%nc_vars)
    else
      call loader%nc_vars%init(component_name)
    end if

Here is the call graph for this function:

Here is the caller graph for this function:

export_cr()

subroutine, public sourceloadmodule::export_cr

Definition at line 300 of file SourceLoad.F90.

    use modelexportmodule, only: modelexports_create, nc_export_active
#if defined(__WITH_NETCDF__)
    use ncexportcreatemodule, only: nc_export_create
#endif
    call modelexports_create(iout)
    ! are netcdf exports elected
    if (nc_export_active()) then
#if defined(__WITH_NETCDF__)
      call nc_export_create()
#else
      write (errmsg, '(a)') &
        'Model namefile NETCDF_STUCTURED or NETCDF_MESH2D option configured &
        &but NetCDF libraries are not available.'
      call store_error(errmsg, .true.)
#endif
    end if

Here is the call graph for this function:

Here is the caller graph for this function:

export_da()

subroutine, public sourceloadmodule::export_da

Definition at line 335 of file SourceLoad.F90.

    use modelexportmodule, only: modelexports_destroy
    call modelexports_destroy()

Here is the call graph for this function:

Here is the caller graph for this function:

export_post_prepare()

subroutine, public sourceloadmodule::export_post_prepare

Definition at line 321 of file SourceLoad.F90.

    use modelexportmodule, only: modelexports_post_prepare
    call modelexports_post_prepare()

Here is the call graph for this function:

Here is the caller graph for this function:

export_post_step()

subroutine, public sourceloadmodule::export_post_step

Definition at line 328 of file SourceLoad.F90.

    use modelexportmodule, only: modelexports_post_step
    call modelexports_post_step()

Here is the call graph for this function:

Here is the caller graph for this function:

load_modelnam()

subroutine, public sourceloadmodule::load_modelnam	(	character(len=*), intent(in)	mtype,
		character(len=*), intent(in)	mfname,
		character(len=*), intent(in)	mname,
		integer(i4b), intent(in)	iout
	)

Definition at line 124 of file SourceLoad.F90.

    use simvariablesmodule, only: simfile
    use sourcecommonmodule, only: package_source_type, idm_component_type
    use idmmf6filemodule, only: input_load
    character(len=*), intent(in) :: mtype
    character(len=*), intent(in) :: mfname
    character(len=*), intent(in) :: mname
    integer(I4B), intent(in) :: iout
    type(ModflowInputType) :: mf6_input
    character(len=LENPACKAGENAME) :: source_type
 
    ! set source type
    source_type = package_source_type(mfname)
 
    ! create description of input
    mf6_input = getmodflowinput(mtype, idm_component_type(mtype), 'NAM', &
                                mname, 'NAM', mfname)
    select case (source_type)
    case ('MF6FILE')
      call input_load(mfname, mf6_input, simfile, iout)
    case default
    end select

Here is the call graph for this function:

Here is the caller graph for this function:

load_simnam()

subroutine, public sourceloadmodule::load_simnam

Definition at line 148 of file SourceLoad.F90.

    use simvariablesmodule, only: simfile, iout
    use memorymanagermodule, only: mem_setptr
    use messagemodule, only: write_message
    use idmmf6filemodule, only: input_load
    use sourcecommonmodule, only: filein_fname
    type(ModflowInputType) :: mf6_input, hpc_input
    character(len=LINELENGTH) :: hpc6_filename
    character(len=LINELENGTH) :: line
    logical(LGP) :: lexist
 
    ! load mfsim.nam if it exists
    inquire (file=trim(adjustl(simfile)), exist=lexist)
 
    if (lexist) then
      ! write name of namfile to stdout
      write (line, '(2(1x,a))') 'Using Simulation name file:', &
        trim(adjustl(simfile))
      call write_message(line, skipafter=1)
      ! create description of input
      mf6_input = getmodflowinput('NAM6', 'SIM', 'NAM', 'SIM', 'NAM', simfile)
      ! open namfile and load to input context
      call input_load(simfile, mf6_input, simfile, iout)
      ! load optional HPC configuration file
      if (filein_fname(hpc6_filename, 'HPC6_FILENAME', mf6_input%mempath, &
                       simfile)) then
        hpc_input = getmodflowinput('HPC6', 'UTL', 'HPC', 'UTL', 'HPC')
        call input_load(hpc6_filename, hpc_input, simfile, iout)
      end if
    end if

Here is the call graph for this function:

Here is the caller graph for this function:

load_simtdis()

subroutine, public sourceloadmodule::load_simtdis

Definition at line 180 of file SourceLoad.F90.

    use simvariablesmodule, only: simfile, iout
    use memoryhelpermodule, only: create_mem_path
    use memorymanagermodule, only: mem_setptr
    use simvariablesmodule, only: idm_context
    use sourcecommonmodule, only: package_source_type
    use idmmf6filemodule, only: input_load
    character(len=LENMEMPATH) :: input_mempath
    type(ModflowInputType) :: mf6_input
    character(len=LENPACKAGENAME) :: source_type
    character(len=:), pointer :: tdis6
    logical(LGP) :: lexist
 
    ! set input memory path
    input_mempath = create_mem_path('SIM', 'NAM', idm_context)
 
    ! set pointers to input context timing params
    call mem_setptr(tdis6, 'TDIS6', input_mempath)
 
    ! create timing
    if (tdis6 /= '') then
      ! set source type
      source_type = package_source_type(tdis6)
      select case (source_type)
      case ('MF6FILE')
        inquire (file=trim(adjustl(tdis6)), exist=lexist)
        if (lexist) then
          ! create description of input
          mf6_input = getmodflowinput('TDIS6', 'SIM', 'TDIS', &
                                      'SIM', 'TDIS', simfile)
          ! open namfile and load to input context
          call input_load(tdis6, mf6_input, simfile, iout)
        else
          write (errmsg, '(a)') &
            'Simulation TIMING input file "'//trim(tdis6)// &
            '" does not exist.'
          call store_error(errmsg)
          call store_error_filename(simfile)
        end if
      case default
      end select
    end if

Here is the call graph for this function:

Here is the caller graph for this function:

nc_close()

subroutine, public sourceloadmodule::nc_close	(	integer(i4b), intent(in)	ncid,
		character(len=*), intent(in)	nc_fname
	)

Definition at line 342 of file SourceLoad.F90.

#if defined(__WITH_NETCDF__)
    use netcdfcommonmodule, only: nc_fclose
#endif
    integer(I4B), intent(in) :: ncid
    character(len=*), intent(in) :: nc_fname
    if (ncid > 0) then
#if defined(__WITH_NETCDF__)
      call nc_fclose(ncid, nc_fname)
#endif
    end if

Here is the call graph for this function:

Here is the caller graph for this function:

netcdf_context()

type(ncfilevarstype) function, pointer, public sourceloadmodule::netcdf_context	(	character(len=*), intent(in)	modeltype,
		character(len=*), intent(in)	component_type,
		character(len=*), intent(in)	modelname,
		character(len=*), intent(in)	modelfname,
		integer(i4b), intent(in)	iout
	)

Definition at line 357 of file SourceLoad.F90.

    use memoryhelpermodule, only: create_mem_path
    use simvariablesmodule, only: idm_context
    use sourcecommonmodule, only: filein_fname
#if defined(__WITH_NETCDF__)
    use nccontextbuildmodule, only: open_ncfile, create_netcdf_context
#endif
    character(len=*), intent(in) :: modeltype
    character(len=*), intent(in) :: component_type
    character(len=*), intent(in) :: modelname
    character(len=*), intent(in) :: modelfname
    integer(I4B), intent(in) :: iout
    type(NCFileVarsType), pointer :: nc_vars
    character(len=LENMEMPATH) :: input_mempath
    character(len=LINELENGTH) :: nc_fname
    integer(I4B) :: ncid
 
    ! set input memory path
    input_mempath = create_mem_path(modelname, 'NAM', idm_context)
 
    ! allocate context object
    allocate (nc_vars)
 
    ! check if optional netcdf input file was provided
    if (filein_fname(nc_fname, 'NETCDF_FNAME', input_mempath, modelfname)) then
#if defined(__WITH_NETCDF__)
      ! open nc input file
      ncid = open_ncfile(nc_fname, iout)
      ! read the file and build the context
      call create_netcdf_context(modeltype, modelname, modelfname, &
                                 nc_vars, nc_fname, ncid, iout)
#else
      write (errmsg, '(a)') &
        'Cannot load model packages. NetCDF &
        &keyword specified in input file but &
        &NetCDF libraries are not available.'
      call store_error(errmsg)
      call store_error_filename(modelfname)
#endif
    else
      ncid = 0
      call nc_vars%init(modelname, '', ncid, '')
    end if

Here is the call graph for this function:

Here is the caller graph for this function:

open_source_file()

integer(i4b) function, public sourceloadmodule::open_source_file	(	character(len=*), intent(in)	pkgtype,
		character(len=*), intent(in)	filename,
		character(len=*), intent(in)	modelfname,
		integer(i4b), intent(in)	iout
	)

Definition at line 101 of file SourceLoad.F90.

    use sourcecommonmodule, only: package_source_type
    use idmmf6filemodule, only: open_mf6file
    character(len=*), intent(in) :: pkgtype
    character(len=*), intent(in) :: filename
    character(len=*), intent(in) :: modelfname
    integer(I4B), intent(in) :: iout
    integer(I4B) :: fd
    character(len=LENPACKAGENAME) :: source_type
 
    ! initialize
    fd = 0
 
    ! set source type
    source_type = package_source_type(filename)
    !
    select case (source_type)
    case ('MF6FILE')
      fd = open_mf6file(pkgtype, filename, modelfname, iout)
    case default
    end select

Here is the call graph for this function:

Here is the caller graph for this function:

package_loader()

class(staticpkgloadbasetype) function, pointer sourceloadmodule::package_loader

(

character(len=*), intent(inout)

source_type

)

Definition at line 77 of file SourceLoad.F90.

    use inputloadtypemodule, only: staticpkgloadbasetype
    use idmmf6filemodule, only: mf6filestaticpkgloadtype
    character(len=*), intent(inout) :: source_type
    class(Mf6FileStaticPkgLoadType), pointer :: mf6file_loader
    class(StaticPkgLoadBaseType), pointer :: loader
 
    ! initialize
    nullify (loader)
 
    ! allocate derived object
    select case (source_type)
    case ('MF6FILE')
      allocate (mf6file_loader)
      loader => mf6file_loader
      allocate (loader%nc_vars)
    case default
      write (errmsg, '(a)') &
        'Simulation package input source type "'//trim(source_type)// &
        '" not currently supported.'
      call store_error(errmsg, .true.)
    end select

Here is the call graph for this function:

Here is the caller graph for this function:

remote_model_ndim()

integer(i4b) function, public sourceloadmodule::remote_model_ndim	(	character(len=*), intent(in)	mtype,
		character(len=*), intent(in)	mfname
	)

Definition at line 224 of file SourceLoad.F90.

    use sourcecommonmodule, only: package_source_type
    use constantsmodule, only: linelength
    use inputoutputmodule, only: openfile, getunit
    use blockparsermodule, only: blockparsertype
    character(len=*), intent(in) :: mtype
    character(len=*), intent(in) :: mfname
    integer(I4B) :: ncelldim
    character(len=LENPACKAGENAME) :: source_type
    type(BlockParserType) :: parser
    integer(I4B) :: ierr, inunit
    logical(LGP) :: isfound, endOfBlock
    character(len=LINELENGTH) :: ptype
 
    ! initialize
    ncelldim = 0
 
    ! set source type
    source_type = package_source_type(mfname)
    select case (source_type)
    case ('MF6FILE')
      ! open name file
      inunit = getunit()
      call openfile(inunit, 0, trim(adjustl(mfname)), mtype, &
                    'FORMATTED', 'SEQUENTIAL', 'OLD')
      ! initialize parser
      call parser%Initialize(inunit, 0)
      ! get options block
      call parser%GetBlock('OPTIONS', isfound, ierr, &
                           supportopenclose=.true., blockrequired=.false.)
      ! iterate through options
      if (isfound) then
        do
          call parser%GetNextLine(endofblock)
          if (endofblock) exit
        end do
      end if
      ! get packages block
      call parser%GetBlock('PACKAGES', isfound, ierr, &
                           supportopenclose=.true., blockrequired=.true.)
      if (isfound) then
        ! read through packages
        do
          call parser%GetNextLine(endofblock)
          if (endofblock) exit
          call parser%GetStringCaps(ptype)
 
          select case (ptype)
          case ('DIS6')
            ncelldim = 3
            exit
          case ('DIS2D6')
            ncelldim = 2
            exit
          case ('DISV6')
            ncelldim = 2
            exit
          case ('DISU6')
            ncelldim = 1
            exit
          case default
            write (errmsg, '(a)') &
              'Unknown discretization type "'//trim(ptype)// &
              '" not currently supported.'
            call store_error(errmsg, .true.)
          end select
        end do
      end if
 
      call parser%clear()
    case default
    end select

Here is the call graph for this function:

Here is the caller graph for this function: