Changeset 7089


Ignore:
Timestamp:
May 26, 2009, 3:52:17 PM (11 years ago)
Author:
Leharne
Message:

updates to gold coast inundation scenario

Location:
anuga_work/production/gold_coast_2009
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/gold_coast_2009/Arc_asc2raster_GDA94z56.py

    r7004 r7089  
    2525output_dir = "anuga\\outputs\\"
    2626
    27 time_dir1 = '20090327_135943_run_final_0_51469_lfountai'
    28 ##time_dir2 = '20090327_153206_run_final_0_51348_lfountai'
    29 ##time_dir3 = '20090327_153231_run_final_0_51253_lfountai'
    30 ##time_dir4 = '20090327_153304_run_final_0_50863_lfountai'
    31 ##time_dir5 = '20090327_153403_run_final_0_51423_lfountai'
    32 ##time_dir6 = '20090327_153602_run_final_0_62833_lfountai'
     27time_dir1 = '20090508_150215_run_final_0_51469_lfountai'
     28time_dir2 = '20090511_165539_run_final_0_50863_lfountai'
     29time_dir3 = '20090511_161526_run_final_1.1_51469_kvanputt'
     30time_dir4 = '20090518_154710_run_final_0_50994_lfountai'
     31##time_dir5 = '20090519_160510_run_final_1.1_50994_lfountai'
     32##time_dir6 = '20090521_220101_run_final_1.1_50863_kvanputt'
     33##time_dir7 = '20090522_164526_run_final_0_51392_lfountai'
     34##time_dir8 = '20090522_164640_run_final_1.1_51392_lfountai'
     35##time_dir9 = '20090522_164948_run_final_0_51423_lfountai'
     36##time_dir10 = '20090522_165600_run_final_1.1_51423_lfountai'
    3337
    34 time_dirs = [time_dir1] #[time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
     38time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4] #, time_dir5, time_dir6]
    3539
    3640for time_dir in time_dirs:
     
    4044    raster_gbd = folder + 'raster.gdb'
    4145    #contour = raster_gbd + '\\contour_dep'
    42     land = scenario_dir + "ArcGIS\\gold_coast.gdb\\Coast_250"
     46##    land = scenario_dir + "ArcGIS\\gold_coast.gdb\\Coast_250"
    4347##    ocean = scenario_dir + "map_work\\Perth.gdb\\Outlines\\initial_condition_ocean"
    4448   
     
    5155   
    5256    #replication dictionary
    53     replicate = (('gold_coast', ''), ('_time_21960_0', 'b'), ('_time_43920_0', 'c'),
    54                  ('time_65880_0', 'd'), ('_', ''), ('max','M_'), ('depth','_dep_'),
     57    replicate = (('gold_coast', ''), ('_time_29220_0', 'b'), ('_time_58440_0', 'c'),
     58                 ('_', ''), ('max','M_'), ('depth','_dep_'),
    5559                 ('speed', '_spe_'), ('elevation', '_ele_'), ('stage','_stage'))
    5660
    5761    generate_filename = []
    58     input_ascii = glob.glob(folder + '*depth_max.asc')
     62    input_ascii = glob.glob(folder + '*_max.asc')
    5963    print time_dir
    6064
     
    8084
    8185   
    82         output_extract = output_DEM + 'E'
    83         print 'Output Extract ',output_extract
    84    
    85         print 'Process: Extract by Mask'
    86         gp.ExtractByMask_sa(output_DEM, land, output_extract)
    87 
     86##        output_extract = output_DEM + 'E'
     87##        print 'Output Extract ',output_extract
     88##   
     89##        print 'Process: Extract by Mask'
     90##        gp.ExtractByMask_sa(output_DEM, land, output_extract)
     91##
    8892       
  • anuga_work/production/gold_coast_2009/export_results_max.py

    r7004 r7089  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090327_135943_run_final_0_51469_lfountai'
    25 ##time_dir2 = '20090327_153206_run_final_0_51348_lfountai'
    26 ##time_dir3 = '20090327_153231_run_final_0_51253_lfountai'
    27 ##time_dir4 = '20090327_153304_run_final_0_50863_lfountai'
    28 ##time_dir5 = '20090327_153403_run_final_0_51423_lfountai'
    29 ##time_dir6 = '20090327_153602_run_final_0_62833_lfountai'
    30 time_dirs = [time_dir1] #time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]   
     24time_dir1 = '20090521_220101_run_final_1.1_50863_kvanputt'
     25time_dir2 = '20090522_164526_run_final_0_51392_lfountai'
     26time_dir3 = '20090522_164640_run_final_1.1_51392_lfountai'
     27time_dir4 = '20090522_164948_run_final_0_51423_lfountai'
     28time_dir5 = '20090522_165600_run_final_1.1_51423_lfountai
     29
     30time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]   
    3131
    3232##cellsize = 250
     
    4141
    4242#area = ['Bunbury', 'Busselton']  # strings must match keys in var_equations below
    43 area = ['Dune']      # 'All' means no special areas - the whole thing
     43area = ['All']      # 'All' means no special areas - the whole thing
    4444
    4545######
     
    5454
    5555# one or more key strings from var_equations above
    56 var = ['depth']
     56var = ['speed', 'depth']
    5757
    5858######
     
    8484
    8585            name1 = join(directory, time_dir, project.scenario_name)
    86             name2 = join(directory, time_dir, project.scenario_name)+'_time_21960_0'
    87             name3 = join(directory, time_dir, project.scenario_name)+'_time_43920_0'
    88             #name4 = join(directory, time_dir, project.scenario_name)+'_time_65880_0'
    89 
     86            name2 = join(directory, time_dir, project.scenario_name)+'_time_29220_0'
     87            name3 = join(directory, time_dir, project.scenario_name)+'_time_58440_0'
     88           
    9089            names = [name1, name2, name3]
    9190     
  • anuga_work/production/gold_coast_2009/project.py

    r7004 r7089  
    3434    event_number = int(sys.argv[1])
    3535else:   
    36     event_number = 51469    # the event number or the mux file name
     36    event_number = 51423    # the event number or the mux file name
    3737
    3838event_number_list = [51469, 50863] # To piggy back multiple events
    3939
    40 tide = 0
     40tide = 1.1
    4141alpha = 0.1             # smoothing parameter for mesh
    4242friction=0.01           # manning's friction coefficient
     
    106106building_exposure_filename = 'gold_coast_res_clip.csv' # from NEXIS
    107107
     108# AREA OF IMAGES - Extent of each image to find out highest runup
     109# Headerr: easting, northing, id, value
     110# Used in get_runup.py
     111images_filename = 'images.csv'
     112
    108113# BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively
    109114# NOTE: when files are put together the points must be in sequence
    110115# For ease go clockwise!
    111116# Check the run_model.py for boundary_tags
     117
    112118
    113119# Thinned ordering file from Hazard Map (geographic)
     
    244250    building_exposure = join(gauges_folder, building_exposure_filename)
    245251
     252# The absolute pathname for the image file
     253# Used for get_runup.py
     254if images_filename:
     255    images = join(polygons_folder, images_filename)
     256
    246257# full path to where MUX files (or meta-files) live
    247258mux_input = join(event_folder, mux_input_filename)
Note: See TracChangeset for help on using the changeset viewer.