Changeset 7677


Ignore:
Timestamp:
Apr 9, 2010, 3:16:13 PM (10 years ago)
Author:
fountain
Message:

updates to mandurah storm surge model

Location:
branches/anuga_1_1/anuga_work/production/mandurah_storm_surge_2009
Files:
2 added
8 edited

Legend:

Unmodified
Added
Removed
  • branches/anuga_1_1/anuga_work/production/mandurah_storm_surge_2009/Arc_asc2raster_GDA94z50.py

    r7614 r7677  
    2525output_dir = "anuga\\outputs\\"
    2626
    27 time_dir1 = '20091229_100831_run_storm_surge_final_0_alby_waves_lfountai'
     27time_dir1 = '20100219_094651_run_storm_surge_final_0_alby_waves_lfountai'
     28# time_dir2 = '20100223_153032_run_storm_surge_final_0_alby_waves_GEMS_compare_lfountai'
    2829
    29 
    30 time_dirs = [time_dir1] ##, time_dir2, time_dir3, time_dir4, time_dir5,
    31              ##time_dir6, time_dir7, time_dir8, time_dir9, time_dir10]
     30time_dirs = [time_dir1] #, time_dir2]
    3231
    3332for time_dir in time_dirs:
     
    3635    folder =  scenario_dir + output_dir + time_dir + '\\'
    3736    raster_gbd = folder + 'raster.gdb'
    38     elevation_gdb = folder + 'elevation.gdb'
    39     #contour = raster_gbd + '\\contour_dep'
    40     land = scenario_dir + "ArcGIS\\data.gdb\\initial_conditions_10m"
    41 ##    ocean = scenario_dir + "map_work\\Perth.gdb\\Outlines\\initial_condition_ocean"
     37    land = scenario_dir + "ArcGIS\\data.gdb\\initial_conditions_10m_clip"
     38    ocean = scenario_dir + "ArcGIS\\data.gdb\\initial_conditions_ocean"
    4239   
    4340    print 'Process: Create File GDB'
    44 ##    gp.createFileGDB_management(folder, "elevation")
    45 ##    gp.CreateFileGDB_management(folder, "raster")
    46 
     41    # gp.CreateFileGDB_management(folder, "raster")
    4742    gp.Workspace = raster_gbd
    48 ##    gp.Workspace = elevation_gdb
    4943
    5044    print gp.Workspace
    5145   
    5246    #replication dictionary
    53     replicate = (('gold_coast', ''), ('_time_29220_0', 'b'), ('_time_58440_0', 'c'),
     47    replicate = (('alby_waves', 'man'), ('_time_29220_0', 'b'), ('_time_58440_0', 'c'),
    5448                 ('_time_28860_0', 'b'), ('_time_57720_0', 'c'),
    5549                 ('_', ''), ('max','M_'), ('depth','_dep_'),
     
    5751
    5852    generate_filename = []
    59  ##   input_ascii = glob.glob(folder + '*elevation*.asc')
    60     input_ascii = glob.glob(folder + '*depth*.asc')
     53    input_ascii = glob.glob(folder + '*elevation*.asc')
     54    # input_ascii = glob.glob(folder + '*stage.asc')
    6155    print time_dir
    6256
     
    8680   
    8781        print 'Process: Extract by Mask'
    88         gp.ExtractByMask_sa(output_DEM, land, output_extract)
     82        # gp.ExtractByMask_sa(output_DEM, land, output_extract) # use for depth to get onshore inundation depths
     83        # gp.ExtractByMask_sa(output_DEM, ocean, output_extract) # use for stage to get offshore wave heights
    8984
    9085       
  • branches/anuga_1_1/anuga_work/production/mandurah_storm_surge_2009/build_elevation.py

    r7647 r7677  
    8585# G2 = geospatial_data['ph10m_ss'].clip_outside(project.elevation_clip_box)
    8686
    87 G1 = geospatial_data['m_peel_aoi'].clip(project.elevation_clip_box)
    88 G2 = geospatial_data['ph10m_ss'].clip_outside(project.elevation_clip_box)
    89 G3 = geospatial_data['DpiU1A03_ss.csv']
    90 G4 = geospatial_data['MA-46893-SNDS_AHD.csv']
    91 G5 = geospatial_data['MS0205HY_AHD.csv']
    92 G6 = geospatial_data['MS0404_AHD.csv']
    93 G7 = geospatial_data['YU0403HY_AHD.csv']
    94 G8 = geospatial_data['original_data_ss.csv']
     87G1 = geospatial_data[keylist[0]].clip(project.elevation_clip_box)
     88G2 = geospatial_data[keylist[1]].clip_outside(project.elevation_clip_box)
     89G3 = geospatial_data[keylist[2]]
     90G4 = geospatial_data[keylist[3]]
     91G5 = geospatial_data[keylist[4]]
     92G6 = geospatial_data[keylist[5]]
     93G7 = geospatial_data[keylist[6]]
     94G8 = geospatial_data[keylist[7]]
    9595
    9696G = G1 + G2 + G3 + G4 + G5 + G6 + G7 + G8
     
    102102# Use for comparision in ARC
    103103G.export_points_file(project.combined_elevation + '.txt')
    104 G1.export_points_file(join(project.topographies_folder, (str(project.ascii_grid_filenames[0]) +'_export.txt')))
    105 G2.export_points_file(join(project.topographies_folder, (str(project.ascii_grid_filenames[1]) + '_export.txt')))
     104G1.export_points_file(join(project.topographies_folder, keylist[0] +'_export.txt'))
     105G2.export_points_file(join(project.topographies_folder, keylist[1] +'_export.txt'))
  • branches/anuga_1_1/anuga_work/production/mandurah_storm_surge_2009/export_results_max.py

    r7614 r7677  
    2121
    2222directory = project.output_folder
    23 time_dir1 = '20091229_100831_run_storm_surge_final_0_alby_waves_lfountai'
     23#time_dir1 = '20100219_094651_run_storm_surge_final_0_alby_waves_lfountai'
     24time_dir1 = '20100223_153032_run_storm_surge_final_0_alby_waves_GEMS_compare_lfountai'
    2425
    25 time_dirs = [time_dir1]
     26time_dirs = [time_dir1] #, time_dir2]
    2627
    2728# sww filename extensions ie. hobart_time_17640_0.sww, input into list 17640
     
    4546
    4647# one or more key strings from var_equations below
    47 var = ['stage']
     48var = ['elevation', 'speed'] #'stage','depth','speed','elevation']
    4849
    4950######
     
    100101                sww2dem(name, basename_out = outname,
    101102                            quantity = quantityname,
    102                             timestep = timestep,
     103                            reduction = timestep,
    103104                            cellsize = cellsize,     
    104105                            easting_min = easting_min,
     
    106107                            northing_min = northing_min,
    107108                            northing_max = northing_max,       
    108                             reduction = max,
    109109                            verbose = True,
    110110                            format = 'asc')
  • branches/anuga_1_1/anuga_work/production/mandurah_storm_surge_2009/gems_comparison/build_elevation_GEMS.py

    r7647 r7677  
    8282print 'project_GEMS.elevation_clip_box', project_GEMS.elevation_clip_box
    8383
    84 G1 = geospatial_data['man_25m_prj'].clip(project_GEMS.elevation_clip_box)
    85 G2 = geospatial_data['man10m_ss'].clip_outside(project_GEMS.elevation_clip_box)
     84G1 = geospatial_data[keylist[0]].clip(project.elevation_clip_box)
     85G2 = geospatial_data[keylist[1]].clip_outside(project.elevation_clip_box)
    8686
    8787G = G1 + G2
     
    9393# Use for comparision in ARC
    9494G.export_points_file(project_GEMS.combined_elevation + '.txt')
    95 G1.export_points_file(join(project_GEMS.topographies_folder, (str(project_GEMS.ascii_grid_filenames[0]) +'_export.txt')))
    96 G2.export_points_file(join(project_GEMS.topographies_folder, (str(project_GEMS.ascii_grid_filenames[1]) + '_export.txt')))
     95# G1.export_points_file(join(project.topographies_folder, keylist[0] +'_export.txt'))
     96# G2.export_points_file(join(project.topographies_folder, keylist[1] +'_export.txt'))
  • branches/anuga_1_1/anuga_work/production/mandurah_storm_surge_2009/gems_comparison/project_GEMS.py

    r7647 r7677  
    180180# Output filename for elevation
    181181# this is a combination of all the data generated in build_elevation.py
    182 combined_elevation_basename = scenario_name + 'GEMS_combined_elevation'
     182combined_elevation_basename = scenario_name + '_GEMS_combined_elevation'
    183183
    184184#-------------------------------------------------------------------------------
  • branches/anuga_1_1/anuga_work/production/mandurah_storm_surge_2009/gems_comparison/run_model_GEMS.py

    r7647 r7677  
    4545##from anuga.shallow_water.data_manager import urs2sts
    4646from anuga.utilities.polygon import read_polygon, Polygon_function
     47from anuga.utilities.system_tools import get_revision_number
    4748
    4849# Application specific imports
     
    6061                             project_GEMS.__name__+'.py'))
    6162start_screen_catcher(project_GEMS.output_run, 0, 1)
     63
     64print 'SVN revision number: ', get_revision_number()
    6265
    6366#-------------------------------------------------------------------------------
  • branches/anuga_1_1/anuga_work/production/mandurah_storm_surge_2009/project.py

    r7647 r7677  
    128128
    129129# GEMS order filename
    130 # Format is index,northing, easting, elevation (without header)
     130# Format is index, longitude, latitude, elevation (without header)
    131131gems_order_filename = 'gems_boundary_order_thinned.csv'
    132132
  • branches/anuga_1_1/anuga_work/production/mandurah_storm_surge_2009/run_model.py

    r7647 r7677  
    4545##from anuga.shallow_water.data_manager import urs2sts
    4646from anuga.utilities.polygon import read_polygon, Polygon_function
     47from anuga.utilities.system_tools import get_revision_number
    4748
    4849# Application specific imports
     
    6061                             project.__name__+'.py'))
    6162start_screen_catcher(project.output_run, 0, 1)
     63
     64print 'SVN revision number: ', get_revision_number()
    6265
    6366#-------------------------------------------------------------------------------
Note: See TracChangeset for help on using the changeset viewer.