Changeset 7057


Ignore:
Timestamp:
May 19, 2009, 11:02:41 AM (11 years ago)
Author:
myall
Message:
 
Location:
anuga_work/production/australia_ph2
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/australia_ph2/darwin/Arc_asc2raster_GDA94.py

    r7023 r7057  
    2525output_dir = "anuga\\outputs\\"
    2626
    27 time_dir1 = '20090415_123212_run_final_0_17898_1829_Tb__kvanputt'
    28 time_dir2 = '20090428_170820_run_final_0_64977_1829_Tb__mhingee'
     27##time_dir1 = '20090415_123212_run_final_0_17898_1829_Tb__kvanputt'
     28##time_dir2 = '20090428_170820_run_final_0_64977_1829_Tb__mhingee'
    2929
    30 ##time_dir3 = '20090428_170944_run_final_0_70906_1829_Tb__mhingee'
    31 ##time_dir4 = '20090428_171421_run_final_0_71261_1829_Tb__mhingee'
     30time_dir3 = '20090428_170944_run_final_0_70906_1829_Tb__mhingee'
     31time_dir4 = '20090428_171421_run_final_0_71261_1829_Tb__mhingee'
    3232
    33 time_dir3 = '20090428_172336_run_final_0_17898_1829_Tb_internal_mhingee'
     33##time_dir3 = '20090428_172336_run_final_0_17898_1829_Tb_internal_mhingee'
    3434
    35 time_dirs = [time_dir1, time_dir2, time_dir3]
     35##time_dirs = [time_dir1, time_dir2, time_dir3]
    3636
    37 ##time_dirs = [time_dir4, time_dir5]
     37time_dirs = [time_dir3, time_dir4]
    3838
    3939for time_dir in time_dirs:
     
    4646    #ocean = scenario_dir + "map_work\\Perth.gdb\\Outlines\\initial_conditions_ocean1"
    4747   
    48 ##    print 'Process: Create File GDB'
    49 ##    gp.CreateFileGDB_management(folder, "raster")
     48    print 'Process: Create File GDB'
     49    gp.CreateFileGDB_management(folder, "raster")
    5050
    5151    gp.Workspace = raster_gbd
     
    6060
    6161    generate_filename = []
    62     input_ascii = glob.glob(folder + '*elevation.asc')
     62    input_ascii = glob.glob(folder + '*max.asc')
    6363    print time_dir
    6464
  • anuga_work/production/australia_ph2/melbourne_east/single_urs/export_results_max.py

    r7023 r7057  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090511_145103_run_trial_0_31853_2649_Tb_both__mhingee'
    25 time_dir2 = ''
    26 time_dir3 = ''
    27 time_dir4 = ''
    28 time_dir5 = ''
    29 time_dir6 = ''
     24time_dir1 = '20090511_152312_run_final_0_31853_2649_Tb_both__mhingee'
     25time_dir2 = '20090511_165621_run_final_0_51449_2649_Tb_both__mhingee'
     26time_dir3 = '20090511_165711_run_final_0_58115_2649_Tb_both__mhingee'
     27time_dir4 = '20090511_165725_run_final_0_64288_2649_Tb_both__mhingee'
     28time_dir5 = '20090511_165801_run_final_0_68793_2649_Tb_both__mhingee'
     29time_dir6 = '20090511_165818_run_final_0_58348_2649_Tb_both__mhingee'
    3030
    3131##time_dir1 = '20090421_152359_run_final_0_58115_2649_Tb_east_internal_mhingee'
     
    8888
    8989            name1 = join(directory, time_dir, project.scenario_name)
    90 ##            name2 = join(directory, time_dir, project.scenario_name)+'_time_10740_0'
    91 ##            name3 = join(directory, time_dir, project.scenario_name)+'_time_21480_0'
    92 ##            name4 = join(directory, time_dir, project.scenario_name)+'_time_32220_0'
    93 ##            name5 = join(directory, time_dir, project.scenario_name)+'_time_42960_0'
    94 ##            name6 = join(directory, time_dir, project.scenario_name)+'_time_53700_0'
     90            name2 = join(directory, time_dir, project.scenario_name)+'_time_10740_0'
     91            name3 = join(directory, time_dir, project.scenario_name)+'_time_21480_0'
     92            name4 = join(directory, time_dir, project.scenario_name)+'_time_32220_0'
     93            name5 = join(directory, time_dir, project.scenario_name)+'_time_42960_0'
     94            name6 = join(directory, time_dir, project.scenario_name)+'_time_53700_0'
    9595
    96             names = [name1]#, name2, name3, name4, name5, name6]
     96            names = [name1, name2, name3, name4, name5, name6]
    9797     
    98 ##            asc_name = []   
     98            asc_name = []   
    9999               
    100100            for name in names:
     
    117117                            format = 'asc')
    118118
    119 ##                asc_name.append(outname + '.asc')
    120 ##   
    121 ##            maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
    122 ##
    123 ##            print 'max asc outname ', maxasc_outname
    124 ##            print 'asc_name ', str(asc_name)
    125 ##           
    126 ##            MaxAsc(maxasc_outname, asc_name)
     119                asc_name.append(outname + '.asc')
     120   
     121            maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
     122
     123            print 'max asc outname ', maxasc_outname
     124            print 'asc_name ', str(asc_name)
     125           
     126            MaxAsc(maxasc_outname, asc_name)
Note: See TracChangeset for help on using the changeset viewer.