Changeset 6603 for anuga_work/production/pt_hedland_2008

Timestamp:

Mar 24, 2009, 3:36:38 PM (16 years ago)

Author:

jgriffin

Message:

 

Location:

anuga_work/production/pt_hedland_2008

Files:

• export_results.py (modified) (2 diffs)
• export_results_all.py (modified) (3 diffs)
• export_results_max.py (modified) (1 diff)
• project.py (modified) (3 diffs)

anuga_work/production/pt_hedland_2008/export_results.py

@@ -5 +5 @@
 from os import sep
 
-time_dir = '20081014_163150_run_final_0_27283_alpha0.1_jgriffin'
+time_dir = '20081210_163235_run_final_3.6_27283_alpha0.1_jgriffin'
 
-cellsize = 30
+cellsize = 20
 #cellsize = 150
 timestep = 0
@@ -24 +24 @@
 #var = [3,4] # depth and speed
 #var = [2] # depth
-var = [0] #stage
+#var = [0] #stage
 #var = [0,4] # stage and elevation
+var = [4] #elevation
 
 for which_var in var:

anuga_work/production/pt_hedland_2008/export_results_all.py

@@ -20 +20 @@
 directory = project.output_dir
 
-time_dir1 = '20081015_140641_run_final_0_27283_alpha0.1_jgriffin'
+time_dir1 = '20081211_095614_run_final_0.0_27283_alpha0.1_jgriffin'
 time_dirs = [time_dir1]
 
@@ -53 +53 @@
 for time_dir in time_dirs:
     name1 = directory+time_dir+sep+project.scenario_name
-#    name2 = directory+time_dir+sep+project.scenario_name+'_time_39900_0' #need to get assistance on how to make this into anything
-    names = [name1]
+    name2 = directory+time_dir+sep+project.scenario_name+'_time_24720_0' #need to get assistance on how to make this into anything
+    name3 = directory+time_dir+sep+project.scenario_name+'_time_49440_0'
+    name4 = directory+time_dir+sep+project.scenario_name+'_time_74160_0' 
+    names = [name1,name2,name3,name4]
     
     for name in names:
@@ -94 +96 @@
                             verbose = True,
                             format = 'asc')
+
+                asc_name.append(outname + '.asc')
+    
+            maxasc_outname = directory+time_dir+sep+project_250m.scenario_name+'_'+which_area+'_'+which_var+'_max.asc'
+
+            print 'max asc outname ', maxasc_outname
+            print 'asc_name ', str(asc_name)
+            
+            MaxAsc(maxasc_outname, asc_name)

anuga_work/production/pt_hedland_2008/export_results_max.py

@@ -20 +20 @@
 directory = project_250m.output_dir
 
-time_dir1 = '20081125_135917_run_final_0_27283_250m_none_lfountai'
+time_dir1 = '20090108_120325_run_final_3.6_27283_alpha0.1_jgriffin'
 
 time_dirs = [time_dir1]#, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]

anuga_work/production/pt_hedland_2008/project.py

@@ -43 +43 @@
 # Model specific parameters. One or all can be changed each time the
 # run_scenario script is executed
-tide = 0                #0.6
+tide = 3.6              #3.6
 #event_number = 27255   # linked to hazard map
 event_number = 27283
@@ -176 +176 @@
 poly_ocean = read_polygon(polygons_dir+'initial_conditions_ocean.csv')
 
+#Island, to set the initial stage/water to be offcoast only
+#poly_island1 = read_polygon(polygons_dir+'initial_conditions_island1.csv')
+
+#Island, to set the initial stage/water to be offcoast only
+poly_island2 = read_polygon(polygons_dir+'initial_conditions_island2.csv')
+
 # Initial bounding polygon for data clipping 
 poly_all = read_polygon(polygons_dir+'poly_all.csv')
-res_poly_all = 100000*res_factor
+#res_poly_all = 100000*res_factor
+res_poly_all = 200000*res_factor
 
 # Area of Interest 1 (pt_hedland)
@@ -186 +193 @@
 # Area of Significance 1 (pt_hedland)
 poly_aos1 = read_polygon(polygons_dir+'area_of_significance.csv')
-res_aos1 = 1000*res_factor
+res_aos1 = 2500*res_factor
 
 # Shallow water 1
 poly_sw1 = read_polygon(polygons_dir+'shallow_water.csv')
-res_sw1 = 25000*res_factor
+res_sw1 = 50000*res_factor