gwfobsmodule Module Reference

Data Types
type	gwfobstype

Functions/Subroutines
subroutine, public	gwf_obs_cr (obs, inobs)
	Create a new GwfObsType object. More...
subroutine	gwf_obs_ar (this, ic, x, flowja)
	Allocate and read. More...
subroutine	gwf_obs_df (this, iout, pkgname, filtyp, dis)
	Define. More...
subroutine	gwf_obs_bd (this)
	Save obs. More...
subroutine	gwf_obs_rp (this)
	Do GWF observations need any checking? If so, add checks here. More...
subroutine	gwf_obs_da (this)
	Deallocate memory. More...
subroutine	set_pointers (this, ic, x, flowja)
	Set pointers. More...
subroutine	gwf_process_head_drawdown_obs_id (obsrv, dis, inunitobs, iout)
	Calculate drawdown observation when requested. More...
subroutine	gwf_process_intercell_obs_id (obsrv, dis, inunitobs, iout)
	Process flow between two cells when requested. More...

Function/Subroutine Documentation

gwf_obs_ar()

subroutine gwfobsmodule::gwf_obs_ar ( class(gwfobstype), intent(inout)  this, type(gwfictype), intent(in), pointer  ic, real(dp), dimension(:), intent(in), pointer, contiguous  x, real(dp), dimension(:), intent(in), pointer, contiguous  flowja  )

private

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

    ! -- dummy
    class(GwfObsType), intent(inout) :: this
    type(GwfIcType), pointer, intent(in) :: ic
    real(DP), dimension(:), pointer, contiguous, intent(in) :: x
    real(DP), dimension(:), pointer, contiguous, intent(in) :: flowja
    !
    ! Call ar method of parent class
    call this%obs_ar()
    !
    ! set pointers
    call this%set_pointers(ic, x, flowja)

gwf_obs_bd()

subroutine gwfobsmodule::gwf_obs_bd ( class(gwfobstype), intent(inout)  this )

private

Definition at line 101 of file gwf-obs.f90.

    ! -- dummy
    class(GwfObsType), intent(inout) :: this
    ! -- local
    integer(I4B) :: i, jaindex, nodenumber
    real(DP) :: v
    character(len=100) :: msg
    class(ObserveType), pointer :: obsrv => null()
    !
    call this%obs_bd_clear()
    !
    ! -- iterate through all GWF observations
    if (this%npakobs > 0) then
      do i = 1, this%npakobs
        obsrv => this%pakobs(i)%obsrv
        nodenumber = obsrv%NodeNumber
        jaindex = obsrv%JaIndex
        select case (obsrv%ObsTypeId)
        case ('HEAD')
          call this%SaveOneSimval(obsrv, this%x(nodenumber))
        case ('DRAWDOWN')
          v = this%ic%strt(nodenumber) - this%x(nodenumber)
          call this%SaveOneSimval(obsrv, v)
        case ('FLOW-JA-FACE')
          call this%SaveOneSimval(obsrv, this%flowja(jaindex))
        case default
          msg = 'Error: Unrecognized observation type: '//trim(obsrv%ObsTypeId)
          call store_error(msg)
        end select
      end do
      !
      ! -- write summary of error messages
      if (count_errors() > 0) then
        call store_error_unit(this%inUnitObs)
      end if
    end if

Here is the call graph for this function:

gwf_obs_cr()

subroutine, public gwfobsmodule::gwf_obs_cr	(	type(gwfobstype), intent(out), pointer	obs,
		integer(i4b), intent(in), pointer	inobs
	)

Create observation object, allocate pointers, initialize values

Definition at line 40 of file gwf-obs.f90.

    ! -- dummy
    type(GwfObsType), pointer, intent(out) :: obs
    integer(I4B), pointer, intent(in) :: inobs
    !
    allocate (obs)
    call obs%allocate_scalars()
    obs%active = .false.
    obs%inputFilename = ''
    obs%inUnitObs => inobs

Here is the caller graph for this function:

gwf_obs_da()

subroutine gwfobsmodule::gwf_obs_da ( class(gwfobstype), intent(inout)  this )

private

Definition at line 150 of file gwf-obs.f90.

    ! -- dummy
    class(GwfObsType), intent(inout) :: this
    !
    nullify (this%ic)
    nullify (this%x)
    nullify (this%flowja)
    call this%ObsType%obs_da()

gwf_obs_df()

subroutine gwfobsmodule::gwf_obs_df ( class(gwfobstype), intent(inout)  this, integer(i4b), intent(in)  iout, character(len=*), intent(in)  pkgname, character(len=*), intent(in)  filtyp, class(disbasetype), pointer  dis  )

private

