Changeset 5274

Timestamp:

May 5, 2008, 12:56:24 PM (17 years ago)

Author:

ole

Message:

Work in progress regarding flow through a cross section, ticket:175

Location:

anuga_core/source/anuga/shallow_water

Files:

• data_manager.py (modified) (1 diff)
• test_data_manager.py (modified) (1 diff)

anuga_core/source/anuga/shallow_water/data_manager.py

@@ -5634 +5634 @@
 
     # Get mesh and quantities from sww file
-    X = get_mesh_and_quantities_from_sww_file(filename, ['elevation',
-                                                         'stage',
-                                                         'xmomentum',
-                                                         'ymomentum'], verbose=verbose)
+    X = get_mesh_and_quantities_from_sww_file(filename, 
+                                              ['elevation',
+                                               'stage',
+                                               'xmomentum',
+                                               'ymomentum'], 
+                                              verbose=verbose)
     interpolator, quantities, geo_reference, time = X
 

anuga_core/source/anuga/shallow_water/test_data_manager.py

@@ -7419 +7419 @@
 
 
+    def test_get_mesh_and_quantities_from_sww_file(self):
+        """test_get_mesh_and_quantities_from_sww_file(self):
+        """     
+        
+        # Generate a test sww file
+        
+        import time, os
+        from Numeric import array, zeros, allclose, Float, concatenate
+        from Scientific.IO.NetCDF import NetCDFFile
+
+        # Setup
+        from mesh_factory import rectangular
+
+        # Create basic mesh (100m x 5m)
+        width = 5
+        len = 100
+        points, vertices, boundary = rectangular(len, width, 100, 5)
+
+        # Create shallow water domain
+        domain = Domain(points, vertices, boundary)
+
+        domain.set_name('flowtest')
+        swwfile = domain.get_name() + '.sww'
+        domain.set_datadir('.')
+
+        Br = Reflective_boundary(domain)    # Side walls
+        Bd = Dirichlet_boundary([1, 0, 0])  # inflow
+
+        domain.set_boundary( {'left': Bd, 'right': Bd, 'top': Br, 'bottom': Br})
+
+        for t in domain.evolve(yieldstep=1, finaltime = 50):
+            pass
+
+        
+        # Read it
+        
+        # FIXME (Ole): TODO     
+        
+        
     def NOtest_get_flow_through_cross_section(self):
         """test_get_flow_through_cross_section(self):