Changeset 6427

Timestamp:

Feb 27, 2009, 9:02:48 AM (16 years ago)

Author:

jgriffin

Message:

updated batemans bay scripts

Location:

anuga_work/production/new_south_wales/batemans_bay

Files:

• build_elevation.py (modified) (3 diffs)
• export_results_max.py (modified) (2 diffs)
• project.py (modified) (5 diffs)
• run_model.py (modified) (1 diff)

anuga_work/production/new_south_wales/batemans_bay/build_elevation.py

@@ -59 +59 @@
         dem2pts(absolute_filename, use_cache=True, verbose=True)
 
-        geospatial_data[filename] = Geospatial_data(file_name=absolute_filename+'.pts',
+        G_grid = Geospatial_data(file_name=absolute_filename+'.pts',
                                                     verbose=True)
+        print 'Clip geospatial object'
+        geospatial_data[filename] = G_grid.clip(project.bounding_polygon)
 
 # Create Geospatial data from TXT files
@@ -66 +68 @@
     for filename in project.point_filenames:
         absolute_filename = join(project.topographies_folder, filename)
-        geospatial_data[filename] = Geospatial_data(file_name=absolute_filename,
+        G_points = Geospatial_data(file_name=absolute_filename,
                                                     verbose=True)
+        print 'Clip geospatial object'
+        geospatial_data[filename] = G_points.clip(project.bounding_polygon)
 
 #-------------------------------------------------------------------------------
@@ -82 +86 @@
     G += geospatial_data[key]
         
-print 'Clip combined geospatial data'
+#print 'Clip combined geospatial data'
 ##G_clip = G.clip_outside(project.poly_aoi1)
 ##G_all = G_clip + geospatial_data[project.point_filenames[4]] 
 #G_clipped = G_all.clip(project.poly_all)
-G_clip = G.clip(project.bounding_polygon)
+#G_clip = G.clip(project.bounding_polygon)
 
 
 print 'Export combined DEM file'
-G_clip.export_points_file(project.combined_elevation + '.pts')
+G.export_points_file(project.combined_elevation + '.pts')
 print 'Do txt version too'
 # Use for comparision in ARC
-G_clip.export_points_file(project.combined_elevation + '.txt')
+G.export_points_file(project.combined_elevation + '.txt')
 

anuga_work/production/new_south_wales/batemans_bay/export_results_max.py

@@ -29 +29 @@
 ##time_dirs = [time_dir1, time_dir2]#, time_dir3]
 
-time_dir1 = '20090130_165532_run_final_0.6_27255_extend_dt_kvanputt'
+time_dir1 = '20090225_153333_run_final_0_31821_jgriffin'
 ##time_dir2 = '20081211_162311_run_final_0_27255_alpha0.1_kvanputt'
 ##time_dir3 = '20081211_162346_run_final_0_68693_alpha0.1_kvanputt'
@@ -39 +39 @@
   
 
-cellsize = 250
+cellsize = 20
 ##cellsize = 5
 

anuga_work/production/new_south_wales/batemans_bay/project.py

@@ -16 +16 @@
 # this section needs to be updated to reflect the modelled community.
 # Note, the user needs to set up the directory system accordingly
-state = 'western_australia'
-scenario_name = 'busselton'
-scenario_folder = 'busselton_tsunami_scenario'
+state = 'new_south_wales'
+scenario_name = 'batemans_bay'
+scenario_folder = 'batemans_bay_tsunami_scenario_2009'
 
 #-------------------------------------------------------------------------------
@@ -27 +27 @@
 # One or all can be changed each time the run_model script is executed
 tide = 0                # difference between MSL and HAT
-event_number = 27197    # the event number or the mux file name
+event_number = 31821    # 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=80000         # final time for simulation
+finaltime=1000         # final time for simulation
 
 setup = 'final'         # This can be one of three values
@@ -58 +58 @@
 # Used in build_elevation.py
 # Format for ascii grids, as produced in ArcGIS + a projection file
-ascii_grid_filenames = ['busselton_v2',   # Topo
-                        'grid_250m_clip'] # Busselton Topo
-
+ascii_grid_filenames = ['1a',# Topographic data
+                        '1b',
+                        '2b',
+                        '2a_3',
+                        '3b',
+                        '3a',
+                        '4a_2',
+                        '4b',]
+                    
 # Format for point is x,y,elevation (with header)
-point_filenames = ['Busselton_Contour0.txt',    # Coastline
-                   'Busselton_BeachSurvey.txt', # Beach survey
-                   'Busselton_NavyFinal.txt',   # Bathymetry
-                   'Busselton_Chart.txt',       # Bathymetry Charts
-                   'Busselton_Digitised.txt',   # Digitised Fairsheet
-                   'Busselton_250m.txt',        # 250m
-                   'Bunbury_TIN.txt',           # Bunbury aoi TIN'd in ArcGIS
-                   'Busselton_TIN.txt',         # Busselton aoi TIN'd in ArcGIS
-                   'XYAHD_clip.txt']            # To extend boundary
+point_filenames = ['Batemans_BBHD_MGA_1995.csv',
+                   'SD100015374_MGA.csv',
+                   'SD100016502_MGA.csv',
+                   'SD100016517_MGA.csv',
+                   'SD100031996_batemans_mga_clip_written.csv',
+                   'tomaga_offshore_AHD_MGA_1997.csv']
+          
 
 ### Add csv header list to all files in point_filenames
@@ -85 +89 @@
 # INTERIOR REGIONS -  for designing the mesh
 # Used in run_model.py
-# Format for points easting,northing (no header)                    
-interior_regions_data = [['busselton_1km.csv', 500],
-                         ['bunbury_1km.csv', 500],
-                         ['busselton_2km.csv', 10000],
-                         ['bunbury_2km.csv', 10000],
-                         ['island1.csv', 10000],
-                         ['island2.csv', 10000],
-                         ['coast_5km_d20m.csv', 40000]]
+# Format for points easting,northing (no header)
+
+interior_regions_data = [['area_of_interest.csv', 500],
+                        ['area_of_significance.csv', 2500],
+                        ['shallow_water.csv', 10000]]
+
 
 # LAND - used to set the initial stage/water to be offcoast only
 # Used in run_model.py.  Format for points easting,northing (no header)
-land_initial_conditions_filename = [['initial_condition_extend.csv', 0],
-                                    ['initial_condition_marina.csv', 0]]
+land_initial_conditions_filename = [['initial_condition_extend.csv', 0]]
+                                  
 
 # GAUGES - for creating timeseries at a specific point
 # Used in get_timeseries.py.  
 # Format easting,northing,name,elevation (with header)
-gauges_filename = 'gauges.csv'
+##gauges_filename = 'gauges.csv'
 
 # BUILDINGS EXPOSURE - for identifying inundated houses
 # Used in run_building_inundation.py
 # Format latitude,longitude etc (geographic)
-building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS
+##building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS
 
 # BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively
@@ -237 +239 @@
 # The absolute pathname for the gauges file
 # Used for get_timeseries.py
-gauges = join(gauges_folder, gauges_filename)       
+##gauges = join(gauges_folder, gauges_filename)       
 
 # The absolute pathname for the building file
 # Used for run_building_inundation.py
-building_exposure = join(gauges_folder, building_exposure_filename)
+##building_exposure = join(gauges_folder, building_exposure_filename)
 
 # full path to where MUX files (or meta-files) live

anuga_work/production/new_south_wales/batemans_bay/run_model.py

@@ -100 +100 @@
                                     interior_regions=project.interior_regions,
                                     mesh_filename=project.meshes,
-                                    use_cache=True,
+                                    use_cache=False,
                                     verbose=True)
 print domain.statistics()