Changeset 6929
- Timestamp:
- Apr 28, 2009, 6:06:14 PM (16 years ago)
- Location:
- anuga_work/production/australia_ph2
- Files:
-
- 2 added
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/australia_ph2/cooktown/export_results_max.py
r6835 r6929 22 22 directory = project.output_folder 23 23 24 time_dir1 = '20090408_153032_run_final_0_31977_2065_Tb__mhingee'25 time_dir2 = '20090408_214931_run_final_0_51469_2065_Tb__mhingee'26 time_dir3 = '20090409_050811_run_final_0_63735_2065_Tb__mhingee'24 ##time_dir1 = '20090408_153032_run_final_0_31977_2065_Tb__mhingee' 25 ##time_dir2 = '20090408_214931_run_final_0_51469_2065_Tb__mhingee' 26 ##time_dir3 = '20090409_050811_run_final_0_63735_2065_Tb__mhingee' 27 27 28 time_dirs = [time_dir1, time_dir2, time_dir3] 28 time_dir1 = '20090421_135741_run_final_0_31977_2065_Tb_internal_lfountai' 29 30 time_dirs = [time_dir1]#, time_dir2, time_dir3] 29 31 30 32 cellsize = 250 -
anuga_work/production/australia_ph2/darwin/project.py
r6814 r6929 50 50 51 51 # index is only used when wave = Tb 52 index = 1829 # index from the PTHA - Y2000 0. 257m52 index = 1829 # index from the PTHA - Y2000 0.592m; using events picked with PTHA point 1833 (0.642m) 53 53 wave = 'Tb' # Bf (sts wave) Tb (index wave) 54 54 55 55 56 internal_polygon = False56 internal_polygon = True 57 57 58 58 #------------------------------------------------------------------------------- -
anuga_work/production/australia_ph2/lockhart_river/Arc_asc2raster_GDA94z50.py
r6926 r6929 25 25 output_dir="anuga\\outputs\\" 26 26 27 time_dir1 = '20090408_142152_run_final_0_51316_1994_Tb__mhingee' 28 time_dir2 = '20090408_152731_run_final_0_51316_1994_Tb__mhingee' 29 time_dir3 = '20090408_200614_run_final_0_62786_1994_Tb__mhingee' 27 time_dir1 = '20090408_152731_run_final_0_51316_1994_Tb__mhingee' 28 ##time_dir2 = '20090408_200614_run_final_0_62786_1994_Tb__mhingee' 30 29 31 time_dirs = [time_dir1, time_dir2, time_dir3] 30 ##time_dir1 = '20090421_135544_run_final_0_51316_1994_Tb_internal_lfountai' 31 32 time_dirs = [time_dir1]#, time_dir2] 32 33 33 34 for time_dir in time_dirs: … … 38 39 ## ocean = scenario_dir + "map_work\\\port_hedland.gdb\\outlines\\initial_conditions_ocean" 39 40 40 print 'Process: Create File GDB'41 gp.CreateFileGDB_management(folder, "raster")41 ## print 'Process: Create File GDB' 42 ## gp.CreateFileGDB_management(folder, "raster") 42 43 43 44 gp.Workspace = raster_gbd … … 52 53 53 54 generate_filename = [] 54 input_ascii = glob.glob(folder + '* stage.asc')55 input_ascii = glob.glob(folder + '*elevation.asc') 55 56 56 57 for infile in input_ascii: -
anuga_work/production/australia_ph2/lockhart_river/export_results_max.py
r6835 r6929 22 22 directory = project.output_folder 23 23 24 time_dir1 = '20090408_142152_run_final_0_51316_1994_Tb__mhingee'25 time_dir2 = '20090408_152731_run_final_0_51316_1994_Tb__mhingee'26 time_dir3 = '20090408_200614_run_final_0_62786_1994_Tb__mhingee'24 ##time_dir1 = '20090408_142152_run_final_0_51316_1994_Tb__mhingee' 25 ##time_dir2 = '20090408_152731_run_final_0_51316_1994_Tb__mhingee' 26 ##time_dir3 = '20090408_200614_run_final_0_62786_1994_Tb__mhingee' 27 27 28 time_dirs = [time_dir1, time_dir2, time_dir3] 28 time_dir1 = '20090421_135544_run_final_0_51316_1994_Tb_internal_lfountai' 29 30 time_dirs = [time_dir1]#, time_dir2, time_dir3] 29 31 30 32 cellsize = 250 -
anuga_work/production/australia_ph2/wyndham/export_results_max.py
r6539 r6929 22 22 directory = project.output_folder 23 23 24 time_dir1 = '20090 316_120226_run_final_0_17898_1833_Bf_mhingee'25 ##time_dir2 = '20090 316_120344_run_final_0_64977_1833_Bf_mhingee'26 ##time_dir3 = '20090 316_120457_run_final_0_70906_1833_Bf_mhingee'27 ##time_dir4 = ' 20090316_120558_run_final_0_71261_1833_Bf_mhingee'24 time_dir1 = '20090415_094009_run_final_0_17918_1823_Tb__kvanputt' 25 ##time_dir2 = '20090418_182759_run_final_0_64962_1823_Tb__kvanputt' 26 ##time_dir3 = '20090421_145541_run_final_0_70920_1823_Tb__kvanputt' 27 ##time_dir4 = '' 28 28 ##time_dir5 = '' 29 29 ##time_dir6 = '' 30 30 ##time_dir7 = '' 31 31 ##time_dir8 = '' 32 time_dirs = [time_dir1] #5, time_dir4]# , time_dir4, time_dir5, time_dir6] 32 33 time_dirs = [time_dir1]#, time_dir2, time_dir3]# , time_dir4, time_dir5, time_dir6] 33 34 34 35 … … 36 37 ##cellsize = 5 37 38 38 #timestep = None # None means no timestep!39 timestep = 039 timestep = None # None means no timestep! 40 #timestep = 0 40 41 41 42 ###### … … 57 58 58 59 # one or more key strings from var_equations above 59 var = [' elevation']60 var = ['stage'] 60 61 61 62 ###### … … 87 88 88 89 name1 = join(directory, time_dir, project.scenario_name) 89 #name2 = join(directory, time_dir, project.scenario_name)+'_time_39600_0' 90 #name3 = join(directory, time_dir, project.scenario_name)+'_time_79200_0' 91 92 names = [name1] #, name2]#, name3] 90 name2 = join(directory, time_dir, project.scenario_name)+'_time_8460_0' 91 name3 = join(directory, time_dir, project.scenario_name)+'_time_16920_0' 92 name4 = join(directory, time_dir, project.scenario_name)+'_time_25380_0' 93 name5 = join(directory, time_dir, project.scenario_name)+'_time_33840_0' 94 name6 = join(directory, time_dir, project.scenario_name)+'_time_42300_0' 95 name7 = join(directory, time_dir, project.scenario_name)+'_time_50760_0' 96 name8 = join(directory, time_dir, project.scenario_name)+'_time_59220_0' 97 98 names = [name1, name2, name3, name4, name5, name6, name7, name8] 93 99 94 #asc_name = []100 asc_name = [] 95 101 96 102 for name in names: … … 113 119 format = 'asc') 114 120 115 ##asc_name.append(outname + '.asc')116 ##117 ##maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'118 ## 119 ##print 'max asc outname ', maxasc_outname120 ##print 'asc_name ', str(asc_name)121 ##122 ##MaxAsc(maxasc_outname, asc_name)121 asc_name.append(outname + '.asc') 122 123 maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc' 124 125 print 'max asc outname ', maxasc_outname 126 print 'asc_name ', str(asc_name) 127 128 MaxAsc(maxasc_outname, asc_name)
Note: See TracChangeset
for help on using the changeset viewer.