Changeset 6315

Timestamp:

Feb 12, 2009, 12:34:41 PM (15 years ago)

Author:

myall

Message:

in run_model.py, removed transmissive boundary and used Dirichlet instead

File:

• anuga_work/production/australia_ph2/eucla_motel/run_model.py (modified) (5 diffs)

anuga_work/production/australia_ph2/eucla_motel/run_model.py

@@ -26 +26 @@
 # Related major packages
 from anuga.interface import create_domain_from_regions
-from anuga.interface import Transmissive_stage_zero_momentum_boundary
 from anuga.interface import Dirichlet_boundary
 from anuga.interface import Reflective_boundary
@@ -32 +31 @@
 from anuga.interface import create_sts_boundary
 from anuga.interface import csv2building_polygons
+from file_length import file_length
 
 from anuga.shallow_water.data_manager import start_screen_catcher
@@ -70 +70 @@
 # Number of boundary segments
 N = len(event_sts)-1
+# Number of landward_boundary points
+M = file_length(project.landward_boundary)
 
 # Boundary tags refer to project.landward_boundary
 # 4 points equals 5 segments start at N
-boundary_tags={'back': [N+1, N+2, N+3, N+4, N+5, N+6, N+7],
-               'side': [N, N+8],
+boundary_tags={'back': range(N+1,N+M),
+               'side': [N,N+M],
                'ocean': range(N)}
 
@@ -116 +118 @@
 
 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,
@@ -125 +127 @@
                     verbose=True)
 
+
 domain.set_boundary({'back': Br,
-                     'side': Bt,
+                     'side': Bd,
                      'ocean': Bf})