1. Richard Mills
  2. pflotran-dev_surfaceice

Commits

Richard Mills  committed 09ce8eb

Added missing parentheses in calculation of surface velocity in SurfaceTHFlux().

  • Participants
  • Parent commits 99bbb26
  • Branches default

Comments (0)

Files changed (1)

File src/pflotran/surface_th.F90

View file
     ! is often taken to be the flow depth, so I believe this makes sense. (?)
     ! The actual total head term ('hw_half' here) is NOT modified by the 
     ! unfrozen fraction.
-    vel = (unfrozen_fraction_half * hw_half**(2.d0/3.d0))/mannings_half* &
+    vel = ((unfrozen_fraction_half * hw_half)**(2.d0/3.d0))/mannings_half* &
           dhead/(abs(dhead)**(1.d0/2.d0))* &
           1.d0/(dist**0.5d0)