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/dampier
Files:
1 added
1 edited

Legend:

Unmodified
Added
Removed
  • 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]
Note: See TracChangeset for help on using the changeset viewer.