Changeset 4301

Timestamp:

Mar 12, 2007, 5:00:07 PM (18 years ago)

Author:

duncan

Message:

starting to add urs2txt stuff

File:

• anuga_core/source/anuga/shallow_water/data_manager.py (modified) (1 diff)

anuga_core/source/anuga/shallow_water/data_manager.py

@@ -4866 +4866 @@
     ymomentum[slice_index] = va*h
 
+def URS2txt(basename_in):
+    """
+    Not finished or tested
+    """
+    files_in = [basename_in+'-z-mux',
+                basename_in+'-e-mux',
+                basename_in+'-n-mux']
+    quantities = ['HA','UA','VA']
+
+    # instanciate urs_points of the three mux files.
+    mux = {}
+    for quantity, file in map(None, quantities, files_in):
+        mux[quantity] = Urs_points(file)
         
+    # Could check that the depth is the same. (hashing)
+
+    # handle to a mux file to do depth stuff
+    a_mux = mux[quantities[0]]
+    
+    # Convert to utm
+    lat = a_mux.lonlatdep[:,1]
+    long = a_mux.lonlatdep[:,0]
+    points_utm, zone = convert_from_latlon_to_utm( \
+        latitudes=lat, longitudes=long)
+    #print "points_utm", points_utm
+    #print "zone", zone
+    elevations = a_mux.lonlatdep[:,2] * -1 #
+    
+    fid = open(basename_in + '.txt', 'w')
+
+    fid.write("zone" + str(zone) + "\n")
+
+    for elevation, point_utm in map(None, elevations, points_utm):
+        pass
+    
 class Urs_points:
     """