MODFLOW 6  version 6.5.0.dev2
MODFLOW 6 Code Documentation
TernarySolveTrack.f90 File Reference

Go to the source code of this file.

Modules

module  ternarysolvetrack
 

Functions/Subroutines

subroutine, public ternarysolvetrack::traverse_triangle (isolv, tol, step, texit, alpexit, betexit, itrifaceenter, itrifaceexit, rxx, rxy, ryx, ryy, alp0, bet0, alp1, bet1, alp2, bet2, alpi, beti, vziodz, az, bary)
 Traverse triangular cell. More...
 
subroutine, public ternarysolvetrack::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, bary)
 Set coordinates to "canonical" configuration. More...
 
subroutine, public ternarysolvetrack::get_w (alp1, bet1, alp2, bet2, waa, wab, wba, wbb, bary)
 Compute elements of W matrix. More...
 
subroutine, public ternarysolvetrack::solve_coefs (alpi, beti)
 Compute analytical solution coefficients depending on case. More...
 
subroutine, public ternarysolvetrack::step_analytical (t, alp, bet)
 Step (evaluate) analytically depending on case. More...
 
subroutine, public ternarysolvetrack::find_exit_bary (isolv, itriface, itrifaceenter, alpi, beti, tol, step, vziodz, az, texit, alpexit, betexit)
 Find the exit time and location in barycentric coordinates. More...
 
real(dp) function ternarysolvetrack::fbary1 (bet)
 Brent's method applied to canonical face 1 (gamma = 0) More...
 
real(dp) function ternarysolvetrack::fbary2 (bet)
 Brent's method applied to canonical face 2 (alpha = 0) More...
 
subroutine, public ternarysolvetrack::get_t_alpt (bet, t, alp)
 Given beta evaluate t and alpha depending on case. More...
 
subroutine, public ternarysolvetrack::get_bet_outflow_bary (vn1, vn2, betoutlo, betouthi)
 Find outflow interval. More...
 
subroutine, public ternarysolvetrack::get_bet_soln_limits (beti, betsollo, betsolhi, ibettrend)
 Find trend of and limits on beta from beta{t} solution. More...
 
subroutine, public ternarysolvetrack::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 ternarysolvetrack::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) ternarysolvetrack::ca1
 
real(dp) ternarysolvetrack::ca2
 
real(dp) ternarysolvetrack::ca3
 
real(dp) ternarysolvetrack::cb1
 
real(dp) ternarysolvetrack::cb2
 Analytical solution coefficients. More...
 
real(dp) ternarysolvetrack::waa
 
real(dp) ternarysolvetrack::wab
 
real(dp) ternarysolvetrack::wba
 
real(dp) ternarysolvetrack::wbb
 Elements of the "velocity matrix," W. More...
 
real(dp), dimension(2) ternarysolvetrack::cv0
 
real(dp), dimension(2) ternarysolvetrack::cv1
 
real(dp), dimension(2) ternarysolvetrack::cv2
 "Canonical" velocity components at corners of triangular subcell More...
 
integer(i4b) ternarysolvetrack::icase
 Case index for analytical solution. More...