ghbmodule Module Reference

Data Types
type	ghbtype

Functions/Subroutines
subroutine, public	ghb_create (packobj, id, ibcnum, inunit, iout, namemodel, pakname, mempath)
	Create a New Ghb Package and point bndobj to the new package. More...
subroutine	ghb_da (this)
	Deallocate memory. More...
subroutine	ghb_options (this)
	Set options specific to GhbType. More...
subroutine	log_ghb_options (this, found_mover)
	Log options specific to GhbType. More...
subroutine	ghb_allocate_arrays (this, nodelist, auxvar)
	Allocate arrays. More...
subroutine	ghb_rp (this)
	Read and prepare. More...
subroutine	ghb_ck (this)
	Check ghb boundary condition data. More...
subroutine	ghb_cf (this)
	Formulate the HCOF and RHS terms. More...
subroutine	ghb_fc (this, rhs, ia, idxglo, matrix_sln)
	Copy rhs and hcof into solution rhs and amat. More...
subroutine	define_listlabel (this)
	Define the list heading that is written to iout when PRINT_INPUT option is used. More...
logical function	ghb_obs_supported (this)
	Return true because GHB package supports observations. More...
subroutine	ghb_df_obs (this)
	Store observation type supported by GHB package. More...
subroutine	ghb_store_user_cond (this)
	Store user-specified conductance for GHB boundary type. More...
real(dp) function	cond_mult (this, row)
	Apply multiplier to GHB conductance if option AUXMULTCOL is used. More...
real(dp) function	ghb_bound_value (this, col, row)
	Return requested boundary value. More...

Variables
character(len=lenftype)	ftype = 'GHB'
character(len=lenpackagename)	text = ' GHB'

Function/Subroutine Documentation

cond_mult()

real(dp) function ghbmodule::cond_mult ( class(ghbtype), intent(inout)  this, integer(i4b), intent(in)  row  )

private

Definition at line 376 of file gwf-ghb.f90.

    ! -- modules
    use constantsmodule, only: dzero
    ! -- dummy variables
    class(GhbType), intent(inout) :: this
    integer(I4B), intent(in) :: row
    ! -- result
    real(DP) :: cond
    !
    if (this%iauxmultcol > 0) then
      cond = this%cond(row) * this%auxvar(this%iauxmultcol, row)
    else
      cond = this%cond(row)
    end if

define_listlabel()

subroutine ghbmodule::define_listlabel ( class(ghbtype), intent(inout)  this )

private

Definition at line 303 of file gwf-ghb.f90.

    ! -- dummy
    class(GhbType), 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), 'STAGE'
    write (this%listlabel, '(a, a16)') trim(this%listlabel), 'CONDUCTANCE'
    if (this%inamedbound == 1) then
      write (this%listlabel, '(a, a16)') trim(this%listlabel), 'BOUNDARY NAME'
    end if

ghb_allocate_arrays()

subroutine ghbmodule::ghb_allocate_arrays ( class(ghbtype)  this, integer(i4b), dimension(:), optional, pointer, contiguous  nodelist, real(dp), dimension(:, :), optional, pointer, contiguous  auxvar  )

private

Definition at line 142 of file gwf-ghb.f90.

    ! -- modules
    use memorymanagermodule, only: mem_allocate, mem_setptr, mem_checkin
    ! -- dummy
    class(GhbType) :: this
    integer(I4B), dimension(:), pointer, contiguous, optional :: nodelist
    real(DP), dimension(:, :), pointer, contiguous, optional :: auxvar
    !
    ! -- call base type allocate arrays
    call this%BndExtType%allocate_arrays(nodelist, auxvar)
    !
    ! -- set ghb input context pointers
    call mem_setptr(this%bhead, 'BHEAD', this%input_mempath)
    call mem_setptr(this%cond, 'COND', this%input_mempath)
    !
    ! --checkin ghb input context pointers
    call mem_checkin(this%bhead, 'BHEAD', this%memoryPath, &
                     'BHEAD', this%input_mempath)
    call mem_checkin(this%cond, 'COND', this%memoryPath, &
                     'COND', this%input_mempath)

ghb_bound_value()

real(dp) function ghbmodule::ghb_bound_value	(	class(ghbtype), intent(inout)	this,
		integer(i4b), intent(in)	col,
		integer(i4b), intent(in)	row
	)

Definition at line 394 of file gwf-ghb.f90.

    ! -- modules
    use constantsmodule, only: dzero
    ! -- dummy
    class(GhbType), intent(inout) :: this
    integer(I4B), intent(in) :: col
    integer(I4B), intent(in) :: row
    ! -- result
    real(DP) :: bndval
    !
    select case (col)
    case (1)
      bndval = this%bhead(row)
    case (2)
      bndval = this%cond_mult(row)
    case default
      errmsg = 'Programming error. GHB bound value requested column '&
               &'outside range of ncolbnd (2).'
      call store_error(errmsg)
      call store_error_filename(this%input_fname)
    end select

