Changeset 6958


Ignore:
Timestamp:
May 5, 2009, 11:42:46 AM (11 years ago)
Author:
jgriffin
Message:

added functionality to read list of point filenames

File:
1 edited

Legend:

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

    r6700 r6958  
    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                # 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 = 58242    #1 in 10000 yr Puysegur
    37 event_number = 58368    #1 in 10000 yr Puysegur
    38 ##event_number = 51436    #1 in 10000 yr New Hebrides
    39 
     30tide = 0                # difference between MSL and HAT in metres
     31zone = 56               # specify zone of model
     32event_number = 31891    # the event number or the mux file name
    4033alpha = 0.1             # smoothing parameter for mesh
    4134friction=0.01           # manning's friction coefficient
    4235starttime=0             # start time for simulation
    43 finaltime=10000         # final time for simulation
    44 
    45 setup = 'final'         # This can be one of three values
     36finaltime=1000         # final time for simulation
     37
     38setup = 'trial'         # This can be one of three values
    4639                        #    trial - coarsest mesh, fast
    4740                        #    basic - coarse mesh
     
    5952#-------------------------------------------------------------------------------
    6053
    61 output_comment = [setup, tide, event_number, 'viz']
     54output_comment = [setup, tide, event_number, 'lidar']
    6255
    6356#-------------------------------------------------------------------------------
     
    6861# Used in build_elevation.py
    6962# Format for ascii grids, as produced in ArcGIS + a projection file
    70 ascii_grid_filenames = ['1a',# Topographic data
    71                         '1b',
    72                         '2b',
    73                         '2a_3',
    74                         '3b',
    75                         '3a',
    76                         '4a_2',
    77                         '4b',
    78                         'off1',
    79                         'off2',
    80                         'off3',
    81                         'bbhd',
    82                         'sd100031996_p',
    83                         'sd100031996_p2',
    84                         'sd100031996_p3',
    85                         'sd100031996_p4']
     63ascii_grid_filenames = []
    8664                   
     65
     66
     67#Get long list of Lidar filenames.
     68ENV_INUNDATIONHOME = 'INUNDATIONHOME'
     69home = join(os.getenv(ENV_INUNDATIONHOME), 'data')
     70anuga_folder = join(home, state, scenario_folder, 'anuga')
     71topographies_folder = join(anuga_folder, 'topographies/original')
     72file_list = join(topographies_folder, 'file_list.csv')
     73fid = csv.reader(open(file_list))
     74lidar_filenames = []
     75for row in fid:
     76    lidar_filenames.append(','.join(row))
     77lidar_filenames = list(lidar_filenames)   
     78
    8779# Format for point is x,y,elevation (with header)
    88 point_filenames = ['SD100031996_jgriffin_clip.csv',
    89                    'tomaga_offshore_AHD_MGA_1997.csv',
    90                    'Batemans_BBHD_MGA_1995.csv',
    91                    'moruya_AHD_MGA_2000.csv']
     80point_filenames = lidar_filenames    #
    9281         
    9382
     
    10796# Format for points easting,northing (no header)
    10897
    109 interior_regions_data = [['area_of_interest.csv', 500],
    110                         ['area_of_significance.csv', 2500],
    111                         ['shallow_water.csv', 10000]]
     98interior_regions_data =[]#[['area_of_interest.csv', 500],
     99                        #['area_of_significance.csv', 2500],
     100                        #['shallow_water.csv', 10000]]
    112101
    113102
     
    210199# check various directories/files that must exist
    211200anuga_folder = join(home, state, scenario_folder, 'anuga')
    212 topographies_folder = join(anuga_folder, 'topographies')
     201topographies_folder = join(anuga_folder, 'topographies/original')
    213202polygons_folder = join(anuga_folder, 'polygons')
    214203boundaries_folder = join(anuga_folder, 'boundaries')
Note: See TracChangeset for help on using the changeset viewer.