Changeset 5415 for anuga_work/production/busselton/run_busselton.py

Timestamp:

Jun 24, 2008, 10:24:21 AM (15 years ago)

Author:

kristy

Message:

 

File:

• anuga_work/production/busselton/run_busselton.py (modified) (4 diffs)

anuga_work/production/busselton/run_busselton.py

@@ -173 +173 @@
     Br = Reflective_boundary(domain)
     Bd = Dirichlet_boundary([kwargs['tide'],0,0])
-   
+    Bo = Dirichlet_boundary([kwargs['tide']+10.0,0,0]) # To be deleted for FESA runs
 
     Bf = Field_boundary(kwargs['boundary_file'],
@@ -181 +181 @@
     domain.set_boundary({'back': Bd,
                              'side': Bd,
-                             'ocean': Bf})
+                             'ocean': Bd}) #changed from Bf to Bd for large wave
 
     kwargs['input_start_time']=domain.starttime
@@ -192 +192 @@
     t0 = time.time()
 
-    for t in domain.evolve(yieldstep = 240, finaltime = kwargs['finaltime']): 
+    for t in domain.evolve(yieldstep = project.yieldstep, finaltime = kwargs['finaltime']
+                            ,skip_initial_step = False ): 
         domain.write_time()
-        domain.write_boundary_statistics(tags = 'ocean')     
+        domain.write_boundary_statistics(tags = 'ocean')
+
+        # To be deleted for FESA runs
+        if allclose(t, 240):
+            domain.set_boundary({'back': Br, 'side': Bd, 'ocean': Bo})
+
+        if allclose(t, 1440):
+            domain.set_boundary({'back': Br, 'side': Bd, 'ocean': Bd})
+
+
 
     x, y = domain.get_maximum_inundation_location()
@@ -222 +232 @@
     kwargs['starttime']=project.starttime
     kwargs['yieldstep']=project.yieldstep
-    kwargs['midtime']=project.midtime
     kwargs['finaltime']=project.finaltime