Changeset 7070


Ignore:
Timestamp:
May 25, 2009, 11:28:13 AM (11 years ago)
Author:
jgriffin
Message:
 
File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/new_south_wales/batemans_bay/project.py

    r7067 r7070  
    88from time import localtime, strftime, gmtime
    99from os.path import join, exists
    10 import csv
    1110
    1211
     
    1817# Note, the user needs to set up the directory system accordingly
    1918state = 'new_south_wales'
    20 scenario_name = 'gosford'
    21 scenario_folder = 'gosford_tsunami_scenario_2009'
     19scenario_name = 'batemans_bay'
     20scenario_folder = 'batemans_bay_tsunami_scenario_2009'
    2221
    2322#-------------------------------------------------------------------------------
    2423# Initial Conditions
    2524#-------------------------------------------------------------------------------
    26 
    2725
    2826# Model specific parameters.
    2927# One or all can be changed each time the run_model script is executed
    30 tide = 1.0                # difference between MSL and HAT in metres
    31 zone = 56               # specify zone of model
    32 event_number = 58242    # Puysegur 1 in 10 000 # the event number or the mux file name
    33 #event_number = 51436    # New Hebrides 1 in 10 000
    34 #event_number = 58349    # Puysegur 1 in 5000
    35 #event_number = 58284    # Puysegur 1 in 2000
    36 #event_number = 58187    # Puysegur 1 in 1000
    37 #event_number = 58113    # Puysegur 1 in 500
    38 ##event_number = 58025    # Puysegur 1 in 200
    39 
     28tide = 1.0                # difference between MSL and HAT (1.0)
     29
     30# the event number or the mux file name
     31##event_number = 58129  #1 in 200 yr Puyesgur
     32##event_number = 58115    #1 in 500 yr
     33##event_number = 58226    #1 in 1000 yr
     34##event_number = 58284    #1 in 2000 yr
     35##event_number = 58286    #1 in 5000 yr
     36##event_number = 58272      #1 in 10000 yr Puysegur
     37event_number = 51445     #1 in 10000 yr New Hebrides
     38##event_number = 58368    #1 in 100 000 yr Puysegur!!!
     39##event_number = 51436    #1 in 100 000 yr New Hebrides!!!
    4040
    4141alpha = 0.1             # smoothing parameter for mesh
     
    6969# Used in build_elevation.py
    7070# Format for ascii grids, as produced in ArcGIS + a projection file
    71 ascii_grid_filenames = ['e151_s34_clip2',
    72                         'aus197_topo2raster_clip',
    73                         'hydro_topo2raster_clip']
    74                         #'aus197_topo2raster_terrigal']   #GA data??
     71ascii_grid_filenames = ['1a',# Topographic data
     72                        '1b',
     73                        '2b',
     74                        '2a_3',
     75                        '3b',
     76                        '3a',
     77                        '4a_2',
     78                        '4b',
     79                        'off1',
     80                        'off2',
     81                        'off3',
     82                        'bbhd',
     83                        'sd100031996_p',
     84                        'sd100031996_p2',
     85                        'sd100031996_p3',
     86                        'sd100031996_p4']
    7587                   
    76 
    77 
    78 #Get long list of Lidar filenames.
    79 ENV_INUNDATIONHOME = 'INUNDATIONHOME'
    80 home = join(os.getenv(ENV_INUNDATIONHOME), 'data')
    81 anuga_folder = join(home, state, scenario_folder, 'anuga')
    82 topographies_folder = join(anuga_folder, 'topographies/original')
    83 
    84 file_list = join(topographies_folder, 'file_list.csv')
    85 ###file_list2 = join(topographies_folder, 'file_list2.csv')
    86 fid = csv.reader(open(file_list))
    87 ###fid2 = csv.reader(open(file_list2))
    88 lidar_filenames = []
    89 for row in fid:
    90     lidar_filenames.append(','.join(row))
    91 ##for row in fid2:
    92 ##    lidar_filenames.append(','.join(row))
    93 ###print lidar_filenames
    94 #lidar_filenames = list(lidar_filenames)
    95 #print lidar_filenames
    96 
    9788# Format for point is x,y,elevation (with header)
    98 point_filenames = ['estuaries.txt',
    99 ##                   'XYHI341_HSDB_SD100035103',
    100 ##                   'XYAHD',
    101 ##                   'XYSCHOOL_11_99_HSDB_SD100021002_MGA',
    102                    'ENT5C1S08_03_AHD_prepared.txt',
    103                    'ENT5C1S09_03_AHD_prepared.txt',
    104                    'ENT5C1S10_03_AHD_prepared.txt',
    105 ##                   'XYAHD',
    106 ##                   'XYHI341_HSDB_SD100035103_MGA',
    107 ##                   'XYHI339_F_HSDB_SD100031669_MGA',
    108 ##                   'AUS197_MGA_AHD_1972.txt',
    109                    'hydro_neg.txt'
    110                    ]
    111 point_filenames.extend(lidar_filenames)
    112 #print point_filenames
     89point_filenames = ['SD100031996_jgriffin_clip.csv',
     90                   'tomaga_offshore_AHD_MGA_1997.csv',
     91                   'Batemans_BBHD_MGA_1995.csv',
     92                   'moruya_AHD_MGA_2000.csv']
    11393         
    11494
     
    128108# Format for points easting,northing (no header)
    129109
    130 interior_regions_data =[['aoi_umina.csv', 500],
    131                         ['aoi_terrigal.csv', 500],
    132                         ['aos_umina.csv', 2500],
    133                         ['aos_terrigal.csv', 2500]]
     110interior_regions_data = [['area_of_interest.csv', 500],
     111                        ['area_of_significance.csv', 2500],
     112                        ['shallow_water.csv', 10000]]
    134113
    135114
     
    156135# Thinned ordering file from Hazard Map (geographic)
    157136# Format is index,latitude,longitude (with header)
    158 urs_order_filename = 'urs_order.csv'
     137urs_order_filename = 'thinned_boundary_ordering_extend.csv'
    159138
    160139# Landward bounding points
    161140# Format easting,northing (no header)
    162 landward_boundary_filename = 'landward_boundary.csv'
     141landward_boundary_filename = 'landward_boundary_extend.csv'
    163142
    164143# MUX input filename.
     
    232211# check various directories/files that must exist
    233212anuga_folder = join(home, state, scenario_folder, 'anuga')
    234 topographies_folder = join(anuga_folder, 'topographies/original')
     213topographies_folder = join(anuga_folder, 'topographies')
    235214polygons_folder = join(anuga_folder, 'polygons')
    236215boundaries_folder = join(anuga_folder, 'boundaries')
Note: See TracChangeset for help on using the changeset viewer.