Changeset 7799

Timestamp:

Jun 7, 2010, 11:41:49 AM (15 years ago)

Author:

habili

Message:

changed the order of assigning "no values" to stage, current_x and current_y.

File:

• branches/anuga_1_1/anuga_work/production/bunbury_storm_surge_2009/process_gems_grids.py (modified) (3 diffs)

branches/anuga_1_1/anuga_work/production/bunbury_storm_surge_2009/process_gems_grids.py

@@ -122 +122 @@
         #    print 'third block'
 
-        if i > 13 + 3*nrows and i < (14 + 4*nrows):   
+        if i > (13 + 3*nrows) and i < (14 + 4*nrows):   
             d += L.strip().split()
             
@@ -158 +158 @@
 no_value_index = numpy.where(((depth < -9000) + (current_x < -9000) + (current_y < -9000)) == True)[0]
 
-numpy.put(stage, no_value_index, -9999)
-numpy.put(current_x, no_value_index, -9999)
-numpy.put(current_y, no_value_index, -9999)
-numpy.put(depth, no_value_index, 0)
+#numpy.put(stage, no_value_index, -9999)
+#numpy.put(current_x, no_value_index, -9999)
+#numpy.put(current_y, no_value_index, -9999)
+#numpy.put(depth, no_value_index, 0)
     
 # Taking absolute value is to account for -ve depths obtained when depth-elevation 
@@ -174 +174 @@
 ymomentum = current_y*depth #momentum*numpy.cos(numpy.radians(current_y))
 
+numpy.put(stage, no_value_index, 0)
 numpy.put(xmomentum, no_value_index, 0)
 numpy.put(ymomentum, no_value_index, 0)