uzfetutilmodule Module Reference

Functions/Subroutines
real(dp) function, public	etfunc_lin (efflndsrf, extdp, resid_pet, deriv_et, trhs, thcof, hgwf, celtop, celbot)
	Calculate gwf ET using linear decay ET function from mf-2005. More...
real(dp) function, public	etfunc_nlin (efflndsrf, extdp, resid_pet, deriv_et, trhs, thcof, hgwf)
	Calculate gwf ET using a square decay ET function with smoothing at the specified extinction depth. More...
pure real(dp) function, public	calc_lin_scaling_fac (hgwf, lndsrf, extdp)
	Calculate the linear scaling factor. More...

Function/Subroutine Documentation

calc_lin_scaling_fac()

pure real(dp) function, public uzfetutilmodule::calc_lin_scaling_fac	(	real(dp), intent(in)	hgwf,
		real(dp), intent(in)	lndsrf,
		real(dp), intent(in)	extdp
	)

Parameters

[in]	hgwf	groundwater head
[in]	lndsrf	effective land surface (after applying surfdep to land surface)
[in]	extdp	extinction depth

Definition at line 107 of file UzfEtUtil.f90.

    ! dummy
    real(DP), intent(in) :: hgwf !< groundwater head
    real(DP), intent(in) :: lndsrf !< effective land surface (after applying surfdep to land surface)
    real(DP), intent(in) :: extdp !< extinction depth
    ! return
    real(DP) :: sclfac
    !
    sclfac = (hgwf - (lndsrf - extdp)) / extdp
    ! the calculated scaling factor cannot exceed 1.0
    if (sclfac > done) sclfac = done

Here is the caller graph for this function:

etfunc_lin()

real(dp) function, public uzfetutilmodule::etfunc_lin	(	real(dp), intent(in)	efflndsrf,
		real(dp), intent(in)	extdp,
		real(dp), intent(in)	resid_pet,
		real(dp), intent(inout)	deriv_et,
		real(dp), intent(inout)	trhs,
		real(dp), intent(inout)	thcof,
		real(dp), intent(in)	hgwf,
		real(dp), intent(in)	celtop,
		real(dp), intent(in)	celbot
	)

Parameters

[in]	efflndsrf	effective land surface elevation after subtracting off 0.5*surfdep
[in]	extdp	extinction depth
[in]	resid_pet	residual pET remaining after subtracting simulated ET in the unsaturated zone
[in,out]	deriv_et	derivative of gw ET for Newton addition to equations in _fn()
[in,out]	trhs	total uzf rhs contribution to GWF model
[in,out]	thcof	total uzf hcof contribution to GWF model
[in]	hgwf	calculated groundwater head
[in]	celtop	elevation of the top of the cell
[in]	celbot	elevation of the bottom of the cell

Definition at line 16 of file UzfEtUtil.f90.

    ! return
    real(DP) :: etfunc_lin
    ! dummy
    real(DP), intent(in) :: efflndsrf !< effective land surface elevation after subtracting off 0.5*surfdep
    real(DP), intent(in) :: extdp !< extinction depth
    real(DP), intent(in) :: resid_pet !< residual pET remaining after subtracting simulated ET in the unsaturated zone
    real(DP), intent(inout) :: deriv_et !< derivative of gw ET for Newton addition to equations in _fn()
    real(DP), intent(inout) :: trhs !< total uzf rhs contribution to GWF model
    real(DP), intent(inout) :: thcof !< total uzf hcof contribution to GWF model
    real(DP), intent(in) :: hgwf !< calculated groundwater head
    real(DP), intent(in) :: celtop !< elevation of the top of the cell
    real(DP), intent(in) :: celbot !< elevation of the bottom of the cell
    ! local
    real(DP) :: etgw
    real(DP) :: range
    real(DP) :: depth, scale, thick, lin_scaling_fac
    !
    ! initialize
    etgw = dzero
    !
    ! if water table between ET surface and extinction depth
    ! (extdp starts at the bottom of the effective land surface elevation,
    !  where the effective land surface elevation accounts for surfdep)
    if (hgwf > (efflndsrf - extdp) .and. hgwf < efflndsrf) THEN
      lin_scaling_fac = calc_lin_scaling_fac(hgwf, efflndsrf, extdp)
      etgw = resid_pet * lin_scaling_fac
      !
      trhs = resid_pet - resid_pet * efflndsrf / extdp
      thcof = -resid_pet / extdp
      etgw = trhs - (thcof * hgwf)
      !
      ! water table above land surface
    else if (hgwf >= efflndsrf) then
      trhs = resid_pet
      etgw = resid_pet
      !
      ! water table below extinction depth
    else
      etfunc_lin = dzero
      return
    end if
    !
    ! calculate rate
    depth = hgwf - (efflndsrf - extdp)
    thick = celtop - celbot
    if (depth > thick) depth = thick
    if (depth < dzero) depth = dzero
    range = dem4 * extdp
    call scubic(depth, range, deriv_et, scale)
    trhs = scale * trhs
    thcof = scale * thcof
    etgw = trhs - (thcof * hgwf)
    deriv_et = -deriv_et * etgw
    etfunc_lin = etgw
    !

Here is the call graph for this function:

Here is the caller graph for this function:

etfunc_nlin()

real(dp) function, public uzfetutilmodule::etfunc_nlin	(	real(dp), intent(in)	efflndsrf,
		real(dp), intent(in)	extdp,
		real(dp), intent(in)	resid_pet,
		real(dp), intent(inout)	deriv_et,
		real(dp), intent(inout)	trhs,
		real(dp), intent(inout)	thcof,
		real(dp), intent(in)	hgwf
	)

Parameters

[in]	efflndsrf	effective land surface elevation after subtracting off 0.5*surfdep
[in]	extdp	extinction depth
[in]	resid_pet	residual pET remaining after subtracting simulated ET in the unsaturated zone
[in,out]	deriv_et	derivative of gw ET for Newton addition to equations in _fn()
[in,out]	trhs	total uzf rhs contribution to GWF model
[in,out]	thcof	total uzf hcof contribution to GWF model
[in]	hgwf	calculated groundwater head

Definition at line 78 of file UzfEtUtil.f90.

    ! -- return
    real(DP) :: etfunc_nlin
    ! -- dummy
    real(DP), intent(in) :: efflndsrf !< effective land surface elevation after subtracting off 0.5*surfdep
    real(DP), intent(in) :: extdp !< extinction depth
    real(DP), intent(in) :: resid_pet !< residual pET remaining after subtracting simulated ET in the unsaturated zone
    real(DP), intent(inout) :: deriv_et !< derivative of gw ET for Newton addition to equations in _fn()
    real(DP), intent(inout) :: trhs !< total uzf rhs contribution to GWF model
    real(DP), intent(inout) :: thcof !< total uzf hcof contribution to GWF model
    real(DP), intent(in) :: hgwf !< calculated groundwater head
    ! -- local
    real(DP) :: etgw
    real(DP) :: range
    real(DP) :: depth, scale
    !
    depth = hgwf - (efflndsrf - extdp)
    if (depth < dzero) depth = dzero
    etgw = resid_pet
    range = dem3 * extdp
    call scubic(depth, range, deriv_et, scale)
    etgw = etgw * scale
    trhs = etgw
    deriv_et = -deriv_et * etgw
    etfunc_nlin = etgw

Here is the call graph for this function:

Here is the caller graph for this function: