Changeset 7060


Ignore:
Timestamp:
May 20, 2009, 2:40:47 PM (15 years ago)
Author:
jgriffin
Message:
 
File:
1 edited

Legend:

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

    r7051 r7060  
    88from time import localtime, strftime, gmtime
    99from os.path import join, exists
     10import csv
    1011
    1112
     
    1718# Note, the user needs to set up the directory system accordingly
    1819state = 'new_south_wales'
    19 scenario_name = 'batemans_bay'
    20 scenario_folder = 'batemans_bay_tsunami_scenario_2009'
     20scenario_name = 'gosford'
     21scenario_folder = 'gosford_tsunami_scenario_2009'
    2122
    2223#-------------------------------------------------------------------------------
    2324# Initial Conditions
    2425#-------------------------------------------------------------------------------
     26
    2527
    2628# Model specific parameters.
    2729# One or all can be changed each time the run_model script is executed
    28 tide = 0.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
    37 ##event_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!!!
     30tide = 1.0                # difference between MSL and HAT in metres
     31zone = 56               # specify zone of model
     32event_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
     38event_number = 58025    # Puysegur 1 in 200
     39
    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 = ['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']
     71ascii_grid_filenames = ['e151_s34_clip2',
     72                        'aus197_topo2raster_clip',
     73                        'hydro_topo2raster_clip']
     74                        #'aus197_topo2raster_terrigal']   #GA data??
    8775                   
     76
     77
     78#Get long list of Lidar filenames.
     79ENV_INUNDATIONHOME = 'INUNDATIONHOME'
     80home = join(os.getenv(ENV_INUNDATIONHOME), 'data')
     81anuga_folder = join(home, state, scenario_folder, 'anuga')
     82topographies_folder = join(anuga_folder, 'topographies/original')
     83
     84file_list = join(topographies_folder, 'file_list.csv')
     85###file_list2 = join(topographies_folder, 'file_list2.csv')
     86fid = csv.reader(open(file_list))
     87###fid2 = csv.reader(open(file_list2))
     88lidar_filenames = []
     89for 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
    8897# Format for point is x,y,elevation (with header)
    89 point_filenames = ['SD100031996_jgriffin_clip.csv',
    90                    'tomaga_offshore_AHD_MGA_1997.csv',
    91                    'Batemans_BBHD_MGA_1995.csv',
    92                    'moruya_AHD_MGA_2000.csv']
     98point_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                   ]
     111point_filenames.extend(lidar_filenames)
     112#print point_filenames
    93113         
    94114
     
    108128# Format for points easting,northing (no header)
    109129
    110 interior_regions_data = [['area_of_interest.csv', 500],
    111                         ['area_of_significance.csv', 2500],
    112                         ['shallow_water.csv', 10000]]
     130interior_regions_data =[['aoi_umina.csv', 500],
     131                        ['aoi_terrigal.csv', 500],
     132                        ['aos_umina.csv', 2500],
     133                        ['aos_terrigal.csv', 2500]]
    113134
    114135
     
    135156# Thinned ordering file from Hazard Map (geographic)
    136157# Format is index,latitude,longitude (with header)
    137 urs_order_filename = 'thinned_boundary_ordering_extend.csv'
     158urs_order_filename = 'urs_order.csv'
    138159
    139160# Landward bounding points
    140161# Format easting,northing (no header)
    141 landward_boundary_filename = 'landward_boundary_extend.csv'
     162landward_boundary_filename = 'landward_boundary.csv'
    142163
    143164# MUX input filename.
     
    211232# check various directories/files that must exist
    212233anuga_folder = join(home, state, scenario_folder, 'anuga')
    213 topographies_folder = join(anuga_folder, 'topographies')
     234topographies_folder = join(anuga_folder, 'topographies/original')
    214235polygons_folder = join(anuga_folder, 'polygons')
    215236boundaries_folder = join(anuga_folder, 'boundaries')
Note: See TracChangeset for help on using the changeset viewer.