Changeset 6958

Timestamp:

May 5, 2009, 11:42:46 AM (16 years ago)

Author:

jgriffin

Message:

added functionality to read list of point filenames

File:

• anuga_work/production/new_south_wales/batemans_bay/project.py (modified) (6 diffs)

anuga_work/production/new_south_wales/batemans_bay/project.py

@@ -8 +8 @@
 from time import localtime, strftime, gmtime
 from os.path import join, exists
+import csv
 
 
@@ -17 +18 @@
 # Note, the user needs to set up the directory system accordingly
 state = 'new_south_wales'
-scenario_name = 'batemans_bay'
-scenario_folder = 'batemans_bay_tsunami_scenario_2009'
+scenario_name = 'gosford'
+scenario_folder = 'gosford_tsunami_scenario_2009'
 
 #-------------------------------------------------------------------------------
 # Initial Conditions
 #-------------------------------------------------------------------------------
+
 
 # Model specific parameters.
 # One or all can be changed each time the run_model script is executed
-tide = 0                # difference between MSL and HAT (1.0)
-
-# the event number or the mux file name
-##event_number = 58129  #1 in 200 yr Puyesgur
-##event_number = 58115    #1 in 500 yr
-##event_number = 58226    #1 in 1000 yr
-##event_number = 58284    #1 in 2000 yr
-##event_number = 58286    #1 in 5000 yr
-##event_number = 58242    #1 in 10000 yr Puysegur
-event_number = 58368    #1 in 10000 yr Puysegur
-##event_number = 51436    #1 in 10000 yr New Hebrides
-
+tide = 0                # difference between MSL and HAT in metres
+zone = 56               # specify zone of model
+event_number = 31891    # the event number or the mux file name
 alpha = 0.1             # smoothing parameter for mesh
 friction=0.01           # manning's friction coefficient
 starttime=0             # start time for simulation
-finaltime=10000         # final time for simulation
-
-setup = 'final'         # This can be one of three values
+finaltime=1000         # final time for simulation
+
+setup = 'trial'         # This can be one of three values
                         #    trial - coarsest mesh, fast
                         #    basic - coarse mesh
@@ -59 +52 @@
 #-------------------------------------------------------------------------------
 
-output_comment = [setup, tide, event_number, 'viz']
+output_comment = [setup, tide, event_number, 'lidar']
 
 #-------------------------------------------------------------------------------
@@ -68 +61 @@
 # Used in build_elevation.py
 # Format for ascii grids, as produced in ArcGIS + a projection file
-ascii_grid_filenames = ['1a',# Topographic data
-                        '1b',
-                        '2b',
-                        '2a_3',
-                        '3b',
-                        '3a',
-                        '4a_2',
-                        '4b',
-                        'off1',
-                        'off2',
-                        'off3',
-                        'bbhd',
-                        'sd100031996_p',
-                        'sd100031996_p2',
-                        'sd100031996_p3',
-                        'sd100031996_p4']
+ascii_grid_filenames = [] 
                     
+
+
+#Get long list of Lidar filenames.
+ENV_INUNDATIONHOME = 'INUNDATIONHOME'
+home = join(os.getenv(ENV_INUNDATIONHOME), 'data')
+anuga_folder = join(home, state, scenario_folder, 'anuga')
+topographies_folder = join(anuga_folder, 'topographies/original')
+file_list = join(topographies_folder, 'file_list.csv')
+fid = csv.reader(open(file_list))
+lidar_filenames = []
+for row in fid:
+    lidar_filenames.append(','.join(row))
+lidar_filenames = list(lidar_filenames)    
+
 # Format for point is x,y,elevation (with header)
-point_filenames = ['SD100031996_jgriffin_clip.csv',
-                   'tomaga_offshore_AHD_MGA_1997.csv',
-                   'Batemans_BBHD_MGA_1995.csv',
-                   'moruya_AHD_MGA_2000.csv']
+point_filenames = lidar_filenames    # 
           
 
@@ -107 +96 @@
 # Format for points easting,northing (no header)
 
-interior_regions_data = [['area_of_interest.csv', 500],
-                        ['area_of_significance.csv', 2500],
-                        ['shallow_water.csv', 10000]]
+interior_regions_data =[]#[['area_of_interest.csv', 500],
+                        #['area_of_significance.csv', 2500],
+                        #['shallow_water.csv', 10000]]
 
 
@@ -210 +199 @@
 # check various directories/files that must exist
 anuga_folder = join(home, state, scenario_folder, 'anuga')
-topographies_folder = join(anuga_folder, 'topographies')
+topographies_folder = join(anuga_folder, 'topographies/original')
 polygons_folder = join(anuga_folder, 'polygons')
 boundaries_folder = join(anuga_folder, 'boundaries')