Changeset 7014


Ignore:
Timestamp:
May 11, 2009, 4:47:04 PM (11 years ago)
Author:
myall
Message:

exporting results and making rasters

Location:
anuga_work/production/australia_ph2
Files:
12 edited

Legend:

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

    r6942 r7014  
    2525output_dir="anuga\\outputs\\"
    2626
    27 time_dir1 = '20090408_152157_run_final_0_27319_1330_Tb__mhingee'
    28 time_dir2 = '20090408_201701_run_final_0_64344_1330_Tb__mhingee'
    29 time_dir3 = '20090409_010309_run_final_0_46697_1330_Tb__mhingee'
     27time_dir1 = '20090507_115824_run_final_30_0_27319_1330_Tb__kvanputt'
     28time_dir2 = '20090507_164822_run_final_30_0_64344_1330_Tb__kvanputt'
     29time_dir3 = '20090507_214813_run_final_30_0_46697_1330_Tb__kvanputt'
    3030
    31 time_dir4 = '20090421_165204_run_final_0_64344_1330_Tb_internal_mhingee'
     31time_dir4 = '20090507_115302_run_final_30_0_64344_1330_Tb_internal_kvanputt'
    3232
    33 ##time_dirs = [time_dir1, time_dir2]#, time_dir3]
    34 time_dirs = [time_dir2]
     33time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4]
     34##time_dirs = [time_dir2]
    3535
    3636for time_dir in time_dirs:
     
    4141##    ocean = scenario_dir + "map_work\\\port_hedland.gdb\\outlines\\initial_conditions_ocean"
    4242
    43 ##    print 'Process: Create File GDB'
    44 ##    gp.CreateFileGDB_management(folder, "raster")
     43    print 'Process: Create File GDB'
     44    gp.CreateFileGDB_management(folder, "raster")
    4545
    4646    gp.Workspace = raster_gbd
     
    5555
    5656    generate_filename = []
    57     input_ascii = glob.glob(folder + '*elevation.asc')
     57    input_ascii = glob.glob(folder + '*.asc')
    5858
    5959    for infile in input_ascii:
  • anuga_work/production/australia_ph2/bamaga/Arc_asc2raster_GDA94cm143_5.py

    r6985 r7014  
    2525output_dir="anuga\\outputs\\"
    2626
    27 time_dir1 = '20090408_152655_run_final_0_31975_2000_Tb__mhingee'
    28 time_dir2 = '20090409_123239_run_final_0_51390_2000_Tb__mhingee'
    29 time_dir3 = '20090410_110658_run_final_0_63735_2000_Tb__mhingee'
    30 time_dir4 = '20090421_135412_run_final_0_31975_2000_Tb_internal_lfountai'
     27time_dir1 = '20090508_132715_run_final_0_31975_2000_Tb__mhingee'
     28time_dir2 = '20090508_140920_run_final_0_51390_2000_Tb__mhingee'
     29time_dir3 = '20090508_140936_run_final_0_63735_2000_Tb__mhingee'
     30time_dir4 = '20090508_133139_run_final_0_31975_2000_Tb_internal_mhingee'
    3131
    32 ##time_dirs = [time_dir2, time_dir3, time_dir4]
    33 time_dirs = [time_dir1, time_dir4]
     32time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4]
     33##time_dirs = [time_dir1, time_dir4]
    3434
    3535for time_dir in time_dirs:
  • anuga_work/production/australia_ph2/bamaga/export_results.py

    r6977 r7014  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090408_152655_run_final_0_31975_2000_Tb__mhingee'
    25 time_dir2 = '20090421_135412_run_final_0_31975_2000_Tb_internal_lfountai'
     24time_dir1 = '20090508_132715_run_final_0_31975_2000_Tb__mhingee'
     25time_dir2 = '20090508_133139_run_final_0_31975_2000_Tb_internal_mhingee'
    2626##time_dir3 = ''
    2727##time_dir4 = ''
  • anuga_work/production/australia_ph2/darwin/export_results_max.py

    r6959 r7014  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090415_123212_run_final_0_17898_1829_Tb__kvanputt'
    25 ##time_dir2 = '20090428_170820_run_final_0_64977_1829_Tb__mhingee'
     24##time_dir1 = '20090415_123212_run_final_0_17898_1829_Tb__kvanputt'
     25time_dir2 = '20090428_170820_run_final_0_64977_1829_Tb__mhingee'
    2626##time_dir3 = '20090428_170944_run_final_0_70906_1829_Tb__mhingee'
    2727##time_dir4 = '20090428_171421_run_final_0_71261_1829_Tb__mhingee'
     
    3131##time_dir8 = ''
    3232
    33 time_dir2 = '20090428_172336_run_final_0_17898_1829_Tb_internal_mhingee'
     33##time_dir2 = '20090428_172336_run_final_0_17898_1829_Tb_internal_mhingee'
    3434
    35 time_dirs = [time_dir1, time_dir2]#, time_dir3]# , time_dir4, time_dir5, time_dir6]
     35time_dirs = [time_dir2]#, time_dir2]#, time_dir3]# , time_dir4, time_dir5, time_dir6]
    3636 
    3737
  • anuga_work/production/australia_ph2/hobart/Arc_asc2raster_GDA94z50.py

    r6959 r7014  
    2525output_dir="anuga\\outputs\\"
    2626
    27 time_dir1 = '20090424_121952_run_final_0_51204_2462_Tb__kvanputt'
    28 time_dir2 = '20090424_122951_run_final_0_58212_2462_Tb__kvanputt'
    29 time_dir3 = '20090424_123104_run_final_0_64425_2462_Tb__kvanputt'
     27time_dir1 = '20090507_111303_run_final_30_0_51204_2462_Tb__kvanputt'
     28time_dir2 = '20090507_150745_run_final_30_0_58212_2462_Tb__kvanputt'
     29time_dir3 = '20090507_185058_run_final_30_0_64425_2462_Tb__kvanputt'
    3030
    31 time_dir4 = '20090424_123250_run_final_0_51204_2462_Tb_internal_kvanputt'
     31time_dir4 = '20090507_111251_run_final_30_0_51204_2462_Tb_internal_kvanputt'
    3232
    3333time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4]
     
    4141##    ocean = scenario_dir + "map_work\\\port_hedland.gdb\\outlines\\initial_conditions_ocean"
    4242
    43     print 'Process: Create File GDB'
    44     gp.CreateFileGDB_management(folder, "raster")
     43##    print 'Process: Create File GDB'
     44##    gp.CreateFileGDB_management(folder, "raster")
    4545
    4646    gp.Workspace = raster_gbd
     
    5555
    5656    generate_filename = []
    57     input_ascii = glob.glob(folder + '*.asc')
     57    input_ascii = glob.glob(folder + '*elevation.asc')
    5858
    5959    for infile in input_ascii:
  • anuga_work/production/australia_ph2/lockhart_river/Arc_asc2raster_GDA94z50.py

    r6929 r7014  
    2525output_dir="anuga\\outputs\\"
    2626
    27 time_dir1 = '20090408_152731_run_final_0_51316_1994_Tb__mhingee'
     27##time_dir1 = '20090408_152731_run_final_0_51316_1994_Tb__mhingee'
    2828##time_dir2 = '20090408_200614_run_final_0_62786_1994_Tb__mhingee'
     29time_dir3 = '20090508_141119_run_final_0_72030_1994_Tb__kvanputt'
    2930
    3031##time_dir1 = '20090421_135544_run_final_0_51316_1994_Tb_internal_lfountai'
    3132
    32 time_dirs = [time_dir1]#, time_dir2]
     33time_dirs = [time_dir3]#, time_dir2]
    3334
    3435for time_dir in time_dirs:
     
    3940##    ocean = scenario_dir + "map_work\\\port_hedland.gdb\\outlines\\initial_conditions_ocean"
    4041
    41 ##    print 'Process: Create File GDB'
    42 ##    gp.CreateFileGDB_management(folder, "raster")
     42    print 'Process: Create File GDB'
     43    gp.CreateFileGDB_management(folder, "raster")
    4344
    4445    gp.Workspace = raster_gbd
     
    5354
    5455    generate_filename = []
    55     input_ascii = glob.glob(folder + '*elevation.asc')
     56    input_ascii = glob.glob(folder + '*.asc')
    5657
    5758    for infile in input_ascii:
  • anuga_work/production/australia_ph2/lockhart_river/export_results_max.py

    r6929 r7014  
    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'
     24time_dir1 = '20090408_152731_run_final_0_51316_1994_Tb__mhingee'
     25time_dir2 = '20090408_200614_run_final_0_62786_1994_Tb__mhingee'
     26time_dir3 = '20090508_141119_run_final_0_72030_1994_Tb__kvanputt'
    2727
    28 time_dir1 = '20090421_135544_run_final_0_51316_1994_Tb_internal_lfountai'
     28##time_dir1 = '20090421_135544_run_final_0_51316_1994_Tb_internal_lfountai'
    2929
    30 time_dirs = [time_dir1]#, time_dir2, time_dir3]
     30time_dirs = [time_dir3]#, time_dir2, time_dir3]
    3131
    3232cellsize = 250
  • anuga_work/production/australia_ph2/perth/export_results.py

    r6926 r7014  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090415_090717_run_final_0_64479_803_Tb__kvanputt'
    25 time_dir2 = '20090421_170526_run_final_0_64479_803_Tb_internal_mhingee'
     24time_dir1 = '20090508_134620_run_final_30_0_64479_803_Tb__kvanputt'
     25time_dir2 = '20090508_134358_run_final_30_0_64479_803_Tb_internal_kvanputt'
    2626##time_dir3 = ''
    2727##time_dir4 = ''
  • anuga_work/production/australia_ph2/perth/export_results_max.py

    r6926 r7014  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090421_170526_run_final_0_64479_803_Tb_internal_mhingee'
    25 ##time_dir2 = ''
    26 ##time_dir3 = ''
    27 ##time_dir4 = ''
     24time_dir1 = '20090508_134620_run_final_30_0_64479_803_Tb__kvanputt'
     25time_dir2 = '20090508_210649_run_final_30_0_70921_803_Tb__kvanputt'
     26time_dir3 = '20090509_052217_run_final_30_0_27181_803_Tb__kvanputt'
     27time_dir4 = '20090508_134358_run_final_30_0_64479_803_Tb_internal_kvanputt'
    2828##time_dir5 = ''
    2929##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
     
    8686
    8787            name1 = join(directory, time_dir, project.scenario_name)
    88             name2 = join(directory, time_dir, project.scenario_name)+'_time_54720_0'
    89 ##            name3 = join(directory, time_dir, project.scenario_name)+'_time_56760_0'
     88            name2 = join(directory, time_dir, project.scenario_name)+'_time_27330_0'
     89            name3 = join(directory, time_dir, project.scenario_name)+'_time_54660_0'
    9090##            name4 = join(directory, time_dir, project.scenario_name)+'_time_32220_0'
    9191##            name5 = join(directory, time_dir, project.scenario_name)+'_time_42960_0'
    9292##            name6 = join(directory, time_dir, project.scenario_name)+'_time_53700_0'
    9393
    94             names = [name1, name2]#, name3]#, name4, name5, name6]
     94            names = [name1, name2, name3]#, name4, name5, name6]
    9595     
    9696            asc_name = []   
  • anuga_work/production/australia_ph2/strahan/Arc_asc2raster_GDA94z50.py

    r6971 r7014  
    2525output_dir="anuga\\outputs\\"
    2626
    27 time_dir1 = '20090408_152513_run_final_0_58337_2044_Tb__mhingee'
    28 time_dir2 = '20090408_163848_run_final_0_64214_2044_Tb__mhingee'
    29 time_dir3 = '20090408_175114_run_final_0_68779_2044_Tb__mhingee'
     27time_dir1 = '20090507_111501_run_final_30_0_58337_2044_Tb_internal_kvanputt'
     28time_dir2 = '20090507_122805_run_final_30_0_64214_2044_Tb_internal_kvanputt'
     29time_dir3 = '20090507_134426_run_final_30_0_68779_2044_Tb_internal_kvanputt'
    3030
    31 time_dir4 = '20090421_151652_run_final_0_58337_2044_Tb_internal_mhingee'
     31time_dir4 = ''
    3232
    33 ##time_dirs = [time_dir1, time_dir2, time_dir3]
    34 time_dirs = [time_dir1]
     33time_dirs = [time_dir1, time_dir2, time_dir3]
     34##time_dirs = [time_dir1]
    3535   
    3636for time_dir in time_dirs:
     
    4141##    ocean = scenario_dir + "map_work\\\port_hedland.gdb\\outlines\\initial_conditions_ocean"
    4242
    43 ##    print 'Process: Create File GDB'
    44 ##    gp.CreateFileGDB_management(folder, "raster")
     43    print 'Process: Create File GDB'
     44    gp.CreateFileGDB_management(folder, "raster")
    4545
    4646    gp.Workspace = raster_gbd
     
    5555
    5656    generate_filename = []
    57     input_ascii = glob.glob(folder + '*elevation.asc')
     57    input_ascii = glob.glob(folder + '*.asc')
    5858
    5959    for infile in input_ascii:
  • anuga_work/production/australia_ph2/sydney/export_results.py

    r6982 r7014  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090408_171539_run_final_0_58152_2938_Tb__kvanputt'
    25 time_dir2 = '20090421_150708_run_final_0_58152_2938_Tb_internal_mhingee'
     24time_dir1 = '20090508_015217_run_final_30_0_58152_2938_Tb__kvanputt'
     25time_dir2 = '20090428_131539_run_final_30_0_58152_2938_Tb_internal_kvanputt'
    2626time_dirs = [time_dir1, time_dir2]#,time_dir3,time_dir4] 
    2727
     
    2929##cellsize = 5
    3030
    31 timestep = None    # None means no timestep!
    32 #timestep = 0
     31##timestep = None    # None means no timestep!
     32timestep = 0
    3333
    3434######
     
    5050
    5151# one or more key strings from var_equations above
    52 var = ['stage']
     52var = ['elevation']
    5353
    5454######
  • anuga_work/production/australia_ph2/townsville/export_results_max.py

    r6985 r7014  
    2929##time_dir6 = ''
    3030##
    31 time_dirs = [time_dir1]#, time_dir2, time_dir3] #time_dir2]# , time_dir4, time_dir5, time_dir6]
     31time_dirs = [time_dir1, time_dir2, time_dir3] #time_dir2]# , time_dir4, time_dir5, time_dir6]
    3232 
    3333
     
    5656
    5757# one or more key strings from var_equations above
    58 var = ['stage']
     58var = ['speed']
    5959
    6060######
Note: See TracChangeset for help on using the changeset viewer.