Changeset 8204 for trunk/anuga_work/development/2010-projects/anuga_1d/sww/sww_domain.py

Timestamp:

Sep 1, 2011, 7:28:58 PM (14 years ago)

Author:

steve

Message:

Getting the 1d tsunami model working

File:

• trunk/anuga_work/development/2010-projects/anuga_1d/sww/sww_domain.py (modified) (2 diffs)

trunk/anuga_work/development/2010-projects/anuga_1d/sww/sww_domain.py

@@ -235 +235 @@
         
     #Calculate height (and fix negatives)better be non-negative!
+    w_C[:] = numpy.maximum(w_C, z_C)
     h_C[:] = w_C - z_C
+
+
     u_C[:]  = uh_C/(h_C + h0/h_C)
 
@@ -282 +285 @@
     #print 'any' ,numpy.any( h_V[:,1] < 0.0)
 
-    h00 = 1e-12
+    #h00 = 1e-12
     #print h00
 
-    h_V[:,:] = numpy.where (h_V <= h00, 0.0, h_V)
-    u_V[:,:] = numpy.where (h_V <= h00, 0.0, u_V)
+    #h_V[:,:] = numpy.where (h_V <= h00, 0.0, h_V)
+    #u_V[:,:] = numpy.where (h_V <= h00, 0.0, u_V)
 
 #    # protect from edge values going negative