Changeset 7611


Ignore:
Timestamp:
Jan 22, 2010, 3:52:44 PM (9 years ago)
Author:
monica
Message:

updated for mandurah

Files:
5 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/wa/mandurah_2009/build_elevation.py

    r7558 r7611  
    5252if not project.ascii_grid_filenames == []:
    5353    for filename in project.ascii_grid_filenames:
     54        print 'hello', filename
     55        print 'hello', project.topographies_folder
    5456        absolute_filename = join(project.topographies_folder, filename)
    5557        convert_dem_from_ascii2netcdf(absolute_filename,
  • anuga_work/production/wa/mandurah_2009/export_results_max.py

    r7545 r7611  
    2323
    2424
    25 time_dir1 = '20090324_134328_run_final_0.6_27255_Bt_kvanputt'
    26 time_dir2 = '20090324_140747_run_final_0.6_68693_Bt_kvanputt'
    27 time_dir3 = '20090327_152717_run_final_0.6_27283_Bt_kvanputt'
    28 time_dir4 = '20090327_160646_run_final_0_27255_Bt_kvanputt'
    29 time_dir5 = '20090327_161520_run_final_0_68693_Bt_kvanputt'
    30 time_dir6 = '20090330_120458_run_final_0_27283_Bt_kvanputt'
     25time_dir1 = '20100121_170639_run_final_0.6_27255_mosuchow'
     26
    3127
    3228time_dirs = [time_dir1] #1,time_dir2, time_dir3, time_dir4, time_dir5]
     
    3430
    3531##cellsize = 250
    36 cellsize = 20
     32cellsize = 25
    3733
    3834#timestep = None    # None means no timestep!
     
    4339######
    4440
    45 area = ['Bunbury', 'Busselton']  # strings must match keys in var_equations below
     41area = ['AoiR']  # strings must match keys in var_equations below
    4642#area = ['All']      # 'All' means no special areas - the whole thing
    4743
     
    5753
    5854# one or more key strings from var_equations above
    59 var = ['elevation']
     55var = ['stage']
    6056
    6157######
  • anuga_work/production/wa/mandurah_2009/project.py

    r7558 r7611  
    3232friction=0.01           # manning's friction coefficient
    3333starttime=0             # start time for simulation
    34 finaltime=80000        # final time for simulation
     34finaltime=10         #80000 # final time for simulation
    3535yieldstep=60            # timestep to save output
    3636
    37 setup = 'trial'         # This can be one of three values
     37setup = 'final'         # This can be one of three values
    3838                        #    trial - coarsest mesh, fast
    3939                        #    basic - coarse mesh
     
    6161# Format for ascii grids, as produced in ArcGIS + a projection file
    6262ascii_grid_filenames = ['m_harvey_10m',   # Topo
    63                         'm_peel_10m']
     63                        'm_peel_10m',
     64                        'peel_aoi_1',
     65                        'peel_aoi_2',
     66                        'peel_aoi_3']
    6467                        #'extract_m_peel' # mandurah Topo
    6568
    6669# Format for point is x,y,elevation (with header)
    67 point_filenames = ['bathymetry_all_data.txt',
     70point_filenames = ['bathy_update.txt',
    6871                   'DpiU1A03.txt',
    6972                   'MA-46893-SNDS_AHD.txt',
     
    8790# Used in run_model.py
    8891# Format for points easting,northing (no header)                   
    89 interior_regions_data = [['area_of_interest1.csv', 500],
    90                          ['area_of_interest2.csv', 500],
    91                          ['area_of_interest3.csv', 500],
     92interior_regions_data = [['aoi_1.csv', 500],
     93                         ['aoi_2.csv', 500],
     94                         ['aoi_3.csv', 500],
     95                         ['aoi_4.csv', 500],
     96                         ['area_of_interest_2.csv', 500],
     97                         ['area_of_interest_3.csv', 500],
     98                         ['aoi_explode_2.csv', 500],
     99                         ['aoi_explode_3.csv', 500],
     100                         ['aoi_explode_4.csv', 500],
     101                         ['aoi_explode_5.csv', 500],
     102                         ['aoi_explode_6.csv', 500],
     103                         ['aoi_explode_7.csv', 500],
     104                         ['aoi_explode_8.csv', 500],
     105                         ['aoi_explode_9.csv', 500],
     106                         ['aoi_explode_10.csv', 500],
     107                         ['aoi_explode_11.csv', 500],
     108                         ['aoi_explode_12.csv', 500],
     109                         ['aoi_clip2a_explode1.csv', 500],
     110                         ['aoi_clip2a_explode2.csv', 500],
     111                         ['aoi_clip2a_explode3.csv', 500],
     112                         ['aoi_clip2a_explode4.csv', 500],
     113                         ['aoi_clip2a_explode5.csv', 500],
     114                         ['aoi_clip2a_explode6.csv', 500],
    92115                         ['bounding_area_of_interest.csv', 2500]]
    93116
     
    141164clip_box = 'clip_box.csv'
    142165
    143 # ASCII export grid for mandurah
    144 xminmandurah = 377708
    145 xmaxmandurah = 380186
    146 yminmandurah = 6399699
    147 ymaxmandurah = 6401544
     166# ASCII export grid for mandurah for inundation maps
     167#xminRegion = 369013
     168#xmaxRegion = 385914
     169#yminRegion = 6390118
     170#ymaxRegion = 6415260
     171
     172xminAoiR = 370407
     173xmaxAoiR = 385745
     174yminAoiR = 6390330
     175ymaxAoiR = 6405626
     176
     177#xminAoi1 = 380524
     178#xmaxAoi1 = 382123
     179#yminAoi1 = 6401707
     180#ymaxAoi1 = 6405240
     181
     182#xminAoi2 = 377764
     183#xmaxAoi2 = 381493
     184#yminAoi2 = 6396471
     185#ymaxAoi2 = 6401692
     186
     187#xminAoi3 = 371074
     188#xmaxAoi3 = 374219
     189#yminAoi3 = 6390639
     190#ymaxAoi3 = 6394731
     191
    148192
    149193
  • misc/tools/event_selection/find_fault/find_fault.py

    r7378 r7611  
    77
    88
    9 search_dist = 0.75 # radius in degrees to search
     9search_dist = 2.0 # radius in degrees to search
    1010
    1111# List coordinates to search for
    12 MaxLon = 120.2605
    13 MinLon = 113.7402
    14 MaxLat = -9.4653
    15 MinLat = -11.2732
     12MaxLon = 99.8
     13MinLon = 97.9
     14MaxLat = -0.6
     15MinLat = -2.7
    1616
    1717# Specify input Folder
    18 inDir = r'/nas/gemd/georisk_models/inundation/data/western_australia/broome_tsunami_scenario_2009/anuga/boundaries/Results_Australia_1629_0.01_15.00/'
    19 inFile = 'fault.xy'
    20 #inFile = 'fault.txt'
     18inDir = r'/nas/gemd/georisk/presentations/tsunamis/FESA_Dec09_Workshop/Results_Australia_1032_0.01_15.00/'
     19#inFile = 'fault.xy'
     20inFile = 'fault.txt'
    2121inFilePath = file(inDir + inFile)
    2222
Note: See TracChangeset for help on using the changeset viewer.