Changeset 6929


Ignore:
Timestamp:
Apr 28, 2009, 6:06:14 PM (15 years ago)
Author:
myall
Message:
 
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  
    2222directory = project.output_folder
    2323
    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'
    2727
    28 time_dirs = [time_dir1, time_dir2, time_dir3]
     28time_dir1 = '20090421_135741_run_final_0_31977_2065_Tb_internal_lfountai'
     29
     30time_dirs = [time_dir1]#, time_dir2, time_dir3]
    2931
    3032cellsize = 250
  • anuga_work/production/australia_ph2/darwin/project.py

    r6814 r6929  
    5050
    5151# index is only used when wave = Tb
    52 index = 1829            # index from the PTHA - Y2000 0.257m
     52index = 1829            # index from the PTHA - Y2000 0.592m; using events picked with PTHA point 1833 (0.642m)
    5353wave = 'Tb'             # Bf (sts wave) Tb (index wave)
    5454
    5555
    56 internal_polygon = False
     56internal_polygon = True
    5757
    5858#-------------------------------------------------------------------------------
  • anuga_work/production/australia_ph2/lockhart_river/Arc_asc2raster_GDA94z50.py

    r6926 r6929  
    2525output_dir="anuga\\outputs\\"
    2626
    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'
     27time_dir1 = '20090408_152731_run_final_0_51316_1994_Tb__mhingee'
     28##time_dir2 = '20090408_200614_run_final_0_62786_1994_Tb__mhingee'
    3029
    31 time_dirs = [time_dir1, time_dir2, time_dir3]
     30##time_dir1 = '20090421_135544_run_final_0_51316_1994_Tb_internal_lfountai'
     31
     32time_dirs = [time_dir1]#, time_dir2]
    3233
    3334for time_dir in time_dirs:
     
    3839##    ocean = scenario_dir + "map_work\\\port_hedland.gdb\\outlines\\initial_conditions_ocean"
    3940
    40     print 'Process: Create File GDB'
    41     gp.CreateFileGDB_management(folder, "raster")
     41##    print 'Process: Create File GDB'
     42##    gp.CreateFileGDB_management(folder, "raster")
    4243
    4344    gp.Workspace = raster_gbd
     
    5253
    5354    generate_filename = []
    54     input_ascii = glob.glob(folder + '*stage.asc')
     55    input_ascii = glob.glob(folder + '*elevation.asc')
    5556
    5657    for infile in input_ascii:
  • anuga_work/production/australia_ph2/lockhart_river/export_results_max.py

    r6835 r6929  
    2222directory = project.output_folder
    2323
    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'
    2727
    28 time_dirs = [time_dir1, time_dir2, time_dir3]
     28time_dir1 = '20090421_135544_run_final_0_51316_1994_Tb_internal_lfountai'
     29
     30time_dirs = [time_dir1]#, time_dir2, time_dir3]
    2931
    3032cellsize = 250
  • anuga_work/production/australia_ph2/wyndham/export_results_max.py

    r6539 r6929  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090316_120226_run_final_0_17898_1833_Bf_mhingee'
    25 ##time_dir2 = '20090316_120344_run_final_0_64977_1833_Bf_mhingee'
    26 ##time_dir3 = '20090316_120457_run_final_0_70906_1833_Bf_mhingee'
    27 ##time_dir4 = '20090316_120558_run_final_0_71261_1833_Bf_mhingee'
     24time_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 = ''
    2828##time_dir5 = ''
    2929##time_dir6 = ''
    3030##time_dir7 = ''
    3131##time_dir8 = ''
    32 time_dirs = [time_dir1] #5, time_dir4]# , time_dir4, time_dir5, time_dir6]
     32
     33time_dirs = [time_dir1]#, time_dir2, time_dir3]# , time_dir4, time_dir5, time_dir6]
    3334 
    3435
     
    3637##cellsize = 5
    3738
    38 #timestep = None    # None means no timestep!
    39 timestep = 0
     39timestep = None    # None means no timestep!
     40#timestep = 0
    4041
    4142######
     
    5758
    5859# one or more key strings from var_equations above
    59 var = ['elevation']
     60var = ['stage']
    6061
    6162######
     
    8788
    8889            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]
    9399     
    94            # asc_name = []   
     100            asc_name = []   
    95101               
    96102            for name in names:
     
    113119                            format = 'asc')
    114120
    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_outname
    120 ##            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.