Changeset 7364


Ignore:
Timestamp:
Aug 14, 2009, 10:30:16 AM (11 years ago)
Author:
Leharne
Message:
 
Location:
anuga_work/production/gold_coast_2009
Files:
1 added
7 edited

Legend:

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

    r7089 r7364  
    2929time_dir3 = '20090511_161526_run_final_1.1_51469_kvanputt'
    3030time_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'
     31time_dir5 = '20090519_160510_run_final_1.1_50994_lfountai'
     32time_dir6 = '20090521_220101_run_final_1.1_50863_kvanputt'
     33time_dir7 = '20090522_164526_run_final_0_51392_lfountai'
     34time_dir8 = '20090522_164640_run_final_1.1_51392_lfountai'
     35time_dir9 = '20090522_164948_run_final_0_51423_lfountai'
     36time_dir10 = '20090522_165600_run_final_1.1_51423_lfountai'
    3737
    38 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,
     39             time_dir6, time_dir7, time_dir8, time_dir9, time_dir10]
    3940
    4041for time_dir in time_dirs:
     
    4445    raster_gbd = folder + 'raster.gdb'
    4546    #contour = raster_gbd + '\\contour_dep'
    46 ##    land = scenario_dir + "ArcGIS\\gold_coast.gdb\\Coast_250"
     47    land = scenario_dir + "ArcGIS\\gold_coast.gdb\\Land_initial_conditions"
    4748##    ocean = scenario_dir + "map_work\\Perth.gdb\\Outlines\\initial_condition_ocean"
    4849   
     
    5657    #replication dictionary
    5758    replicate = (('gold_coast', ''), ('_time_29220_0', 'b'), ('_time_58440_0', 'c'),
     59                 ('_time_28860_0', 'b'), ('_time_57720_0', 'c'),
    5860                 ('_', ''), ('max','M_'), ('depth','_dep_'),
    5961                 ('speed', '_spe_'), ('elevation', '_ele_'), ('stage','_stage'))
     
    8486
    8587   
    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 ##
     88        output_extract = output_DEM + 'E'
     89        print 'Output Extract ',output_extract
     90   
     91        print 'Process: Extract by Mask'
     92        gp.ExtractByMask_sa(output_DEM, land, output_extract)
     93
    9294       
  • anuga_work/production/gold_coast_2009/For_DVD/project.py

    r7306 r7364  
    9191interior_regions_data = [['area_of_interest.csv', 500],
    9292                         ['intermediate.csv', 25000]]
    93 PriorityArea_filename = 'PriorityAreas.csv'
     93# If there are several priority areas of interest these can be defined
     94# in a single .csv file with headers 'easting, northing, id, value' where
     95# id defines each polygon, and the value is the mesh resolution
     96PriorityArea_filename = None
    9497
    9598# LAND - used to set the initial stage/water to be offcoast only
     
    137140# Environment variable names.
    138141# The inundation directory, not the data directory.
    139 ENV_INUNDATIONHOME = 'INUNDATIONHOME'
     142ENV_INUNDATIONHOME = 'ANUGADATA'
    140143
    141144#-------------------------------------------------------------------------------
     
    159162# create paths generated from environment variables.
    160163home = join(os.getenv(ENV_INUNDATIONHOME), 'data') # Absolute path for data folder
    161 muxhome = os.getenv(ENV_MUXHOME)
    162164   
    163165# check various directories/files that must exist
  • anuga_work/production/gold_coast_2009/For_DVD/run_model.py

    r7306 r7364  
    1 """Run a tsunami inundation scenario for Busselton, WA, Australia.
     1"""Run a tsunami inundation scenario for Gold Coast, Queensland, Australia.
    22
    33The scenario is defined by a triangular mesh created from project.polygon, the
  • anuga_work/production/gold_coast_2009/For_DVD/setup_model.py

    r7306 r7364  
    1111from os.path import join, exists
    1212from anuga.utilities.polygon import read_polygon, number_mesh_triangles
    13 
     13from anuga.shallow_water.data_manager import csv2polygons
    1414import project
    1515
  • anuga_work/production/gold_coast_2009/export_results_max.py

    r7089 r7364  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090521_220101_run_final_1.1_50863_kvanputt'
    25 time_dir2 = '20090522_164526_run_final_0_51392_lfountai'
    26 time_dir3 = '20090522_164640_run_final_1.1_51392_lfountai'
    27 time_dir4 = '20090522_164948_run_final_0_51423_lfountai'
    28 time_dir5 = '20090522_165600_run_final_1.1_51423_lfountai
     24time_dir1 = '20090605_120328_run_final_0_250m_51469_lfountai'
     25##time_dir1 = '20090522_164640_run_final_1.1_51392_lfountai'
    2926
    30 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]   
    3127
    32 ##cellsize = 250
    33 cellsize = 20
     28time_dirs = [time_dir1] #, time_dir2] #, time_dir3, time_dir4, time_dir5, time_dir6]   
     29
     30cellsize = 250
     31##cellsize = 20
    3432
    3533#timestep = None    # None means no timestep!
     
    8482
    8583            name1 = join(directory, time_dir, project.scenario_name)
    86             name2 = join(directory, time_dir, project.scenario_name)+'_time_29220_0'
    87             name3 = join(directory, time_dir, project.scenario_name)+'_time_58440_0'
     84##            name2 = join(directory, time_dir, project.scenario_name)+'_time_29220_0'
     85##            name3 = join(directory, time_dir, project.scenario_name)+'_time_58440_0'
    8886           
    89             names = [name1, name2, name3]
     87            names = [name1] #, name2, name3]
    9088     
    9189            asc_name = []   
     
    112110                asc_name.append(outname + '.asc')
    113111   
    114             maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
    115 
    116             print 'max asc outname ', maxasc_outname
    117             print 'asc_name ', str(asc_name)
    118            
    119             MaxAsc(maxasc_outname, asc_name)
     112##            maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
     113##
     114##            print 'max asc outname ', maxasc_outname
     115##            print 'asc_name ', str(asc_name)
     116##           
     117##            MaxAsc(maxasc_outname, asc_name)
  • anuga_work/production/gold_coast_2009/get_timeseries.py

    r7004 r7364  
    1919directory = project.output_folder
    2020
    21 time_dir1 = '20090327_135943_run_final_0_51469_lfountai'
    22 ##time_dir2 = '20090327_153206_run_final_0_51348_lfountai'
    23 ##time_dir3 = '20090327_153231_run_final_0_51253_lfountai'
    24 ##time_dir4 = '20090327_153304_run_final_0_50863_lfountai'
    25 time_dir5 = '20090327_153403_run_final_0_51423_lfountai'
    26 time_dir6 = '20090327_153602_run_final_0_62833_lfountai'
    27 time_dirs = [time_dir1, time_dir5, time_dir6] #time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
     21time_dir1 = '20090508_150215_run_final_0_51469_lfountai'
     22time_dir2 = '20090511_161526_run_final_1.1_51469_kvanputt'
     23time_dir3 = '20090511_165539_run_final_0_50863_lfountai'
     24time_dir4 = '20090521_220101_run_final_1.1_50863_kvanputt'
     25time_dir5 = '20090522_164526_run_final_0_51392_lfountai'
     26time_dir6 = '20090522_164640_run_final_1.1_51392_lfountai'
     27time_dir7 = '20090518_154710_run_final_0_50994_lfountai'
     28time_dir8 = '20090519_160510_run_final_1.1_50994_lfountai'
     29time_dir9 = '20090522_164948_run_final_0_51423_lfountai'
     30time_dir10 = '20090522_165600_run_final_1.1_51423_lfountai'
     31
     32time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5,
     33             time_dir6, time_dir7, time_dir8, time_dir9, time_dir10]
    2834
    2935for time_dir in time_dirs:
  • anuga_work/production/gold_coast_2009/project.py

    r7247 r7364  
    3434    event_number = int(sys.argv[1])
    3535else:   
    36     event_number = 51469    # the event number or the mux file name
     36    event_number = 51253    # the event number or the mux file name
    3737
    3838event_number_list = [51469, 51392, 50863] # To piggy back multiple events
Note: See TracChangeset for help on using the changeset viewer.