Changeset 6522

Timestamp:

Mar 16, 2009, 1:55:40 PM (15 years ago)

Author:

ole

Message:

Added test that zone, xllcorner and yllcorner are stored correctly in the sww file.

File:

• anuga_core/source/anuga/shallow_water/test_data_manager.py (modified) (6 diffs)

anuga_core/source/anuga/shallow_water/test_data_manager.py

@@ -46 +46 @@
         
         self.verbose = Test_Data_Manager.verbose
-        #Create basic mesh
+        # Create basic mesh
         points, vertices, boundary = rectangular(2, 2)
 
-        #Create shallow water domain
+        # Create shallow water domain
         domain = Domain(points, vertices, boundary)
         domain.default_order = 2
 
-        #Set some field values
+        # Set some field values
         domain.set_quantity('elevation', lambda x,y: -x)
         domain.set_quantity('friction', 0.03)
@@ -1220 +1220 @@
 
     def test_sww2dem_asc_elevation_depth(self):
-        """
-        test_sww2dem_asc_elevation_depth(self):
+        """test_sww2dem_asc_elevation_depth
+        
         Test that sww information can be converted correctly to asc/prj
         format readable by e.g. ArcView
+        
+        Also check geo_reference is correct
         """
 
@@ -1229 +1231 @@
         from Scientific.IO.NetCDF import NetCDFFile
 
-        #Setup
+        # Setup
         self.domain.set_name('datatest')
 
@@ -1242 +1244 @@
         self.domain.set_quantity('stage', 1.0)
 
-        self.domain.geo_reference = Geo_reference(56,308500,6189000)
+        self.domain.geo_reference = Geo_reference(56, 308500, 6189000)
 
         sww = get_dataobject(self.domain)
@@ -1248 +1250 @@
         sww.store_timestep()
 
-        #self.domain.tight_slope_limiters = 1
 
         self.domain.evolve_to_end(finaltime = 0.01)
@@ -1265 +1266 @@
         time = fid.variables['time'][:]
         stage = fid.variables['stage'][:]
+        
+        # Check georeferencig: zone, xllcorner and yllcorner
+        assert fid.zone[0] == 56
+        assert fid.xllcorner[0] == 308500
+        assert fid.yllcorner[0] == 6189000
+                
 
         fid.close()