Changeset 1836

Timestamp:

Sep 15, 2005, 5:37:00 PM (20 years ago)

Author:

ole

Message:

Added points where timeseries is plottet

File:

• production/karratha_2005/get_timeseries.py (modified) (2 diffs)

production/karratha_2005/get_timeseries.py

@@ -13 +13 @@
 #timestep = 0.05
 
+gauges = [[60000, 83000],
+          [60000, 80000],
+          [60000, 75000],                    
+          [60000, 70000],          
+          [60000, 60000],
+          [60000, 50000],          
+          [60000, 40000],
+          [60000, 30000],          
+          [60000, 20000],
+          [60000, 10000],          
+          [60000, 0]]
+#gauge_names = ['SW', 'N', 'N2', ]
+ 
 
 #Read model output
-f = file_function(swwfile,
-                  quantities = 'stage',
-                  #interpolation_points = gauges,
-                  verbose = True)
+from caching import cache
+#f = file_function(swwfile,
+#                  quantities = 'stage',
+#                  interpolation_points = gauges,
+#                  verbose = True)
+
+f = cache(file_function, swwfile,
+          {'quantities': 'stage',
+           'interpolation_points': gauges,
+           'verbose': True},
+          dependencies = [swwfile],
+          verbose = True)
 
 print f
+model_time = f.T
 
 
-for k, name in enumerate(gauge_names):
+for k, g in enumerate(gauges):
     model_values = []
-    for i, t in enumerate(reference_time):
+    for i, t in enumerate(f.T):
         val = f(t, point_id = k)[0]
-        model.append(val)
+        model_values.append(val)
 
     from pylab import *
@@ -33 +55 @@
     hold(False)
     plot(model_time, model_values, '-b')
+    name = str(g)
     title('Gauge %s' %name)
     savefig(name, dpi = 300)