Changeset 6959


Ignore:
Timestamp:
May 5, 2009, 11:50:53 AM (15 years ago)
Author:
myall
Message:
 
Location:
anuga_work/production/australia_ph2
Files:
8 edited

Legend:

Unmodified
Added
Removed

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

    r6730 r6959  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090327_103428_run_final_0_70906_1833_Bf_None_mhingee'
    25 time_dir2 = '20090327_121728_run_final_0_70906_1833_Tb_None_mhingee'
    26 ##time_dir3 = ''
    27 ##time_dir4 = ''
     24time_dir1 = '20090415_123212_run_final_0_17898_1829_Tb__kvanputt'
     25##time_dir2 = '20090428_170820_run_final_0_64977_1829_Tb__mhingee'
     26##time_dir3 = '20090428_170944_run_final_0_70906_1829_Tb__mhingee'
     27##time_dir4 = '20090428_171421_run_final_0_71261_1829_Tb__mhingee'
    2828##time_dir5 = ''
    2929##time_dir6 = ''
    3030##time_dir7 = ''
    3131##time_dir8 = ''
    32 time_dirs = [time_dir1, time_dir2] # time_dir4]# , time_dir4, time_dir5, time_dir6]
     32
     33time_dir2 = '20090428_172336_run_final_0_17898_1829_Tb_internal_mhingee'
     34
     35time_dirs = [time_dir1, time_dir2]#, time_dir3]# , time_dir4, time_dir5, time_dir6]
    3336 
    3437
    … …  
    3639##cellsize = 5
    3740
    38 #timestep = None    # None means no timestep!
    39 timestep = 0
     41timestep = None    # None means no timestep!
     42#timestep = 0
    4043
    4144######
    … …  
    5760
    5861# one or more key strings from var_equations above
    59 var = ['stage','speed','depth']
     62var = ['stage']
    6063
    6164######
    … …  
    8790
    8891            name1 = join(directory, time_dir, project.scenario_name)
    89             name2 = join(directory, time_dir, project.scenario_name)+'_time_7380_0'
    90             name3 = join(directory, time_dir, project.scenario_name)+'_time_14760_0'
    91             name4 = join(directory, time_dir, project.scenario_name)+'_time_22140_0'
    92             name5 = join(directory, time_dir, project.scenario_name)+'_time_29520_0'
    93             name6 = join(directory, time_dir, project.scenario_name)+'_time_36900_0'
    94             name7 = join(directory, time_dir, project.scenario_name)+'_time_44280_0'
    95             name8 = join(directory, time_dir, project.scenario_name)+'_time_51660_0'
    96             name9 = join(directory, time_dir, project.scenario_name)+'_time_59040_0'
    97             name10 = join(directory, time_dir, project.scenario_name)+'_time_66420_0'
    98             name11 = join(directory, time_dir, project.scenario_name)+'_time_73800_0'
    99 
    100             names = [name1, name2, name3, name4, name5, name6, name7, name8, name9, name10, name11]
     92            name2 = join(directory, time_dir, project.scenario_name)+'_time_8520_0'
     93            name3 = join(directory, time_dir, project.scenario_name)+'_time_17040_0'
     94            name4 = join(directory, time_dir, project.scenario_name)+'_time_25560_0'
     95            name5 = join(directory, time_dir, project.scenario_name)+'_time_34080_0'
     96            name6 = join(directory, time_dir, project.scenario_name)+'_time_42600_0'
     97            name7 = join(directory, time_dir, project.scenario_name)+'_time_51120_0'
     98            name8 = join(directory, time_dir, project.scenario_name)+'_time_59640_0'
    10199           
     100            names = [name1, name2, name3, name4, name5, name6, name7, name8]
     101     
    102102            asc_name = []   
    103103               

  • anuga_work/production/australia_ph2/hobart/Arc_asc2raster_GDA94z50.py

    r6842 r6959  
    2525output_dir="anuga\\outputs\\"
    2626
    27 time_dir1 = '20090408_152504_run_final_0_51268_2528_Tb__mhingee'
    28 time_dir2 = '20090408_191713_run_final_0_58082_2528_Tb__mhingee'
    29 time_dir3 = '20090408_225208_run_final_0_64309_2528_Tb__mhingee'
     27time_dir1 = '20090424_121952_run_final_0_51204_2462_Tb__kvanputt'
     28time_dir2 = '20090424_122951_run_final_0_58212_2462_Tb__kvanputt'
     29time_dir3 = '20090424_123104_run_final_0_64425_2462_Tb__kvanputt'
    3030
    31 time_dirs = [time_dir1, time_dir2, time_dir3]
    32 events = [[time_dir1,],[time_dir2,],[time_dir3,]]
     31time_dir4 = '20090424_123250_run_final_0_51204_2462_Tb_internal_kvanputt'
    3332
    34 for event in events:
    35 ##for time_dir in time_dirs:
    36     time_dir = event[0]
    37     max_wave = event[1]
    38     print time_dir
    39     print max_wave
     33time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4]
     34
     35for time_dir in time_dirs:
     36
    4037    # Local variables...
    4138    folder = scenario_dir + output_dir +  time_dir +'\\'
    … …  
    4441##    ocean = scenario_dir + "map_work\\\port_hedland.gdb\\outlines\\initial_conditions_ocean"
    4542
    46 ##    print 'Process: Create File GDB'
    47 ##    gp.CreateFileGDB_management(folder, "raster")
     43    print 'Process: Create File GDB'
     44    gp.CreateFileGDB_management(folder, "raster")
    4845
    4946    gp.Workspace = raster_gbd
    … …  
    8380##        print 'Process: Extract by Mask'
    8481##        gp.ExtractByMask_sa(output_DEM, land, output_extract)
    85         # do this bit only if there are only stage asc files
    86         div_file = output_DEM.replace('stage','div_stage')
    87         print 'divide', output_DEM,' by', max_wave,' and call',div_file
    88         gp.Divide_sa(output_DEM,max_wave,div_file)
    8982
    9083

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

    r6842 r6959  
    2525output_dir="anuga\\outputs\\"
    2626
    27 time_dir1 = ''
    28 time_dir2 = ''
    29 time_dir3 = ''
     27time_dir1 = '20090417_151345_run_final_0_58115_2649_Tb_east__kvanputt'
     28time_dir2 = '20090416_103113_run_final_0_51449_2649_Tb_east__kvanputt'
     29time_dir3 = '20090421_152343_run_final_0_31853_2649_Tb_east__kvanputt'
     30time_dir4 = '20090416_102418_run_final_0_58348_1979_Tb_west__kvanputt'
     31time_dir5 = '20090421_151342_run_final_0_64288_1979_Tb_west__kvanputt'
     32##time_dir6 = '20090426_042056_run_final_0_68793_1979_Tb_west__kvanputt'
    3033
    31 events = [[time_dir1,],[time_dir2,],[time_dir3,]]
     34time_dir7 = '20090421_152359_run_final_0_58115_2649_Tb_east_internal_mhingee'
    3235
    33 for event in events:
    34 ##for time_dir in time_dirs:
    35     time_dir = event[0]
    36     max_wave = event[1]
    37     print time_dir
    38     print max_wave
     36time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir7]
     37##time_dirs = [time_dir4]
     38   
     39for time_dir in time_dirs:
    3940    # Local variables...
    4041    folder = scenario_dir + output_dir +  time_dir +'\\'
    … …  
    5152   
    5253    #replication dictionary
    53     replicate = (('strahan', ''),('_', ''),('max','_M'),
     54    replicate = (('melbourne', ''),('_', ''),('max','_M'),
    5455                 ('CBD', 'CDB'),('All',''),
    5556                 ('depth','_depth'),('speed', '_speed'),
    … …  
    5758
    5859    generate_filename = []
    59     input_ascii = glob.glob(folder + '*.asc')
     60    input_ascii = glob.glob(folder + '*max.asc')
    6061
    6162    for infile in input_ascii:
    … …  
    7475        gp.ASCIIToRaster_conversion(infile, output_DEM, "FLOAT")
    7576
    76         print 'Process: Define Projection'
     77        print 'Process: Define Projection'
     78## GDA_1994_MGA_Zone_55       
    7779        gp.DefineProjection_management(output_DEM, "PROJCS['CM_146',GEOGCS['GCS_GDA_1994',DATUM['D_GDA_1994',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]"
    7880                                                   ",PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',146.0],PARAMETER['Scale_Factor',0.9996]"
    … …  
    8284##        print 'Process: Extract by Mask'
    8385##        gp.ExtractByMask_sa(output_DEM, land, output_extract)
    84         # do this bit only if there are only stage asc files
    85         div_file = output_DEM.replace('stage','div_stage')
    86         print 'divide', output_DEM,' by', max_wave,' and call',div_file
    87         gp.Divide_sa(output_DEM,max_wave,div_file)
    8886
     87

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

    r6926 r6959  
    2727time_dir4 = '20090417_151345_run_final_0_58115_2649_Tb_east__kvanputt'
    2828time_dir5 = '20090421_152343_run_final_0_31853_2649_Tb_east__kvanputt'
    29 time_dir6 = '20090421_152359_run_final_0_58115_2649_Tb_east_internal_mhingee'
     29time_dir6 = '20090426_042056_run_final_0_68793_1979_Tb_west__kvanputt'
    3030
    31 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
    32  
     31##time_dir1 = '20090421_152359_run_final_0_58115_2649_Tb_east_internal_mhingee'
     32
     33##time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
     34time_dirs = [time_dir6]
    3335
    3436cellsize = 250

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

    r6842 r6959  
    2929time_dir3 = '20090408_175114_run_final_0_68779_2044_Tb__mhingee'
    3030
    31 time_dirs = [time_dir1, time_dir2, time_dir3]
    32 events = [[time_dir1,0.593696],[time_dir2,0.566183],[time_dir3,0.561357]]
     31time_dir4 = '20090421_151652_run_final_0_58337_2044_Tb_internal_mhingee'
    3332
    34 for event in events:
    35 ##for time_dir in time_dirs:
    36     time_dir = event[0]
    37     max_wave = event[1]
    38     print time_dir
    39     print max_wave
     33##time_dirs = [time_dir1, time_dir2, time_dir3]
     34time_dirs = [time_dir4]
     35   
     36for time_dir in time_dirs:
    4037    # Local variables...
    4138    folder = scenario_dir + output_dir +  time_dir +'\\'
    … …  
    8380##        print 'Process: Extract by Mask'
    8481##        gp.ExtractByMask_sa(output_DEM, land, output_extract)
    85         # do this bit only if there are only stage asc files
    86         div_file = output_DEM.replace('stage','div_stage')
    87         print 'divide', output_DEM,' by', max_wave,' and call',div_file
    88         gp.Divide_sa(output_DEM,max_wave,div_file)
    8982
     83

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

    r6842 r6959  
    2828time_dir2 = '20090409_011331_run_final_0_58222_2938_Tb__kvanputt'
    2929time_dir3 = '20090409_085101_run_final_0_58326_2938_Tb__kvanputt'
    30 time_dirs = [time_dir1, time_dir2, time_dir3]
    3130
    32 ##time_dir1 = '20090311_145740_run_final_0_58152_2956_Tb_jgriffin'
    33 ##time_dir2 = '20090311_145828_run_final_0_58152_2956_Bf_jgriffin'
    34 ##time_dir3 = '20090313_154712_run_final_0_58222_2956_Bf_jgriffin'
    35 ##time_dir4 = '20090313_154810_run_final_0_58326_2956_Bf_jgriffin'
    36 ##time_dir5 = '20090316_142137_run_final_0_58326_2956_Bf_jgriffin'
    37 ##time_dir6 = '20081031_133925_run_final_0.6_27283_alpha0.1_kvanputt'
     31##time_dir1 = '20090421_150708_run_final_0_58152_2938_Tb_internal_mhingee'
    3832
    39 
    40 ##time_dirs = [time_dir1]#, time_dir4 , time_dir3, time_dir4, time_dir5]#, time_dir6]
     33time_dirs = [time_dir1]#, time_dir2, time_dir3]
    4134
    4235for time_dir in time_dirs:
    … …  
    5750   
    5851    #replication dictionary
    59     replicate = (('sydney', ''),
     52    replicate = (('sydney', ''),('All',''),
    6053                 ('_', ''),('Geordie', 'Geo'),('Sorrento', 'Sor'), ('max','M_'),
    6154                 ('Fremantle', 'Fre'),('Rockingham', 'Roc'),('depth','_dep_'),
    … …  
    6356
    6457    generate_filename = []
    65     input_ascii = glob.glob(folder + '*.asc')
     58    input_ascii = glob.glob(folder + '*elevation.asc')
    6659    print time_dir
    6760

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

    r6864 r6959  
    3232    event_number = int(sys.argv[1])
    3333else:   
    34     event_number = 58152    # the event number or the mux file name
    35 
    36 event_number_list = [58152, 58222, 58326] # To piggy back multiple events
     34    event_number = 51418    # the event number or the mux file name
     35
     36##event_number_list = [58152, 58222, 58326] # To piggy back multiple events
     37event_number_list = [31865, 51418]# 58152
    3738
    3839# events: 27346 (java 2.1e-6), 58334 (puysegur 8.8e-7), 64322 (sandwich 8.2e-6),
    … …  
    5455
    5556
    56 internal_polygon = True
     57internal_polygon = False
    5758
    5859#-------------------------------------------------------------------------------

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

    r6929 r6959  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090415_094009_run_final_0_17918_1823_Tb__kvanputt'
     24##time_dir1 = '20090415_094009_run_final_0_17918_1823_Tb__kvanputt'
    2525##time_dir2 = '20090418_182759_run_final_0_64962_1823_Tb__kvanputt'
    2626##time_dir3 = '20090421_145541_run_final_0_70920_1823_Tb__kvanputt'
    … …  
    3030##time_dir7 = ''
    3131##time_dir8 = ''
     32
     33time_dir1 = '20090428_145948_run_final_0_64962_1823_Tb_internal_mhingee'
    3234
    3335time_dirs = [time_dir1]#, time_dir2, time_dir3]# , time_dir4, time_dir5, time_dir6]
Note: See TracChangeset for help on using the changeset viewer.