Changeset 7011


Ignore:
Timestamp:
May 11, 2009, 3:43:32 PM (11 years ago)
Author:
kristy
Message:
 
Location:
anuga_work/production/australia_ph2
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/australia_ph2/carnarvon/Arc_asc2raster_GDA94z50.py

    r6949 r7011  
    2525output_dir="anuga\\outputs\\"
    2626
    27 time_dir1 = '20090429_170742_run_trial_0_27285_84_Tb__kvanputt'
    28 time_dir2 = ''
     27time_dir1 = '20090504_114635_run_final_0_68754_84_Tb__mhingee'
     28time_dir2 = '20090504_114715_run_final_0_70881_84_Tb__mhingee'
     29time_dir3 = '20090505_100250_run_final_0_27285_84_Tb__mhingee'
     30time_dir4 = '20090505_173023_run_final_0_27285_84_Tb_internal_mhingee'
    2931
    30 time_dirs = [time_dir1] #, time_dir2] #4, time_dir5] #, time_dir3, time_dir4, time_dir5, time_dir6]
     32time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4]
    3133 
    3234for time_dir in time_dirs:
     
    3537    folder = scenario_dir + output_dir + time_dir +'\\'
    3638    raster_gbd = folder + 'raster.gdb'
    37     land = scenario_dir + "map_work\\Busselton.gdb\\Internal_polygons\\initial_conditions_extend_Di"
    38     ocean = scenario_dir + "map_work\\Busselton.gdb\\input_boundaries\\Ocean"
    39 
    40    
    41     print 'Process: Create File GDB'
    42     gp.CreateFileGDB_management(folder, "raster")
     39       
     40##    print 'Process: Create File GDB'
     41##    gp.CreateFileGDB_management(folder, "raster")
    4342
    4443    gp.Workspace = raster_gbd
     
    7372                                                   ",PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',113.5],PARAMETER['Scale_Factor',0.9996]"
    7473                                                   ",PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]")
    75 ##        output_extract = output_DEM + '_E'
    76 ##        print 'Output Extract ',output_extract
    77 ##        print 'Process: Extract by Mask'
    78 ##        gp.ExtractByMask_sa(output_DEM, land, output_extract)
    7974
    80 
  • anuga_work/production/australia_ph2/gold_coast/export_results_max.py

    r6926 r7011  
    2222directory = project.output_folder
    2323
    24 ##time_dir1 = '20090311_152848_run_final_0_70385_1203_Bf_jgriffin'
    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 = '20090317_144100_run_final_0_51314_3727_Tb_jgriffin'
    28 time_dir1 = '20090421_150448_run_final_0_50883_3727_Tb_internal_mhingee'
    29 time_dirs = [time_dir1]#,time_dir2,time_dir3,time_dir4] 
     24time_dir1 = '20090505_145250_run_final_0_31854_3727_Tb__lfountai'
     25time_dir2 = '20090506_045043_run_final_0_64479_3727_Tb__lfountai'
     26time_dir3 = '20090506_190156_run_final_0_50883_3727_Tb__lfountai'
     27time_dirs = [time_dir1,time_dir2,time_dir3] #,time_dir4] 
    3028
    3129cellsize = 250
     
    8381
    8482            name1 = join(directory, time_dir, project.scenario_name)
    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'
    89 ##            name5 = join(directory, time_dir, project.scenario_name)+'_time_53280_0'
    90 ##            name6 = join(directory, time_dir, project.scenario_name)+'_time_66600_0'
    91 ##            name7 = join(directory, time_dir, project.scenario_name)+'_time_79920_0'
    92            
    93             names = [name1, name2]#, name3,name4]#,name5,name6,name7]
     83            name2 = join(directory, time_dir, project.scenario_name)+'_time_36480_0'
     84
     85            names = [name1, name2]
    9486     
    9587            asc_name = []   
  • anuga_work/production/australia_ph2/townsville/Arc_asc2raster_GDA94z50.py

    r6939 r7011  
    2727time_dir1 = '20090414_104809_run_final_0_31965_2740_Tb__kvanputt'
    2828time_dir2 = '20090420_064036_run_final_0_51469_2740_Tb__kvanputt'
    29 ##time_dir3 = '20090429_145415_run_final_0_62890_2740_Tb__mhingee'
     29time_dir3 = '20090421_141635_run_final_0_31965_2740_Tb_internal_lfountai'
     30time_dir4 = '20090429_145415_run_final_0_62890_2740_Tb__mhingee'
    3031
    31 time_dir3 = '20090421_141635_run_final_0_31965_2740_Tb_internal_lfountai'
    32 
    33 ##time_dirs = [time_dir1, time_dir2, time_dir3]
    34 time_dirs = [time_dir1, time_dir3]
     32time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4]
    3533
    3634for time_dir in time_dirs:
  • anuga_work/production/australia_ph2/wyndham/Arc_asc2raster_GDA94z50.py

    r6971 r7011  
    2525output_dir="anuga\\outputs\\"
    2626
    27 time_dir1 = '20090415_094009_run_final_0_17918_1823_Tb__kvanputt'
    28 time_dir2 = '20090418_182759_run_final_0_64962_1823_Tb__kvanputt'
    29 time_dir3 = '20090421_145541_run_final_0_70920_1823_Tb__kvanputt'
     27##time_dir1 = '20090415_094009_run_final_0_17918_1823_Tb__kvanputt'
     28##time_dir2 = '20090418_182759_run_final_0_64962_1823_Tb__kvanputt'
     29##time_dir3 = '20090421_145541_run_final_0_70920_1823_Tb__kvanputt'
    3030
    31 ##time_dir1 = '20090428_145948_run_final_0_64962_1823_Tb_internal_mhingee'
     31time_dir1 = '20090428_145948_run_final_0_64962_1823_Tb_internal_mhingee'
    3232
    33 time_dirs = [time_dir2]#, time_dir2, time_dir3]
     33time_dirs = [time_dir1]#, time_dir2, time_dir3]
    3434
    3535for time_dir in time_dirs:
     
    5454
    5555    generate_filename = []
    56     input_ascii = glob.glob(folder + '*elevation.asc')
     56    input_ascii = glob.glob(folder + '*max.asc')
    5757
    5858    for infile in input_ascii:
Note: See TracChangeset for help on using the changeset viewer.