Changeset 1992 for inundation/pyvolution/test_data_manager.py

Timestamp:

Nov 2, 2005, 10:20:13 AM (18 years ago)

Author:

duncan

Message:

adding code to convert gippsland to sww

File:

• inundation/pyvolution/test_data_manager.py (modified) (2 diffs)

inundation/pyvolution/test_data_manager.py

@@ -6 +6 @@
 from Numeric import zeros, array, allclose, Float
 from util import mean
+import tempfile
 
 from data_manager import *
@@ -2653 +2654 @@
                 verbose = True)
 
-
-
+    def test_read_asc(self):
+        """Test conversion from dem in ascii format to native NetCDF xya format
+        """
+
+        import time, os
+        from Numeric import array, zeros, allclose, Float, concatenate
+        from Scientific.IO.NetCDF import NetCDFFile
+
+        import data_manager
+        #Write test asc file
+        filename = tempfile.mktemp(".000")
+        fid = open(filename, 'w')
+        fid.write("""ncols         7
+nrows         4
+xllcorner     2000.5
+yllcorner     3000.5
+cellsize      25
+NODATA_value  -9999
+    97.921    99.285   125.588   180.830   258.645   342.872   415.836
+   473.157   514.391   553.893   607.120   678.125   777.283   883.038
+   984.494  1040.349  1008.161   900.738   730.882   581.430   514.980
+   502.645   516.230   504.739   450.604   388.500   338.097   514.980
+""")
+        fid.close()
+        xllcorner, yllcorner, cellsize, NODATA_value, grid = \
+                   data_manager._read_asc(filename, verbose=False)
+        self.failUnless(xllcorner  == 2000.5,  'Failed')
+        self.failUnless(yllcorner  == 3000.5,  'Failed')
+        self.failUnless(cellsize  == 25,  'Failed')
+        self.failUnless(NODATA_value  == -9999,  'Failed')
+        self.failUnless(grid[0][0]  == 97.921,  'Failed')
+        self.failUnless(grid[3][6]  == 514.980,  'Failed')
+        
 
 #-------------------------------------------------------------