Changeset 7002

Timestamp:

May 8, 2009, 2:22:03 PM (16 years ago)

Author:

myall

Message:

 

Location:

anuga_work/production/australia_ph2

Files:

• albany/export_results.py (modified) (1 diff)
• albany/export_results_max.py (modified) (1 diff)
• bamaga/project.py (modified) (1 diff)
• bamaga/run_model.py (modified) (1 diff)
• gold_coast/Arc_asc2raster_GDA94z50.py (modified) (3 diffs)
• gold_coast/project.py (modified) (1 diff)
• hobart/export_results.py (modified) (1 diff)
• hobart/export_results_max.py (modified) (2 diffs)
• portland/export_results.py (modified) (2 diffs)
• portland/export_results_max.py (modified) (4 diffs)
• strahan/export_results.py (modified) (1 diff)
• strahan/export_results_max.py (modified) (1 diff)

anuga_work/production/australia_ph2/albany/export_results.py

@@ -22 +22 @@
 directory = project.output_folder
 
-time_dir1 = '20090408_201701_run_final_0_64344_1330_Tb__mhingee'
-time_dir2 = '20090421_165204_run_final_0_64344_1330_Tb_internal_mhingee'
+time_dir1 = '20090507_164822_run_final_30_0_64344_1330_Tb__kvanputt'
+time_dir2 = '20090507_115302_run_final_30_0_64344_1330_Tb_internal_kvanputt'
 ##time_dir3 = ''
 ##time_dir4 = ''

anuga_work/production/australia_ph2/albany/export_results_max.py

@@ -22 +22 @@
 directory = project.output_folder
 
-time_dir1 = '20090421_165204_run_final_0_64344_1330_Tb_internal_mhingee'
-##time_dir2 = ''
-##time_dir3 = ''
-##time_dir4 = ''
+time_dir1 = '20090507_115302_run_final_30_0_64344_1330_Tb_internal_kvanputt'
+time_dir2 = '20090507_115824_run_final_30_0_27319_1330_Tb__kvanputt'
+time_dir3 = '20090507_164822_run_final_30_0_64344_1330_Tb__kvanputt'
+time_dir4 = '20090507_214813_run_final_30_0_46697_1330_Tb__kvanputt'
 ##time_dir5 = ''
 ##time_dir6 = ''
 
-time_dirs = [time_dir1]#, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
+time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4]#, time_dir5, time_dir6]
   
 

anuga_work/production/australia_ph2/bamaga/project.py

@@ -32 +32 @@
     event_number = int(sys.argv[1])
 else:    
-    event_number = 31975    # the event number or the mux file name
+    event_number = 63735    # the event number or the mux file name
 
 event_number_list = [31975, 51390, 63735] # To piggy back multiple events

anuga_work/production/australia_ph2/bamaga/run_model.py

@@ -92 +92 @@
 # Boundary tags refer to project.landward_boundary
 # 4 points equals 5 segments start at N
