Changeset 1113

Timestamp:

Mar 22, 2005, 11:00:33 AM (20 years ago)

Author:

prow

Message:

making ferret2sww read the elevations

File:

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

inundation/ga/storm_surge/pyvolution/data_manager.py

@@ -1230 +1230 @@
     file_u = NetCDFFile(basename_in + '_ua.nc', 'r') #Velocity (x) (cm/s)
     file_v = NetCDFFile(basename_in + '_va.nc', 'r') #Velocity (y) (cm/s)   
+    file_e = NetCDFFile(basename_in + '_va.nc', 'r') #Elevation (z) (m)   
 
     if basename_out is None:
@@ -1279 +1280 @@
     uspeed = file_u.variables['UA'][jmin:jmax, kmin:kmax, lmin:lmax] #Lon
     vspeed = file_v.variables['VA'][jmin:jmax, kmin:kmax, lmin:lmax] #Lat
-
+    elevations = file_e.variables['ELEVATION'][kmin:kmax, lmin:lmax]
     #Get missing values
     nan_ha = file_h.variables['HA'].missing_value[0]
     nan_ua = file_u.variables['UA'].missing_value[0]
     nan_va = file_v.variables['VA'].missing_value[0]
+    nan_e  = file_e.variables['ELEVATION'].missing_value[0]
 
     
@@ -1317 +1319 @@
 
 
+    missing = (elevations == nan_e)
+    if sometrue (missing):
+        if fail_on_NaN:
+            msg = 'NetCDFFile %s contains missing values'\
+                  %(basename_in+'_e.nc')
+            raise msg
+        else:
+            elevations = elevations*(missing==0) + missing*NaN_filler
 
     #######
@@ -1355 +1365 @@
         print '  %s in [%f, %f]' %(name, min(q), max(q))                
 
+        q = elevations.flat
+        name = 'Elevations (e) [m]'                
+        print '  %s in [%f, %f]' %(name, min(q), max(q))                
+
 
     #print number_of_latitudes, number_of_longitudes
@@ -1364 +1378 @@
     file_u.close()
     file_v.close()    
+    file_e.close()    
 
 
@@ -1371 +1386 @@
     #Create new file
     outfile.institution = 'Geoscience Australia'
-    outfile.description = 'Converted from Ferret files: %s, %s, %s'\
+    outfile.description = 'Converted from Ferret files: %s, %s, %s, %s'\
                           %(basename_in + '_ha.nc',
                             basename_in + '_ua.nc',
-                            basename_in + '_va.nc')
+                            basename_in + '_va.nc'
+                            basename_in + '_e.nc')
 
 
@@ -1479 +1495 @@
 
 
-    if elevation is not None:
-        z = elevation
+    if elevations is not None:
+        z = elevations
     else:
         pass