Changeset 2443

Timestamp:

Feb 23, 2006, 11:02:06 AM (19 years ago)

Author:

nick

Message:

 

Location:

production/onslow_2006

Files:

• project.py (modified) (2 diffs)
• project.pyc (modified) (previous)
• run_onslow.py (modified) (4 diffs)

production/onslow_2006/project.py

@@ -13 +13 @@
 scenario_dir_name = 'onslow_tsunami_scenario_2006'
 
+
+# 250m data to be provided
+coarsename = 'onsl_bathydem250' # get from Neil/Ingo (DEM or topo data)
 """
-# 250m data to be provided
-coarsename = 'onshore_250' # get from Neil/Ingo (DEM or topo data)
 # 30m data to be provided
 finename = 'onshore_30' # get from Neil/Ingo (DEM or topo data) 
@@ -41 +42 @@
 
 meshname = meshdir + basename
+
+coarsedemname = datadir + coarsename
 """
-coarsedemname = datadir + coarsename
 finedemname = datadir + finename
 combineddemname = datadir + 'onslow_combined_elevation'

production/onslow_2006/run_onslow.py

@@ -38 +38 @@
 # Fine pts file to be clipped to area of interest
 #-------------------------------------------------------------------------------
-'''
+
 # filenames
 coarsedemname = project.coarsedemname
+'''
 #finedemname = project.finedemname
 '''
 meshname = project.meshname+'.msh'
-'''
+
 # coarse data
 convert_dem_from_ascii2netcdf(coarsedemname, use_cache=True, verbose=True)
 dem2pts(coarsedemname, use_cache=True, verbose=True)
-
+'''
 # fine data (clipping the points file to smaller area)
 convert_dem_from_ascii2netcdf(finedemname, use_cache=True, verbose=True)
@@ -133 +134 @@
 domain.set_quantity('stage', 0.)
 domain.set_quantity('friction', 0.0) 
-domain.set_quantity('elevation', 0.
+domain.set_quantity('elevation', 
 #                    filename = project.combineddemname + '.pts',
-#                    filename = project.coarsedemname + '.pts',
-#                    use_cache = True,
-#                    verbose = True
+                    filename = project.coarsedemname + '.pts',
+                    use_cache = True,
+                    verbose = True
                     )
 
@@ -150 +151 @@
 # 10 min square wave starting at 1 min, 6m high
 Bw = Time_boundary(domain = domain,
-                   f=lambda t: [(20<t<200)*6, 0, 0])
+                   f=lambda t: [(50<t<800)*6, 0, 0])
 
-domain.set_boundary( {'top': Br, 'topleft': Br,
+domain.set_boundary( {'top': Bw, 'topleft': Br,
                              'left': Br, 'bottom': Br,
                              'bottomright': Br, 'topright': Br} )
@@ -164 +165 @@
 t0 = time.time()
 
-for t in domain.evolve(yieldstep = 10, finaltime = 20): 
+for t in domain.evolve(yieldstep = 50, finaltime = 5000): 
     domain.write_time()
     domain.write_boundary_statistics(tags = 'top')