rchmodule Module Reference

Data Types
type	rchtype

Functions/Subroutines
subroutine, public	rch_create (packobj, id, ibcnum, inunit, iout, namemodel, pakname, mempath)
	Create a New Recharge Package. More...
subroutine	rch_allocate_scalars (this)
	Allocate scalar members. More...
subroutine	rch_allocate_arrays (this, nodelist, auxvar)
	Allocate package arrays. More...
subroutine	rch_source_options (this)
	Source options specific to RchType. More...
subroutine	log_rch_options (this, found_fixed_cell)
	Log options specific to RchType. More...
subroutine	rch_read_initial_attr (this)
	Part of allocate and read. More...
subroutine	rch_rp (this)
	Read and Prepare. More...
subroutine	set_nodesontop (this)
	Store nodelist in nodesontop. More...
subroutine	rch_cf (this)
	Formulate the HCOF and RHS terms. More...
subroutine	rch_fc (this, rhs, ia, idxglo, matrix_sln)
	Copy rhs and hcof into solution rhs and amat. More...
subroutine	rch_da (this)
	Deallocate memory. More...
subroutine	rch_define_listlabel (this)
	Define the list heading that is written to iout when PRINT_INPUT option is used. More...
logical function	rch_obs_supported (this)
	Overrides BndTypebnd_obs_supported() More...
subroutine	rch_df_obs (this)
	Implements bnd_df_obs. More...
real(dp) function	rch_bound_value (this, col, row)
	Return requested boundary value. More...

Variables
character(len=lenftype)	ftype = 'RCH'
character(len=lenpackagename)	text = ' RCH'
character(len=lenpackagename)	texta = ' RCHA'

Function/Subroutine Documentation

log_rch_options()

subroutine rchmodule::log_rch_options	(	class(rchtype), intent(inout)	this,
		logical(lgp), intent(in)	found_fixed_cell
	)

Definition at line 154 of file gwf-rch.f90.

    implicit none
    ! -- dummy
    class(RchType), intent(inout) :: this
    logical(LGP), intent(in) :: found_fixed_cell
    ! -- formats
    character(len=*), parameter :: fmtfixedcell = &
      &"(4x, 'RECHARGE WILL BE APPLIED TO SPECIFIED CELL.')"
    !
    ! -- log found options
    write (this%iout, '(/1x,a)') 'PROCESSING '//trim(adjustl(this%text)) &
      //' OPTIONS'
    !
    if (found_fixed_cell) then
      write (this%iout, fmtfixedcell)
    end if
    !
    ! -- close logging block
    write (this%iout, '(1x,a)') &
      'END OF '//trim(adjustl(this%text))//' OPTIONS'

rch_allocate_arrays()

subroutine rchmodule::rch_allocate_arrays ( class(rchtype)  this, integer(i4b), dimension(:), optional, pointer, contiguous  nodelist, real(dp), dimension(:, :), optional, pointer, contiguous  auxvar  )

private

Definition at line 107 of file gwf-rch.f90.

    ! -- modules
    use memorymanagermodule, only: mem_setptr, mem_checkin
    ! -- dummy
    class(RchType) :: this
    integer(I4B), dimension(:), pointer, contiguous, optional :: nodelist
    real(DP), dimension(:, :), pointer, contiguous, optional :: auxvar
    !
    ! -- allocate base arrays
    call this%BndExtType%allocate_arrays(nodelist, auxvar)
    !
    ! -- set recharge input context pointer
    call mem_setptr(this%recharge, 'RECHARGE', this%input_mempath)
    !
    ! -- checkin recharge input context pointer
    call mem_checkin(this%recharge, 'RECHARGE', this%memoryPath, &
                     'RECHARGE', this%input_mempath)

rch_allocate_scalars()

subroutine rchmodule::rch_allocate_scalars ( class(rchtype), intent(inout)  this )

private

