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

Legend:

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