Changeset 6678

Timestamp:

Apr 1, 2009, 2:46:58 PM (15 years ago)

Author:

kristy

Message:

 

Location:

anuga_work/production/busselton/busselton_rerun

Files:

• Arc_asc2raster_GDA94z50.py (added)
• Arc_txt2feature_class.py (added)
• Arc_txt2feature_class_res.py (added)
• export_results_max.py (modified) (5 diffs)
• get_runup.py (added)
• project.py (modified) (2 diffs)
• run_building_inundation.py (added)

anuga_work/production/busselton/busselton_rerun/export_results_max.py

@@ -22 +22 @@
 directory = project.output_folder
 
-##time_dir1 = '20081217_101310_run_final_0_27283_250m_none_dt_kvanputt'
-##time_dir2 = '20081217_115336_run_final_0_27283_250m_none_dp_kvanputt'
-###time_dir1 = '20080924_123626_run_final_0_27283_250m_all_kvanputt' # This uses the 250m bathymetry with all interior polygons
-###time_dir2 = '20080912_154716_run_final_0_27283_alpha0.1_kvanputt' # This uses original bathyemetry data
-###time_dir3 = '20080924_123601_run_final_0_27283_250m_none_kvanputt' # This uses the 250m bathymetry without any interior polygons
-##time_dirs = [time_dir1, time_dir2]#, time_dir3]
 
-time_dir1 = '20090130_165532_run_final_0.6_27255_extend_dt_kvanputt'
-##time_dir2 = '20081211_162311_run_final_0_27255_alpha0.1_kvanputt'
-##time_dir3 = '20081211_162346_run_final_0_68693_alpha0.1_kvanputt'
-##time_dir4 = '20081211_162433_run_final_0.6_68693_alpha0.1_kvanputt'
-##time_dir5 = '20081211_162656_run_final_0.6_27283_alpha0.1_kvanputt'
-##time_dir6 = '20081211_162744_run_final_0_27283_alpha0.1_kvanputt'
-##
-time_dirs = [time_dir1] #5, time_dir4]# , time_dir4, time_dir5, time_dir6]
+time_dir1 = '20090324_134328_run_final_0.6_27255_Bt_kvanputt'
+time_dir2 = '20090324_140747_run_final_0.6_68693_Bt_kvanputt'
+time_dir3 = '20090327_152717_run_final_0.6_27283_Bt_kvanputt'
+time_dir4 = '20090327_160646_run_final_0_27255_Bt_kvanputt'
+time_dir5 = '20090327_161520_run_final_0_68693_Bt_kvanputt'
+time_dir6 = '20090330_120458_run_final_0_27283_Bt_kvanputt'
+
+time_dirs = [time_dir1,time_dir2, time_dir3, time_dir4, time_dir5]
   
 
-cellsize = 250
-##cellsize = 5
+##cellsize = 250
+cellsize = 20
 
 #timestep = None    # None means no timestep!
@@ -49 +43 @@
 ######
 
-#area = ['Bunbury', 'Busselton']  # strings must match keys in var_equations below
-area = ['All']      # 'All' means no special areas - the whole thing
+area = ['Bunbury', 'Busselton']  # strings must match keys in var_equations below
+#area = ['All']      # 'All' means no special areas - the whole thing
 
 ######
@@ -63 +57 @@
 
 # one or more key strings from var_equations above
-var = ['elevation']
+var = ['depth', 'speed']
 
 ######
@@ -93 +87 @@
 
             name1 = join(directory, time_dir, project.scenario_name)
-            #name2 = join(directory, time_dir, project.scenario_name)+'_time_39600_0'
-            #name3 = join(directory, time_dir, project.scenario_name)+'_time_79200_0'
+            name2 = join(directory, time_dir, project.scenario_name)+'_time_35520_0'
+            name3 = join(directory, time_dir, project.scenario_name)+'_time_71040_0'
 
-            names = [name1] #, name2]#, name3]
+            names = [name1, name2, name3]
       
-           # asc_name = []   
+            asc_name = []   
                
             for name in names:
@@ -119 +113 @@
                             format = 'asc')
 
-##                asc_name.append(outname + '.asc')
-##    
-##            maxasc_outname = join(directory, time_dir, project.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)
+                asc_name.append(outname + '.asc')
+    
+            maxasc_outname = join(directory, time_dir, project.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/busselton/busselton_rerun/project.py

@@ -26 +26 @@
 # Model specific parameters.
 # One or all can be changed each time the run_model script is executed
-tide = 0.6              # difference between MSL and HAT in metres
+tide = 0              # difference between MSL and HAT in metres
 zone = 50               # specify zone of model
-event_number = 68693    # the event number or the mux file name
+event_number = 27283    # the event number or the mux file name
 alpha = 0.1             # smoothing parameter for mesh
 friction=0.01           # manning's friction coefficient
@@ -50 +50 @@
 #-------------------------------------------------------------------------------
 
-output_comment = [setup, tide, event_number]
+output_comment = [setup, tide, event_number, 'Bt']
 
 #-------------------------------------------------------------------------------