Changeset 1178 for inundation/ga/storm_surge/pyvolution/data_manager.py

Timestamp:

Mar 31, 2005, 4:51:54 PM (20 years ago)

Author:

duncan

Message:

incorporating geo_ref object into domain class

File:

• inundation/ga/storm_surge/pyvolution/data_manager.py (modified) (8 diffs)

inundation/ga/storm_surge/pyvolution/data_manager.py

@@ -50 +50 @@
 from Numeric import concatenate
 
+from coordinate_transforms.geo_reference import Geo_reference, DEFAULT_ZONE
 
 def make_filename(s):
@@ -248 +249 @@
             #FIXME: Use Georef
             fid.starttime = domain.starttime
-            fid.xllcorner = domain.xllcorner
-            fid.yllcorner = domain.yllcorner
-            fid.zone = domain.zone
 
             # dimension definitions
@@ -267 +265 @@
             fid.createVariable('y', self.precision, ('number_of_points',))
             fid.createVariable('elevation', self.precision, ('number_of_points',))
-
+            if domain.geo_reference is not None:
+                domain.geo_reference.write_NetCDF(fid)
             #FIXME: Backwards compatibility
             fid.createVariable('z', self.precision, ('number_of_points',))
@@ -1077 +1076 @@
     false_northing = 10000000
     NODATA_value = -9999
-    #
-
     
     if quantity is None:
@@ -1110 +1107 @@
     number_of_points = fid.dimensions['number_of_points']
     if origin is None:
-        xllcorner = fid.xllcorner[0]
-        yllcorner = fid.yllcorner[0]
-        zone = fid.zone[0]
+        
+        # get geo_reference
+        #sww files don't have to have a geo_ref
+        try:
+            geo_reference = Geo_reference(NetCDFObject=fid)
+        except AttributeError, e:
+            geo_reference = Geo_reference(DEFAULT_ZONE,0,0)
+        
+        xllcorner = geo_reference.get_xllcorner()
+        yllcorner = geo_reference.get_yllcorner()
+        zone = geo_reference.get_zone()
     else:
         zone = origin[0]
@@ -1872 +1877 @@
     ymomentum = fid.variables['ymomentum']   #Momentum in the y-direction
 
-    xllcorner = fid.xllcorner[0]
-    yllcorner = fid.yllcorner[0]
     starttime = fid.starttime[0]
-    zone = fid.zone
     volumes = fid.variables['volumes'][:] #Connectivity
     coordinates=transpose(asarray([x.tolist(),y.tolist()]))
@@ -1883 +1885 @@
     other_quantities = []
 
+    # get geo_reference
+    #sww files don't have to have a geo_ref
+    try:
+        geo_reference = Geo_reference(NetCDFObject=fid)
+    except AttributeError, e:
+        geo_reference = None
+        
     if verbose: print '    interpolating quantities'
     for quantity in fid.variables.keys():
@@ -1916 +1925 @@
     if not boundary is None:
         domain.boundary = boundary
-    domain.zone=zone
-    domain.xllcorner=float(xllcorner)
-    domain.yllcorner=float(yllcorner)
+    
+    domain.geo_reference = geo_reference
 
     domain.starttime=float(starttime)+float(t)