Here is the call graph for this function:

ghb_cf()

subroutine ghbmodule::ghb_cf

(

class(ghbtype)

this

)

Skip if no GHBs

Definition at line 240 of file gwf-ghb.f90.

    ! -- dummy
    class(GhbType) :: this
    ! -- local
    integer(I4B) :: i, node
    !
    ! -- Return if no ghbs
    if (this%nbound .eq. 0) return
    !
    ! -- Calculate hcof and rhs for each ghb entry
    do i = 1, this%nbound
      node = this%nodelist(i)
      if (node == 0) cycle
      if (this%ibound(node) .le. 0) then
        this%hcof(i) = dzero
        this%rhs(i) = dzero
        cycle
      end if
      this%hcof(i) = -this%cond_mult(i)
      this%rhs(i) = -this%cond_mult(i) * this%bhead(i)
    end do

ghb_ck()

subroutine ghbmodule::ghb_ck

(

class(ghbtype), intent(inout)

this

)

Definition at line 185 of file gwf-ghb.f90.

    ! -- modules
    use constantsmodule, only: linelength
    use simmodule, only: store_error, count_errors, store_error_unit
    ! -- dummy
    class(GhbType), intent(inout) :: this
    ! -- local
    character(len=LINELENGTH) :: errmsg
    integer(I4B) :: i
    integer(I4B) :: node
    real(DP) :: bt
    ! -- formats
    character(len=*), parameter :: fmtghberr = &
      "('GHB BOUNDARY (',i0,') HEAD (',f10.3,') IS LESS THAN CELL &
      &BOTTOM (',f10.3,')')"
    character(len=*), parameter :: fmtcondmulterr = &
      "('GHB BOUNDARY (',i0,') CONDUCTANCE MULTIPLIER (',g10.3,') IS &
      &LESS THAN ZERO')"
    character(len=*), parameter :: fmtconderr = &
      "('GHB BOUNDARY (',i0,') CONDUCTANCE (',g10.3,') IS LESS THAN &
      &ZERO')"
    !
    ! -- check stress period data
    do i = 1, this%nbound
      node = this%nodelist(i)
      if (node == 0) cycle
      bt = this%dis%bot(node)
      ! -- accumulate errors
      if (this%bhead(i) < bt .and. this%icelltype(node) /= 0) then
        write (errmsg, fmt=fmtghberr) i, this%bhead(i), bt
        call store_error(errmsg)
      end if
      if (this%iauxmultcol > 0) then
        if (this%auxvar(this%iauxmultcol, i) < dzero) then
          write (errmsg, fmt=fmtcondmulterr) &
            i, this%auxvar(this%iauxmultcol, i)
          call store_error(errmsg)
        end if
      end if
      if (this%cond(i) < dzero) then
        write (errmsg, fmt=fmtconderr) i, this%cond(i)
        call store_error(errmsg)
      end if
    end do
    !
    !write summary of ghb package error messages
    if (count_errors() > 0) then
      call store_error_unit(this%inunit)
    end if

Here is the call graph for this function:

ghb_create()

subroutine, public ghbmodule::ghb_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
	)

Definition at line 46 of file gwf-ghb.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(GhbType), pointer :: ghbobj
    !
    ! -- allocate the object and assign values to object variables
    allocate (ghbobj)
    packobj => ghbobj
    !
    ! -- create name and memory path
    call packobj%set_names(ibcnum, namemodel, pakname, ftype, mempath)
    packobj%text = text
    !
    ! -- allocate scalars
    call packobj%allocate_scalars()
    !
    ! -- initialize package
    call packobj%pack_initialize()
    !
    packobj%inunit = inunit
    packobj%iout = iout
    packobj%id = id
    packobj%ibcnum = ibcnum
    packobj%ncolbnd = 2
    packobj%ictMemPath = create_mem_path(namemodel, 'NPF')

Here is the call graph for this function:

Here is the caller graph for this function:

ghb_da()

subroutine ghbmodule::ghb_da ( class(ghbtype)  this )

private

Definition at line 84 of file gwf-ghb.f90.

    ! -- modules
    use memorymanagermodule, only: mem_deallocate
    ! -- dummy
    class(GhbType) :: this
    !
    ! -- Deallocate parent package
    call this%BndExtType%bnd_da()
    !
    ! -- arrays
    call mem_deallocate(this%bhead, 'BHEAD', this%memoryPath)
    call mem_deallocate(this%cond, 'COND', this%memoryPath)

ghb_df_obs()

subroutine ghbmodule::ghb_df_obs ( class(ghbtype)  this )

private

Overrides BndTypebnd_df_obs

