Changeset 6926


Ignore:
Timestamp:
Apr 28, 2009, 4:00:49 PM (16 years ago)
Author:
myall
Message:

exporting max stage for runs with internal polygons, and exporting elevation using export_results.py for this run, and the equivalent normal run

Location:
anuga_work/production/australia_ph2
Files:
17 added
18 edited

Legend:

Unmodified
Added
Removed

  • anuga_work/production/australia_ph2/adelaide/export_results_max.py

    r6845 r6926  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090414_104456_run_final_0_27346_1936_Tb__kvanputt'
    25 time_dir2 = '20090416_012314_run_final_0_58334_1936_Tb__kvanputt'
    26 time_dir3 = '20090417_082435_run_final_0_64322_1936_Tb__kvanputt'
     24time_dir1 = '20090421_153519_run_final_0_64322_1936_Tb_internal_mhingee'
     25##time_dir2 = '20090416_012314_run_final_0_58334_1936_Tb__kvanputt'
     26##time_dir3 = '20090417_082435_run_final_0_64322_1936_Tb__kvanputt'
    2727
    28 time_dirs = [time_dir1, time_dir2, time_dir3]
     28time_dirs = [time_dir1]#, time_dir2, time_dir3]
    2929
    3030cellsize = 250

  • anuga_work/production/australia_ph2/albany/export_results_max.py

    r6835 r6926  
    2222directory = project.output_folder
    2323
    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'
     24time_dir1 = '20090421_165204_run_final_0_64344_1330_Tb_internal_mhingee'
     25##time_dir2 = ''
     26##time_dir3 = ''
     27##time_dir4 = ''
     28##time_dir5 = ''
     29##time_dir6 = ''
    2730
    28 time_dirs = [time_dir1, time_dir2, time_dir3]
     31time_dirs = [time_dir1]#, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
     32 
    2933
    3034cellsize = 250
    … …  
    8185        for time_dir in time_dirs:
    8286
    83             name = join(directory, time_dir, project.scenario_name)
     87            name1 = join(directory, time_dir, project.scenario_name)
     88##            name2 = join(directory, time_dir, project.scenario_name)+'_time_28380_0'
     89##            name3 = join(directory, time_dir, project.scenario_name)+'_time_56760_0'
     90##            name4 = join(directory, time_dir, project.scenario_name)+'_time_32220_0'
     91##            name5 = join(directory, time_dir, project.scenario_name)+'_time_42960_0'
     92##            name6 = join(directory, time_dir, project.scenario_name)+'_time_53700_0'
    8493
    85             outname = name + '_' + which_area + '_' + which_var
    86             quantityname = var_equations[which_var]
     94            names = [name1]#, name2, name3]#, name4, name5, name6]
     95     
     96##            asc_name = []   
     97               
     98            for name in names:
     99                         
     100                outname = name + '_' + which_area + '_' + which_var
     101                quantityname = var_equations[which_var]
    87102
    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')
     103                print 'start sww2dem: time_dir=%s' % time_dir
     104               
     105                sww2dem(name, basename_out = outname,
     106                            quantity = quantityname,
     107                            timestep = timestep,
     108                            cellsize = cellsize,     
     109                            easting_min = easting_min,
     110                            easting_max = easting_max,
     111                            northing_min = northing_min,
     112                            northing_max = northing_max,       
     113                            reduction = max,
     114                            verbose = True,
     115                            format = 'asc')
    101116
    102          
     117##                asc_name.append(outname + '.asc')
     118##   
     119##            maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
     120##
     121##            print 'max asc outname ', maxasc_outname
     122##            print 'asc_name ', str(asc_name)
     123##           
     124##            MaxAsc(maxasc_outname, asc_name)

  • anuga_work/production/australia_ph2/broome/export_results_max.py

    r6845 r6926  
    2222directory = project.output_folder
    2323
    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'
     24time_dir1 = '20090421_174625_run_final_0_70869_1661_Tb_internal_mhingee'
    2725
    28 time_dirs = [time_dir1, time_dir2, time_dir3]
     26time_dirs = [time_dir1]#, time_dir2, time_dir3]
    2927
    3028cellsize = 250
    … …  
    8280
    8381            name1 = join(directory, time_dir, project.scenario_name)
    84             name2 = join(directory, time_dir, project.scenario_name)+'_time_31380_0'
     82            name2 = join(directory, time_dir, project.scenario_name)+'_time_31320_0'
     83##            name3 = join(directory, time_dir, project.scenario_name)+'_time_45480_0'
    8584
    86             names = [name1, name2]
     85            names = [name1, name2]#, name3]
    8786
    8887            asc_name = []
    … …  
    115114           
    116115            MaxAsc(maxasc_outname, asc_name)
    117              

  • anuga_work/production/australia_ph2/ceduna/export_results_max.py

    r6540 r6926  
    2222directory = project.output_folder
    2323
    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'
     24time_dir1 = '20090421_155143_run_final_0_27347_1884_Tb_internal_mhingee'
     25##time_dir2 = ''
     26##time_dir3 = ''
    2727##time_dir4 = ''
    2828##time_dir5 = ''
    2929##time_dir6 = ''
    30 ##
    31 time_dirs = [time_dir1, time_dir2, time_dir3] #time_dir2]# , time_dir4, time_dir5, time_dir6]
     30
     31time_dirs = [time_dir1]#, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
    3232 
    3333
    … …  
    3535##cellsize = 5
    3636
    37 #timestep = None    # None means no timestep!
    38 timestep = 0
     37timestep = None    # None means no timestep!
     38#timestep = 0
    3939
    4040######
    … …  
    5656
    5757# one or more key strings from var_equations above
    58 var = ['elevation','stage']
     58var = ['stage']
    5959
    6060######
    … …  
    8686
    8787            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'
     88            name2 = join(directory, time_dir, project.scenario_name)+'_time_30780_0'
     89##            name3 = join(directory, time_dir, project.scenario_name)+'_time_21480_0'
     90##            name4 = join(directory, time_dir, project.scenario_name)+'_time_32220_0'
     91##            name5 = join(directory, time_dir, project.scenario_name)+'_time_42960_0'
     92##            name6 = join(directory, time_dir, project.scenario_name)+'_time_53700_0'
    9093
    91             names = [name1] #, name2]#, name3]
     94            names = [name1, name2]#, name3, name4, name5, name6]
    9295     
    93            # asc_name = []   
     96            asc_name = []   
    9497               
    9598            for name in names:

  • anuga_work/production/australia_ph2/dampier/export_results_max.py

    r6797 r6926  
    2222directory = project.output_folder
    2323
    24 ##time_dir1 = '20090324_110442_run_final_0_70385_Bf_30km_kvanputt'
    25 ##time_dir2 = '20090324_123534_run_final_0_70385_Bf_60km_kvanputt'
    26 ##time_dir3 = '20090311_153009_run_final_0_70454_1203_Bf_jgriffin'
    27 time_dir4 = '20090323_165436_run_final_0_70385_1203_Tb_jgriffin'
    28 time_dirs = [time_dir4]#,time_dir2,time_dir3,time_dir4] 
     24time_dir1 = '20090422_102457_run_final_0_70844_1307_Tb_internal_mhingee'
     25
     26time_dirs = [time_dir1]#, time_dir2, time_dir3]
    2927
    3028cellsize = 250
    … …  
    5250
    5351# one or more key strings from var_equations above
    54 var = ['elevation']
     52var = ['stage']
    5553
    5654######
    … …  
    6967            northing_min = None
    7068            northing_max = None
    71 ##            easting_min = 600000
    72 ##            easting_max = 660000
    73 ##            northing_min = 7880000
    74 ##            northing_max = 7910000
    7569        else:
    7670            try:
    … …  
    8680
    8781            name1 = join(directory, time_dir, project.scenario_name)
    88             name2 = join(directory, time_dir, project.scenario_name)+'_time_13320_0'
    89 #            name3 = join(directory, time_dir, project.scenario_name)+'_time_26640_0'
    90 ##            name4 = join(directory, time_dir, project.scenario_name)+'_time_6660_0'
    91 ##            name5 = join(directory, time_dir, project.scenario_name)+'_time_8880_0'
    92 ####            name6 = join(directory, time_dir, project.scenario_name)+'_time_66600_0'
    93 ##            name7 = join(directory, time_dir, project.scenario_name)+'_time_79920_0'
     82            name2 = join(directory, time_dir, project.scenario_name)+'_time_22740_0'
     83            name3 = join(directory, time_dir, project.scenario_name)+'_time_45480_0'
     84
     85            names = [name1, name2, name3]
     86
     87            asc_name = []
    9488           
    95             names = [name1, name2,]#name3,name4,name5]#,name6,name7]
    96      
    97             asc_name = []   
    98                
    9989            for name in names:
    100                          
     90               
    10191                outname = name + '_' + which_area + '_' + which_var
    10292                quantityname = var_equations[which_var]

  • anuga_work/production/australia_ph2/derby/export_results_max.py

    r6845 r6926  
    2222directory = project.output_folder
    2323
    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'
     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'
    2727
    28 time_dirs = [time_dir1, time_dir2, time_dir3]
     28time_dir1 = '20090421_175503_run_final_0_27268_1708_Tb_internal_mhingee'
     29
     30time_dirs = [time_dir1]#, time_dir2, time_dir3]
    2931
    3032cellsize = 250

  • anuga_work/production/australia_ph2/esperance/export_results_max.py

    r6835 r6926  
    2222directory = project.output_folder
    2323
    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'
     24time_dir1 = '20090421_164728_run_final_0_64449_1699_Tb_internal_mhingee'
     25##time_dir2 = ''
     26##time_dir3 = ''
     27##time_dir4 = ''
     28##time_dir5 = ''
     29##time_dir6 = ''
    2730
    28 time_dirs = [time_dir1, time_dir2, time_dir3]
     31time_dirs = [time_dir1]#, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
     32 
    2933
    3034cellsize = 250
    … …  
    8185        for time_dir in time_dirs:
    8286
    83             name = join(directory, time_dir, project.scenario_name)
     87            name1 = join(directory, time_dir, project.scenario_name)
     88##            name2 = join(directory, time_dir, project.scenario_name)+'_time_28380_0'
     89##            name3 = join(directory, time_dir, project.scenario_name)+'_time_56760_0'
     90##            name4 = join(directory, time_dir, project.scenario_name)+'_time_32220_0'
     91##            name5 = join(directory, time_dir, project.scenario_name)+'_time_42960_0'
     92##            name6 = join(directory, time_dir, project.scenario_name)+'_time_53700_0'
    8493
    85             outname = name + '_' + which_area + '_' + which_var
    86             quantityname = var_equations[which_var]
     94            names = [name1]#, name2, name3]#, name4, name5, name6]
     95     
     96##            asc_name = []   
     97               
     98            for name in names:
     99                         
     100                outname = name + '_' + which_area + '_' + which_var
     101                quantityname = var_equations[which_var]
    87102
    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')
     103                print 'start sww2dem: time_dir=%s' % time_dir
     104               
     105                sww2dem(name, basename_out = outname,
     106                            quantity = quantityname,
     107                            timestep = timestep,
     108                            cellsize = cellsize,     
     109                            easting_min = easting_min,
     110                            easting_max = easting_max,
     111                            northing_min = northing_min,
     112                            northing_max = northing_max,       
     113                            reduction = max,
     114                            verbose = True,
     115                            format = 'asc')
    101116
    102          
     117##                asc_name.append(outname + '.asc')
     118##   
     119##            maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
     120##
     121##            print 'max asc outname ', maxasc_outname
     122##            print 'asc_name ', str(asc_name)
     123##           
     124##            MaxAsc(maxasc_outname, asc_name)

  • anuga_work/production/australia_ph2/eucla_motel/export_results_max.py

    r6455 r6926  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090225_153734_run_final_0_65371_mhingee'
     24time_dir1 = '20090421_164640_run_final_0_64449_1837_Tb_internal_mhingee'
    2525##time_dir2 = ''
    2626##time_dir3 = ''
    … …  
    2828##time_dir5 = ''
    2929##time_dir6 = ''
    30 ##
    31 time_dirs = [time_dir1] #time_dir2]# , time_dir4, time_dir5, time_dir6]
     30
     31time_dirs = [time_dir1]#, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
    3232 
    3333
    … …  
    3535##cellsize = 5
    3636
    37 #timestep = None    # None means no timestep!
    38 timestep = 0
     37timestep = None    # None means no timestep!
     38#timestep = 0
    3939
    4040######
    … …  
    5656
    5757# one or more key strings from var_equations above
    58 var = ['elevation','stage']
     58var = ['stage']
    5959
    6060######
    … …  
    8686
    8787            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'
     88            name2 = join(directory, time_dir, project.scenario_name)+'_time_28380_0'
     89            name3 = join(directory, time_dir, project.scenario_name)+'_time_56760_0'
     90##            name4 = join(directory, time_dir, project.scenario_name)+'_time_32220_0'
     91##            name5 = join(directory, time_dir, project.scenario_name)+'_time_42960_0'
     92##            name6 = join(directory, time_dir, project.scenario_name)+'_time_53700_0'
    9093
    91             names = [name1] #, name2]#, name3]
     94            names = [name1, name2, name3]#, name4, name5, name6]
    9295     
    93            # asc_name = []   
     96            asc_name = []   
    9497               
    9598            for name in names:
    … …  
    112115                            format = 'asc')
    113116
    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)
     117                asc_name.append(outname + '.asc')
     118   
     119            maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
     120
     121            print 'max asc outname ', maxasc_outname
     122            print 'asc_name ', str(asc_name)
     123           
     124            MaxAsc(maxasc_outname, asc_name)

  • anuga_work/production/australia_ph2/geraldton/export_results_max.py

    r6845 r6926  
    2222directory = project.output_folder
    2323
    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'
     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'
    2727
    28 time_dirs = [time_dir1, time_dir2, time_dir3]
     28time_dir1 = '20090421_172030_run_final_0_64475_576_Tb_internal_mhingee'
     29
     30time_dirs = [time_dir1]#, time_dir2, time_dir3]
    2931
    3032cellsize = 250

  • anuga_work/production/australia_ph2/gold_coast/export_results_max.py

    r6589 r6926  
    2525##time_dir2 = '20090311_152929_run_final_0_70385_1203_Tb_jgriffin'
    2626##time_dir3 = '20090311_153009_run_final_0_70454_1203_Bf_jgriffin'
    27 time_dir4 = '20090317_144100_run_final_0_51314_3727_Tb_jgriffin'
    28 time_dirs = [time_dir4]#,time_dir2,time_dir3,time_dir4] 
     27##time_dir4 = '20090317_144100_run_final_0_51314_3727_Tb_jgriffin'
     28time_dir1 = '20090421_150448_run_final_0_50883_3727_Tb_internal_mhingee'
     29time_dirs = [time_dir1]#,time_dir2,time_dir3,time_dir4] 
    2930
    3031cellsize = 250
    … …  
    8283
    8384            name1 = join(directory, time_dir, project.scenario_name)
    84             name2 = join(directory, time_dir, project.scenario_name)+'_time_21060_0'
    85             name3 = join(directory, time_dir, project.scenario_name)+'_time_42120_0'
    86             name4 = join(directory, time_dir, project.scenario_name)+'_time_63180_0'
     85            name2 = join(directory, time_dir, project.scenario_name)+'_time_36420_0'
     86##            name2 = join(directory, time_dir, project.scenario_name)+'_time_21060_0'
     87##            name3 = join(directory, time_dir, project.scenario_name)+'_time_42120_0'
     88##            name4 = join(directory, time_dir, project.scenario_name)+'_time_63180_0'
    8789##            name5 = join(directory, time_dir, project.scenario_name)+'_time_53280_0'
    8890##            name6 = join(directory, time_dir, project.scenario_name)+'_time_66600_0'
    8991##            name7 = join(directory, time_dir, project.scenario_name)+'_time_79920_0'
    9092           
    91             names = [name1, name2, name3,name4]#,name5,name6,name7]
     93            names = [name1, name2]#, name3,name4]#,name5,name6,name7]
    9294     
    9395            asc_name = []   

  • anuga_work/production/australia_ph2/melbourne_east/export_results_max.py

    r6692 r6926  
    2222directory = project.output_folder
    2323
    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]
     24time_dir1 = '20090416_102418_run_final_0_58348_1979_Tb_west__kvanputt'
     25time_dir2 = '20090421_151342_run_final_0_64288_1979_Tb_west__kvanputt'
     26time_dir3 = '20090416_103113_run_final_0_51449_2649_Tb_east__kvanputt'
     27time_dir4 = '20090417_151345_run_final_0_58115_2649_Tb_east__kvanputt'
     28time_dir5 = '20090421_152343_run_final_0_31853_2649_Tb_east__kvanputt'
     29time_dir6 = '20090421_152359_run_final_0_58115_2649_Tb_east_internal_mhingee'
     30
     31time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
    3232 
    3333
    … …  
    3535##cellsize = 5
    3636
    37 #timestep = None    # None means no timestep!
    38 timestep = 0
     37timestep = None    # None means no timestep!
     38#timestep = 0
    3939
    4040######
    … …  
    5656
    5757# one or more key strings from var_equations above
    58 var = ['elevation']
     58var = ['stage']
    5959
    6060######
    … …  
    8686
    8787            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'
     88            name2 = join(directory, time_dir, project.scenario_name)+'_time_10740_0'
     89            name3 = join(directory, time_dir, project.scenario_name)+'_time_21480_0'
     90            name4 = join(directory, time_dir, project.scenario_name)+'_time_32220_0'
     91            name5 = join(directory, time_dir, project.scenario_name)+'_time_42960_0'
     92            name6 = join(directory, time_dir, project.scenario_name)+'_time_53700_0'
    9093
    91             names = [name1] #, name2]#, name3]
     94            names = [name1, name2, name3, name4, name5, name6]
    9295     
    93            # asc_name = []   
     96            asc_name = []   
    9497               
    9598            for name in names:
    … …  
    112115                            format = 'asc')
    113116
    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)
     117                asc_name.append(outname + '.asc')
     118   
     119            maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
     120
     121            print 'max asc outname ', maxasc_outname
     122            print 'asc_name ', str(asc_name)
     123           
     124            MaxAsc(maxasc_outname, asc_name)

  • anuga_work/production/australia_ph2/onslow/export_results_max.py

    r6845 r6926  
    2222directory = project.output_folder
    2323
    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'
     24time_dir1 = '20090421_172702_run_final_0_27279_711_Tb_internal_mhingee'
    2725
    28 time_dirs = [time_dir1, time_dir2, time_dir3]
     26time_dirs = [time_dir1]#, time_dir2, time_dir3]
    2927
    3028cellsize = 250
    … …  
    8280
    8381            name1 = join(directory, time_dir, project.scenario_name)
    84             name2 = join(directory, time_dir, project.scenario_name)+'_time_44400_0'
     82            name2 = join(directory, time_dir, project.scenario_name)+'_time_44340_0'
    8583
    8684            names = [name1, name2]
    8785
    8886            asc_name = []
    89 
     87           
    9088            for name in names:
    9189               
    … …  
    115113           
    116114            MaxAsc(maxasc_outname, asc_name)
    117              
    118              

  • anuga_work/production/australia_ph2/perth/export_results_max.py

    r6845 r6926  
    2222directory = project.output_folder
    2323
    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'
     24time_dir1 = '20090421_170526_run_final_0_64479_803_Tb_internal_mhingee'
     25##time_dir2 = ''
     26##time_dir3 = ''
     27##time_dir4 = ''
     28##time_dir5 = ''
     29##time_dir6 = ''
    2730
    28 time_dirs = [time_dir1, time_dir2, time_dir3]
     31time_dirs = [time_dir1]#, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
     32 
    2933
    3034cellsize = 250
    … …  
    8286
    8387            name1 = join(directory, time_dir, project.scenario_name)
    84             name2 = join(directory, time_dir, project.scenario_name)+'_time_54660_0'
     88            name2 = join(directory, time_dir, project.scenario_name)+'_time_54720_0'
     89##            name3 = join(directory, time_dir, project.scenario_name)+'_time_56760_0'
     90##            name4 = join(directory, time_dir, project.scenario_name)+'_time_32220_0'
     91##            name5 = join(directory, time_dir, project.scenario_name)+'_time_42960_0'
     92##            name6 = join(directory, time_dir, project.scenario_name)+'_time_53700_0'
    8593
    86             names = [name1, name2]
    87 
    88             asc_name = []
    89            
     94            names = [name1, name2]#, name3]#, name4, name5, name6]
     95     
     96            asc_name = []   
     97               
    9098            for name in names:
    91                
     99                         
    92100                outname = name + '_' + which_area + '_' + which_var
    93101                quantityname = var_equations[which_var]
    … …  
    115123           
    116124            MaxAsc(maxasc_outname, asc_name)
    117              

  • anuga_work/production/australia_ph2/portland/export_results_max.py

    r6716 r6926  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090327_143655_run_final_0_64214_2044_Tb_None_mhingee'
    25 time_dir2 = '20090327_143836_run_final_0_58337_2044_Tb_None_mhingee'
    26 time_dir3 = '20090330_095007_run_final_0_68779_2044_Tb_None_mhingee'
    27 time_dir4 = '20090330_095719_run_final_0_64214_2044_Bf_None_mhingee'
    28 time_dir5 = '20090330_095845_run_final_0_58337_2044_Bf_None_mhingee'
    29 time_dir6 = '20090330_143221_run_final_0_68779_2044_Bf_None_mhingee'
     24time_dir1 = '20090421_153138_run_final_0_64322_1969_Tb_internal_mhingee'
     25##time_dir2 = ''
     26##time_dir3 = ''
     27##time_dir4 = ''
     28##time_dir5 = ''
     29##time_dir6 = ''
    3030
    31 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
     31time_dirs = [time_dir1]#, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
    3232 
    3333
    … …  
    3535##cellsize = 5
    3636
    37 #timestep = None    # None means no timestep!
    38 timestep = None
     37timestep = None    # None means no timestep!
     38#timestep = 0
    3939
    4040######
    … …  
    5656
    5757# one or more key strings from var_equations above
    58 var = ['stage','speed']
     58var = ['stage']
    5959
    6060######
    … …  
    8686
    8787            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'
     88##            name2 = join(directory, time_dir, project.scenario_name)+'_time_10740_0'
     89##            name3 = join(directory, time_dir, project.scenario_name)+'_time_21480_0'
     90##            name4 = join(directory, time_dir, project.scenario_name)+'_time_32220_0'
     91##            name5 = join(directory, time_dir, project.scenario_name)+'_time_42960_0'
     92##            name6 = join(directory, time_dir, project.scenario_name)+'_time_53700_0'
    9093
    91             names = [name1] #, name2]#, name3]
     94            names = [name1]#, name2, name3, name4, name5, name6]
    9295     
    93             asc_name = []   
     96##            asc_name = []   
    9497               
    9598            for name in names:
    … …  
    112115                            format = 'asc')
    113116
    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)
     117##                asc_name.append(outname + '.asc')
     118##   
     119##            maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
     120##
     121##            print 'max asc outname ', maxasc_outname
     122##            print 'asc_name ', str(asc_name)
     123##           
     124##            MaxAsc(maxasc_outname, asc_name)

  • anuga_work/production/australia_ph2/strahan/export_results_max.py

    r6672 r6926  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090327_143655_run_final_0_64214_2044_Tb_None_mhingee'
    25 time_dir2 = '20090327_143836_run_final_0_58337_2044_Tb_None_mhingee'
    26 time_dir3 = '20090330_095007_run_final_0_68779_2044_Tb_None_mhingee'
    27 time_dir4 = '20090330_095719_run_final_0_64214_2044_Bf_None_mhingee'
    28 time_dir5 = '20090330_095845_run_final_0_58337_2044_Bf_None_mhingee'
    29 time_dir6 = '20090330_143221_run_final_0_68779_2044_Bf_None_mhingee'
     24time_dir1 = '20090421_151652_run_final_0_58337_2044_Tb_internal_mhingee'
     25##time_dir2 = ''
     26##time_dir3 = ''
     27##time_dir4 = ''
     28##time_dir5 = ''
     29##time_dir6 = ''
    3030
    31 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
     31time_dirs = [time_dir1]#, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
    3232 
    3333
    … …  
    3535##cellsize = 5
    3636
    37 #timestep = None    # None means no timestep!
    38 timestep = None
     37timestep = None    # None means no timestep!
     38#timestep = 0
    3939
    4040######
    … …  
    5656
    5757# one or more key strings from var_equations above
    58 var = ['stage','speed']
     58var = ['stage']
    5959
    6060######
    … …  
    8686
    8787            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'
     88##            name2 = join(directory, time_dir, project.scenario_name)+'_time_10740_0'
     89##            name3 = join(directory, time_dir, project.scenario_name)+'_time_21480_0'
     90##            name4 = join(directory, time_dir, project.scenario_name)+'_time_32220_0'
     91##            name5 = join(directory, time_dir, project.scenario_name)+'_time_42960_0'
     92##            name6 = join(directory, time_dir, project.scenario_name)+'_time_53700_0'
    9093
    91             names = [name1] #, name2]#, name3]
     94            names = [name1]#, name2, name3, name4, name5, name6]
    9295     
    93             asc_name = []   
     96##            asc_name = []   
    9497               
    9598            for name in names:
    … …  
    112115                            format = 'asc')
    113116
    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)
     117##                asc_name.append(outname + '.asc')
     118##   
     119##            maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
     120##
     121##            print 'max asc outname ', maxasc_outname
     122##            print 'asc_name ', str(asc_name)
     123##           
     124##            MaxAsc(maxasc_outname, asc_name)

  • anuga_work/production/australia_ph2/sydney/export_results.py

    r6731 r6926  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090331_150428_run_final_0_58152_2956_Tb_lfountai'
    25 ##time_dir2 = '20090311_152929_run_final_0_70385_1203_Tb_jgriffin'
    26 ##time_dir3 = '20090311_153009_run_final_0_70454_1203_Bf_jgriffin'
    27 ##time_dir4 = '20090311_153054_run_final_0_70339_1203_Bf_jgriffin'
    28 time_dirs = [time_dir1]#,time_dir2,time_dir3,time_dir4] 
     24time_dir1 = '20090408_171539_run_final_0_58152_2938_Tb__kvanputt'
     25time_dir2 = '20090421_150708_run_final_0_58152_2938_Tb_internal_mhingee'
     26time_dirs = [time_dir1, time_dir2]#,time_dir3,time_dir4] 
    2927
    3028cellsize = 250
    … …  
    5250
    5351# one or more key strings from var_equations above
    54 var = ['elevation', 'stage']
     52var = ['elevation']
    5553
    5654######
    … …  
    6967            northing_min = None
    7068            northing_max = None
    71 ##            easting_min = 290000
    72 ##            easting_max = 316000
    73 ##            northing_min = 6117000
    74 ##            northing_max = 6170000
    7569        else:
    7670            try:
    … …  
    8680
    8781            name1 = join(directory, time_dir, project.scenario_name)
    88             name2 = join(directory, time_dir, project.scenario_name)+'_time_66360_0'
    89 ##            name3 = join(directory, time_dir, project.scenario_name)+'_time_26640_0'
    90 ##            name4 = join(directory, time_dir, project.scenario_name)+'_time_39960_0'
     82##            name2 = join(directory, time_dir, project.scenario_name)+'_time_21060_0'
     83##            name3 = join(directory, time_dir, project.scenario_name)+'_time_42120_0'
     84##            name4 = join(directory, time_dir, project.scenario_name)+'_time_63180_0'
    9185##            name5 = join(directory, time_dir, project.scenario_name)+'_time_53280_0'
    9286##            name6 = join(directory, time_dir, project.scenario_name)+'_time_66600_0'
    9387##            name7 = join(directory, time_dir, project.scenario_name)+'_time_79920_0'
    9488           
    95             names = [name1]#, name2]#, name3,name4,name5,name6,name7]
     89            names = [name1]#, name2, name3,name4]#,name5,name6,name7]
    9690     
    97 ##            asc_name = []   
     91            asc_name = []   
    9892               
    9993            for name in names:

  • anuga_work/production/australia_ph2/sydney/export_results_max.py

    r6596 r6926  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090311_145740_run_final_0_58152_2956_Tb_jgriffin'
     24time_dir1 = '20090421_150708_run_final_0_58152_2938_Tb_internal_mhingee'
    2525##time_dir2 = '20090311_152929_run_final_0_70385_1203_Tb_jgriffin'
    2626##time_dir3 = '20090311_153009_run_final_0_70454_1203_Bf_jgriffin'
    … …  
    2828time_dirs = [time_dir1]#,time_dir2,time_dir3,time_dir4] 
    2929
    30 cellsize = 20
     30cellsize = 250
    3131##cellsize = 5
    3232
    … …  
    8686
    8787            name1 = join(directory, time_dir, project.scenario_name)
    88             name2 = join(directory, time_dir, project.scenario_name)+'_time_66360_0'
     88##            name2 = join(directory, time_dir, project.scenario_name)+'_time_66360_0'
    8989##            name3 = join(directory, time_dir, project.scenario_name)+'_time_26640_0'
    9090##            name4 = join(directory, time_dir, project.scenario_name)+'_time_39960_0'
    … …  
    116116                            format = 'asc')
    117117
    118                 asc_name.append(outname + '.asc')
    119    
    120             maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
    121 
    122             print 'max asc outname ', maxasc_outname
    123             print 'asc_name ', str(asc_name)
    124            
    125             MaxAsc(maxasc_outname, asc_name)
     118##                asc_name.append(outname + '.asc')
     119##   
     120##            maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
     121##
     122##            print 'max asc outname ', maxasc_outname
     123##            print 'asc_name ', str(asc_name)
     124##           
     125##            MaxAsc(maxasc_outname, asc_name)

  • anuga_work/production/australia_ph2/wyndham/project.py

    r6807 r6926  
    5353
    5454
    55 internal_polygon = False
     55internal_polygon = True
    5656
    5757#-------------------------------------------------------------------------------
Note: See TracChangeset for help on using the changeset viewer.