Changeset 341 in lmdz_wrf

Timestamp:

Mar 6, 2015, 12:54:12 PM (10 years ago)

Author:

lfita

Message:

Too many after 5MetMed fast changes to kept track..

Location:

trunk/tools

Files:

• drawing_tools.py (modified) (2 diffs)
• validation_sim.py (modified) (2 diffs)

trunk/tools/drawing_tools.py

@@ -9 +9 @@
 import os
 from netCDF4 import Dataset as NetCDFFile
+import nc_var_tools as ncvar
 
 errormsg = 'ERROR -- error -- ERROR -- error'
@@ -1490 +1491 @@
 
     potinterval = np.log10(interval)
-
+    print 'LLuis minv, maxv:',minv,maxv,'potinterval:',potinterval
     Ipotint = int(potinterval)
     intvalue = np.float(interval / np.float(Nint))

trunk/tools/validation_sim.py

@@ -707 +707 @@
                 varNOcheckv = 360.*(1. + theta/(2.*np.pi))
 
-                dimensions = ncobj.variables['U'].dimensions
-                shape = ncobj.variables['U'].shape
-
             elif varn == 'WRFz':
                 grav = 9.81
@@ -1272 +1269 @@
 
 # Around values
-    if dims.has_key('Z'):
-        if not onewnc.dimensions.has_key('zaround'):
-            newdim = onewnc.createDimension('zaround',Ngrid*2+1)
-
-        newvar = onewnc.createVariable(valvars[ivar][0] + 'around', 'f',            \
-          ('time','zaround','yaround','xaround'), fill_value=fillValueF)
-    else:
-        newvar = onewnc.createVariable(valvars[ivar][0] + 'around', 'f',            \
-          ('time','yaround','xaround'), fill_value=fillValueF)
-
-    descvar = 'around simulated values +/- grid values: ' + valvars[ivar][0]
-    basicvardef(newvar, varsimobs + 'around', descvar, ovobs.getncattr('units'))
-    newvar[:] = simobsSvalues
+    if not onewnc.variables.has_key(valvars[ivar][0] + 'around'):
+        if dims.has_key('Z'):
+            if not onewnc.dimensions.has_key('zaround'):
+                newdim = onewnc.createDimension('zaround',Ngrid*2+1)
+                newvar = onewnc.createVariable(valvars[ivar][0] + 'around', 'f',     \
+                  ('time','zaround','yaround','xaround'), fill_value=fillValueF)
+        else:
+            newvar = onewnc.createVariable(valvars[ivar][0] + 'around', 'f',         \
+              ('time','yaround','xaround'), fill_value=fillValueF)
+
+        descvar = 'around simulated values +/- grid values: ' + valvars[ivar][0]
+        basicvardef(newvar, varsimobs + 'around', descvar, ovobs.getncattr('units'))
+        newvar[:] = simobsSvalues
 
 # Statistics
-    newvar = onewnc.createVariable(valvars[ivar][0] + 'staround', 'f',               \
-      ('time','stats'), fill_value=fillValueF)
-    descvar = 'around simulated statisitcs: ' + valvars[ivar][0]
-    basicvardef(newvar, varsimobs + 'staround', descvar, ovobs.getncattr('units'))
-    newvar[:] = aroundstats.transpose()
-
-    onewnc.sync()
+    if not searchInlist(onewnc.variables,valvars[ivar][0]):
+        newvar = onewnc.createVariable(valvars[ivar][0] + 'staround', 'f',           \
+          ('time','stats'), fill_value=fillValueF)
+        descvar = 'around simulated statisitcs: ' + valvars[ivar][0]
+        basicvardef(newvar, varsimobs + 'staround', descvar, ovobs.getncattr('units'))
+        newvar[:] = aroundstats.transpose()
+
+        onewnc.sync()
 
 newvar = onewnc.createVariable('simtrj','i',('time','simtrj'))