Definition at line 91 of file gwf-rch.f90.

    ! -- dummy
    class(RchType), intent(inout) :: this
    !
    ! -- allocate base scalars
    call this%BndExtType%allocate_scalars()
    !
    ! -- allocate internal members
    allocate (this%fixed_cell)
    !
    ! -- Set values
    this%fixed_cell = .false.

rch_bound_value()

real(dp) function rchmodule::rch_bound_value ( class(rchtype), intent(inout)  this, integer(i4b), intent(in)  col, integer(i4b), intent(in)  row  )

private

Parameters

[in,out]

this

BndExtType object

Definition at line 401 of file gwf-rch.f90.

    ! -- modules
    use constantsmodule, only: dzero
    ! -- dummy
    class(RchType), intent(inout) :: this !< BndExtType object
    integer(I4B), intent(in) :: col
    integer(I4B), intent(in) :: row
    ! -- result
    real(DP) :: bndval
    !
    select case (col)
    case (1)
      if (this%iauxmultcol > 0) then
        bndval = this%recharge(row) * this%auxvar(this%iauxmultcol, row)
      else
        bndval = this%recharge(row)
      end if
    case default
      errmsg = 'Programming error. RCH bound value requested column '&
               &'outside range of ncolbnd (1).'
      call store_error(errmsg)
      call store_error_filename(this%input_fname)
    end select

Here is the call graph for this function:

rch_cf()

subroutine rchmodule::rch_cf ( class(rchtype)  this )

private

Skip if no recharge. Otherwise, calculate hcof and rhs

Definition at line 243 of file gwf-rch.f90.

    ! -- dummy
    class(rchtype) :: this
    ! -- local
    integer(I4B) :: i, node
    !
    ! -- Return if no recharge
    if (this%nbound == 0) return
    !
    ! -- Calculate hcof and rhs for each recharge entry
    do i = 1, this%nbound
      !
      ! -- Find the node number
      if (this%fixed_cell) then
        node = this%nodelist(i)
      else
        node = this%nodesontop(i)
      end if
      !
      ! -- cycle if nonexistent bound
      if (node <= 0) then
        this%hcof(i) = dzero
        this%rhs(i) = dzero
        cycle
      end if
      !
      ! -- reset nodelist to highest active
      if (.not. this%fixed_cell) then
        if (this%ibound(node) == 0) &
          call this%dis%highest_active(node, this%ibound)
        this%nodelist(i) = node
      end if
      !
      ! -- Set rhs and hcof
      this%hcof(i) = dzero
      if (this%iauxmultcol > 0) then
        this%rhs(i) = -this%recharge(i) * this%dis%get_area(node) * &
                      this%auxvar(this%iauxmultcol, i)
      else
        this%rhs(i) = -this%recharge(i) * this%dis%get_area(node)
      end if
      if (this%ibound(node) <= 0) then
        this%rhs(i) = dzero
        cycle
      end if
      if (this%ibound(node) == iwetlake) then
        this%rhs(i) = dzero
        cycle
      end if
    end do

rch_create()

subroutine, public rchmodule::rch_create	(	class(bndtype), pointer	packobj,
		integer(i4b), intent(in)	id,
		integer(i4b), intent(in)	ibcnum,
		integer(i4b), intent(in)	inunit,
		integer(i4b), intent(in)	iout,
		character(len=*), intent(in)	namemodel,
		character(len=*), intent(in)	pakname,
		character(len=*), intent(in)	mempath
	)

Create new RCH package and point packobj to the new package

Definition at line 54 of file gwf-rch.f90.

    ! -- dummy
    class(BndType), pointer :: packobj
    integer(I4B), intent(in) :: id
    integer(I4B), intent(in) :: ibcnum
    integer(I4B), intent(in) :: inunit
    integer(I4B), intent(in) :: iout
    character(len=*), intent(in) :: namemodel
    character(len=*), intent(in) :: pakname
    character(len=*), intent(in) :: mempath
    ! -- local
    type(rchtype), pointer :: rchobj
    !
    ! -- allocate recharge object and scalar variables
    allocate (rchobj)
    packobj => rchobj
    !
    ! -- create name and memory path
    call packobj%set_names(ibcnum, namemodel, pakname, ftype, mempath)
    packobj%text = text
    !
    ! -- allocate scalars
    call rchobj%rch_allocate_scalars()
    !
    ! -- initialize package
    call packobj%pack_initialize()
    !
    packobj%inunit = inunit
    packobj%iout = iout
    packobj%id = id
    packobj%ibcnum = ibcnum
    packobj%ictMemPath = create_mem_path(namemodel, 'NPF')

