Changeset 4856 for anuga_work/production/onslow_2006/export_results.py

Timestamp:

Nov 27, 2007, 4:18:39 PM (17 years ago)

Author:

sexton

Message:

rename plot_polygons_points to original name of plot_polygons

File:

• anuga_work/production/onslow_2006/export_results.py (modified) (6 diffs)

anuga_work/production/onslow_2006/export_results.py

@@ -1 @@
-import project, os
+import project_grad, os
 import sys
 
@@ -11 +11 @@
 #time_dir = '' # HAT 1 in Pt Hedland 10000 yr
 #time_dir = '' # MSL 1 in Pt Hedland 10000 yr
-#time_dir = '20070530_082920_run_final_1.5_onslow_nbartzis' # HAT 1 in Onslow 10000 yr
+#time_dir = '20070608_060316_run_final_1.5_onslow_nbartzis' # HAT 1 in Onslow 10000 yr
 #time_dir = '20070531_002753_run_final_0.0_onslow_nbartzis' # MSL 1 in Onslow 10000 yr
-time_dir = '20070608_062811_run_final_1.5_exmouth_nbartzis' # HAT 1 in Exmouth 10000 yr
+#time_dir = '20070608_062811_run_final_1.5_exmouth_nbartzis' # HAT 1 in Exmouth 10000 yr
 #time_dir = '20070605_234905_run_final_0.0_exmouth_nbartzis' # MSL 1 in Exmouth 10000 yr
 #time_dir = '20060704_063005' # 2006 simulation HAT Mw9
+#time_dir = '20070619_042140_run_final_1.5_exmouth_nbartzis' # exmouth HAT longer time
+#time_dir = '20070904_235118_run_final_1.5_jsexton_exmouth_original' # revised model original data
+#time_dir = '20070905_052311_run_final_1.5_jsexton_exmouth_revised' # revised model original + survey data
+time_dir = '20070906_061418_run_final_1.5_jsexton_exmouth_original' # refined again
+#time_dir = '20070907_055936_run_final_1.5_jsexton_exmouth_revised' # refined again
 cellsize = 25
 #cellsize = 150
 timestep = None
-directory = project.outputdir
-#name = directory + time_dir + sep + project.scenario_name
-name = directory+time_dir+sep+project.scenario_name
+directory = project_grad.output_dir
+#name = directory + time_dir + sep + project_grad.scenario_name
+name = directory+time_dir+sep+project_grad.scenario_name
+#name = directory+time_dir+sep+'onslow_time_35760_0'
 #name = directory+time_dir+sep+'source'
 
 from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf, dem2pts
 
-elevgrid = directory+time_dir+sep+'onslow_elevation_town'
+#elevgrid = directory+time_dir+sep+'onslow_elevation_town'
 
-convert_dem_from_ascii2netcdf(elevgrid, use_cache=True, verbose=True)
+#convert_dem_from_ascii2netcdf(elevgrid, use_cache=True, verbose=True)
 
-dem2pts(elevgrid,run_up_elev=True,run_up_loc=True)
+#dem2pts(elevgrid,run_up_elev=True,run_up_loc=True)
 
-"""
 is_parallel = False
 #is_parallel = True
@@ -40 +45 @@
 #var = [2,3,4] # depth and speed
 #var = [2] # depth
-var = [0] 
+var = [4] 
 
 for which_var in var:
@@ -48 +53 @@
 
     if which_var == 1:  # Absolute Momentum
-        outname = name + '_momentum_i1'
+        outname = name + '_momentum'
         quantityname = '(xmomentum**2 + ymomentum**2)**0.5'  
 
@@ -73 +78 @@
                         timestep = timestep,
                         cellsize = cellsize,      
-                        easting_min = project.e_min_area,
-                        easting_max = project.e_max_area,
-                        northing_min = project.n_min_area,
-                        northing_max = project.n_max_area,        
+                        easting_min = project_grad.e_min_area,
+                        easting_max = project_grad.e_max_area,
+                        northing_min = project_grad.n_min_area,
+                        northing_max = project_grad.n_max_area,        
                         reduction = max, 
                         verbose = True,
@@ -86 +91 @@
                     timestep = timestep,
                     cellsize = cellsize,      
-                    easting_min = project.e_min_area,
-                    easting_max = project.e_max_area,
-                    northing_min = project.n_min_area,
-                    northing_max = project.n_max_area,        
+                    easting_min = project_grad.e_min_area,
+                    easting_max = project_grad.e_max_area,
+                    northing_min = project_grad.n_min_area,
+                    northing_max = project_grad.n_max_area,        
                     reduction = max, 
                     verbose = True,
                     format = 'asc')
 
-"""