Definition at line 70 of file gwf-obs.f90.

    ! -- dummy
    class(GwfObsType), intent(inout) :: this
    integer(I4B), intent(in) :: iout
    character(len=*), intent(in) :: pkgname
    character(len=*), intent(in) :: filtyp
    class(DisBaseType), pointer :: dis
    ! -- local
    integer(I4B) :: indx
    !
    ! Call overridden method of parent class
    call this%ObsType%obs_df(iout, pkgname, filtyp, dis)
    !
    ! -- StoreObsType arguments are: (ObserveType, cumulative, indx);
    !    indx is returned.
    !
    ! -- Store obs type and assign procedure pointer for head observation type
    call this%StoreObsType('head', .false., indx)
    this%obsData(indx)%ProcessIdPtr => gwf_process_head_drawdown_obs_id
    !
    ! -- Store obs type and assign procedure pointer for drawdown observation type
    call this%StoreObsType('drawdown', .false., indx)
    this%obsData(indx)%ProcessIdPtr => gwf_process_head_drawdown_obs_id
    !
    ! -- Store obs type and assign procedure pointer for flow-ja-face observation type
    call this%StoreObsType('flow-ja-face', .true., indx)
    this%obsData(indx)%ProcessIdPtr => gwf_process_intercell_obs_id

Here is the call graph for this function:

gwf_obs_rp()

subroutine gwfobsmodule::gwf_obs_rp ( class(gwfobstype), intent(inout)  this )

private

Definition at line 141 of file gwf-obs.f90.

    ! -- dummy
    class(GwfObsType), intent(inout) :: this
    !
    ! Do GWF observations need any checking? If so, add checks here

gwf_process_head_drawdown_obs_id()

subroutine gwfobsmodule::gwf_process_head_drawdown_obs_id ( type(observetype), intent(inout)  obsrv, class(disbasetype), intent(in)  dis, integer(i4b), intent(in)  inunitobs, integer(i4b), intent(in)  iout  )

private

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

    ! -- dummy
    type(ObserveType), intent(inout) :: obsrv
    class(DisBaseType), intent(in) :: dis
    integer(I4B), intent(in) :: inunitobs
    integer(I4B), intent(in) :: iout
    ! -- local
    integer(I4B) :: nn1
    integer(I4B) :: icol, istart, istop
    character(len=LINELENGTH) :: ermsg, string
    !
    ! -- Initialize variables
    string = obsrv%IDstring
    icol = 1
    !
    ! Get node number, with option for ID string to be either node
    ! number or lay, row, column (when dis is structured).
    nn1 = dis%noder_from_string(icol, istart, istop, inunitobs, &
                                iout, string, .false.)
    !
    if (nn1 > 0) then
      obsrv%NodeNumber = nn1
    else
      ermsg = 'Error reading data from ID string'
      call store_error(ermsg)
      call store_error_unit(inunitobs)
    end if

Here is the call graph for this function:

Here is the caller graph for this function:

gwf_process_intercell_obs_id()

subroutine gwfobsmodule::gwf_process_intercell_obs_id ( type(observetype), intent(inout)  obsrv, class(disbasetype), intent(in)  dis, integer(i4b), intent(in)  inunitobs, integer(i4b), intent(in)  iout  )

private

Definition at line 209 of file gwf-obs.f90.

    ! -- dummy
    type(ObserveType), intent(inout) :: obsrv
    class(DisBaseType), intent(in) :: dis
    integer(I4B), intent(in) :: inunitobs
    integer(I4B), intent(in) :: iout
    ! -- local
    integer(I4B) :: nn1, nn2
    integer(I4B) :: icol, istart, istop, jaidx
    character(len=LINELENGTH) :: ermsg, string
    ! formats
70  format('Error: No connection exists between cells identified in text: ', a)
    !
    ! -- Initialize variables
    string = obsrv%IDstring
    icol = 1
    !
    ! Get node number, with option for ID string to be either node
    ! number or lay, row, column (when dis is structured).
    nn1 = dis%noder_from_string(icol, istart, istop, inunitobs, &
                                iout, string, .false.)
    !
    if (nn1 > 0) then
      obsrv%NodeNumber = nn1
    else
      ermsg = 'Error reading data from ID string: '//string(istart:istop)
      call store_error(ermsg)
    end if
    !
    ! Get node number, with option for ID string to be either node
    ! number or lay, row, column (when dis is structured).
    nn2 = dis%noder_from_string(icol, istart, istop, inunitobs, &
                                iout, string, .false.)
    if (nn2 > 0) then
      obsrv%NodeNumber2 = nn2
    else
      ermsg = 'Error reading data from ID string: '//string(istart:istop)
      call store_error(ermsg)
    end if
    !
    ! -- store JA index
    jaidx = dis%con%getjaindex(nn1, nn2)
    if (jaidx == 0) then
      write (ermsg, 70) trim(string)
      call store_error(ermsg)
    end if
    obsrv%JaIndex = jaidx
    !
    if (count_errors() > 0) then
      call store_error_unit(inunitobs)
    end if

Here is the call graph for this function:

Here is the caller graph for this function:

set_pointers()

subroutine gwfobsmodule::set_pointers ( class(gwfobstype), intent(inout)  this, type(gwfictype), intent(in), pointer  ic, real(dp), dimension(:), intent(in), pointer, contiguous  x, real(dp), dimension(:), intent(in), pointer, contiguous  flowja  )

private

Definition at line 162 of file gwf-obs.f90.

    ! -- dummy
    class(GwfObsType), intent(inout) :: this
    type(GwfIcType), pointer, intent(in) :: ic
    real(DP), dimension(:), pointer, contiguous, intent(in) :: x
    real(DP), dimension(:), pointer, contiguous, intent(in) :: flowja
    !
    this%ic => ic
    this%x => x
    this%flowja => flowja