Here is the call graph for this function:

Here is the caller graph for this function:

rch_da()

subroutine rchmodule::rch_da ( class(rchtype)  this )

private

Definition at line 325 of file gwf-rch.f90.

    ! -- modules
    use memorymanagermodule, only: mem_deallocate
    ! -- dummy
    class(RchType) :: this
    !
    ! -- Deallocate parent package
    call this%BndExtType%bnd_da()
    !
    ! -- scalars
    deallocate (this%fixed_cell)
    !
    ! -- arrays
    if (associated(this%nodesontop)) deallocate (this%nodesontop)
    call mem_deallocate(this%recharge, 'RECHARGE', this%memoryPath)

rch_define_listlabel()

subroutine rchmodule::rch_define_listlabel

(

class(rchtype), intent(inout)

this

)

Definition at line 345 of file gwf-rch.f90.

    ! -- dummy
    class(RchType), intent(inout) :: this
    !
    ! -- create the header list label
    this%listlabel = trim(this%filtyp)//' NO.'
    if (this%dis%ndim == 3) then
      write (this%listlabel, '(a, a7)') trim(this%listlabel), 'LAYER'
      write (this%listlabel, '(a, a7)') trim(this%listlabel), 'ROW'
      write (this%listlabel, '(a, a7)') trim(this%listlabel), 'COL'
    elseif (this%dis%ndim == 2) then
      write (this%listlabel, '(a, a7)') trim(this%listlabel), 'LAYER'
      write (this%listlabel, '(a, a7)') trim(this%listlabel), 'CELL2D'
    else
      write (this%listlabel, '(a, a7)') trim(this%listlabel), 'NODE'
    end if
    write (this%listlabel, '(a, a16)') trim(this%listlabel), 'RECHARGE'
!    if(this%multindex > 0) &
!      write(this%listlabel, '(a, a16)') trim(this%listlabel), 'MULTIPLIER'
    if (this%inamedbound == 1) then
      write (this%listlabel, '(a, a16)') trim(this%listlabel), 'BOUNDARY NAME'
    end if

rch_df_obs()

subroutine rchmodule::rch_df_obs ( class(rchtype)  this )

private

Store observation type supported by RCH package. Overrides BndTypebnd_df_obs

Definition at line 388 of file gwf-rch.f90.

    implicit none
    ! -- dummy
    class(RchType) :: this
    ! -- local
    integer(I4B) :: indx
    !
    call this%obs%StoreObsType('rch', .true., indx)
    this%obs%obsData(indx)%ProcessIdPtr => defaultobsidprocessor

Here is the call graph for this function:

rch_fc()

subroutine rchmodule::rch_fc ( class(rchtype)  this, real(dp), dimension(:), intent(inout)  rhs, integer(i4b), dimension(:), intent(in)  ia, integer(i4b), dimension(:), intent(in)  idxglo, class(matrixbasetype), pointer  matrix_sln  )

private

Definition at line 297 of file gwf-rch.f90.

    ! -- dummy
    class(RchType) :: this
    real(DP), dimension(:), intent(inout) :: rhs
    integer(I4B), dimension(:), intent(in) :: ia
    integer(I4B), dimension(:), intent(in) :: idxglo
    class(MatrixBaseType), pointer :: matrix_sln
    ! -- local
    integer(I4B) :: i, n, ipos
    !
    ! -- Copy package rhs and hcof into solution rhs and amat
    do i = 1, this%nbound
      n = this%nodelist(i)
      if (n <= 0) cycle
      ! -- reset hcof and rhs for excluded cells
      if (this%ibound(n) == iwetlake) then
        this%hcof(i) = dzero
        this%rhs(i) = dzero
        cycle
      end if
      rhs(n) = rhs(n) + this%rhs(i)
      ipos = ia(n)
      call matrix_sln%add_value_pos(idxglo(ipos), this%hcof(i))
    end do

