Changeset 3300 for production/onslow_2006/get_building_inundation.py

Timestamp:

Jul 11, 2006, 11:30:36 AM (18 years ago)

Author:

sexton

Message:

updates

File:

• production/onslow_2006/get_building_inundation.py (modified) (2 diffs)

production/onslow_2006/get_building_inundation.py

@@ -30 +30 @@
 # Inputs
 #timestampdir = '20060704_063005' # HAT
-#timestampdir = '20060704_063112' # LAT
+#timestampdir = '20060706_235246' # LAT
 timestampdir = '20060704_063234' # MSL
+#timestampdir = '20060515_001733' # MSL DTED
 file_loc = project.outputdir + timestampdir + sep 
 swwfile = file_loc + project.basename + '.sww'
@@ -101 +102 @@
     zero_depth = 0.0
     zero_momentum = 0.0
-    if inside_polygon(g,project.polyAll) == True:   
-        for i, t in enumerate(T):
-            w = f(t, point_id = k)[0]
-            z = f(t, point_id = k)[1]
-            uh = f(t, point_id = k)[2]
-            vh = f(t, point_id = k)[3]
-            depth = w-z
-            m = sqrt(uh*uh + vh*vh)   
-            #vel = m / (depth + 1.e-30) 
-            if depth > max_depth: max_depth = depth
-            if m > max_momentum: max_momentum = m
-    else:
-        max_depth = 0.0
-        max_momentum = 0.0
-
+    #if inside_polygon(g,project.polyAll) == True:   
+    for i, t in enumerate(T):
+        w = f(t, point_id = k)[0]
+        z = f(t, point_id = k)[1]
+        uh = f(t, point_id = k)[2]
+        vh = f(t, point_id = k)[3]
+        depth = w-z
+        m = sqrt(uh*uh + vh*vh)   
+        #vel = m / (depth + 1.e-30) 
+        if depth > max_depth: max_depth = depth
+        if m > max_momentum: max_momentum = m
+    #else:
+    #    max_depth = 0.0
+    #    max_momentum = 0.0
+    print 'gauge ', g
+    print 'max depth', max_depth
     lines[k+1] = lines[k+1].strip() +\
                  ',%f' %max_depth +\