GwfSfrCommon.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
real(dp) function	kinematic_residual (qa, qb, qc, qd, aa, ab, ac, ad, qsrc, length, weight, delt, courant)
	Kinematic routing equation residual. More...
real(dp) function	kinematic_storage (aa, ab, ac, ad, length, delt, courant)
	Kinematic routing equation storage term. More...

Function/Subroutine Documentation

kinematic_residual()

real(dp) function kinematic_residual	(	real(dp), intent(in)	qa,
		real(dp), intent(in)	qb,
		real(dp), intent(in)	qc,
		real(dp), intent(in)	qd,
		real(dp), intent(in)	aa,
		real(dp), intent(in)	ab,
		real(dp), intent(in)	ac,
		real(dp), intent(in)	ad,
		real(dp), intent(in)	qsrc,
		real(dp), intent(in)	length,
		real(dp), intent(in)	weight,
		real(dp), intent(in)	delt,
		real(dp), intent(in)	courant
	)

Method to calculate the kinematic-wave routing residual.

Returns

kinematic-wave residual

Parameters

[in]	qa	upstream flow at previous time
[in]	qb	downstream flow at previous time
[in]	qc	upstream flow
[in]	qd	downstream flow
[in]	aa	upstream area at previous time
[in]	ab	downstream area at previous time
[in]	ac	upstream area
[in]	ad	downstream area
[in]	qsrc	lateral flow term (L3T-1L-1)
[in]	length	reach length
[in]	weight	temporal weight
[in]	delt	time step length
[in]	courant	courant number

Definition at line 7 of file GwfSfrCommon.f90.

    use kindmodule, only: dp
    use constantsmodule, only: dzero, dhalf, done
    ! -- return variable
    real(DP) :: kinematic_residual !< kinematic-wave residual
    ! -- dummy variables
    real(DP), intent(in) :: qa !< upstream flow at previous time
    real(DP), intent(in) :: qb !< downstream flow at previous time
    real(DP), intent(in) :: qc !< upstream flow
    real(DP), intent(in) :: qd !< downstream flow
    real(DP), intent(in) :: aa !< upstream area at previous time
    real(DP), intent(in) :: ab !< downstream area at previous time
    real(DP), intent(in) :: ac !< upstream area
    real(DP), intent(in) :: ad !< downstream area
    real(DP), intent(in) :: qsrc !< lateral flow term (L3T-1L-1)
    real(DP), intent(in) :: length !< reach length
    real(DP), intent(in) :: weight !< temporal weight
    real(DP), intent(in) :: delt !< time step length
    real(DP), intent(in) :: courant !< courant number
    ! --local variables
    real(DP) :: f11
    real(DP) :: f12
 
    if (courant > done) then
      f11 = (qd - qc) / length
      f12 = (ac - aa) / delt
    else
      f11 = (weight * (qd - qc) + (done - weight) * (qb - qa)) / length
      f12 = dhalf * ((ad - ab) + (ac - aa)) / delt
    end if
    kinematic_residual = f11 + f12 - qsrc
 

kinematic_storage()

real(dp) function kinematic_storage	(	real(dp), intent(in)	aa,
		real(dp), intent(in)	ab,
		real(dp), intent(in)	ac,
		real(dp), intent(in)	ad,
		real(dp), intent(in)	length,
		real(dp), intent(in)	delt,
		real(dp), intent(in)	courant
	)

Method to calculate the kinematic-wave routing storage term.

Returns

kinematic-wave storage change

Parameters

[in]	aa	upstream area at previous time
[in]	ab	downstream area at previous time
[in]	ac	upstream area
[in]	ad	downstream area
[in]	length	reach length
[in]	delt	time step length
[in]	courant	courant number

Definition at line 50 of file GwfSfrCommon.f90.

    use kindmodule, only: dp
    use constantsmodule, only: dhalf
 
    real(DP) :: kinematic_storage !< kinematic-wave storage change
 
    real(DP), intent(in) :: aa !< upstream area at previous time
    real(DP), intent(in) :: ab !< downstream area at previous time
    real(DP), intent(in) :: ac !< upstream area
    real(DP), intent(in) :: ad !< downstream area
    real(DP), intent(in) :: length !< reach length
    real(DP), intent(in) :: delt !< time step length
    real(DP), intent(in) :: courant !< courant number
 
    if (courant > done) then
      kinematic_storage = (aa - ac) * length / delt
    else
      kinematic_storage = dhalf * ((ac + ad) - (aa + ab)) * length / delt
    end if