Changeset 623 for inundation/ga/storm_surge/pyvolution/test_generic_boundary_conditions.py

Timestamp:

Nov 24, 2004, 4:43:40 PM (20 years ago)

Author:

ole

Message:

Created and tested general File_function (reading and interpolating time series from file) and used it to refactor and simplify File_boundary

File:

• inundation/ga/storm_surge/pyvolution/test_generic_boundary_conditions.py (modified) (2 diffs)

inundation/ga/storm_surge/pyvolution/test_generic_boundary_conditions.py

@@ -144 +144 @@
             t_string = time.strftime(time_format, time.gmtime(t))
 
-            
             fid.write('%s,%f %f\n' %(t_string, 1.0*i, sin(i*2*pi/10))) 
         fid.close()
                       
         
-        F = File_boundary(domain, filename)
+        F = File_boundary(filename, domain)
 
         from config import default_boundary_tag
@@ -164 +163 @@
 
         os.remove(filename)
+
+
+
+    def test_fileboundary_exception(self):
+        """Test that boundary object complians if number of
+        conserved quantities are wrong
+        """
+
+        from domain import Domain
+        from quantity import Quantity
+        import time, os
+        from math import sin, pi
+        from config import time_format
+        
+        a = [0.0, 0.0]
+        b = [0.0, 2.0]
+        c = [2.0,0.0]
+        d = [0.0, 4.0]
+        e = [2.0, 2.0]
+        f = [4.0,0.0]
+
+        points = [a, b, c, d, e, f]
+
+        #bac, bce, ecf, dbe
+        elements = [ [1,0,2], [1,2,4], [4,2,5], [3,1,4] ]
+        
+        domain = Domain(points, elements)
+        domain.conserved_quantities = ['level', 'xmomentum', 'ymomentum']
+        domain.quantities['level'] =\
+                                   Quantity(domain, [[1,2,3], [5,5,5],
+                                                     [0,0,9], [-6, 3, 3]])
+
+        domain.quantities['xmomentum'] =\
+                                   Quantity(domain, [[2,3,4], [5,5,5],
+                                                     [0,0,9], [-6, 3, 3]])
+        domain.quantities['ymomentum'] =\
+                                   Quantity(domain, [[2,3,4], [5,5,5],
+                                                     [0,0,9], [-6, 3, 3]])
+
+        domain.check_integrity()
+
+        #Write file (with only two values)
+        filename = 'boundarytest' + str(time.time())
+        fid = open(filename, 'w')
+        start = time.mktime(time.strptime('2000', '%Y'))
+        dt = 5*60  #Five minute intervals
+        for i in range(10):
+            t = start + i*dt
+            t_string = time.strftime(time_format, time.gmtime(t))
+
+            fid.write('%s,%f %f\n' %(t_string, 1.0*i, sin(i*2*pi/10))) 
+        fid.close()
+                      
+
+        try:
+            F = File_boundary(filename, domain)
+        except AssertionError:
+            pass
+        else:
+            raise 'Should have raised an exception'
+
+        os.remove(filename)
+
+
         
 #-------------------------------------------------------------