Changeset 845

Timestamp:

Feb 8, 2005, 1:00:25 PM (20 years ago)

Author:

ole

Message:

Started working on spatio-temporal boundary

Location:

inundation/ga/storm_surge/pyvolution

Files:

• shallow_water.py (modified) (2 diffs)
• test_shallow_water.py (modified) (1 diff)

inundation/ga/storm_surge/pyvolution/shallow_water.py

@@ -716 +716 @@
 
 ###############################################
-#Boundary - specific to the shallow water wave equation
+#Boundaries - specific to the shallow water wave equation
 class Reflective_boundary(Boundary):
     """Reflective boundary returns same conserved quantities as
@@ -765 +765 @@
         return q
 
+        
+class Spatio_temporal_boundary(Boundary):
+    """The spatio-temporal boundary, reads values for the conserved 
+    quantities from an sww NetCDF file, and returns interpolated values 
+    at the midpoints of each associated boundaty segment.
+    Time dependency is interpolated linearly as in util.File_function.
+      
+    """
+
+    pass
+        
 
 #########################

inundation/ga/storm_surge/pyvolution/test_shallow_water.py

@@ -2372 +2372 @@
 
 
-        
+
+    def test_spatio_temporal_boundary(self):
+        """Test that boundary values can be read from file and interpolated
+        in both time and space
+        """
+        import time
+       
+        #Create sww file of simple propagation from left to right 
+        #through rectangular domain: domain1
+        
+        from mesh_factory import rectangular
+
+        #Create basic mesh
+        points, vertices, boundary = rectangular(3, 3)
+
+        #Create shallow water domain
+        domain = Domain(points, vertices, boundary)
+        domain.smooth = False
+
+        domain.visualise = False
+        domain.default_order = 2
+
+        domain.set_datadir('.')
+        domain.store = True     
+        domain.set_name('test_spatio_temporal_boundary_' + str(time.time()))
+        
+        #Bed-slope and friction at vertices (and interpolated elsewhere)
+        domain.set_quantity('elevation', 0)
+        domain.set_quantity('friction', 0)      
+
+        # Boundary conditions
+        Br = Reflective_boundary(domain)
+        Bd = Dirichlet_boundary([0.3,0,0])      
+        Bt = Transmissive_boundary(domain)
+        domain.set_boundary({'left': Bd, 'right': Bt, 'top': Br, 'bottom': Br})
+
+        #Initial condition
+        domain.set_quantity('stage', 0)
+        domain.check_integrity()
+        
+        #Evolution
+        for t in domain.evolve(yieldstep = 0.05, finaltime = 4):
+            import time
+            #time.sleep(0.01)
+            #domain.write_time()            
+
+
+        
+        
+        
+        #Create an triangle shaped domain, formed from the 
+        #lower and right hand  boundaries and the sw-ne diagonal 
+        #from domain 1. Called it domain2
+        
+        #Use output from domain1 as spatio-temporal boundary for domain2
+        #and verify that results at right hand side are close.  
+
+        
+        
     
 #-------------------------------------------------------------