Changeset 6929

Timestamp:

Apr 28, 2009, 6:06:14 PM (14 years ago)

Author:

myall

Message:

 

Location:

anuga_work/production/australia_ph2

Files:

• Communities_model_tracking.xls (modified) (previous)
• cooktown/Arc_asc2raster_GDA94z50.py (added)
• cooktown/export_results_max.py (modified) (1 diff)
• darwin/project.py (modified) (1 diff)
• lockhart_river/Arc_asc2raster_GDA94z50.py (modified) (3 diffs)
• lockhart_river/export_results.py (added)
• lockhart_river/export_results_max.py (modified) (1 diff)
• wyndham/export_results_max.py (modified) (5 diffs)

anuga_work/production/australia_ph2/cooktown/export_results_max.py

@@ -22 +22 @@
 directory = project.output_folder
 
-time_dir1 = '20090408_153032_run_final_0_31977_2065_Tb__mhingee'
-time_dir2 = '20090408_214931_run_final_0_51469_2065_Tb__mhingee'
-time_dir3 = '20090409_050811_run_final_0_63735_2065_Tb__mhingee'
+##time_dir1 = '20090408_153032_run_final_0_31977_2065_Tb__mhingee'
+##time_dir2 = '20090408_214931_run_final_0_51469_2065_Tb__mhingee'
+##time_dir3 = '20090409_050811_run_final_0_63735_2065_Tb__mhingee'
 
-time_dirs = [time_dir1, time_dir2, time_dir3]
+time_dir1 = '20090421_135741_run_final_0_31977_2065_Tb_internal_lfountai'
+
+time_dirs = [time_dir1]#, time_dir2, time_dir3]
 
 cellsize = 250

anuga_work/production/australia_ph2/darwin/project.py

@@ -50 +50 @@
 
 # index is only used when wave = Tb
-index = 1829            # index from the PTHA - Y2000 0.257m
+index = 1829            # index from the PTHA - Y2000 0.592m; using events picked with PTHA point 1833 (0.642m)
 wave = 'Tb'             # Bf (sts wave) Tb (index wave)
 
 
-internal_polygon = False
+internal_polygon = True
 
 #-------------------------------------------------------------------------------

anuga_work/production/australia_ph2/lockhart_river/Arc_asc2raster_GDA94z50.py

@@ -25 +25 @@
 output_dir="anuga\\outputs\\"
 
-time_dir1 = '20090408_142152_run_final_0_51316_1994_Tb__mhingee'
-time_dir2 = '20090408_152731_run_final_0_51316_1994_Tb__mhingee'
-time_dir3 = '20090408_200614_run_final_0_62786_1994_Tb__mhingee'
+time_dir1 = '20090408_152731_run_final_0_51316_1994_Tb__mhingee'
+##time_dir2 = '20090408_200614_run_final_0_62786_1994_Tb__mhingee'
 
-time_dirs = [time_dir1, time_dir2, time_dir3]
+##time_dir1 = '20090421_135544_run_final_0_51316_1994_Tb_internal_lfountai'
+
+time_dirs = [time_dir1]#, time_dir2]
 
 for time_dir in time_dirs:
@@ -38 +39 @@
 ##    ocean = scenario_dir + "map_work\\\port_hedland.gdb\\outlines\\initial_conditions_ocean"
 
-    print 'Process: Create File GDB'
-    gp.CreateFileGDB_management(folder, "raster")
+##    print 'Process: Create File GDB'
+##    gp.CreateFileGDB_management(folder, "raster")
 
     gp.Workspace = raster_gbd
@@ -52 +53 @@
 
     generate_filename = []
-    input_ascii = glob.glob(folder + '*stage.asc')
+    input_ascii = glob.glob(folder + '*elevation.asc')
 
     for infile in input_ascii:

anuga_work/production/australia_ph2/lockhart_river/export_results_max.py

@@ -22 +22 @@
 directory = project.output_folder
 
-time_dir1 = '20090408_142152_run_final_0_51316_1994_Tb__mhingee'
-time_dir2 = '20090408_152731_run_final_0_51316_1994_Tb__mhingee'
-time_dir3 = '20090408_200614_run_final_0_62786_1994_Tb__mhingee'
+##time_dir1 = '20090408_142152_run_final_0_51316_1994_Tb__mhingee'
+##time_dir2 = '20090408_152731_run_final_0_51316_1994_Tb__mhingee'
+##time_dir3 = '20090408_200614_run_final_0_62786_1994_Tb__mhingee'
 
-time_dirs = [time_dir1, time_dir2, time_dir3]
+time_dir1 = '20090421_135544_run_final_0_51316_1994_Tb_internal_lfountai'
+
+time_dirs = [time_dir1]#, time_dir2, time_dir3]
 
 cellsize = 250

anuga_work/production/australia_ph2/wyndham/export_results_max.py

@@ -22 +22 @@
 directory = project.output_folder
 
-time_dir1 = '20090316_120226_run_final_0_17898_1833_Bf_mhingee'
-##time_dir2 = '20090316_120344_run_final_0_64977_1833_Bf_mhingee'
-##time_dir3 = '20090316_120457_run_final_0_70906_1833_Bf_mhingee'
-##time_dir4 = '20090316_120558_run_final_0_71261_1833_Bf_mhingee'
+time_dir1 = '20090415_094009_run_final_0_17918_1823_Tb__kvanputt'
+##time_dir2 = '20090418_182759_run_final_0_64962_1823_Tb__kvanputt'
+##time_dir3 = '20090421_145541_run_final_0_70920_1823_Tb__kvanputt'
+##time_dir4 = ''
 ##time_dir5 = ''
 ##time_dir6 = ''
 ##time_dir7 = ''
 ##time_dir8 = ''
-time_dirs = [time_dir1] #5, time_dir4]# , time_dir4, time_dir5, time_dir6]
+
+time_dirs = [time_dir1]#, time_dir2, time_dir3]# , time_dir4, time_dir5, time_dir6]
   
 
@@ -36 +37 @@
 ##cellsize = 5
 
-#timestep = None    # None means no timestep!
-timestep = 0
+timestep = None    # None means no timestep!
+#timestep = 0
 
 ######
@@ -57 +58 @@
 
 # one or more key strings from var_equations above
-var = ['elevation']
+var = ['stage']
 
 ######
@@ -87 +88 @@
 
             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'
-
-            names = [name1] #, name2]#, name3]
+            name2 = join(directory, time_dir, project.scenario_name)+'_time_8460_0'
+            name3 = join(directory, time_dir, project.scenario_name)+'_time_16920_0'
+            name4 = join(directory, time_dir, project.scenario_name)+'_time_25380_0'
+            name5 = join(directory, time_dir, project.scenario_name)+'_time_33840_0'
+            name6 = join(directory, time_dir, project.scenario_name)+'_time_42300_0'
+            name7 = join(directory, time_dir, project.scenario_name)+'_time_50760_0'
+            name8 = join(directory, time_dir, project.scenario_name)+'_time_59220_0'
+            
+            names = [name1, name2, name3, name4, name5, name6, name7, name8]
       
-           # asc_name = []   
+            asc_name = []   
                
             for name in names:
@@ -113 +119 @@
                             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)