Changeset 500 in lmdz_wrf for trunk/tools/validation_sim.py

Timestamp:

Jun 15, 2015, 4:48:34 PM (9 years ago)

Author:

lfita

Message:

Working version with the function and the addition of the station description

File:

• trunk/tools/validation_sim.py (modified) (9 diffs)

trunk/tools/validation_sim.py

@@ -992 +992 @@
     fname = 'adding_station_desc'
 
-    newvar = objfile.createVariable( 'station', 'c', ('StrLength'))
+    newvar = onc.createVariable( 'station', 'c', ('StrLength'))
     newvar[0:len(stdesc[0])] = stdesc[0]
 
@@ -1005 +1005 @@
         lonname = 'lon'
 
-    newvar = objfile.createVariable( lonname, 'f4', ('nst'))
+    newvar = onc.createVariable( lonname, 'f4', ('nst'))
     basicvardef(newvar, lonname, 'longitude', 'degrees_West' )
     newvar[:] = stdesc[1]
@@ -1017 +1017 @@
         latname = 'lat'
 
-    newvar = objfile.createVariable( latname, 'f4', ('nst'))
+    newvar = onc.createVariable( latname, 'f4', ('nst'))
     basicvardef(newvar, lonname, 'latitude', 'degrees_North' )
     newvar[:] = stdesc[2]
@@ -1029 +1029 @@
         heightname = 'height'
 
-    newvar = objfile.createVariable( heightname, 'f4', ('nst'))
+    newvar = onc.createVariable( heightname, 'f4', ('nst'))
     basicvardef(newvar, heightname, 'height above sea level', 'm' )
     newvar[:] = stdesc[3]
@@ -1293 +1293 @@
             print '    [DIM]@[dimnsim]@[dimnobs]'
             quit(-1)
-        dims[dsecs[0]] = [dsecs[1], dsecs[2]]
-        simdims[dsecs[0]] = dsecs[1]
-        obsdims[dsecs[0]] = dsecs[2]
-
-        print '  ',dsecs[0],':',dsecs[1],',',dsecs[2]
+        if dsecs[1] != 'None':
+            dims[dsecs[0]] = [dsecs[1], dsecs[2]]
+            simdims[dsecs[0]] = dsecs[1]
+            obsdims[dsecs[0]] = dsecs[2]
+
+            print '  ',dsecs[0],':',dsecs[1],',',dsecs[2]
         
 if opts.vardims is None:
@@ -1316 +1317 @@
             print '    [DIM]@[vardimnsim]@[vardimnobs]'
             quit(-1)
-        vardims[dsecs[0]] = [dsecs[1], dsecs[2]]
-        print '  ',dsecs[0],':',dsecs[1],',',dsecs[2]
+        if dsecs[1] != 'None':
+            vardims[dsecs[0]] = [dsecs[1], dsecs[2]]
+            print '  ',dsecs[0],':',dsecs[1],',',dsecs[2]
 
 if opts.obskind is None:
@@ -1551 +1553 @@
 
 print main + ': found',Nexactt,'Temporal same values in simulation and observations'
-print '  ',exacttvalues
-
 # Sim Coincident times
 ##
@@ -1571 +1571 @@
 print main + ': found',Ncoindt,'Simulation time-interval (within consecutive ' +     \
   'observed times) coincident times between simulation and observations'
-print '  ',coindtvalues
 
 if Ncoindt == 0:
@@ -1840 +1839 @@
         onewnc.sync()
 
-newvar = onewnc.createVariable('simtrj','i',('time','simtrj'))
-basicvardef(newvar, 'simtrj', 'coordinates [X,Y,Z,T] of the coincident trajectory ' +\
-  'in sim', obstunits)
-newvar[:] = trjsim.transpose()
+if trjsim is not None:
+    newvar = onewnc.createVariable('simtrj','i',('time','simtrj'))
+    basicvardef(newvar,'simtrj','coordinates [X,Y,Z,T] of the coincident ' +         \
+      'trajectory in sim', obstunits)
+    newvar[:] = trjsim.transpose()
 
 # Adding three variables with the station location, longitude, latitude and height
 if obskind == 'single-station':
-    adding_station_desc(objfile,stationdesc)
+    adding_station_desc(onewnc,stationdesc)
 
 # Global attributes