Changeset 6835
- Timestamp:
- Apr 17, 2009, 4:06:53 PM (15 years ago)
- Location:
- anuga_work/production/australia_ph2
- Files:
-
- 15 edited
-
Communities_model_tracking.xls (modified) (previous)
-
albany/export_results_max.py (modified) (3 diffs)
-
bamaga/export_results_max.py (modified) (4 diffs)
-
broome/export_results_max.py (modified) (4 diffs)
-
bundaberg/export_results_max.py (modified) (3 diffs)
-
cairns/export_results_max.py (modified) (3 diffs)
-
cooktown/export_results_max.py (modified) (3 diffs)
-
derby/export_results_max.py (modified) (4 diffs)
-
esperance/export_results_max.py (modified) (3 diffs)
-
geraldton/export_results_max.py (modified) (4 diffs)
-
lockhart_river/export_results_max.py (modified) (3 diffs)
-
lucinda/export_results_max.py (modified) (3 diffs)
-
onslow/export_results_max.py (modified) (4 diffs)
-
perth/export_results_max.py (modified) (4 diffs)
-
rockhampton/export_results_max.py (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/australia_ph2/albany/export_results_max.py
r6579 r6835 22 22 directory = project.output_folder 23 23 24 time_dir1 = '20090318_171220_run_final_0_27319_1330_Tb_mhingee' 25 time_dir2 = '20090318_172003_run_final_0_64344_1330_Tb_mhingee' 26 time_dir3 = '20090318_172103_run_final_0_46697_1330_Tb_mhingee' 27 time_dir4 = '20090318_175318_run_final_1e-05_46697_1330_Tb_mhingee' 28 time_dir5 = '20090320_144313_run_final_1e-05_27319_1330_Tb_mhingee' 29 time_dir6 = '20090320_144404_run_final_1e-05_64344_1330_Tb_mhingee' 24 time_dir1 = '20090408_152157_run_final_0_27319_1330_Tb__mhingee' 25 time_dir2 = '20090408_201701_run_final_0_64344_1330_Tb__mhingee' 26 time_dir3 = '20090409_010309_run_final_0_46697_1330_Tb__mhingee' 30 27 31 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6] 32 28 time_dirs = [time_dir1, time_dir2, time_dir3] 33 29 34 30 cellsize = 250 35 31 ##cellsize = 5 36 32 37 #timestep = None # None means no timestep!38 timestep = 033 timestep = None # None means no timestep! 34 #timestep = 0 39 35 40 36 ###### … … 56 52 57 53 # one or more key strings from var_equations above 58 var = [' elevation','stage','speed']54 var = ['stage'] 59 55 60 56 ###### … … 85 81 for time_dir in time_dirs: 86 82 87 name1 = join(directory, time_dir, project.scenario_name) 88 #name2 = join(directory, time_dir, project.scenario_name)+'_time_39600_0' 89 #name3 = join(directory, time_dir, project.scenario_name)+'_time_79200_0' 83 name = join(directory, time_dir, project.scenario_name) 90 84 91 names = [name1] #, name2]#, name3] 92 93 # asc_name = [] 94 95 for name in names: 96 97 outname = name + '_' + which_area + '_' + which_var 98 quantityname = var_equations[which_var] 85 outname = name + '_' + which_area + '_' + which_var 86 quantityname = var_equations[which_var] 99 87 100 101 102 103 104 105 106 107 108 109 110 111 112 88 print 'start sww2dem: time_dir=%s' % time_dir 89 90 sww2dem(name, basename_out = outname, 91 quantity = quantityname, 92 timestep = timestep, 93 cellsize = cellsize, 94 easting_min = easting_min, 95 easting_max = easting_max, 96 northing_min = northing_min, 97 northing_max = northing_max, 98 reduction = max, 99 verbose = True, 100 format = 'asc') 113 101 114 ## asc_name.append(outname + '.asc') 115 ## 116 ## maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc' 117 ## 118 ## print 'max asc outname ', maxasc_outname 119 ## print 'asc_name ', str(asc_name) 120 ## 121 ## MaxAsc(maxasc_outname, asc_name) 102 -
anuga_work/production/australia_ph2/bamaga/export_results_max.py
r6742 r6835 22 22 directory = project.output_folder 23 23 24 time_dir1 = '20090306_131720_run_final_0_64448_1889_Tb_mhingee' 25 time_dir2 = '20090306_165354_run_final_0_58331_1889_Tb_mhingee' 26 time_dir3 = '20090306_165528_run_final_0_27347_1889_Tb_mhingee' 27 ##time_dir4 = '' 28 ##time_dir5 = '' 29 ##time_dir6 = '' 30 ## 31 time_dirs = [time_dir1, time_dir2, time_dir3] #time_dir2]# , time_dir4, time_dir5, time_dir6] 32 24 time_dir1 = '20090408_152655_run_final_0_31975_2000_Tb__mhingee' 25 time_dir2 = '20090409_123239_run_final_0_51390_2000_Tb__mhingee' 26 time_dir3 = '20090410_110658_run_final_0_63735_2000_Tb__mhingee' 27 28 time_dirs = [time_dir1, time_dir2, time_dir3] 33 29 34 30 cellsize = 250 35 31 ##cellsize = 5 36 32 37 #timestep = None # None means no timestep!38 timestep = 033 timestep = None # None means no timestep! 34 #timestep = 0 39 35 40 36 ###### … … 56 52 57 53 # one or more key strings from var_equations above 58 var = [' elevation','stage']54 var = ['stage'] 59 55 60 56 ###### … … 86 82 87 83 name1 = join(directory, time_dir, project.scenario_name) 88 #name2 = join(directory, time_dir, project.scenario_name)+'_time_39600_0' 89 #name3 = join(directory, time_dir, project.scenario_name)+'_time_79200_0' 84 name2 = join(directory, time_dir, project.scenario_name)+'_time_37920_0' 90 85 91 names = [name1] #, name2]#, name3] 92 93 # asc_name = [] 94 86 names = [name1, name2] 87 95 88 for name in names: 96 89 97 90 outname = name + '_' + which_area + '_' + which_var 98 91 quantityname = var_equations[which_var] … … 112 105 format = 'asc') 113 106 114 asc_name.append(outname + '.asc') 115 116 maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc' 117 118 print 'max asc outname ', maxasc_outname 119 print 'asc_name ', str(asc_name) 120 121 MaxAsc(maxasc_outname, asc_name) 107 -
anuga_work/production/australia_ph2/broome/export_results_max.py
r6339 r6835 22 22 directory = project.output_folder 23 23 24 ##time_dir1 = '20081217_101310_run_final_0_27283_250m_none_dt_kvanputt' 25 ##time_dir2 = '20081217_115336_run_final_0_27283_250m_none_dp_kvanputt' 26 ###time_dir1 = '20080924_123626_run_final_0_27283_250m_all_kvanputt' # This uses the 250m bathymetry with all interior polygons 27 ###time_dir2 = '20080912_154716_run_final_0_27283_alpha0.1_kvanputt' # This uses original bathyemetry data 28 ###time_dir3 = '20080924_123601_run_final_0_27283_250m_none_kvanputt' # This uses the 250m bathymetry without any interior polygons 29 ##time_dirs = [time_dir1, time_dir2]#, time_dir3] 24 time_dir1 = '20090408_152311_run_final_0_27337_1661_Tb__mhingee' 25 time_dir2 = '20090409_070404_run_final_0_70336_1661_Tb__mhingee' 26 time_dir3 = '20090409_224602_run_final_0_70869_1661_Tb__mhingee' 30 27 31 time_dir1 = '20090130_165532_run_final_0.6_27255_extend_dt_kvanputt' 32 ##time_dir2 = '20081211_162311_run_final_0_27255_alpha0.1_kvanputt' 33 ##time_dir3 = '20081211_162346_run_final_0_68693_alpha0.1_kvanputt' 34 ##time_dir4 = '20081211_162433_run_final_0.6_68693_alpha0.1_kvanputt' 35 ##time_dir5 = '20081211_162656_run_final_0.6_27283_alpha0.1_kvanputt' 36 ##time_dir6 = '20081211_162744_run_final_0_27283_alpha0.1_kvanputt' 37 ## 38 time_dirs = [time_dir1] #5, time_dir4]# , time_dir4, time_dir5, time_dir6] 39 28 time_dirs = [time_dir1, time_dir2, time_dir3] 40 29 41 30 cellsize = 250 42 31 ##cellsize = 5 43 32 44 #timestep = None # None means no timestep!45 timestep = 033 timestep = None # None means no timestep! 34 #timestep = 0 46 35 47 36 ###### … … 63 52 64 53 # one or more key strings from var_equations above 65 var = [' elevation']54 var = ['stage'] 66 55 67 56 ###### … … 93 82 94 83 name1 = join(directory, time_dir, project.scenario_name) 95 #name2 = join(directory, time_dir, project.scenario_name)+'_time_39600_0' 96 #name3 = join(directory, time_dir, project.scenario_name)+'_time_79200_0' 84 name2 = join(directory, time_dir, project.scenario_name)+'_time_31380_0' 97 85 98 names = [name1] #, name2]#, name3] 99 100 # asc_name = [] 101 86 names = [name1, name2] 87 102 88 for name in names: 103 89 104 90 outname = name + '_' + which_area + '_' + which_var 105 91 quantityname = var_equations[which_var] … … 119 105 format = 'asc') 120 106 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) 107 -
anuga_work/production/australia_ph2/bundaberg/export_results_max.py
r6743 r6835 22 22 directory = project.output_folder 23 23 24 time_dir1 = '20090306_131720_run_final_0_64448_1889_Tb_mhingee' 25 time_dir2 = '20090306_165354_run_final_0_58331_1889_Tb_mhingee' 26 time_dir3 = '20090306_165528_run_final_0_27347_1889_Tb_mhingee' 27 ##time_dir4 = '' 28 ##time_dir5 = '' 29 ##time_dir6 = '' 30 ## 31 time_dirs = [time_dir1, time_dir2, time_dir3] #time_dir2]# , time_dir4, time_dir5, time_dir6] 32 24 time_dir1 = '20090408_153233_run_final_0_31756_3274_Tb__mhingee' 25 time_dir2 = '20090409_015134_run_final_0_51469_3274_Tb__mhingee' 26 time_dir3 = '20090409_112941_run_final_0_62840_3274_Tb__mhingee' 27 28 time_dirs = [time_dir1, time_dir2, time_dir3] 33 29 34 30 cellsize = 250 35 31 ##cellsize = 5 36 32 37 #timestep = None # None means no timestep!38 timestep = 033 timestep = None # None means no timestep! 34 #timestep = 0 39 35 40 36 ###### … … 56 52 57 53 # one or more key strings from var_equations above 58 var = [' elevation','stage']54 var = ['stage'] 59 55 60 56 ###### … … 85 81 for time_dir in time_dirs: 86 82 87 name1 = join(directory, time_dir, project.scenario_name) 88 #name2 = join(directory, time_dir, project.scenario_name)+'_time_39600_0' 89 #name3 = join(directory, time_dir, project.scenario_name)+'_time_79200_0' 83 name = join(directory, time_dir, project.scenario_name) 90 84 91 names = [name1] #, name2]#, name3] 92 93 # asc_name = [] 94 95 for name in names: 96 97 outname = name + '_' + which_area + '_' + which_var 98 quantityname = var_equations[which_var] 85 outname = name + '_' + which_area + '_' + which_var 86 quantityname = var_equations[which_var] 99 87 100 101 102 103 104 105 106 107 108 109 110 111 112 88 print 'start sww2dem: time_dir=%s' % time_dir 89 90 sww2dem(name, basename_out = outname, 91 quantity = quantityname, 92 timestep = timestep, 93 cellsize = cellsize, 94 easting_min = easting_min, 95 easting_max = easting_max, 96 northing_min = northing_min, 97 northing_max = northing_max, 98 reduction = max, 99 verbose = True, 100 format = 'asc') 113 101 114 asc_name.append(outname + '.asc') 115 116 maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc' 117 118 print 'max asc outname ', maxasc_outname 119 print 'asc_name ', str(asc_name) 120 121 MaxAsc(maxasc_outname, asc_name) 102 -
anuga_work/production/australia_ph2/cairns/export_results_max.py
r6744 r6835 22 22 directory = project.output_folder 23 23 24 time_dir1 = '20090306_131720_run_final_0_64448_1889_Tb_mhingee' 25 time_dir2 = '20090306_165354_run_final_0_58331_1889_Tb_mhingee' 26 time_dir3 = '20090306_165528_run_final_0_27347_1889_Tb_mhingee' 27 ##time_dir4 = '' 28 ##time_dir5 = '' 29 ##time_dir6 = '' 30 ## 31 time_dirs = [time_dir1, time_dir2, time_dir3] #time_dir2]# , time_dir4, time_dir5, time_dir6] 32 24 time_dir1 = '20090408_153044_run_final_0_31975_2172_Tb__mhingee' 25 time_dir2 = '20090408_223355_run_final_0_51469_2172_Tb__mhingee' 26 time_dir3 = '20090409_055521_run_final_0_63737_2172_Tb__mhingee' 27 28 time_dirs = [time_dir1, time_dir2, time_dir3] 33 29 34 30 cellsize = 250 35 31 ##cellsize = 5 36 32 37 #timestep = None # None means no timestep!38 timestep = 033 timestep = None # None means no timestep! 34 #timestep = 0 39 35 40 36 ###### … … 56 52 57 53 # one or more key strings from var_equations above 58 var = [' elevation','stage']54 var = ['stage'] 59 55 60 56 ###### … … 85 81 for time_dir in time_dirs: 86 82 87 name1 = join(directory, time_dir, project.scenario_name) 88 #name2 = join(directory, time_dir, project.scenario_name)+'_time_39600_0' 89 #name3 = join(directory, time_dir, project.scenario_name)+'_time_79200_0' 83 name = join(directory, time_dir, project.scenario_name) 90 84 91 names = [name1] #, name2]#, name3] 92 93 # asc_name = [] 94 95 for name in names: 96 97 outname = name + '_' + which_area + '_' + which_var 98 quantityname = var_equations[which_var] 85 outname = name + '_' + which_area + '_' + which_var 86 quantityname = var_equations[which_var] 99 87 100 101 102 103 104 105 106 107 108 109 110 111 112 88 print 'start sww2dem: time_dir=%s' % time_dir 89 90 sww2dem(name, basename_out = outname, 91 quantity = quantityname, 92 timestep = timestep, 93 cellsize = cellsize, 94 easting_min = easting_min, 95 easting_max = easting_max, 96 northing_min = northing_min, 97 northing_max = northing_max, 98 reduction = max, 99 verbose = True, 100 format = 'asc') 113 101 114 asc_name.append(outname + '.asc') 115 116 maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc' 117 118 print 'max asc outname ', maxasc_outname 119 print 'asc_name ', str(asc_name) 120 121 MaxAsc(maxasc_outname, asc_name) 102 -
anuga_work/production/australia_ph2/cooktown/export_results_max.py
r6745 r6835 22 22 directory = project.output_folder 23 23 24 time_dir1 = '20090306_131720_run_final_0_64448_1889_Tb_mhingee' 25 time_dir2 = '20090306_165354_run_final_0_58331_1889_Tb_mhingee' 26 time_dir3 = '20090306_165528_run_final_0_27347_1889_Tb_mhingee' 27 ##time_dir4 = '' 28 ##time_dir5 = '' 29 ##time_dir6 = '' 30 ## 31 time_dirs = [time_dir1, time_dir2, time_dir3] #time_dir2]# , time_dir4, time_dir5, time_dir6] 32 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 28 time_dirs = [time_dir1, time_dir2, time_dir3] 33 29 34 30 cellsize = 250 35 31 ##cellsize = 5 36 32 37 #timestep = None # None means no timestep!38 timestep = 033 timestep = None # None means no timestep! 34 #timestep = 0 39 35 40 36 ###### … … 56 52 57 53 # one or more key strings from var_equations above 58 var = [' elevation','stage']54 var = ['stage'] 59 55 60 56 ###### … … 85 81 for time_dir in time_dirs: 86 82 87 name1 = join(directory, time_dir, project.scenario_name) 88 #name2 = join(directory, time_dir, project.scenario_name)+'_time_39600_0' 89 #name3 = join(directory, time_dir, project.scenario_name)+'_time_79200_0' 83 name = join(directory, time_dir, project.scenario_name) 90 84 91 names = [name1] #, name2]#, name3] 92 93 # asc_name = [] 94 95 for name in names: 96 97 outname = name + '_' + which_area + '_' + which_var 98 quantityname = var_equations[which_var] 85 outname = name + '_' + which_area + '_' + which_var 86 quantityname = var_equations[which_var] 99 87 100 101 102 103 104 105 106 107 108 109 110 111 112 88 print 'start sww2dem: time_dir=%s' % time_dir 89 90 sww2dem(name, basename_out = outname, 91 quantity = quantityname, 92 timestep = timestep, 93 cellsize = cellsize, 94 easting_min = easting_min, 95 easting_max = easting_max, 96 northing_min = northing_min, 97 northing_max = northing_max, 98 reduction = max, 99 verbose = True, 100 format = 'asc') 113 101 114 asc_name.append(outname + '.asc') 115 116 maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc' 117 118 print 'max asc outname ', maxasc_outname 119 print 'asc_name ', str(asc_name) 120 121 MaxAsc(maxasc_outname, asc_name) 102 -
anuga_work/production/australia_ph2/derby/export_results_max.py
r6767 r6835 22 22 directory = project.output_folder 23 23 24 ##time_dir1 = '20081217_101310_run_final_0_27283_250m_none_dt_kvanputt' 25 ##time_dir2 = '20081217_115336_run_final_0_27283_250m_none_dp_kvanputt' 26 ###time_dir1 = '20080924_123626_run_final_0_27283_250m_all_kvanputt' # This uses the 250m bathymetry with all interior polygons 27 ###time_dir2 = '20080912_154716_run_final_0_27283_alpha0.1_kvanputt' # This uses original bathyemetry data 28 ###time_dir3 = '20080924_123601_run_final_0_27283_250m_none_kvanputt' # This uses the 250m bathymetry without any interior polygons 29 ##time_dirs = [time_dir1, time_dir2]#, time_dir3] 24 time_dir1 = '20090408_171910_run_final_0_27268_1708_Tb__kvanputt' 25 time_dir2 = '20090410_070302_run_final_0_70437_1708_Tb__kvanputt' 26 time_dir3 = '20090411_153927_run_final_0_70862_1708_Tb__kvanputt' 30 27 31 time_dir1 = '20090130_165532_run_final_0.6_27255_extend_dt_kvanputt' 32 ##time_dir2 = '20081211_162311_run_final_0_27255_alpha0.1_kvanputt' 33 ##time_dir3 = '20081211_162346_run_final_0_68693_alpha0.1_kvanputt' 34 ##time_dir4 = '20081211_162433_run_final_0.6_68693_alpha0.1_kvanputt' 35 ##time_dir5 = '20081211_162656_run_final_0.6_27283_alpha0.1_kvanputt' 36 ##time_dir6 = '20081211_162744_run_final_0_27283_alpha0.1_kvanputt' 37 ## 38 time_dirs = [time_dir1] #5, time_dir4]# , time_dir4, time_dir5, time_dir6] 39 28 time_dirs = [time_dir1, time_dir2, time_dir3] 40 29 41 30 cellsize = 250 42 31 ##cellsize = 5 43 32 44 #timestep = None # None means no timestep!45 timestep = 033 timestep = None # None means no timestep! 34 #timestep = 0 46 35 47 36 ###### … … 63 52 64 53 # one or more key strings from var_equations above 65 var = [' elevation']54 var = ['stage'] 66 55 67 56 ###### … … 93 82 94 83 name1 = join(directory, time_dir, project.scenario_name) 95 #name2 = join(directory, time_dir, project.scenario_name)+'_time_39600_0' 96 #name3 = join(directory, time_dir, project.scenario_name)+'_time_79200_0' 84 name2 = join(directory, time_dir, project.scenario_name)+'_time_22020_0' 85 name3 = join(directory, time_dir, project.scenario_name)+'_time_44040_0' 86 87 names = [name1, name2, name3] 97 88 98 names = [name1] #, name2]#, name3]99 100 # asc_name = []101 102 89 for name in names: 103 90 104 91 outname = name + '_' + which_area + '_' + which_var 105 92 quantityname = var_equations[which_var] … … 119 106 format = 'asc') 120 107 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) 108 -
anuga_work/production/australia_ph2/esperance/export_results_max.py
r6455 r6835 22 22 directory = project.output_folder 23 23 24 time_dir1 = '20090302_104336_run_final_0_27255_mhingee' 25 ##time_dir2 = '' 26 ##time_dir3 = '' 27 ##time_dir4 = '' 28 ##time_dir5 = '' 29 ##time_dir6 = '' 30 ## 31 time_dirs = [time_dir1] #time_dir2]# , time_dir4, time_dir5, time_dir6] 32 24 time_dir1 = '20090408_152320_run_final_0_27329_1699_Tb__mhingee' 25 time_dir2 = '20090408_205317_run_final_0_58367_1699_Tb__mhingee' 26 time_dir3 = '20090409_024308_run_final_0_64449_1699_Tb__mhingee' 27 28 time_dirs = [time_dir1, time_dir2, time_dir3] 33 29 34 30 cellsize = 250 35 31 ##cellsize = 5 36 32 37 #timestep = None # None means no timestep!38 timestep = 033 timestep = None # None means no timestep! 34 #timestep = 0 39 35 40 36 ###### … … 56 52 57 53 # one or more key strings from var_equations above 58 var = [' elevation','stage']54 var = ['stage'] 59 55 60 56 ###### … … 85 81 for time_dir in time_dirs: 86 82 87 name1 = join(directory, time_dir, project.scenario_name) 88 #name2 = join(directory, time_dir, project.scenario_name)+'_time_39600_0' 89 #name3 = join(directory, time_dir, project.scenario_name)+'_time_79200_0' 83 name = join(directory, time_dir, project.scenario_name) 90 84 91 names = [name1] #, name2]#, name3] 92 93 # asc_name = [] 94 95 for name in names: 96 97 outname = name + '_' + which_area + '_' + which_var 98 quantityname = var_equations[which_var] 85 outname = name + '_' + which_area + '_' + which_var 86 quantityname = var_equations[which_var] 99 87 100 101 102 103 104 105 106 107 108 109 110 111 112 88 print 'start sww2dem: time_dir=%s' % time_dir 89 90 sww2dem(name, basename_out = outname, 91 quantity = quantityname, 92 timestep = timestep, 93 cellsize = cellsize, 94 easting_min = easting_min, 95 easting_max = easting_max, 96 northing_min = northing_min, 97 northing_max = northing_max, 98 reduction = max, 99 verbose = True, 100 format = 'asc') 113 101 114 ## asc_name.append(outname + '.asc') 115 ## 116 ## maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc' 117 ## 118 ## print 'max asc outname ', maxasc_outname 119 ## print 'asc_name ', str(asc_name) 120 ## 121 ## MaxAsc(maxasc_outname, asc_name) 102 -
anuga_work/production/australia_ph2/geraldton/export_results_max.py
r6801 r6835 22 22 directory = project.output_folder 23 23 24 ##time_dir1 = '20081217_101310_run_final_0_27283_250m_none_dt_kvanputt' 25 ##time_dir2 = '20081217_115336_run_final_0_27283_250m_none_dp_kvanputt' 26 ###time_dir1 = '20080924_123626_run_final_0_27283_250m_all_kvanputt' # This uses the 250m bathymetry with all interior polygons 27 ###time_dir2 = '20080912_154716_run_final_0_27283_alpha0.1_kvanputt' # This uses original bathyemetry data 28 ###time_dir3 = '20080924_123601_run_final_0_27283_250m_none_kvanputt' # This uses the 250m bathymetry without any interior polygons 29 ##time_dirs = [time_dir1, time_dir2]#, time_dir3] 24 time_dir1 = '20090414_114540_run_final_0_64475_576_Tb__kvanputt' 25 time_dir2 = '20090415_014342_run_final_0_70921_576_Tb__kvanputt' 26 time_dir3 = '20090415_153319_run_final_0_26985_576_Tb__kvanputt' 30 27 31 time_dir1 = '20090130_165532_run_final_0.6_27255_extend_dt_kvanputt' 32 ##time_dir2 = '20081211_162311_run_final_0_27255_alpha0.1_kvanputt' 33 ##time_dir3 = '20081211_162346_run_final_0_68693_alpha0.1_kvanputt' 34 ##time_dir4 = '20081211_162433_run_final_0.6_68693_alpha0.1_kvanputt' 35 ##time_dir5 = '20081211_162656_run_final_0.6_27283_alpha0.1_kvanputt' 36 ##time_dir6 = '20081211_162744_run_final_0_27283_alpha0.1_kvanputt' 37 ## 38 time_dirs = [time_dir1] #5, time_dir4]# , time_dir4, time_dir5, time_dir6] 39 28 time_dirs = [time_dir1, time_dir2, time_dir3] 40 29 41 30 cellsize = 250 42 31 ##cellsize = 5 43 32 44 #timestep = None # None means no timestep!45 timestep = 033 timestep = None # None means no timestep! 34 #timestep = 0 46 35 47 36 ###### … … 63 52 64 53 # one or more key strings from var_equations above 65 var = [' elevation']54 var = ['stage'] 66 55 67 56 ###### … … 93 82 94 83 name1 = join(directory, time_dir, project.scenario_name) 95 #name2 = join(directory, time_dir, project.scenario_name)+'_time_39600_0' 96 #name3 = join(directory, time_dir, project.scenario_name)+'_time_79200_0' 84 name2 = join(directory, time_dir, project.scenario_name)+'_time_41700_0' 97 85 98 names = [name1] #, name2]#, name3] 99 100 # asc_name = [] 101 86 names = [name1, name2] 87 102 88 for name in names: 103 89 104 90 outname = name + '_' + which_area + '_' + which_var 105 91 quantityname = var_equations[which_var] … … 119 105 format = 'asc') 120 106 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) 107 -
anuga_work/production/australia_ph2/lockhart_river/export_results_max.py
r6746 r6835 22 22 directory = project.output_folder 23 23 24 time_dir1 = '20090306_131720_run_final_0_64448_1889_Tb_mhingee' 25 time_dir2 = '20090306_165354_run_final_0_58331_1889_Tb_mhingee' 26 time_dir3 = '20090306_165528_run_final_0_27347_1889_Tb_mhingee' 27 ##time_dir4 = '' 28 ##time_dir5 = '' 29 ##time_dir6 = '' 30 ## 31 time_dirs = [time_dir1, time_dir2, time_dir3] #time_dir2]# , time_dir4, time_dir5, time_dir6] 32 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 28 time_dirs = [time_dir1, time_dir2, time_dir3] 33 29 34 30 cellsize = 250 35 31 ##cellsize = 5 36 32 37 #timestep = None # None means no timestep!38 timestep = 033 timestep = None # None means no timestep! 34 #timestep = 0 39 35 40 36 ###### … … 56 52 57 53 # one or more key strings from var_equations above 58 var = [' elevation','stage']54 var = ['stage'] 59 55 60 56 ###### … … 85 81 for time_dir in time_dirs: 86 82 87 name1 = join(directory, time_dir, project.scenario_name) 88 #name2 = join(directory, time_dir, project.scenario_name)+'_time_39600_0' 89 #name3 = join(directory, time_dir, project.scenario_name)+'_time_79200_0' 83 name = join(directory, time_dir, project.scenario_name) 90 84 91 names = [name1] #, name2]#, name3] 92 93 # asc_name = [] 94 95 for name in names: 96 97 outname = name + '_' + which_area + '_' + which_var 98 quantityname = var_equations[which_var] 85 outname = name + '_' + which_area + '_' + which_var 86 quantityname = var_equations[which_var] 99 87 100 101 102 103 104 105 106 107 108 109 110 111 112 88 print 'start sww2dem: time_dir=%s' % time_dir 89 90 sww2dem(name, basename_out = outname, 91 quantity = quantityname, 92 timestep = timestep, 93 cellsize = cellsize, 94 easting_min = easting_min, 95 easting_max = easting_max, 96 northing_min = northing_min, 97 northing_max = northing_max, 98 reduction = max, 99 verbose = True, 100 format = 'asc') 113 101 114 asc_name.append(outname + '.asc') 115 116 maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc' 117 118 print 'max asc outname ', maxasc_outname 119 print 'asc_name ', str(asc_name) 120 121 MaxAsc(maxasc_outname, asc_name) 102 -
anuga_work/production/australia_ph2/lucinda/export_results_max.py
r6747 r6835 22 22 directory = project.output_folder 23 23 24 time_dir1 = '20090306_131720_run_final_0_64448_1889_Tb_mhingee' 25 time_dir2 = '20090306_165354_run_final_0_58331_1889_Tb_mhingee' 26 time_dir3 = '20090306_165528_run_final_0_27347_1889_Tb_mhingee' 27 ##time_dir4 = '' 28 ##time_dir5 = '' 29 ##time_dir6 = '' 30 ## 31 time_dirs = [time_dir1, time_dir2, time_dir3] #time_dir2]# , time_dir4, time_dir5, time_dir6] 32 24 time_dir1 = '20090408_153121_run_final_0_31982_2313_Tb__mhingee' 25 time_dir2 = '20090408_202614_run_final_0_51311_2313_Tb__mhingee' 26 time_dir3 = '20090409_022558_run_final_0_63735_2313_Tb__mhingee' 27 28 time_dirs = [time_dir1, time_dir2, time_dir3] 33 29 34 30 cellsize = 250 35 31 ##cellsize = 5 36 32 37 #timestep = None # None means no timestep!38 timestep = 033 timestep = None # None means no timestep! 34 #timestep = 0 39 35 40 36 ###### … … 56 52 57 53 # one or more key strings from var_equations above 58 var = [' elevation','stage']54 var = ['stage'] 59 55 60 56 ###### … … 85 81 for time_dir in time_dirs: 86 82 87 name1 = join(directory, time_dir, project.scenario_name) 88 #name2 = join(directory, time_dir, project.scenario_name)+'_time_39600_0' 89 #name3 = join(directory, time_dir, project.scenario_name)+'_time_79200_0' 83 name = join(directory, time_dir, project.scenario_name) 90 84 91 names = [name1] #, name2]#, name3] 92 93 # asc_name = [] 94 95 for name in names: 96 97 outname = name + '_' + which_area + '_' + which_var 98 quantityname = var_equations[which_var] 85 outname = name + '_' + which_area + '_' + which_var 86 quantityname = var_equations[which_var] 99 87 100 101 102 103 104 105 106 107 108 109 110 111 112 88 print 'start sww2dem: time_dir=%s' % time_dir 89 90 sww2dem(name, basename_out = outname, 91 quantity = quantityname, 92 timestep = timestep, 93 cellsize = cellsize, 94 easting_min = easting_min, 95 easting_max = easting_max, 96 northing_min = northing_min, 97 northing_max = northing_max, 98 reduction = max, 99 verbose = True, 100 format = 'asc') 113 101 114 asc_name.append(outname + '.asc') 115 116 maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc' 117 118 print 'max asc outname ', maxasc_outname 119 print 'asc_name ', str(asc_name) 120 121 MaxAsc(maxasc_outname, asc_name) 102 -
anuga_work/production/australia_ph2/onslow/export_results_max.py
r6802 r6835 22 22 directory = project.output_folder 23 23 24 ##time_dir1 = '20081217_101310_run_final_0_27283_250m_none_dt_kvanputt' 25 ##time_dir2 = '20081217_115336_run_final_0_27283_250m_none_dp_kvanputt' 26 ###time_dir1 = '20080924_123626_run_final_0_27283_250m_all_kvanputt' # This uses the 250m bathymetry with all interior polygons 27 ###time_dir2 = '20080912_154716_run_final_0_27283_alpha0.1_kvanputt' # This uses original bathyemetry data 28 ###time_dir3 = '20080924_123601_run_final_0_27283_250m_none_kvanputt' # This uses the 250m bathymetry without any interior polygons 29 ##time_dirs = [time_dir1, time_dir2]#, time_dir3] 24 time_dir1 = '20090414_121216_run_final_0_27279_711_Tb__kvanputt' 25 time_dir2 = '20090415_001215_run_final_0_70371_711_Tb__kvanputt' 26 time_dir3 = '20090415_114548_run_final_0_27168_711_Tb__kvanputt' 30 27 31 time_dir1 = '20090130_165532_run_final_0.6_27255_extend_dt_kvanputt' 32 ##time_dir2 = '20081211_162311_run_final_0_27255_alpha0.1_kvanputt' 33 ##time_dir3 = '20081211_162346_run_final_0_68693_alpha0.1_kvanputt' 34 ##time_dir4 = '20081211_162433_run_final_0.6_68693_alpha0.1_kvanputt' 35 ##time_dir5 = '20081211_162656_run_final_0.6_27283_alpha0.1_kvanputt' 36 ##time_dir6 = '20081211_162744_run_final_0_27283_alpha0.1_kvanputt' 37 ## 38 time_dirs = [time_dir1] #5, time_dir4]# , time_dir4, time_dir5, time_dir6] 39 28 time_dirs = [time_dir1, time_dir2, time_dir3] 40 29 41 30 cellsize = 250 42 31 ##cellsize = 5 43 32 44 #timestep = None # None means no timestep!45 timestep = 033 timestep = None # None means no timestep! 34 #timestep = 0 46 35 47 36 ###### … … 63 52 64 53 # one or more key strings from var_equations above 65 var = [' elevation']54 var = ['stage'] 66 55 67 56 ###### … … 93 82 94 83 name1 = join(directory, time_dir, project.scenario_name) 95 #name2 = join(directory, time_dir, project.scenario_name)+'_time_39600_0' 96 #name3 = join(directory, time_dir, project.scenario_name)+'_time_79200_0' 84 name2 = join(directory, time_dir, project.scenario_name)+'_time_44400_0' 97 85 98 names = [name1] #, name2]#, name3] 99 100 # asc_name = [] 101 86 names = [name1, name2] 87 102 88 for name in names: 103 89 104 90 outname = name + '_' + which_area + '_' + which_var 105 91 quantityname = var_equations[which_var] … … 119 105 format = 'asc') 120 106 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) 107 -
anuga_work/production/australia_ph2/perth/export_results_max.py
r6429 r6835 22 22 directory = project.output_folder 23 23 24 time_dir1 = '20090130_165532_run_final_0.6_27255_extend_dt_kvanputt' 25 ##time_dir2 = '' 26 ##time_dir3 = '' 27 ##time_dir4 = '' 28 ##time_dir5 = '' 29 ##time_dir6 = '' 30 ## 31 time_dirs = [time_dir1] #5, time_dir4]# , time_dir4, time_dir5, time_dir6] 32 24 time_dir1 = '20090414_110305_run_final_0_26978_803_Tb__kvanputt' 25 time_dir2 = '20090415_090717_run_final_0_64479_803_Tb__kvanputt' 26 time_dir3 = '20090416_050507_run_final_0_70921_803_Tb__kvanputt' 27 28 time_dirs = [time_dir1, time_dir2, time_dir3] 33 29 34 30 cellsize = 250 35 31 ##cellsize = 5 36 32 37 #timestep = None # None means no timestep!38 timestep = 033 timestep = None # None means no timestep! 34 #timestep = 0 39 35 40 36 ###### … … 56 52 57 53 # one or more key strings from var_equations above 58 var = [' elevation']54 var = ['stage'] 59 55 60 56 ###### … … 86 82 87 83 name1 = join(directory, time_dir, project.scenario_name) 88 #name2 = join(directory, time_dir, project.scenario_name)+'_time_39600_0' 89 #name3 = join(directory, time_dir, project.scenario_name)+'_time_79200_0' 84 name2 = join(directory, time_dir, project.scenario_name)+'_time_546600_0' 90 85 91 names = [name1] #, name2]#, name3] 92 93 # asc_name = [] 94 86 names = [name1, name2] 87 95 88 for name in names: 96 89 97 90 outname = name + '_' + which_area + '_' + which_var 98 91 quantityname = var_equations[which_var] … … 112 105 format = 'asc') 113 106 114 ## asc_name.append(outname + '.asc') 115 ## 116 ## maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc' 117 ## 118 ## print 'max asc outname ', maxasc_outname 119 ## print 'asc_name ', str(asc_name) 120 ## 121 ## MaxAsc(maxasc_outname, asc_name) 107 -
anuga_work/production/australia_ph2/rockhampton/export_results_max.py
r6748 r6835 22 22 directory = project.output_folder 23 23 24 time_dir1 = '20090306_131720_run_final_0_64448_1889_Tb_mhingee' 25 time_dir2 = '20090306_165354_run_final_0_58331_1889_Tb_mhingee' 26 time_dir3 = '20090306_165528_run_final_0_27347_1889_Tb_mhingee' 27 ##time_dir4 = '' 28 ##time_dir5 = '' 29 ##time_dir6 = '' 30 ## 31 time_dirs = [time_dir1, time_dir2, time_dir3] #time_dir2]# , time_dir4, time_dir5, time_dir6] 32 24 time_dir1 = '20090408_153436_run_final_0_31939_3163_Tb__mhingee' 25 time_dir2 = '20090409_043830_run_final_0_51314_3163_Tb__mhingee' 26 time_dir3 = '20090409_192933_run_final_0_63737_3163_Tb__mhingee' 27 28 time_dirs = [time_dir1, time_dir2, time_dir3] 33 29 34 30 cellsize = 250 35 31 ##cellsize = 5 36 32 37 #timestep = None # None means no timestep!38 timestep = 033 timestep = None # None means no timestep! 34 #timestep = 0 39 35 40 36 ###### … … 56 52 57 53 # one or more key strings from var_equations above 58 var = [' elevation','stage']54 var = ['stage'] 59 55 60 56 ###### … … 86 82 87 83 name1 = join(directory, time_dir, project.scenario_name) 88 #name2 = join(directory, time_dir, project.scenario_name)+'_time_39600_0' 89 #name3 = join(directory, time_dir, project.scenario_name)+'_time_79200_0' 84 name2 = join(directory, time_dir, project.scenario_name)+'_time_54300_0' 90 85 91 names = [name1] #, name2]#, name3] 92 93 # asc_name = [] 94 86 names = [name1, name2] 87 95 88 for name in names: 96 89 97 90 outname = name + '_' + which_area + '_' + which_var 98 91 quantityname = var_equations[which_var] … … 112 105 format = 'asc') 113 106 114 asc_name.append(outname + '.asc') 115 116 maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc' 117 118 print 'max asc outname ', maxasc_outname 119 print 'asc_name ', str(asc_name) 120 121 MaxAsc(maxasc_outname, asc_name) 107
Note: See TracChangeset
for help on using the changeset viewer.
