Changeset 2083

Timestamp:

Nov 28, 2005, 7:17:33 PM (19 years ago)

Author:

ole

Message:

Karratha material for talk

Location:

production/karratha_2005

Files:

• get_timeseries.py (modified) (6 diffs)
• process_buildings.py (added)

production/karratha_2005/get_timeseries.py

@@ -1 +1 @@
 """Read in sww file, interpolate at specified locations and plot time series
-Additionally, store augmented building database with max inundation depths
+
 """
 
@@ -9 +9 @@
 
 
-plot = False
     
 swwfile = project.outputname + '.sww'
@@ -41 +40 @@
 gauges, buildings = get_gauges_from_file(project.gauge_filename)
 
+gauges = [[easting, northing-10]]
+
+
 #Read model output
 quantities = ['stage', 'elevation', 'xmomentum', 'ymomentum']
@@ -50 +52 @@
 
 
-buildings[0] = buildings[0].strip() +\
-               ',MAX INUNDATION DEPTH (m)' +\
-               ',MAX MOMENTUM (m^2/s)\n'            
                
 from math import sqrt
@@ -68 +67 @@
     max_momentum = 0    
     for i, t in enumerate(f.T):
+        #if tmin < t < tmax:
         w = f(t, point_id = k)[0]
         z = f(t, point_id = k)[1]
@@ -85 +85 @@
             max_momentum = m    
 
-    #Augment building data
-    buildings[k+1] = buildings[k+1].strip() +\
-                     ',%f' %max_depth +\
-                     ',%f\n' %max_momentum
     
-    #print buildings[k]
-    
-    if plot is True:
-        #Plot only those gauges that have been inundated by more than a threshold
-        if max_depth < 0.2:
-            print 'Skipping gauge %d' %k
-            continue
 
-        #FIXME Reduce to tmin and tmax
-        #i0 = model_timefind
-    
-        ion()
-        hold(False)
+    #Plot only those gauges that have been inundated by more than a threshold
+    #if max_depth < 0.2:
+    #    print 'Skipping gauge %d' %k
+    #    continue
 
-        if elevations[0] < -10:
-            plot(model_time, stages, '-b')
-        else:    
-            plot(model_time, stages, '-b',
-                 model_time, elevations, '-k')
-        name = 'Gauge_%d: (%.1f, %.1f)' %(k, g[0], g[1])
-        title(name)
+    ion()
+    hold(False)
 
-        title('%s (stage)' %name)
-        xlabel('time [s]')
-        ylabel('elevation [m]')    
-        legend(('Stage', 'Bed = %.1f' %elevations[0]),
-               shadow=True,
-               loc='upper right')
-        savefig('Gauge_%d_stage' %k)
+    if elevations[0] < -10:
+        plot(model_time, stages, '-b')
+    else:    
+        plot(model_time, stages, '-b',
+             model_time, elevations, '-k')
+    name = 'Gauge_%d: (%.1f, %.1f)' %(k, g[0], g[1])
+    title(name)
 
-        raw_input('Next')
+    title('%s (stage)' %name)
+    xlabel('time [s]')
+    ylabel('elevation [m]')    
+    legend(('Stage', 'Bed = %.1f' %elevations[0]),
+           shadow=True,
+           loc='upper right')
+    savefig('Gauge_%d_stage' %k)
+
+    raw_input('Next')
 
 
-        #Momentum plot
-        ion()
-        hold(False)
-        plot(model_time, momenta, '-r')
-        title(name)
+    #Momentum plot
+    ion()
+    hold(False)
+    plot(model_time, momenta, '-r')
+    title(name)
 
-        title('%s (momentum)' %name)
-        xlabel('time [s]')
-        ylabel('sqrt( uh^2 + vh^2 ) [m^2/s]')    
-        savefig('Gauge_%d_momentum' %k)
-
-        raw_input('Next')
+    title('%s (momentum)' %name)
+    xlabel('time [s]')
+    ylabel('sqrt( uh^2 + vh^2 ) [m^2/s]')    
+    savefig('Gauge_%d_momentum' %k)
+    
+    raw_input('Next')
     
 
-#Store new building file with mak depths added
+show()
 
-FN = 'inundation_augmented_' + project.gauge_filename
-fid = open(FN, 'w')
-for line in buildings:
-    fid.write(line)
-fid.close()    
-
-
-if plot is True:    
-    show()
-