DynamicPackageParams.f90

Go to the documentation of this file.

!> @brief This module contains the DynamicPackageParamsModule
!!
!!
!<
module dynamicpackageparamsmodule
 
  use kindmodule, only: dp, i4b, lgp
  use constantsmodule, only: linelength, lenboundname, dzero, izero
  use simvariablesmodule, only: errmsg
  use simmodule, only: store_error, store_error_filename
  use memorymanagermodule, only: mem_allocate
  use modflowinputmodule, only: modflowinputtype
  use inputdefinitionmodule, only: inputparamdefinitiontype
  use definitionselectmodule, only: get_param_definition_type, &
                                    get_aggregate_definition_type, &
                                    idt_parse_rectype
  use arrayhandlersmodule, only: expandarray
  use characterstringmodule, only: characterstringtype
 
  implicit none
  private
  public :: dynamicpackageparamstype
  public :: allocate_param_charstr
  public :: allocate_param_int1d, allocate_param_int2d
  public :: allocate_param_dbl1d, allocate_param_dbl2d
 
  !> @brief dynamic parameter filter type
  !!
  !!
  !<
  type :: dynamicpackageparamstype
    character(len=LINELENGTH), dimension(:), allocatable :: params !< in scope param tags
    character(len=LINELENGTH) :: blockname !< name of block
    integer(I4B) :: iauxiliary !< package auxiliary active, 0=inactive, active for values > 0
    integer(I4B) :: inamedbound !< package inamedbound setting
    integer(I4B) :: nparam !< number of in scope params
    type(modflowinputtype) :: mf6_input !< description of input
  contains
    procedure :: init
    procedure :: destroy
    procedure :: set_filtered_list
    procedure :: set_filtered_grid
    procedure :: package_params
  end type dynamicpackageparamstype
 
