MODFLOW 6  version 6.7.0.dev2
USGS Modular Hydrologic Model
TernarySolveUtils.f90 File Reference

Go to the source code of this file.

Modules

module  ternarysolveutils
 

Functions/Subroutines

subroutine, public ternarysolveutils::traverse_triangle (isolv, tol, texit, alpexit, betexit, itrifaceenter, itrifaceexit, alp1, bet1, alp2, bet2, alpi, beti)
 Traverse triangular cell. More...
 
subroutine, public ternarysolveutils::canonical (x0, y0, x1, y1, x2, y2, v0x, v0y, v1x, v1y, v2x, v2y, xi, yi, rxx, rxy, ryx, ryy, sxx, sxy, syy, alp0, bet0, alp1, bet1, alp2, bet2, alpi, beti)
 Set coordinates to "canonical" configuration. More...
 
subroutine, public ternarysolveutils::get_w (alp1, bet1, alp2, bet2, waa, wab, wba, wbb)
 Compute elements of W matrix. More...
 
subroutine, public ternarysolveutils::solve_coefs (alpi, beti)
 Compute analytical solution coefficients depending on case. More...
 
subroutine, public ternarysolveutils::step_analytical (t, alp, bet)
 Step (evaluate) analytically depending on case. More...
 
subroutine, public ternarysolveutils::find_exit_bary (isolv, itriface, itrifaceenter, alpi, beti, tol, texit, alpexit, betexit)
 Find the exit time and location in barycentric coordinates. More...
 
real(dp) function ternarysolveutils::fbary1 (bet)
 Brent's method applied to canonical face 1 (gamma = 0) More...
 
real(dp) function ternarysolveutils::fbary2 (bet)
 Brent's method applied to canonical face 2 (alpha = 0) More...
 
subroutine, public ternarysolveutils::get_t_alpt (bet, t, alp)
 Given beta evaluate t and alpha depending on case. More...
 
subroutine, public ternarysolveutils::get_bet_outflow_bary (vn1, vn2, betoutlo, betouthi)
 Find outflow interval. More...
 
subroutine, public ternarysolveutils::get_bet_soln_limits (beti, betsollo, betsolhi, ibettrend)
 Find trend of and limits on beta from beta{t} solution. More...
 
subroutine, public ternarysolveutils::soln_brent (itriface, betlo, bethi, tol, texit, alpexit, betexit)
 Use Brent's method with initial bounds on beta of betlo and bethi. More...
 
subroutine, public ternarysolveutils::soln_chand (itriface, betlo, bethi, tol, texit, alpexit, betexit)
 Use Chandrupatla's method with initial bounds on beta of betlo and bethi. More...
 

Variables

real(dp) ternarysolveutils::ca1
 
real(dp) ternarysolveutils::ca2
 
real(dp) ternarysolveutils::ca3
 
real(dp) ternarysolveutils::cb1
 
real(dp) ternarysolveutils::cb2
 Analytical solution coefficients. More...
 
real(dp) ternarysolveutils::waa
 
real(dp) ternarysolveutils::wab
 
real(dp) ternarysolveutils::wba
 
real(dp) ternarysolveutils::wbb
 Elements of the "velocity matrix," W. More...
 
real(dp), dimension(2) ternarysolveutils::cv0
 
real(dp), dimension(2) ternarysolveutils::cv1
 
real(dp), dimension(2) ternarysolveutils::cv2
 "Canonical" velocity components at corners of triangular subcell More...
 
integer(i4b) ternarysolveutils::icase
 Case index for analytical solution. More...