Changeset 6674

Timestamp:

Apr 1, 2009, 10:17:59 AM (16 years ago)

Author:

myall

Message:

fixing Tb boundary - default was still Bd (the Dirichlet with tide added), but should have been Dirichlet with tide=0. this isn't a problem for phase 2 work as project.tide is 0 anyway, but just thought I'd fix it in case of future changes.

Location:

anuga_work/production/australia_ph2

Files:

• adelaide/run_model.py (modified) (1 diff)
• albany/run_model.py (modified) (1 diff)
• broome/run_model.py (modified) (2 diffs)
• carnarvon/run_model.py (modified) (1 diff)
• ceduna/run_model.py (modified) (1 diff)
• darwin/run_model.py (modified) (1 diff)
• esperance/run_model.py (modified) (1 diff)
• eucla_motel/run_model.py (modified) (1 diff)
• hobart/run_model.py (modified) (1 diff)
• perth/run_model.py (modified) (1 diff)
• strahan/run_model.py (modified) (1 diff)
• wyndham/run_model.py (modified) (1 diff)

anuga_work/production/australia_ph2/adelaide/run_model.py

@@ -160 +160 @@
     
 elif project.wave == 'Tb':
-    Tb = Time_boundary(domain,f,default_boundary=Bd )
+    Tb = Time_boundary(domain,f,default_boundary=Dirichlet_boundary([0, 0, 0]) )
 
     domain.set_boundary({'back': Br,

anuga_work/production/australia_ph2/albany/run_model.py

@@ -160 +160 @@
     
 elif project.wave == 'Tb':
-    Tb = Time_boundary(domain,f,default_boundary=Bd )
+    Tb = Time_boundary(domain,f,default_boundary=Dirichlet_boundary([0, 0, 0]) )
 
     domain.set_boundary({'back': Br,

anuga_work/production/australia_ph2/broome/run_model.py

@@ -1 @@
-"""Run a tsunami inundation scenario for Busselton, WA, Australia.
+"""Run a tsunami inundation scenario for phase 2 model.
 
 The scenario is defined by a triangular mesh created from project.polygon, the
@@ -156 +156 @@
                         verbose=True)
     domain.set_boundary({'back': Br,
-                     'side': Bd,
+                     'side': Bt,
                      'ocean': Bf}) 
     
 elif project.wave == 'Tb':
-    Tb = Time_boundary(domain,f,default_boundary=Bd )
+    Tb = Time_boundary(domain,f,default_boundary=Dirichlet_boundary([0, 0, 0]))
 
     domain.set_boundary({'back': Br,
-                         'side': Bd,
+                         'side': Bt,
                          'ocean': Tb})
 else:

anuga_work/production/australia_ph2/carnarvon/run_model.py

@@ -160 +160 @@
     
 elif project.wave == 'Tb':
-    Tb = Time_boundary(domain,f,default_boundary=Bd )
+    Tb = Time_boundary(domain,f,default_boundary=Dirichlet_boundary([0, 0, 0]) )
 
     domain.set_boundary({'back': Br,

anuga_work/production/australia_ph2/ceduna/run_model.py

@@ -160 +160 @@
     
 elif project.wave == 'Tb':
-    Tb = Time_boundary(domain,f,default_boundary=Bd )
+    Tb = Time_boundary(domain,f,default_boundary=Dirichlet_boundary([0, 0, 0]) )
 
     domain.set_boundary({'back': Br,

anuga_work/production/australia_ph2/darwin/run_model.py

@@ -161 +161 @@
     
 elif project.wave == 'Tb':
-    Tb = Time_boundary(domain,f,default_boundary=Bd )
+    Tb = Time_boundary(domain,f,default_boundary=Dirichlet_boundary([0, 0, 0]) )
 
     domain.set_boundary({'back': Br,

anuga_work/production/australia_ph2/esperance/run_model.py

@@ -160 +160 @@
     
 elif project.wave == 'Tb':
-    Tb = Time_boundary(domain,f,default_boundary=Bd )
+    Tb = Time_boundary(domain,f,default_boundary=Dirichlet_boundary([0, 0, 0]) )
 
     domain.set_boundary({'back': Br,

anuga_work/production/australia_ph2/eucla_motel/run_model.py

@@ -160 +160 @@
     
 elif project.wave == 'Tb':
-    Tb = Time_boundary(domain,f,default_boundary=Bd )
+    Tb = Time_boundary(domain,f,default_boundary=Dirichlet_boundary([0, 0, 0]) )
 
     domain.set_boundary({'back': Br,

anuga_work/production/australia_ph2/hobart/run_model.py

@@ -160 +160 @@
     
 elif project.wave == 'Tb':
-    Tb = Time_boundary(domain,f,default_boundary=Bd )
+    Tb = Time_boundary(domain,f,default_boundary=Dirichlet_boundary([0, 0, 0]) )
 
     domain.set_boundary({'back': Br,

anuga_work/production/australia_ph2/perth/run_model.py

@@ -160 +160 @@
     
 elif project.wave == 'Tb':
-    Tb = Time_boundary(domain,f,default_boundary=Bd )
+    Tb = Time_boundary(domain,f,default_boundary=Dirichlet_boundary([0, 0, 0]) )
 
     domain.set_boundary({'back': Br,

anuga_work/production/australia_ph2/strahan/run_model.py

@@ -160 +160 @@
     
 elif project.wave == 'Tb':
-    Tb = Time_boundary(domain,f,default_boundary=Bd )
+    Tb = Time_boundary(domain,f,default_boundary=Dirichlet_boundary([0, 0, 0]) )
 
     domain.set_boundary({'back': Br,

anuga_work/production/australia_ph2/wyndham/run_model.py

@@ -160 +160 @@
     
 elif project.wave == 'Tb':
-    Tb = Time_boundary(domain,f,default_boundary=Bd )
+    Tb = Time_boundary(domain,f,default_boundary=Dirichlet_boundary([0, 0, 0]) )
 
     domain.set_boundary({'back': Br,