-boundary_tags={'back': (num_ocean_segments+1),
+boundary_tags={'back': [num_ocean_segments+1],
                'side': [num_ocean_segments,
                         num_ocean_segments+2,

anuga_work/production/australia_ph2/gold_coast/Arc_asc2raster_GDA94z50.py

@@ -25 +25 @@
 output_dir="anuga\\outputs\\"
 
-time_dir1 = '20090408_171713_run_final_0_50883_3727_Tb__kvanputt'
-time_dir2 = '20090409_054747_run_final_0_51318_3727_Tb__kvanputt'
-time_dir3 = '20090409_175154_run_final_0_51453_3727_Tb__kvanputt'
+time_dir1 = '20090505_145250_run_final_0_31854_3727_Tb__lfountai'
+time_dir2 = '20090506_045043_run_final_0_64479_3727_Tb__lfountai'
+time_dir3 = '20090506_190156_run_final_0_50883_3727_Tb__lfountai'
 
-##time_dir1 = '20090421_150448_run_final_0_50883_3727_Tb_internal_mhingee'
+##time_dir1 = '20090508_112905_run_final_0_31854_3727_Tb_internal_mhingee'
 
-time_dirs = [time_dir1]#, time_dir2, time_dir3]
+time_dirs = [time_dir1, time_dir2, time_dir3]
 
 for time_dir in time_dirs:
@@ -41 +41 @@
 ##    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
@@ -55 +55 @@
 
     generate_filename = []
-    input_ascii = glob.glob(folder + '*elevation.asc')
+    input_ascii = glob.glob(folder + '*max.asc')
 
     for infile in input_ascii:

anuga_work/production/australia_ph2/gold_coast/project.py

@@ -54 +54 @@
 
 
-internal_polygon = False
+internal_polygon = True
 
 #-------------------------------------------------------------------------------

anuga_work/production/australia_ph2/hobart/export_results.py

@@ -22 +22 @@
 directory = project.output_folder
 
-time_dir1 = '20090424_121952_run_final_0_51204_2462_Tb__kvanputt'
-time_dir2 = '20090424_123250_run_final_0_51204_2462_Tb_internal_kvanputt'
+time_dir1 = '20090507_111303_run_final_30_0_51204_2462_Tb__kvanputt'
+time_dir2 = '20090507_111251_run_final_30_0_51204_2462_Tb_internal_kvanputt'
 ##time_dir3 = ''
 ##time_dir4 = ''

anuga_work/production/australia_ph2/hobart/export_results_max.py

@@ -22 +22 @@
 directory = project.output_folder
 
-time_dir1 = '20090424_121952_run_final_0_51204_2462_Tb__kvanputt'
-time_dir2 = '20090424_123250_run_final_0_51204_2462_Tb_internal_kvanputt'
-time_dir3 = '20090424_122951_run_final_0_58212_2462_Tb__kvanputt'
-time_dir4 = '20090424_123104_run_final_0_64425_2462_Tb__kvanputt'
+time_dir1 = '20090507_111303_run_final_30_0_51204_2462_Tb__kvanputt'
+time_dir2 = '20090507_150745_run_final_30_0_58212_2462_Tb__kvanputt'
+time_dir3 = '20090507_185058_run_final_30_0_64425_2462_Tb__kvanputt'
+##time_dir4 = '20090507_111251_run_final_30_0_51204_2462_Tb_internal_kvanputt'
 
-time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4]
+time_dirs = [time_dir1, time_dir2, time_dir3]#, time_dir4]
   
 
@@ -82 +82 @@
         for time_dir in time_dirs:
 
-            name = join(directory, time_dir, project.scenario_name)
-                           
-            outname = name + '_' + which_area + '_' + which_var
-            quantityname = var_equations[which_var]
+            name1 = join(directory, time_dir, project.scenario_name)
+            name2 = join(directory, time_dir, project.scenario_name)+'_time_54840_0'
+##            name3 = join(directory, time_dir, project.scenario_name)+'_time_26640_0'
+##            name4 = join(directory, time_dir, project.scenario_name)+'_time_39960_0'
+##            name5 = join(directory, time_dir, project.scenario_name)+'_time_53280_0'
+##            name6 = join(directory, time_dir, project.scenario_name)+'_time_66600_0'
+##            name7 = join(directory, time_dir, project.scenario_name)+'_time_79920_0'
+            
+            names = [name1, name2]#, name3,name4,name5,name6,name7]
+      
+            asc_name = []   
+               
+            for name in names:
+                         
+                outname = name + '_' + which_area + '_' + which_var
+                quantityname = var_equations[which_var]
 
-            print 'start sww2dem: time_dir=%s' % time_dir
+                print 'start sww2dem: time_dir=%s' % time_dir
+                
+                sww2dem(name, basename_out = outname,
+                            quantity = quantityname,
+                            timestep = timestep,
+                            cellsize = cellsize,      
+                            easting_min = easting_min,
+                            easting_max = easting_max,
+                            northing_min = northing_min,
+                            northing_max = northing_max,        
+                            reduction = max, 
+                            verbose = True,
+                            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)
             
-            sww2dem(name, basename_out = outname,
-                        quantity = quantityname,
-                        timestep = timestep,
-                        cellsize = cellsize,      
-                        easting_min = easting_min,
-                        easting_max = easting_max,
-                        northing_min = northing_min,
-                        northing_max = northing_max,        
-                        reduction = max, 
-                        verbose = True,
-                        format = 'asc')
+            MaxAsc(maxasc_outname, asc_name)
 

anuga_work/production/australia_ph2/portland/export_results.py

@@ -22 +22 @@
 directory = project.output_folder
 
-time_dir1 = '20090409_001857_run_final_0_64322_1969_Tb__mhingee'
-time_dir2 = '20090421_153138_run_final_0_64322_1969_Tb_internal_mhingee'
+time_dir1 = '20090507_112926_run_final_30_0_27289_1969_Tb__kvanputt'
+##time_dir2 = ''
 ##time_dir3 = ''
 ##time_dir4 = ''
@@ -29 +29 @@
 ##time_dir6 = ''
 ##
-time_dirs = [time_dir1, time_dir2]# , time_dir4, time_dir5, time_dir6]
+time_dirs = [time_dir1]#, time_dir2]# , time_dir4, time_dir5, time_dir6]
   
 

anuga_work/production/australia_ph2/portland/export_results_max.py

@@ -22 +22 @@
 directory = project.output_folder
 
-time_dir1 = '20090421_153138_run_final_0_64322_1969_Tb_internal_mhingee'
-##time_dir2 = ''
-##time_dir3 = ''
+time_dir1 = '20090507_112926_run_final_30_0_27289_1969_Tb__kvanputt'
+time_dir2 = '20090507_113007_run_final_30_0_64322_1969_Tb_internal_kvanputt'
+time_dir3 = '20090507_155610_run_final_30_0_58358_1969_Tb_internal_kvanputt'
 ##time_dir4 = ''
 ##time_dir5 = ''
 ##time_dir6 = ''
 
-time_dirs = [time_dir1]#, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
+time_dirs = [time_dir1, time_dir2, time_dir3]#, time_dir4]#, time_dir5, time_dir6]
   
 
@@ -86 +86 @@
 
             name1 = join(directory, time_dir, project.scenario_name)
-##            name2 = join(directory, time_dir, project.scenario_name)+'_time_10740_0'
+            name2 = join(directory, time_dir, project.scenario_name)+'_time_44430_0'
 ##            name3 = join(directory, time_dir, project.scenario_name)+'_time_21480_0'
 ##            name4 = join(directory, time_dir, project.scenario_name)+'_time_32220_0'
@@ -92 +92 @@
 ##            name6 = join(directory, time_dir, project.scenario_name)+'_time_53700_0'
 
-            names = [name1]#, name2, name3, name4, name5, name6]
+            names = [name1, name2]#, name3, name4, name5, name6]
       
-##            asc_name = []   
+            asc_name = []   
                
             for name in names:
@@ -115 +115 @@
                             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/australia_ph2/strahan/export_results.py

@@ -22 +22 @@
 directory = project.output_folder
 
-time_dir1 = '20090408_152513_run_final_0_58337_2044_Tb__mhingee'
-time_dir2 = '20090421_151652_run_final_0_58337_2044_Tb_internal_mhingee'
+time_dir1 = '20090507_111151_run_final_30_0_58337_2044_Tb__kvanputt'
+time_dir2 = '20090507_111501_run_final_30_0_58337_2044_Tb_internal_kvanputt'
 ##time_dir3 = ''
 ##time_dir4 = ''

anuga_work/production/australia_ph2/strahan/export_results_max.py

@@ -22 +22 @@
 directory = project.output_folder
 
-time_dir1 = '20090421_151652_run_final_0_58337_2044_Tb_internal_mhingee'
-##time_dir2 = ''
-##time_dir3 = ''
-##time_dir4 = ''
+time_dir1 = '20090507_111151_run_final_30_0_58337_2044_Tb__kvanputt'
+time_dir2 = '20090507_111501_run_final_30_0_58337_2044_Tb_internal_kvanputt'
+time_dir3 = '20090507_122805_run_final_30_0_64214_2044_Tb_internal_kvanputt'
+time_dir4 = '20090507_134426_run_final_30_0_68779_2044_Tb_internal_kvanputt'
 ##time_dir5 = ''
 ##time_dir6 = ''
 
-time_dirs = [time_dir1]#, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
+time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4]#, time_dir5, time_dir6]