gwf-ic.f90

Go to the documentation of this file.

module gwficmodule
 
  use kindmodule, only: dp, i4b, lgp
  use constantsmodule, only: linelength
  use numericalpackagemodule, only: numericalpackagetype
  use blockparsermodule, only: blockparsertype
  use basedismodule, only: disbasetype
 
  implicit none
  private
  public :: gwfictype
  public :: ic_cr
 
  type, extends(numericalpackagetype) :: gwfictype
 
    real(dp), dimension(:), pointer, contiguous :: strt => null() ! starting head
 
  contains
 
    procedure :: ic_ar
    procedure :: ic_da
    procedure, private :: ic_load
    procedure, private :: allocate_arrays
    procedure, private :: source_griddata
 
  end type gwfictype
 
contains
 
  !> @brief Create a new initial conditions object
  !<
  subroutine ic_cr(ic, name_model, input_mempath, inunit, iout, dis)
    ! -- modules
    use memorymanagerextmodule, only: mem_set_value
    ! -- dummy
    type(gwfictype), pointer :: ic
    character(len=*), intent(in) :: name_model
    character(len=*), intent(in) :: input_mempath
    integer(I4B), intent(in) :: inunit
    integer(I4B), intent(in) :: iout
    class(disbasetype), pointer, intent(in) :: dis
    ! -- formats
    character(len=*), parameter :: fmtic = &
      "(1x, /1x, 'IC -- Initial Conditions Package, Version 8, 3/28/2015', &
      &' input read from mempath: ', A, //)"
    !
    ! -- create IC object
    allocate (ic)
    !
    ! -- create name and memory path
    call ic%set_names(1, name_model, 'IC', 'IC', input_mempath)
    !
    ! -- allocate scalars
    call ic%allocate_scalars()
    !
    ! -- set variables
    ic%inunit = inunit
    ic%iout = iout
    !
    ! -- set pointers
    ic%dis => dis
    !
    ! -- check if pkg is enabled,
    if (inunit > 0) then
      ! print message identifying pkg
      write (ic%iout, fmtic) input_mempath
    end if
  end subroutine ic_cr
 
  !> @brief Load data from IDM into package
  !<
  subroutine ic_load(this)
    ! -- modules
    use basedismodule, only: disbasetype
    ! -- dummy
    class(gwfictype) :: this
    !
    call this%source_griddata()
  end subroutine ic_load
 
  !> @brief Allocate arrays, load from IDM, and assign head
  !<
  subroutine ic_ar(this, x)
    ! -- dummy
    class(gwfictype) :: this
    real(DP), dimension(:), intent(inout) :: x
    ! -- local
    integer(I4B) :: n
    !
    ! -- allocate arrays
    call this%allocate_arrays(this%dis%nodes)
    !
    ! -- load from IDM
    call this%ic_load()
    !
    ! -- assign starting head
    do n = 1, this%dis%nodes
      x(n) = this%strt(n)
    end do
  end subroutine ic_ar
 
  !> @brief Deallocate
  !<
  subroutine ic_da(this)
    ! -- modules
    use memorymanagermodule, only: mem_deallocate
    use memorymanagerextmodule, only: memorystore_remove
    use simvariablesmodule, only: idm_context
    ! -- dummy
    class(gwfictype) :: this
    !
    ! -- deallocate IDM memory
    call memorystore_remove(this%name_model, 'IC', idm_context)
    !
    ! -- deallocate arrays
    call mem_deallocate(this%strt)
    !
    ! -- deallocate parent
    call this%NumericalPackageType%da()
  end subroutine ic_da
 
  !> @brief Allocate arrays
  !<
  subroutine allocate_arrays(this, nodes)
    ! -- modules
    use memorymanagermodule, only: mem_allocate
    ! -- dummy
    class(gwfictype) :: this
    integer(I4B), intent(in) :: nodes
    !
    ! -- Allocate
    call mem_allocate(this%strt, nodes, 'STRT', this%memoryPath)
  end subroutine allocate_arrays
 
  !> @brief Copy grid data from IDM into package
  !<
  subroutine source_griddata(this)
    ! -- modules
    use simmodule, only: store_error, store_error_filename
    use memorymanagerextmodule, only: mem_set_value
    use gwficinputmodule, only: gwficparamfoundtype
    ! -- dummy
    class(gwfictype) :: this
    ! -- local
    character(len=LINELENGTH) :: errmsg
    type(gwficparamfoundtype) :: found
    integer(I4B), dimension(:), pointer, contiguous :: map
    !
    ! -- set map to convert user to reduced node data
    map => null()
    if (this%dis%nodes < this%dis%nodesuser) map => this%dis%nodeuser
    !
    ! -- set values
    call mem_set_value(this%strt, 'STRT', this%input_mempath, map, found%strt)
    !
    ! -- ensure STRT was found
    if (.not. found%strt) then
      write (errmsg, '(a)') 'Error in GRIDDATA block: STRT not found.'
      call store_error(errmsg, terminate=.false.)
      call store_error_filename(this%input_fname)
    else if (this%iout > 0) then
      write (this%iout, '(4x,a)') 'STRT set from input file'
    end if
  end subroutine source_griddata
 
end module gwficmodule