Changeset 3190 for production/sydney_2006/get_timeseries.py

Timestamp:

Jun 21, 2006, 9:31:57 AM (18 years ago)

Author:

sexton

Message:

MOST and ANUGA comparisons and updates

File:

• production/sydney_2006/get_timeseries.py (modified) (10 diffs)

production/sydney_2006/get_timeseries.py

@@ -5 +5 @@
 import project
 from pyvolution.util import file_function
-from coordinate_transforms.redfearn import degminsec2decimal_degrees, redfearn
+#from coordinate_transforms.redfearn import degminsec2decimal_degrees, redfearn
 from pylab import *
 from matplotlib.ticker import MultipleLocator, FormatStrFormatter
 
-swwfile = project.outputname + '.sww' 
-
-
+swwfile = project.outputname + '.sww'
+#swwfile = project.outputdir + timestampdir + sep + project.outputname + '.sww'
+# place to store figures
+#graphloc = project.outputdir + timestampdir + sep
+graphloc = project.outputdir
 
 #Time interval to plot
-tmin = 13000
-tmax = 21000
+tmin = 2*60
+tmax = 10*60
 
 def get_gauges_from_file(filename):
@@ -34 +36 @@
         gaugelocation.append(location)
         
-    #Return gauges and raw data for subsequent storage
     return gauges, lines, gaugelocation
 
@@ -40 +41 @@
 gauges, lines, locations = get_gauges_from_file(project.gauge_filename)
 
-print 'number of gauges for Benfield:   ', len(gauges)
+print 'number of gauges:   ', len(gauges)
 
 #Read model output
@@ -50 +51 @@
                   use_cache = True)
 
-
-               
+print 'size f', size(f.quantities['stage'],axis=0), size(f.quantities['stage'],axis=1)
+
 from math import sqrt, atan, degrees
 from Numeric import ones
@@ -85 +86 @@
         uh = f(t, point_id = k)[2]
         vh = f(t, point_id = k)[3]
-        myloc = locations[k]
+        gaugeloc = locations[k]
         depth = w-z
-
+        
         m = sqrt(uh*uh + vh*vh)   #Absolute momentum
         vel = sqrt(uh*uh + vh*vh) / (w-z + 1.e-30) #Absolute velocity
@@ -135 +136 @@
              model_time, elevations, '-k')
     #name = 'Gauge_%d: (%.1f, %.1f)' %(k, g[0], g[1])
-    name = 'Gauge_%d: (%.1f, %.1f) Location: %s' %(k, g[0], g[1], myloc)
+    name = 'Gauge_%d: (%.1f, %.1f) Location: %s' %(k, g[0], g[1], gaugeloc)
     title(name)
 
@@ -144 +145 @@
            shadow=True,
            loc='upper right')
-    #savefig('Gauge_%d_stage' %k)
-    savefig('Gauge_%s_stage' %myloc)
-
-    # raw_input('Next')
+    #('Gauge_%d_stage' %k)
+    savefig('%sGauge_%s_stage' %(graphloc, gaugeloc))
+    #savefig('Gauge_%s_stage.eps' %gaugeloc)
     
     #Momentum plot
@@ -159 +159 @@
     ylabel('sqrt( uh^2 + vh^2 ) [m^2/s]')    
     #savefig('Gauge_%d_momentum' %k)
-    savefig('Gauge_%s_momentum' %myloc)
+    savefig('%sGauge_%s_momentum' %(graphloc, gaugeloc))
     
-    # raw_input('Next')
-
     #Bearing plot
     ion()
@@ -181 +179 @@
     ylabel(' atan(vh/uh) [degrees from North]')    
     #savefig('Gauge_%d_bearing' %k)
-    savefig('Gauge_%s_bearing' %myloc)
-    
-    # raw_input('Next')
+    savefig('%sGauge_%s_bearing' %(graphloc, gaugeloc))
 
     #Speed plot
@@ -195 +191 @@
     ylabel('sqrt( uh^2 + vh^2 ) / depth [m/s]')    
     #savefig('Gauge_%d_speed' %k)
-    savefig('Gauge_%s_speed' %myloc)
-    
-    # raw_input('Next')
-
-    whichone = '_%s' %myloc
+    savefig('%sGauge_%s_speed' %(graphloc, gaugeloc))
+
+    whichone = '_%s' %gaugeloc
     thisfile = project.gaugetimeseries+whichone+'.csv'
     fid = open(thisfile, 'w')