Changeset 6460 for anuga_work/production/australia_ph2/broome/run_model.py

Timestamp:

Mar 5, 2009, 10:20:00 AM (14 years ago)

Author:

kristy

Message:

Changed all files for new wave input

File:

• anuga_work/production/australia_ph2/broome/run_model.py (modified) (3 diffs)

anuga_work/production/australia_ph2/broome/run_model.py

@@ -35 +35 @@
 from anuga.interface import Reflective_boundary
 from anuga.interface import Field_boundary
+from anuga.interface import Time_boundary
+from anuga.interface import file_function
+
 from anuga.interface import create_sts_boundary
 from anuga.interface import csv2building_polygons
@@ -47 +50 @@
 from setup_model import project
 import build_urs_boundary as bub
+import prepare_timeboundary as TB
 
 #-------------------------------------------------------------------------------
@@ -135 +139 @@
 print 'Set boundary - available tags:', domain.get_boundary_tags()
 
+# Prepare time boundary
+TB.prepare_timeboundary(project.boundary_csv)
+f = file_function(project.boundary_csv[:-4] + '.tms')
+
 Br = Reflective_boundary(domain)
 Bt = Transmissive_stage_zero_momentum_boundary(domain)
 Bd = Dirichlet_boundary([project.tide, 0, 0])
-Bf = Field_boundary(project.event_sts+'.sts',
-                    domain, mean_stage=project.tide,
-                    time_thinning=1,
-                    default_boundary=Bd,
-                    boundary_polygon=bounding_polygon_sts,                    
-                    use_cache=True,
-                    verbose=True)
 
-domain.set_boundary({'back': Br,
+if project.wave == 'Bf':
+    Bf = Field_boundary(project.event_sts+'.sts',
+                        domain, mean_stage=project.tide,
+                        time_thinning=1,
+                        default_boundary=Bd,
+                        boundary_polygon=bounding_polygon_sts,                    
+                        use_cache=True,
+                        verbose=True)
+    domain.set_boundary({'back': Br,
                      'side': Bd,
                      'ocean': Bf}) 
+    
+elif project.wave == 'Tb':
+    Tb = Time_boundary(domain,f,default_boundary=Bd )
+
+    domain.set_boundary({'back': Br,
+                         'side': Bd,
+                         'ocean': Tb})
+else:
+    print 'No wave specified in project script (Bf or Tb)'
+    
 
 #-------------------------------------------------------------------------------