Definition at line 344 of file gwf-ghb.f90.

    implicit none
    ! -- dummy
    class(GhbType) :: this
    ! -- local
    integer(I4B) :: indx
    !
    call this%obs%StoreObsType('ghb', .true., indx)
    this%obs%obsData(indx)%ProcessIdPtr => defaultobsidprocessor
    !
    ! -- Store obs type and assign procedure pointer
    !    for to-mvr observation type.
    call this%obs%StoreObsType('to-mvr', .true., indx)
    this%obs%obsData(indx)%ProcessIdPtr => defaultobsidprocessor

Here is the call graph for this function:

ghb_fc()

subroutine ghbmodule::ghb_fc ( class(ghbtype)  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 265 of file gwf-ghb.f90.

    ! -- dummy
    class(GhbType) :: 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
    real(DP) :: cond, bhead, qghb
    !
    ! -- pakmvrobj fc
    if (this%imover == 1) then
      call this%pakmvrobj%fc()
    end if
    !
    ! -- Copy package rhs and hcof into solution rhs and amat
    do i = 1, this%nbound
      n = this%nodelist(i)
      if (n == 0) cycle
      rhs(n) = rhs(n) + this%rhs(i)
      ipos = ia(n)
      call matrix_sln%add_value_pos(idxglo(ipos), this%hcof(i))
      !
      ! -- If mover is active and this boundary is discharging,
      !    store available water (as positive value).
      bhead = this%bhead(i)
      if (this%imover == 1 .and. this%xnew(n) > bhead) then
        cond = this%cond_mult(i)
        qghb = cond * (this%xnew(n) - bhead)
        call this%pakmvrobj%accumulate_qformvr(i, qghb)
      end if
    end do

ghb_obs_supported()

logical function ghbmodule::ghb_obs_supported ( class(ghbtype)  this )

private

Overrides BndTypebnd_obs_supported()

Definition at line 332 of file gwf-ghb.f90.

    implicit none
    ! -- dummy
    class(GhbType) :: this
    !
    ghb_obs_supported = .true.

ghb_options()

subroutine ghbmodule::ghb_options

(

class(ghbtype), intent(inout)

this

)

Definition at line 100 of file gwf-ghb.f90.

    ! -- modules
    use memorymanagerextmodule, only: mem_set_value
    use characterstringmodule, only: characterstringtype
    ! -- dummy
    class(GhbType), intent(inout) :: this
    ! -- local
    logical(LGP) :: found_mover
    !
    ! -- source base class options
    call this%BndExtType%source_options()
    !
    ! -- source options from input context
    call mem_set_value(this%imover, 'MOVER', this%input_mempath, found_mover)
    !
    ! -- log ghb specific options
    call this%log_ghb_options(found_mover)

ghb_rp()

subroutine ghbmodule::ghb_rp

(

class(ghbtype), intent(inout)

this

)

Definition at line 166 of file gwf-ghb.f90.

    ! -- modules
    use tdismodule, only: kper
    ! -- dummy
    class(GhbType), intent(inout) :: this
    !
    if (this%iper /= kper) return
    !
    ! -- Call the parent class read and prepare
    call this%BndExtType%bnd_rp()
    !
    ! -- store user cond
    if (this%ivsc == 1) then
      call this%ghb_store_user_cond()
    end if

ghb_store_user_cond()

subroutine ghbmodule::ghb_store_user_cond ( class(ghbtype), intent(inout)  this )

private

Definition at line 362 of file gwf-ghb.f90.

    ! -- dummy
    class(GhbType), intent(inout) :: this
    ! -- local
    integer(I4B) :: n
    !
    ! -- store backup copy of conductance values
    do n = 1, this%nbound
      this%condinput(n) = this%cond_mult(n)
    end do

log_ghb_options()

subroutine ghbmodule::log_ghb_options	(	class(ghbtype), intent(inout)	this,
		logical(lgp), intent(in)	found_mover
	)

Definition at line 121 of file gwf-ghb.f90.

    ! -- modules
    ! -- dummy
    class(GhbType), intent(inout) :: this
    logical(LGP), intent(in) :: found_mover
    !
    ! -- log found options
    write (this%iout, '(/1x,a)') 'PROCESSING '//trim(adjustl(this%text)) &
      //' OPTIONS'
    !
    if (found_mover) then
      write (this%iout, '(4x,A)') 'MOVER OPTION ENABLED'
    end if
    !
    ! -- close logging block
    write (this%iout, '(1x,a)') &
      'END OF '//trim(adjustl(this%text))//' OPTIONS'

Variable Documentation

ftype

character(len=lenftype) ghbmodule::ftype = 'GHB'

private

Definition at line 18 of file gwf-ghb.f90.

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

text

character(len=lenpackagename) ghbmodule::text = ' GHB'

private

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

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