Changeset 3190 for documentation/user_manual/examples/runsydney.py

Timestamp:

Jun 21, 2006, 9:31:57 AM (19 years ago)

Author:

sexton

Message:

MOST and ANUGA comparisons and updates

File:

• documentation/user_manual/examples/runsydney.py (modified) (4 diffs)

documentation/user_manual/examples/runsydney.py

@@ -57 +57 @@
 interior_regions = [[project.northern_polygon, interior_res],
                     [project.harbour_polygon, interior_res],
-                    [project.southern_polygon, interior_res]]
+                    [project.southern_polygon, interior_res],
+                    [project.manly_polygon, interior_res],
+                    [project.top_polygon, interior_res]]
 
 
 create_mesh_from_regions(project.demopoly,
                          boundary_tags= {'oceannorth': [0], 
-                                         'onshorewest': [1],
-                                         'harbournorth': [2],
-                                         'harbourwest': [3],
-                                         'harboursouth': [4],
-                                         'oceansouth': [5],
-                                         'oceaneast': [6]},
+                                         'onshorenorth1': [1],
+                                         'onshorenorthwest1': [2],
+                                         'onshorenorthwest2': [3],
+                                         'onshorenorth2': [4],
+                                         'onshorewest': [5],
+                                         'onshoresouth': [6],
+                                         'oceansouth': [7],
+                                         'oceaneast': [8]},
                          maximum_triangle_area=100000,
                          filename=meshname,           
@@ -115 +119 @@
                     verbose = True)
 
-
 #------------------------------------------------------------------------------
 # Setup boundary conditions (all Dirichlet)
@@ -123 +126 @@
 
 Bd = Dirichlet_boundary([0.0,0.0,0.0])
-domain.set_boundary( {'oceannorth': Bd, 'onshorewest': Bd,
-                      'harbournorth': Bd, 'harbourwest': Bd,
-                      'harboursouth': Bd, 'oceansouth': Bd,
-                      'oceaneast': Bd})
-
+domain.set_boundary( {'oceannorth': Bd,
+                      'onshorenorth1': Bd,
+                      'onshorenorthwest1': Bd,
+                      'onshorenorthwest2': Bd,
+                      'onshorenorth2': Bd,
+                      'onshorewest': Bd,
+                      'onshoresouth': Bd,
+                      'oceansouth': Bd,
+                      'oceaneast': Bd} )
 
 #------------------------------------------------------------------------------
@@ -135 +142 @@
 import time
 t0 = time.time()
-tend = 7200.0
-print 'Simulate to time = %d secs' %tend
 
-for t in domain.evolve(yieldstep = 60, duration = tend):
+for t in domain.evolve(yieldstep = 60, finaltime = 7200):
     print domain.timestepping_statistics()
     print domain.boundary_statistics(tags = 'oceaneast')