Changeset 5592

Timestamp:

Aug 1, 2008, 1:33:26 PM (16 years ago)

Author:

kristy

Message:

 

Location:

anuga_work/production/perth

Files:

• export_results.py (modified) (1 diff)
• get_gauges.py (modified) (2 diffs)
• get_timeseries.py (modified) (2 diffs)
• project.py (modified) (3 diffs)

anuga_work/production/perth/export_results.py

@@ -10 +10 @@
 #time_dir = '20080714_095738_run_final_0.6_exmouth_kvanputt'
 #time_dir = '20080714_134935_run_final_0_exmouth_kvanputt'
-time_dir = '20080718_112228_run_final_0_exmouth_kvanputt'
+#time_dir = '20080718_112228_run_final_0_exmouth_kvanputt'
+#time_dir = '20080725_173911_run_final_0.6_polyline_alpha0.1_kvanputt'
+#time_dir = '20080728_112519_run_final_0.6_polyline_alpha0.2_kvanputt'
+time_dir = '20080728_154123_run_final_0.6_polyline_alpha0.15_kvanputt'
 
 cellsize = 25

anuga_work/production/perth/get_gauges.py

@@ -36 +36 @@
     for i, name in enumerate(quantity_names):
         quantities[name] = fid.variables[name][:]
-        if inside_points is not None:
-            quantities[name] = take(quantities[name],indices,1)
-        else:
-            msg = 'No gauges found inside polygon'
-            raise msg
-            
+##        if inside_points is not None:
+##            quantities[name] = take(quantities[name],indices,1)
+##        else:
+##            msg = 'No gauges found inside polygon'
+##            raise msg
+##            
     for j in range(len(x)-1):
         locx=int(x[j])
@@ -63 +63 @@
 polygon=project.poly_all
 points,quantities,elevation,time=get_sts_gauge_data(os.path.join(project.boundaries_dir,project.scenario_name),polygon,verbose=False)
-
+print quantities['stage'][0,:]
+print len(points), len(elevation), len(quantities['stage'][0,:])
 assert len(points)==len(elevation)==len(quantities['stage'][0,:])

anuga_work/production/perth/get_timeseries.py

@@ -17 +17 @@
 
 #timestamp='20080724_121200_run_trial_0.6_polyline_alpha0.1_kvanputt'
-timestamp='20080724_161830_run_final_0.6_polyline_alpha0.1_kvanputt'
+#timestamp='20080724_161830_run_final_0.6_polyline_alpha0.1_kvanputt'
+timestamp='20080725_173911_run_final_0.6_polyline_alpha0.1_kvanputt'
 
 filename=project.output_dir+timestamp+sep+project.scenario_name+'.sww'
@@ -23 +24 @@
 
 sww2csv_gauges(filename,
-                project.gauges_dir_name,
+                #project.gauges_dir_name,
+                project.gauges_dir_name2,
                 quantities = ['stage','speed','depth','elevation'],
                 verbose=True)

anuga_work/production/perth/project.py

@@ -33 +33 @@
 #Making assumptions about the location of scenario data
 state = 'western_australia'
+#scenario_name = 'perth_44unitsources'
 scenario_name = 'perth'
 scenario = 'perth_tsunami_scenario'
@@ -38 +39 @@
 tide = 0.0 #0.6
 
-alpha = 0.1
+alpha = 0.15
 friction=0.01
 starttime=0
@@ -152 +153 @@
 #gauges
 gauge_name = 'perth.csv'
+gauge_name2 = 'thinned_MGA50+50.csv'
 
 gauges_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'gauges'+sep
 beach_gauges = gauges_dir + 'beach_gauges.csv'
 gauges_dir_name = gauges_dir + gauge_name
+gauges_dir_name2 = gauges_dir + gauge_name2
 
 buildings_filename = gauges_dir + 'Perth_resA.csv'