Changeset 2470

Timestamp:

Mar 3, 2006, 10:03:34 AM (18 years ago)

Author:

nick

Message:

Update source of boundary conditions for Onslow, however still needs more work

Location:

production/onslow_2006

Files:

• project.py (modified) (1 diff)
• run_onslow.py (modified) (5 diffs)

production/onslow_2006/project.py

@@ -40 +40 @@
 outputdir = home+sep+scenario_dir_name+sep+'output'+sep
 polygondir = home+sep+scenario_dir_name+sep+'polygons'+sep
+boundarydir = home+sep+scenario_dir_name+sep+'boundaries'+sep
 
 meshname = meshdir + basename

production/onslow_2006/run_onslow.py

@@ -20 +20 @@
 
 # Related major packages
-from pyvolution.shallow_water import Domain, Reflective_boundary
-from pyvolution.shallow_water import Domain, Time_boundary
+from pyvolution.shallow_water import Domain, Reflective_boundary, \
+                            Dirichlet_boundary, Time_boundary
 from pyvolution.data_manager import convert_dem_from_ascii2netcdf, dem2pts
 from pyvolution.combine_pts import combine_rectangular_points_files 
@@ -132 +132 @@
 #-------------------------------------------------------------------------------
 
-domain.set_quantity('stage', 0.)
+tide = 0.
+
+domain.set_quantity('stage', tide)
 domain.set_quantity('friction', 0.0) 
 domain.set_quantity('elevation', 
@@ -145 +147 @@
 # Setup boundary conditions (all reflective)
 #-------------------------------------------------------------------------------
+
+from pyvolution.data_manager import ferret2sww
+
+south = project.south 
+north = project.north
+west = project.west
+east = project.east
 
 cache(ferret2sww,
@@ -166 +175 @@
 
 Br = Reflective_boundary(domain)
-# 10 min square wave starting at 1 min, 6m high
+Bd = Dirichlet_boundary([tide,0,0])
+
+
+# 7 min square wave starting at 1 min, 6m high
 Bw = Time_boundary(domain = domain,
-                   f=lambda t: [(50<t<800)*6, 0, 0])
+                   f=lambda t: [(60<t<480)*6, 0, 0])
 
 domain.set_boundary( {'top': Bw, 'topleft': Br,
@@ -182 +194 @@
 t0 = time.time()
 
-for t in domain.evolve(yieldstep = 50, finaltime = 5000): 
+for t in domain.evolve(yieldstep = 50, finaltime = 100): 
     domain.write_time()
     domain.write_boundary_statistics(tags = 'top')