Changeset 947

Timestamp:

Feb 23, 2005, 9:54:56 AM (20 years ago)

Author:

duncan

Message:

added max and min of stage to extent_sww

Location:

inundation/ga/storm_surge/pyvolution

Files:

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

inundation/ga/storm_surge/pyvolution/data_manager.py

@@ -1375 +1375 @@
     x = fid.variables['x'][:]
     y = fid.variables['y'][:]
-    return [min(x),max(x),min(y),max(y)]
+    stage = fid.variables['stage'][:]
+    return [min(x),max(x),min(y),max(y),min(min(stage)),max(max(stage))]
 
 

inundation/ga/storm_surge/pyvolution/test_data_manager.py

@@ -691 +691 @@
 
         assert allclose(stage[time_index, linear_point_index], test_value/100)
-        
-
+
+        
+    def test_sww_DSG(self):
+        """Not a test, rather a look at the sww format
+        """
+
+        import time, os
+        from Numeric import array, zeros, allclose, Float, concatenate
+        from Scientific.IO.NetCDF import NetCDFFile        
+
+        self.domain.filename = 'datatest' + str(time.time())
+        self.domain.format = 'sww'
+        self.domain.smooth = True
+        self.domain.reduction = mean 
+        self.domain.set_datadir('.')
+      
+        sww = get_dataobject(self.domain)
+        sww.store_connectivity()
+        sww.store_timestep('stage')
+        self.domain.time = 2. 
+        sww.store_timestep('stage')
+        file_and_extension_name = self.domain.filename + ".sww"
+        #print "file_and_extension_name",file_and_extension_name
+        [xmin, xmax, ymin, ymax, stagemin, stagemax] = \
+               extent_sww(file_and_extension_name )
+        #print "[xmin, xmax, ymin, ymax, stagemin, stagemax",[xmin, xmax, ymin, ymax, stagemin, stagemax] 
+        assert allclose(xmin, 0.0) 
+        assert allclose(xmax, 1.0) 
+        assert allclose(ymin, 0.0) 
+        assert allclose(ymax, 1.0) 
+        assert allclose(stagemin, -0.88) 
+        assert allclose(stagemax, 0.12)
         
     def test_ferret2sww_nz_origin(self):