Changeset 2326 for production/sydney_2006

Timestamp:

Feb 2, 2006, 5:01:56 PM (19 years ago)

Author:

sexton

Message:

Fixes to velocity calculation

Location:

production/sydney_2006

Files:

• export_results.py (modified) (1 diff)
• get_timeseries.py (modified) (4 diffs)

production/sydney_2006/export_results.py

@@ -27 +27 @@
     #quantityname = '((xmomentum/(stage-elevation))**2 + (ymomentum/(stage-elevation))**2)**0.5'  #Speed
     #quantityname = 'xmomentum/(stage-elevation)'  #Speed
-    quantityname = '(xmomentum**2 + ymomentum**2 + 1.e-30)**0.5/(stage-elevation)'  #Speed
+    quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-30)'  #Speed
 
     

production/sydney_2006/get_timeseries.py

@@ -5 +5 @@
 import project
 from pyvolution.util import file_function
-from pyvolution.coordinate_transforms.redfearn import degminsec2decimal_degrees, redfearn
+#from pyvolution.coordinate_transforms.redfearn import degminsec2decimal_degrees, redfearn
+from coordinate_transforms.redfearn import degminsec2decimal_degrees, redfearn
 from pylab import *
 
@@ -61 +62 @@
     elevations = []
     momenta = []
+    velocity = []
 
     max_depth = 0
-    max_momentum = 0    
+    max_momentum = 0
+    max_velocity = 0
     for i, t in enumerate(f.T): # T is a list of times
         #if tmin < t < tmax:
@@ -73 +76 @@
 
         m = sqrt(uh*uh + vh*vh)   #Absolute momentum
+        vel = sqrt(uh*uh + vh*vh) / (w-z + 1.e-30) #Absolute velocity
+        print vel
+        #dep = w-z
+        #vel = sqrt(uh*uh + vh*vh) / dep #Absolute velocity
         
         model_time.append(t)        
         stages.append(w)
         elevations.append(z)  #Should be constant over time
-        momenta.append(m) 
+        momenta.append(m)
+        velocity.append(vel)
 
         if w-z > max_depth:
             max_depth = w-z
         if m > max_momentum:            
-            max_momentum = m    
+            max_momentum = m
+        if vel > max_velocity:
+            max_velocity = vel
+            print 'max speed', max_velocity
 
     
@@ -127 +138 @@
     
     raw_input('Next')
+
+    #Speed plot
+    ion()
+    hold(False)
+    plot(model_time, velocity, '-r')
+    title(name)
+
+    title('%s (velocity)' %name)
+    xlabel('time [s]')
+    ylabel('sqrt( uh^2 + vh^2 ) / depth [m/s]')    
+    savefig('Gauge_%d_speed' %k)
+    #savefig('Gauge_%s_speed' %myloc)
+    
+    raw_input('Next')