Changeset 6939


Ignore:
Timestamp:
Apr 29, 2009, 5:12:28 PM (11 years ago)
Author:
myall
Message:

exporting stage and elevation and making rasters for GBR models

Location:
anuga_work/production/australia_ph2
Files:
12 added
9 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/australia_ph2/bamaga/export_results_max.py

    r6922 r6939  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090424_103709_run_final_0_64977_1914_Tb__jgriffin'
    25 #time_dir2 = '20090424_165617_run_final_0_71114_1914_Tb__jgriffin'
    26 #time_dir3 = '20090410_110658_run_final_0_63735_2000_Tb__mhingee'
     24##time_dir1 = '20090408_152655_run_final_0_31975_2000_Tb__mhingee'
     25##time_dir2 = '20090409_123239_run_final_0_51390_2000_Tb__mhingee'
     26##time_dir3 = '20090410_110658_run_final_0_63735_2000_Tb__mhingee'
    2727
    28 time_dirs = [time_dir1]#, time_dir2]#, time_dir3]
     28time_dir1 = '20090421_135412_run_final_0_31975_2000_Tb_internal_lfountai'
     29
     30time_dirs = [time_dir1]#, time_dir2, time_dir3]
    2931
    3032cellsize = 250
     
    5254
    5355# one or more key strings from var_equations above
    54 var = ['stage','elevation']
     56var = ['stage']
    5557
    5658######
     
    8284
    8385            name1 = join(directory, time_dir, project.scenario_name)
    84             name2 = join(directory, time_dir, project.scenario_name)+'_time_35040_0'
    85             name3 = join(directory, time_dir, project.scenario_name)+'_time_70080_0'
    86 
    87             names = [name1, name2, name3]
     86##            name2 = join(directory, time_dir, project.scenario_name)+'_time_37920_0'
     87            name2 = join(directory, time_dir, project.scenario_name)+'_time_37860_0'
     88           
     89            names = [name1, name2]
    8890
    8991            asc_name = []
  • anuga_work/production/australia_ph2/bundaberg/export_results_max.py

    r6835 r6939  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090408_153233_run_final_0_31756_3274_Tb__mhingee'
    25 time_dir2 = '20090409_015134_run_final_0_51469_3274_Tb__mhingee'
    26 time_dir3 = '20090409_112941_run_final_0_62840_3274_Tb__mhingee'
     24##time_dir1 = '20090408_153233_run_final_0_31756_3274_Tb__mhingee'
     25##time_dir2 = '20090409_015134_run_final_0_51469_3274_Tb__mhingee'
     26##time_dir3 = '20090409_112941_run_final_0_62840_3274_Tb__mhingee'
    2727
    28 time_dirs = [time_dir1, time_dir2, time_dir3]
     28time_dir1 = '20090421_143037_run_final_0_31756_3274_Tb_internal_lfountai'
     29
     30time_dirs = [time_dir1]#, time_dir2, time_dir3]
    2931
    3032cellsize = 250
  • anuga_work/production/australia_ph2/cairns/export_results_max.py

    r6835 r6939  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090408_153044_run_final_0_31975_2172_Tb__mhingee'
    25 time_dir2 = '20090408_223355_run_final_0_51469_2172_Tb__mhingee'
    26 time_dir3 = '20090409_055521_run_final_0_63737_2172_Tb__mhingee'
     24##time_dir1 = '20090408_153044_run_final_0_31975_2172_Tb__mhingee'
     25##time_dir2 = '20090408_223355_run_final_0_51469_2172_Tb__mhingee'
     26##time_dir3 = '20090409_055521_run_final_0_63737_2172_Tb__mhingee'
    2727
    28 time_dirs = [time_dir1, time_dir2, time_dir3]
     28time_dir1 = '20090421_140811_run_final_0_31975_2172_Tb_internal_lfountai'
     29
     30time_dirs = [time_dir1]#, time_dir2, time_dir3]
    2931
    3032cellsize = 250
  • anuga_work/production/australia_ph2/cooktown/Arc_asc2raster_GDA94z50.py

    r6929 r6939  
    2626
    2727time_dir1 = '20090408_153032_run_final_0_31977_2065_Tb__mhingee'
    28 time_dir2 = '20090408_214931_run_final_0_51469_2065_Tb__mhingee'
    29 time_dir3 = '20090409_050811_run_final_0_63735_2065_Tb__mhingee'
     28##time_dir2 = '20090408_214931_run_final_0_51469_2065_Tb__mhingee'
     29##time_dir3 = '20090409_050811_run_final_0_63735_2065_Tb__mhingee'
    3030
    31 ##time_dir1 = '20090421_135741_run_final_0_31977_2065_Tb_internal_lfountai'
     31##time_dir2 = '20090421_135741_run_final_0_31977_2065_Tb_internal_lfountai'
    3232
    33 time_dirs = [time_dir1, time_dir2, time_dir3]
     33time_dirs = [time_dir1]#, time_dir2]#, time_dir3]
    3434
    3535for time_dir in time_dirs:
     
    4040##    ocean = scenario_dir + "map_work\\\port_hedland.gdb\\outlines\\initial_conditions_ocean"
    4141
    42     print 'Process: Create File GDB'
    43     gp.CreateFileGDB_management(folder, "raster")
     42##    print 'Process: Create File GDB'
     43##    gp.CreateFileGDB_management(folder, "raster")
    4444
    4545    gp.Workspace = raster_gbd
  • anuga_work/production/australia_ph2/lucinda/export_results_max.py

    r6835 r6939  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090408_153121_run_final_0_31982_2313_Tb__mhingee'
    25 time_dir2 = '20090408_202614_run_final_0_51311_2313_Tb__mhingee'
    26 time_dir3 = '20090409_022558_run_final_0_63735_2313_Tb__mhingee'
     24##time_dir1 = '20090408_153121_run_final_0_31982_2313_Tb__mhingee'
     25##time_dir2 = '20090408_202614_run_final_0_51311_2313_Tb__mhingee'
     26##time_dir3 = '20090409_022558_run_final_0_63735_2313_Tb__mhingee'
    2727
    28 time_dirs = [time_dir1, time_dir2, time_dir3]
     28time_dir1 = '20090421_140629_run_final_0_31982_2313_Tb_internal_lfountai'
     29
     30time_dirs = [time_dir1]#, time_dir2, time_dir3]
    2931
    3032cellsize = 250
  • anuga_work/production/australia_ph2/rockhampton/export_results_max.py

    r6845 r6939  
    2626time_dir3 = '20090409_192933_run_final_0_63737_3163_Tb__mhingee'
    2727
    28 time_dirs = [time_dir1, time_dir2, time_dir3]
     28time_dir1 = '20090421_142658_run_final_0_31939_3163_Tb_internal_lfountai'
     29
     30time_dirs = [time_dir1]#, time_dir2, time_dir3]
    2931
    3032cellsize = 250
     
    8284
    8385            name1 = join(directory, time_dir, project.scenario_name)
    84             name2 = join(directory, time_dir, project.scenario_name)+'_time_54300_0'
    85 
     86##            name2 = join(directory, time_dir, project.scenario_name)+'_time_54300_0'
     87            name2 = join(directory, time_dir, project.scenario_name)+'_time_54240_0'
     88           
    8689            names = [name1, name2]
    8790
  • anuga_work/production/australia_ph2/townsville/export_results_max.py

    r6749 r6939  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090306_131720_run_final_0_64448_1889_Tb_mhingee'
    25 time_dir2 = '20090306_165354_run_final_0_58331_1889_Tb_mhingee'
    26 time_dir3 = '20090306_165528_run_final_0_27347_1889_Tb_mhingee'
    27 ##time_dir4 = ''
     24time_dir1 = '20090414_104809_run_final_0_31965_2740_Tb__kvanputt'
     25time_dir2 = '20090420_064036_run_final_0_51469_2740_Tb__kvanputt'
     26time_dir3 = '20090421_141635_run_final_0_31965_2740_Tb_internal_lfountai'
     27time_dir4 = ''
    2828##time_dir5 = ''
    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 = ['elevation','stage']
     58var = ['stage']
    5959
    6060######
     
    8686
    8787            name1 = join(directory, time_dir, project.scenario_name)
    88             #name2 = join(directory, time_dir, project.scenario_name)+'_time_39600_0'
    89             #name3 = join(directory, time_dir, project.scenario_name)+'_time_79200_0'
     88            name2 = join(directory, time_dir, project.scenario_name)+'_time_24060_0'
     89            name3 = join(directory, time_dir, project.scenario_name)+'_time_48120_0'
    9090
    91             names = [name1] #, name2]#, name3]
     91            names = [name1, name2, name3]
    9292     
    93            # asc_name = []   
     93            asc_name = []   
    9494               
    9595            for name in names:
  • anuga_work/production/australia_ph2/townsville/project.py

    r6804 r6939  
    3232    event_number = int(sys.argv[1])
    3333else:   
    34     event_number = 31965    # the event number or the mux file name
    35 
    36 event_number_list = [31965, 51469, 65840] # To piggy back multiple events
     34    event_number = 62890    # the event number or the mux file name
     35
     36event_number_list = [31965, 51469, 62890] # To piggy back multiple events
    3737
    3838tide = 0                # difference between MSL and HAT
Note: See TracChangeset for help on using the changeset viewer.