Changeset 2649 for inundation/pyvolution/shallow_water_ext.c

Timestamp:

Apr 2, 2006, 8:07:36 PM (18 years ago)

Author:

steve

Message:

 

File:

• inundation/pyvolution/shallow_water_ext.c (modified) (1 diff)

inundation/pyvolution/shallow_water_ext.c

@@ -245 +245 @@
       h = w[k]-z[k];
       if (h >= eps) {
-        S = -g * eta[k]*eta[k] * sqrt((uh[k]*uh[k] + vh[k]*vh[k]));
-        S /= pow(h, 7.0/3);      //Expensive (on Ole's home computer)
-        //S /= exp(7.0/3.0*log(h));      //seems to save about 15% over manning_friction
-        //S /= h*h*(1 + h/3.0 - h*h/9.0); //FIXME: Could use a Taylor expansion
-
-
-        //Update momentum
-        xmom[k] += S*uh[k];
-        ymom[k] += S*vh[k];
+        S = -g * eta[k]*eta[k] * sqrt((uh[k]*uh[k] + vh[k]*vh[k]));
+        S /= pow(h, 7.0/3);      //Expensive (on Ole's home computer)
+        //S /= exp(7.0/3.0*log(h));      //seems to save about 15% over manning_friction
+        //S /= h*h*(1 + h/3.0 - h*h/9.0); //FIXME: Could use a Taylor expansion
+
+
+        //Update momentum
+        xmom[k] += S*uh[k];
+        ymom[k] += S*vh[k];
       }
     }