Changeset 941

Timestamp:

Feb 22, 2005, 3:23:43 PM (20 years ago)

Author:

duncan

Message:

find the extent of an sww

File:

• inundation/ga/storm_surge/pyvolution/data_manager.py (modified) (1 diff)

inundation/ga/storm_surge/pyvolution/data_manager.py

@@ -1349 +1349 @@
     outfile.close()                
     
-    
-
+ 
+def extent_sww(file_name):
+    """
+    Read in an sww file.
+    
+    Input;
+    file_name - the sww file
+    
+    Output;
+    z - Vector of bed elevation
+    volumes - Array.  Each row has 3 values, representing
+    the vertices that define the volume
+    time - Vector of the times where there is stage information
+    stage - array with respect to time and vertices (x,y)
+    """
+        
+    
+    from Scientific.IO.NetCDF import NetCDFFile         
+    
+    #Check contents
+    #Get NetCDF
+    fid = NetCDFFile(file_name, 'r') 
+    
+    # Get the variables
+    x = fid.variables['x'][:]
+    y = fid.variables['y'][:]
+    return [min(x),max(x),min(y),max(y)]