Changeset 1703 in lmdz_wrf for trunk/tools

Timestamp:

Dec 8, 2017, 7:42:45 PM (8 years ago)

Author:

lfita

Message:

Fixing the right config files and the right variable/dimensions to use

Location:

trunk/tools

Files:

• config_get_stations.NNRP.py (modified) (1 diff)
• config_get_stations.WRF.py (modified) (1 diff)
• get_stations.py (modified) (2 diffs)

trunk/tools/config_get_stations.NNRP.py

@@ -47 +47 @@
 sfcvariables = {'tas':'T2', 'uas':'WRFuas', 'vas':'WRFvas', 'huss':'Q2', 'ps':'PSFC',\
   'zhgt':'HGT'}
-sfcrefvar = 'tas'
+# Surface reference variable (it must be one variable from file which has 
+#   the right right shape (e.g. un-staggered) and rang as 3D dimensions [T,Y,X])
+sfcrefvar = 'T2'
 
 # Vertical Variables
 sndvariables= {'ta':'air', 'ua':'uwnd', 'va':'vwnd', 'hus':'shum',             \
   'zg':'hgt', 'plev':'level', 'tda': 'C_td', 'wd': 'C_wd', 'ws': 'C_ws'}
+# Sounding reference variable (it must be one variable from file which has 
+#   the right right shape (e.g. un-staggered) and rang as 4D dimensions [T,Z,Y,X])
 sndrefvar = 'air'
 

trunk/tools/config_get_stations.WRF.py

@@ -44 +44 @@
   'T': ['WRFtime']}
 
-# dim Variables ('H' to determine station height from file)
+# dim Variables
 dimvariables = {'west_east': 'XLONG', 'south_north': 'XLAT', 'bottom_top': 'ZNU',    \
   'Time': 'WRFtime', 'west_east_stag': 'XLONG_U', 'south_north_stag': 'XLAT_U',      \
-  'bottom_top_stag': 'ZNW', 'H': 'HGT'}
+  'bottom_top_stag': 'ZNW'}
 
 # SFC Variables
-sfcvariables = {'tas':'T2', 'uas':'WRFuas', 'vas':'WRFvas', 'huss':'Q2', 'ps':'PSFC'}
+sfcvariables = {'tas':'T2', 'uas':'WRFuas', 'vas':'WRFvas', 'huss':'Q2', 'ps':'PSFC',\
+  'zhgt':'HGT'}
 sfcrefvar = 'tas'
 

trunk/tools/get_stations.py

@@ -386 +386 @@
 
     if dimn == Pressdimref: dz = len(onc.dimensions[dimn])
+    dt = onc.variables[sfcrefvar].shape[0]
 
 # Looking for 2D 'lon', 'lat' related variables
-for xn in axesvars['X']:
+for xn in axesdims['X']:
     Xvarvals = dimvarvalues[dimvariables[xn]]
-    Yvarvals = dimvarvalues[CFdims['lat']]
+    Yvarvals = dimvarvalues[dimvariables[CFdims['lat']]]
     if len(Xvarvals.shape) == 1:
         newXvarvals, newYvarvals = np.meshgrid(Xvarvals, Yvarvals) 
         dimvarvalues[dimvariables[xn]] = newXvarvals
         if len(Yvarvals.shape) == 1:
-            dimvarvalues[dimvariables['lat']] = newYvarvals
-
-for yn in axesvars['Y']:
-    Xvarvals = dimvarvalues[CFdims['lon']]
+            dimvarvalues[dimvariables[CFdims['lat']]] = newYvarvals
+
+for yn in axesdims['Y']:
+    Xvarvals = dimvarvalues[dimvariables[CFdims['lon']]]
     Yvarvals = dimvarvalues[dimvariables[yn]]
     if len(Yvarvals.shape) == 1:
@@ -676 +677 @@
         # Writting data and slicing
         onewvar = onewnc.variables[sndv]
-        print 'Lluis sndv:', sndv, ' shapes onewvar:', onewvar.shape, 'slice:', slicevar
         onewvar[:] = ogetvar[tuple(slicevar)]    
         onewnc.sync()