contains
 
  !> @brief initialize dynamic param filter
  !!
  !<
  subroutine init(this, mf6_input, blockname, readasarrays, iauxiliary, &
                  inamedbound)
    ! -- modules
    ! -- dummy
    class(dynamicpackageparamstype) :: this
    type(modflowinputtype), intent(in) :: mf6_input
    character(len=*) :: blockname
    logical(LGP), intent(in) :: readasarrays
    integer(I4B), intent(in) :: iauxiliary
    integer(I4B), intent(in) :: inamedbound
    !integer(I4B) :: iparam
    ! -- local
    !
    this%mf6_input = mf6_input
    this%blockname = blockname
    this%nparam = 0
    this%iauxiliary = iauxiliary
    this%inamedbound = inamedbound
    !
    ! -- determine in scope input params
    if (readasarrays) then
      call this%set_filtered_grid()
    else
      call this%set_filtered_list()
    end if
  end subroutine init
 
  !> @brief destroy
  !!
  !<
  subroutine destroy(this)
    ! -- modules
    ! -- dummy
    class(dynamicpackageparamstype) :: this
    !
    ! -- deallocate
    if (allocated(this%params)) deallocate (this%params)
  end subroutine destroy
 
  !> @brief array based input dynamic param filter
  !!
  !<
  subroutine set_filtered_grid(this)
    ! -- modules
    ! -- dummy
    class(dynamicpackageparamstype) :: this
    ! -- local
    type(inputparamdefinitiontype), pointer :: idt
    integer(I4B), dimension(:), allocatable :: idt_idxs
    type(characterstringtype), dimension(:), pointer, contiguous :: boundname
    real(DP), dimension(:, :), pointer, contiguous :: auxvar
    integer(I4B) :: keepcnt, iparam
    logical(LGP) :: keep
    !
    ! -- initialize
    keepcnt = 0
    !
    ! -- allocate dfn input params
    do iparam = 1, size(this%mf6_input%param_dfns)
      !
      keep = .true.
      !
      ! -- assign param definition pointer
      idt => this%mf6_input%param_dfns(iparam)
      !
      if (idt%blockname /= this%blockname) then
        keep = .false.
      end if
      !
      if (idt%tagname == 'AUX') then
        if (this%iauxiliary == 0) then
          keep = .false.
          call mem_allocate(auxvar, 0, 0, 'AUXVAR', this%mf6_input%mempath)
        end if
        if (this%inamedbound == 0) then
          call mem_allocate(boundname, lenboundname, 0, 'BOUNDNAME', &
                            this%mf6_input%mempath)
        end if
      end if
      !
      if (keep) then
        keepcnt = keepcnt + 1
        call expandarray(idt_idxs)
        idt_idxs(keepcnt) = iparam
      end if
    end do
    !
    ! -- update nparam
    this%nparam = keepcnt
    !
    ! -- allocate filtcols
    allocate (this%params(this%nparam))
    !
    ! -- set filtcols
    do iparam = 1, this%nparam
      idt => this%mf6_input%param_dfns(idt_idxs(iparam))
      this%params(iparam) = trim(idt%tagname)
    end do
    !
    ! -- cleanup
    deallocate (idt_idxs)
  end subroutine set_filtered_grid
 
  !> @brief create array of in scope list input columns
  !!
  !! Filter the recarray description of list input parameters
  !! to determine which columns are to be read in this run.
  !<
  subroutine set_filtered_list(this)
    ! -- modules
    ! -- dummy
    class(dynamicpackageparamstype) :: this
    ! -- local
    type(inputparamdefinitiontype), pointer :: ra_idt, idt
    character(len=LINELENGTH), dimension(:), allocatable :: ra_cols
    type(characterstringtype), dimension(:), pointer, contiguous :: boundname
    real(DP), dimension(:, :), pointer, contiguous :: auxvar
    integer(I4B) :: ra_ncol, icol, keepcnt
    logical(LGP) :: keep
    !
    ! -- initialize
    keepcnt = 0
    !
    ! -- get aggregate param definition for period block
    ra_idt => &
      get_aggregate_definition_type(this%mf6_input%aggregate_dfns, &
                                    this%mf6_input%component_type, &
                                    this%mf6_input%subcomponent_type, &
                                    this%blockname)
    !
    ! -- split recarray definition
    call idt_parse_rectype(ra_idt, ra_cols, ra_ncol)
    !
    ! -- determine which columns are in scope
    do icol = 1, ra_ncol
      !
      keep = .false.
      !
      ! -- set dfn pointer to recarray parameter
      idt => get_param_definition_type(this%mf6_input%param_dfns, &
                                       this%mf6_input%component_type, &
                                       this%mf6_input%subcomponent_type, &
                                       this%blockname, ra_cols(icol), '')
      !
      if (ra_cols(icol) == 'RECARRAY') then
        ! no-op
      else if (ra_cols(icol) == 'AUX') then
        if (this%iauxiliary > 0) then
          keep = .true.
        else
          call mem_allocate(auxvar, 0, 0, 'AUXVAR', this%mf6_input%mempath)
        end if
      else if (ra_cols(icol) == 'BOUNDNAME') then
        if (this%inamedbound /= 0) then
          keep = .true.
        else
          call mem_allocate(boundname, lenboundname, 0, 'BOUNDNAME', &
                            this%mf6_input%mempath)
        end if
      else
        ! -- determine if the param is scope
        keep = pkg_param_in_scope(this%mf6_input, this%blockname, ra_cols(icol))
      end if
      !
      if (keep) then
        keepcnt = keepcnt + 1
        call expandarray(this%params)
        this%params(keepcnt) = trim(ra_cols(icol))
      end if
    end do
    !
    ! -- update nparam
    this%nparam = keepcnt
    !
    ! -- cleanup
    deallocate (ra_cols)
  end subroutine set_filtered_list
 
  !> @brief allocate and set input array to filtered param set
  !!
  !<
  subroutine package_params(this, params, nparam)
    ! -- modules
    ! -- dummy
    class(dynamicpackageparamstype) :: this
    character(len=LINELENGTH), dimension(:), allocatable, &
      intent(inout) :: params
    integer(I4B), intent(inout) :: nparam
    integer(I4B) :: n
    !
    if (allocated(params)) deallocate (params)
    !
    nparam = this%nparam
    !
    allocate (params(nparam))
    !
    do n = 1, nparam
      params(n) = this%params(n)
    end do
  end subroutine package_params
 
  !> @brief allocate character string type array
  !<
  subroutine allocate_param_charstr(strlen, nrow, varname, mempath)
    integer(I4B), intent(in) :: strlen !< string number of characters
    integer(I4B), intent(in) :: nrow !< integer array number of rows
    character(len=*), intent(in) :: varname !< variable name
    character(len=*), intent(in) :: mempath !< variable mempath
    type(characterstringtype), dimension(:), pointer, &
      contiguous :: charstr1d
    integer(I4B) :: n
    !
    call mem_allocate(charstr1d, strlen, nrow, varname, mempath)
    do n = 1, nrow
      charstr1d(n) = ''
    end do
  end subroutine allocate_param_charstr
 
  !> @brief allocate int1d
  !<
  subroutine allocate_param_int1d(nrow, varname, mempath)
    integer(I4B), intent(in) :: nrow !< integer array number of rows
    character(len=*), intent(in) :: varname !< variable name
    character(len=*), intent(in) :: mempath !< variable mempath
    integer(I4B), dimension(:), pointer, contiguous :: int1d
    integer(I4B) :: n
    !
    call mem_allocate(int1d, nrow, varname, mempath)
    do n = 1, nrow
      int1d(n) = izero
    end do
  end subroutine allocate_param_int1d
 
  !> @brief allocate int2d
  !<
  subroutine allocate_param_int2d(ncol, nrow, varname, mempath)
    integer(I4B), intent(in) :: ncol !< integer array number of cols
    integer(I4B), intent(in) :: nrow !< integer array number of rows
    character(len=*), intent(in) :: varname !< variable name
    character(len=*), intent(in) :: mempath !< variable mempath
    integer(I4B), dimension(:, :), pointer, contiguous :: int2d
    integer(I4B) :: n, m
    !
    call mem_allocate(int2d, ncol, nrow, varname, mempath)
    do m = 1, nrow
      do n = 1, ncol
        int2d(n, m) = izero
      end do
    end do
  end subroutine allocate_param_int2d
 
  !> @brief allocate dbl1d
  !<
  subroutine allocate_param_dbl1d(nrow, varname, mempath)
    integer(I4B), intent(in) :: nrow !< integer array number of rows
    character(len=*), intent(in) :: varname !< variable name
    character(len=*), intent(in) :: mempath !< variable mempath
    real(dp), dimension(:), pointer, contiguous :: dbl1d
    integer(I4B) :: n
    !
    call mem_allocate(dbl1d, nrow, varname, mempath)
    do n = 1, nrow
      dbl1d(n) = dzero
    end do
  end subroutine allocate_param_dbl1d
 
  !> @brief allocate dbl2d
  !<
  subroutine allocate_param_dbl2d(ncol, nrow, varname, mempath)
    integer(I4B), intent(in) :: ncol !< integer array number of cols
    integer(I4B), intent(in) :: nrow !< integer array number of rows
    character(len=*), intent(in) :: varname !< variable name
    character(len=*), intent(in) :: mempath !< variable mempath
    real(dp), dimension(:, :), pointer, contiguous :: dbl2d
    integer(I4B) :: n, m
    !
    call mem_allocate(dbl2d, ncol, nrow, varname, mempath)
    do m = 1, nrow
      do n = 1, ncol
        dbl2d(n, m) = dzero
      end do
    end do
  end subroutine allocate_param_dbl2d
 
  !> @brief determine if input param is in scope for a package
  !!
  !<
  function pkg_param_in_scope(mf6_input, blockname, tagname) result(in_scope)
    ! -- modules
    use memorymanagermodule, only: get_isize, mem_setptr
    ! -- dummy
    type(modflowinputtype), intent(in) :: mf6_input
    character(len=*), intent(in) :: blockname
    character(len=*), intent(in) :: tagname
    ! -- return
    logical(LGP) :: in_scope
    ! -- local
    type(inputparamdefinitiontype), pointer :: idt
    integer(I4B) :: pdim_isize, popt_isize
    integer(I4B), pointer :: pdim
    !
    ! -- initialize
    in_scope = .false.
    !
    idt => get_param_definition_type(mf6_input%param_dfns, &
                                     mf6_input%component_type, &
                                     mf6_input%subcomponent_type, &
                                     blockname, tagname, '')
    !
    if (idt%required) then
      ! -- required params always included
      in_scope = .true.
    else
      !
      ! -- package specific logic to determine if input params to be read
      select case (mf6_input%subcomponent_type)
      case ('EVT')
        !
        if (tagname == 'PXDP' .or. tagname == 'PETM') then
          call get_isize('NSEG', mf6_input%mempath, pdim_isize)
          if (pdim_isize > 0) then
            call mem_setptr(pdim, 'NSEG', mf6_input%mempath)
            if (pdim > 1) then
              in_scope = .true.
            end if
          end if
        else if (tagname == 'PETM0') then
          call get_isize('SURFRATESPEC', mf6_input%mempath, popt_isize)
          if (popt_isize > 0) then
            in_scope = .true.
          end if
        end if
        !
      case ('NAM')
        in_scope = .true.
      case default
        errmsg = 'IDM unimplemented. DynamicPackageParamsType::pkg_param_in_scope &
                 &add case tagname='//trim(idt%tagname)
        call store_error(errmsg, .true.)
        !call store_error_filename(sourcename)
      end select
    end if
  end function pkg_param_in_scope
 
end module dynamicpackageparamsmodule