Changeset 1706 in lmdz_wrf for trunk/tools/get_stations.py

Timestamp:

Dec 8, 2017, 10:06:34 PM (7 years ago)

Author:

lfita

Message:

Adding includion of time values!!

File:

• trunk/tools/get_stations.py (modified) (8 diffs)

trunk/tools/get_stations.py

@@ -30 +30 @@
 NONcheckingvars = Cvars + WRFvars + ['None']
 
-def creation_sfcstation_file(filen, lv, Lv, hv, lab, tunits, sfcvars, dimt, ifilen):
+def creation_sfcstation_file(filen, lv, Lv, hv, tv, lab, tunits, sfcvars, dimt,      \
+  ifilen):
     """ Function to create the structure of the surface station file
       filen: name of the file
@@ -53 +54 @@
     # Variable-dimensions
     newvar = onewnc.createVariable('time', 'f8', ('time'))
+    newvar[:] = tv
     ncvar.basicvardef(newvar, 'time', 'Time', tunits)
     ncvar.set_attribute(newvar, 'calendar', 'standard')
@@ -133 +135 @@
     return
 
-def creation_sndstation_file(filen, lv, Lv, hv, lab, tunits, punits, ptime, sndvars, \
-  dimt, dimz, ifilen):
+def creation_sndstation_file(filen, lv, Lv, hv, tv, lab, tunits, punits, ptime,      \
+  sndvars, dimt, dimz, ifilen):
     """ Function to create the structure of the surface station file
       filen: name of the file
@@ -161 +163 @@
     # Variable-dimensions
     newvar = onewnc.createVariable('time', 'f8', ('time'))
+    newvar[:] = tv
     ncvar.basicvardef(newvar, 'time', 'Time', tunits)
     ncvar.set_attribute(newvar, 'calendar', 'standard')
@@ -248 +251 @@
 
     # Global attributes
-    ncvar.add_global_PyNCplot(onewnc, 'get_stations.py', fname, '0.1')
+    ncvar.add_global_PyNCplot(onewnc, 'get_stations.py', fname, '0.1')<
     ncvar.set_attribute(newvar, 'data_origin', 'SMN')
 
@@ -409 +412 @@
         dimvarvalues[dimvariables[yn]] = newYvarvals
 
+tvals = dimvarvalues[dimvariables[CFdims['time']]]
+
 # Retrieving surface data
 ##
@@ -428 +433 @@
     stfilen = 'sfc_station_' + labelv + '.nc'
 
-    creation_sfcstation_file(stfilen, lonv, latv, heightv, labelv, utime,            \
+    creation_sfcstation_file(stfilen, lonv, latv, heightv, tvals, labelv, utime,     \
       sfcvariables.keys(), dt, simfilen)
 
@@ -568 +573 @@
     stfilen = 'snd_station_' + labelv + '.nc'
 
-    creation_sndstation_file(stfilen, lonv, latv, heightv, labelv, utime, UnitsPress,\
-      presstime, sndvariables.keys(), dt, dz, simfilen)
+    creation_sndstation_file(stfilen, lonv, latv, heightv, tvals, labelv, utime,     \
+      UnitsPress, presstime, sndvariables.keys(), dt, dz, simfilen)
 
     onewnc = NetCDFFile(stfilen, 'a')