Changeset 6597

Timestamp:

Mar 24, 2009, 2:15:16 PM (16 years ago)

Author:

jgriffin

Message:

 

Location:

anuga_work/production/new_south_wales/batemans_bay

Files:

• get_timeseries.py (added)
• project.py (modified) (5 diffs)
• run_model.py (modified) (1 diff)

anuga_work/production/new_south_wales/batemans_bay/project.py

@@ -17 +17 @@
 # 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'
 
 #-------------------------------------------------------------------------------
@@ -26 +26 @@
 # Model specific parameters.
 # One or all can be changed each time the run_model script is executed
-tide = 0                # difference between MSL and HAT
-event_number = 58368    # the event number or the mux file name
+tide = 0            # difference between MSL and HAT
+
+# the event number or the mux file name
+##event_number =     #1 in 200 yr Puyesgur
+##event_number =     #1 in 500 yr
+##event_number =     #1 in 1000 yr
+##event_number =     #1 in 2000 yr
+##event_number =     #1 in 5000 yr
+event_number = 58242    #1 in 10000 yr Puysegur
+event_number = 51436    #1 in 10000 yr New Hebrides
+
 alpha = 0.1             # smoothing parameter for mesh
 friction=0.01           # manning's friction coefficient
 starttime=0             # start time for simulation
-finaltime=6000         # final time for simulation
+finaltime=60000         # final time for simulation
 
 setup = 'final'         # This can be one of three values
@@ -58 +67 @@
 # 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 = []
                     
 # 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 = []
           
 
@@ -110 +101 @@
 # 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
@@ -245 +236 @@
 # 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

anuga_work/production/new_south_wales/batemans_bay/run_model.py

@@ -161 +161 @@
     print domain.boundary_statistics(tags='ocean')
 
-##for t in domain.evolve(yieldstep=5, 
-##                       finaltime=project.finaltime,
-##                       skip_initial_step=True): 
-##    print domain.timestepping_statistics()
-##    print domain.boundary_statistics(tags='ocean')
+for t in domain.evolve(yieldstep=20, 
+                       finaltime=6000,
+                       skip_initial_step=True): 
+    print domain.timestepping_statistics()
+    print domain.boundary_statistics(tags='ocean')
+
+
+for t in domain.evolve(yieldstep=project.yieldstep, 
+                       finaltime=project.finaltime,
+                       skip_initial_step=True): 
+    print domain.timestepping_statistics()
+    print domain.boundary_statistics(tags='ocean')
 
 print 'Simulation took %.2f seconds' % (time.time()-t0)