rch_obs_supported()

logical function rchmodule::rch_obs_supported ( class(rchtype)  this )

private

Definition at line 375 of file gwf-rch.f90.

    implicit none
    ! -- dummy
    class(RchType) :: this
    !
    rch_obs_supported = .true.

rch_read_initial_attr()

subroutine rchmodule::rch_read_initial_attr ( class(rchtype), intent(inout)  this )

private

Definition at line 178 of file gwf-rch.f90.

    ! -- dummy
    class(RchType), intent(inout) :: this
    !
    if (this%readasarrays) then
      call this%default_nodelist()
    end if
    !
    ! -- if fixed_cell option not set, then need to store nodelist
    !    in the nodesontop array
    if (.not. this%fixed_cell) call this%set_nodesontop()

rch_rp()

subroutine rchmodule::rch_rp ( class(rchtype), intent(inout)  this )

private

Read itmp and read new boundaries if itmp > 0

Definition at line 195 of file gwf-rch.f90.

    ! -- modules
    use tdismodule, only: kper
    implicit none
    ! -- dummy
    class(RchType), intent(inout) :: this
    !
    if (this%iper /= kper) return
    !
    call this%BndExtType%bnd_rp()
    !
    ! -- copy nodelist to nodesontop if not fixed cell
    if (.not. this%fixed_cell) call this%set_nodesontop()
    !
    if (this%iprpak /= 0) then
      if (this%readasarrays) then
        ! no-op
      else
        ! -- Write the list to iout
        call this%write_list()
      end if
    end if

rch_source_options()

subroutine rchmodule::rch_source_options

(

class(rchtype), intent(inout)

this

)

Definition at line 128 of file gwf-rch.f90.

    ! -- modules
    use memorymanagerextmodule, only: mem_set_value
    implicit none
    ! -- dummy
    class(RchType), intent(inout) :: this
    ! -- local
    logical(LGP) :: found_fixed_cell = .false.
    !
    ! -- source common bound options
    call this%BndExtType%source_options()
    !
    ! -- update defaults with idm sourced values
    call mem_set_value(this%fixed_cell, 'FIXED_CELL', this%input_mempath, &
                       found_fixed_cell)
    !
    if (this%readasarrays) then
      this%text = texta
    end if
    !
    ! -- log rch params
    call this%log_rch_options(found_fixed_cell)

set_nodesontop()

subroutine rchmodule::set_nodesontop

(

class(rchtype), intent(inout)

this

)

Definition at line 221 of file gwf-rch.f90.

    implicit none
    ! -- dummy
    class(RchType), intent(inout) :: this
    ! -- local
    integer(I4B) :: n
    !
    ! -- allocate if necessary
    if (.not. associated(this%nodesontop)) then
      allocate (this%nodesontop(this%maxbound))
    end if
    !
    ! -- copy nodelist into nodesontop
    do n = 1, this%nbound
      this%nodesontop(n) = this%nodelist(n)
    end do

Variable Documentation

ftype

character(len=lenftype) rchmodule::ftype = 'RCH'

private

Definition at line 19 of file gwf-rch.f90.

  character(len=LENFTYPE) :: ftype = 'RCH'

text

character(len=lenpackagename) rchmodule::text = ' RCH'

private

Definition at line 20 of file gwf-rch.f90.

  character(len=LENPACKAGENAME) :: text = '             RCH'

texta

character(len=lenpackagename) rchmodule::texta = ' RCHA'

private

Definition at line 21 of file gwf-rch.f90.

  character(len=LENPACKAGENAME) :: texta = '            RCHA'