Changeset 1125

Timestamp:

Mar 22, 2005, 5:48:05 PM (20 years ago)

Author:

ole

Message:

Working on ferret2sww3 test (Peter)

File:

• inundation/ga/storm_surge/pyvolution/test_data_manager.py (modified) (4 diffs)

inundation/ga/storm_surge/pyvolution/test_data_manager.py

@@ -728 +728 @@
         from coordinate_transforms.redfearn import redfearn
 
-        fid1 = NetCDFFile('small_ha.nc')
-        fid2 = NetCDFFile('small_e.nc')
-        fid3 = NetCDFFile('small_va.nc')
+        fid1 = NetCDFFile('small_ha.nc','a')
+        fid3 = NetCDFFile('small_va.nc','a')
+
+        #Populate with decent values
+        fid1.variables['HA'][:][0,0,3]=130.7
+        fid3.variables['VA'][:][0,0,3]=10.0
+
+
+        fid1.sync(); fid1.close()
+        fid3.sync(); fid3.close()
+
+        fid1 = NetCDFFile('small_ha.nc','r')
+        fid2 = NetCDFFile('small_e.nc','r')
+        fid3 = NetCDFFile('small_va.nc','r')
+
+        print fid3.variables['VA'][:]
+        
 
         first_amp = fid1.variables['HA'][:][0,0,0]
@@ -737 +751 @@
         fourth_elevation= fid2.variables['ELEVATION'][:][0,3]
         first_speed = fid3.variables['VA'][0,0]
-        fourth_speed = fid3.variables['VA'][:][0,3]
+        fourth_speed = fid3.variables['VA'][:][0,0,3]
+        print 'fourth speed', fourth_speed
 
         fid1.close()
@@ -748 +763 @@
 
 
-        #Work out the UTM coordinates for first point
-        zone, e, n = redfearn(-34.5, 150.66667)
-        #print zone, e, n
-
         #Read output file 'small.sww'
         fid = NetCDFFile('small.sww')
 
-        x = fid.variables['x'][:]
-        y = fid.variables['y'][:]
-        
-
-        #Check that first coordinate is correctly represented
-        assert allclose(x[0], e)
-        assert allclose(y[0], n)
 
         #Check first value
@@ -776 +780 @@
 
         assert allclose(ymomentum[0][0],first_momentum)  #Meters
+
+        print 'YM1', ymomentum[0][1]
+        print 'YM4', ymomentum[0][3]        
+        print '4M', fourth_momentum
         assert allclose(ymomentum[0][1],fourth_momentum)  #Meters