Changeset 1893 for production/karratha_2005/get_timeseries.py

Timestamp:

Oct 11, 2005, 3:21:38 PM (19 years ago)

Author:

ole

Message:

New gauge points

File:

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

production/karratha_2005/get_timeseries.py

@@ -4 +4 @@
 import project
 from pyvolution.util import file_function
+from pyvolution.coordinate_transforms.redfearn import degminsec2decimal_degrees, redfearn
 
-swwfile = project.outputdir + project.basename + '.sww'
+#swwfile = project.outputdir + project.basename + '.sww'
+swwfile = project.outputname + '_0.0tide' + '.sww'
+#swwfile = project.outputdir + 'karratha_100m_12m.sww'
 
-#gauges = [[0.000, 1.696]]  #Boundary
-#gauges += [[4.521, 1.196],  [4.521, 1.696],  [4.521, 2.196]] #Ch 5-7-9
-#gauge_names = ['Boundary', 'ch5', 'ch7', 'ch9']
-#finaltime = 22.5
-#timestep = 0.05
+#MOST Boundary
+y = project.north
+x = (project.west + project.east)/2
 
-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', ]
+z, easting, northing  = redfearn(y, x)
+
+
+#From Neil
+gauges = [[easting, northing],
+          [470882, 7717952],
+          [469390, 7715426],
+          [469214, 7714979],
+          [468899, 7715177],
+          [469038, 7715725],
+          [470285, 7717470]]
+          
+
+    
+#gauge_names = ['MOST',
+#               'Gauge_1', 'Gauge_2', 'Gauge_3',
+#               'Gauge_4', 'Gauge_5', 'Gauge_6']
  
 
@@ -45 +51 @@
 
 
+
 for k, g in enumerate(gauges):
     model_values = []
@@ -55 +62 @@
     hold(False)
     plot(model_time, model_values, '-b')
-    name = str(g)
-    title('Gauge %s' %name)
+    name = 'Gauge_%d' %k
+    title(name)
     savefig(name, dpi = 300)