Changeset 1830

Timestamp:

Sep 13, 2005, 2:07:41 PM (18 years ago)

Author:

ole

Message:

Karratha study

Files:

• inundation/pyvolution/domain.py (modified) (2 diffs)
• production/karratha_2005/create_mesh.py (modified) (1 diff)
• production/karratha_2005/project.py (modified) (1 diff)
• production/karratha_2005/run_karratha.py (modified) (2 diffs)

inundation/pyvolution/domain.py

@@ -362 +362 @@
             
     def write_boundary_statistics(self, quantities = None, tags = None):
-        print self.boundary_statistics(quantitites, tags)
+        print self.boundary_statistics(quantities, tags)
         
     def boundary_statistics(self, quantities = None, tags = None):
@@ -376 +376 @@
         quantities and tags can be either None, a string or a list of strings
 
-        If quantitites are specified only report on those. Otherwise take all conserved quantities.
+        If quantities are specified only report on those. Otherwise take all conserved quantities.
         If tags are specified only report on those, otherwise take all tags.
 

production/karratha_2005/create_mesh.py

@@ -77 +77 @@
     
     dict['segments'] = [[0,1], [1,2], [2,3], [3,0]]
-    dict['segment_tags'] = ['wall', 'wall', 'wall', 'wall']
+    dict['segment_tags'] = ['wall', 'wall', 'ocean', 'wall']
         
     m.addVertsSegs(dict)

production/karratha_2005/project.py

@@ -9 +9 @@
 #Making assumptions abot the location of scenario data
 scenario_dir_name = 'karratha_tsunami_scenario_2005'
-basename = 'karratha_250m'
+basename = 'karratha_100m'
 boundary_basename = 'SU-AU'
 

production/karratha_2005/run_karratha.py

@@ -7 +7 @@
 
 
-from pyvolution.shallow_water import Domain, Reflective_boundary, File_boundary
+from pyvolution.shallow_water import Domain, Reflective_boundary, File_boundary, Dirichlet_boundary, Time_boundary
 from pyvolution.data_manager import convert_dem_from_ascii2netcdf,\
      dem2pts, ferret2sww
@@ -80 +80 @@
 
 #Setup BC
-Bf = File_boundary(project.boundary_basename + '.sww', domain, verbose = True)
+
+print domain.get_boundary_tags()
+
+#Bf = File_boundary(project.boundary_basename + '.sww', domain, verbose = True)
 
 Br = Reflective_boundary(domain)
-domain.set_boundary( {'wall': Br} )
+Bd = Dirichlet_boundary([3,0,0])
+Bw = Time_boundary(domain=domain,
+                   f=lambda t: [(60<t<660)*12, 0, 0])
+domain.set_boundary( {'wall': Br, 'ocean': Bw} )
 
 
 #Run
-for t in domain.evolve(yieldstep = 1, finaltime = 100):
+for t in domain.evolve(yieldstep = 15, finaltime = 4800):
     domain.write_time()
+    domain.write_boundary_statistics(quantities = 'stage')