Changeset 1692 in lmdz_wrf

Timestamp:

Dec 6, 2017, 3:52:49 PM (7 years ago)

Author:

lfita

Message:

Adding station's height from file

Location:

trunk/tools

Files:

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

trunk/tools/config_get_stations.WRF.py

@@ -44 +44 @@
   'T': ['WRFtime']}
 
-# dim Variables
+# dim Variables ('H' to determine station height from file)
 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'}
+  'bottom_top_stag': 'ZNW', 'H': 'HGT'}
 
 # SFC Variables
-sfcvariables = {'tas':'T2', 'uas':'WRFuas', 'vas':'WRFvas', 'huss':'Q2', 'ps':'PSFC',\
-  'zhgt':'HGT'}
+sfcvariables = {'tas':'T2', 'uas':'WRFuas', 'vas':'WRFvas', 'huss':'Q2', 'ps':'PSFC'}
 sfcrefvar = 'tas'
 

trunk/tools/get_stations.py

@@ -280 +280 @@
 
 # surface stations 
-if stcstatfile != 'None':
+if sfcstatfile != 'None':
     osfcstatf = open(sfcstatfile, 'r')
 
@@ -315 +315 @@
 
 # sounding stations 
-if stcstatfile != 'None':
+if sndstatfile != 'None':
     osndstatf = open(sndstatfile, 'r')
 
@@ -479 +479 @@
             ojvar = onewnc.variables['ipoint']
             ojvar[:] = minji[1]
+            ohvar = onewnc.variables['fheight']
+            ohvar[:] = dimvarvalues[dimvariables['H']][minji[0],minji[1]]
             onewnc.sync()
 
@@ -616 +618 @@
             ojvar = onewnc.variables['ipoint']
             ojvar[:] = minji[1]
+            ohvar = onewnc.variables['fheight']
+            ohvar[:] = dimvarvalues[dimvariables['H']][minji[0],minji[1]]
             onewnc.sync()