Changeset 1413 in lmdz_wrf

Timestamp:

Jan 17, 2017, 5:17:10 PM (8 years ago)

Author:

lfita

Message:

Adding `WRFghtrel': Relative altitude from surface from WRF

Location:

trunk/tools

Files:

• diagnostics.py (modified) (2 diffs)
• documentation/other/diagnostics.html (modified) (1 diff)
• variables_values.dat (modified) (1 diff)

trunk/tools/diagnostics.py

@@ -880 +880 @@
   'OMEGAw', 'RAINTOT',                                                               \
   'rvors', 'td', 'turbulence', 'WRFgeop', 'WRFp', 'WRFrvors', 'ws', 'wds', 'wss',    \
-  'WRFheight', 'WRFua', 'WRFva']
+  'WRFheight', 'WRFheightrel', 'WRFua', 'WRFva']
 
 methods = ['accum', 'deaccum']
@@ -1703 +1703 @@
         ncvar.insert_variable(ncobj, 'zhgt', diagout, dnames, diagoutvd, newnc)
 
+# WRFheightrel relative-height from WRF geopotential as WRFgeop(PH + PHB)/g-HGT 'WRFheightrel|PH@PHB@HGT
+    elif diag == 'WRFheightrel':
+        var0 = ncobj.variables[depvars[0]][:]
+        var1 = ncobj.variables[depvars[1]][:]
+        var2 = ncobj.variables[depvars[2]][:]
+
+        dimz = var0.shape[1]
+        diagout = np.zeros(tuple(var0.shape), dtype=np.float)
+        for iz in range(dimz):
+            diagout[:,iz,:,:] = var0[:,iz,:,:]+ var1[:,iz,:,:]/grav - var2
+
+        # Removing the nonChecking variable-dimensions from the initial list
+        varsadd = []
+        diagoutvd = list(dvnames)
+        for nonvd in NONchkvardims:
+            if gen.searchInlist(dvnames,nonvd): diagoutvd.remove(nonvd)
+            varsadd.append(nonvd)
+
+        ncvar.insert_variable(ncobj, 'zhgtrel', diagout, dnames, diagoutvd, newnc)
+
     else:
         print errormsg

trunk/tools/documentation/other/diagnostics.html

@@ -56 +56 @@
          <TR><TD>wssturb</TD><TD>turbulence</TD><TD>wss</TD><TD>Wind speed turbulence (Taylor's wss*= < wss > - wss)</TD></TR>
          <TR><TD>zg</TD><TD>WRFght</TD><TD>PH@PHB</TD><TD>Geopotential height from WRF</TD></TR>
+         <TR><TD>zhgt</TD><TD>WRFheight</TD><TD>PH@PHB@WRFgeop</TD>Altitude from WRF</TD></TR>
+         <TR><TD>zhgtrel</TD><TD>WRFheightrel</TD><TD>PH@PHB@HGT</TD>Relative altitude from surface from WRF</TD></TR>
       </TABLE>
     </DIV>

trunk/tools/variables_values.dat

@@ -595 +595 @@
 Z, z, z,  0., 100., z, -, Greens
 zhgt, zhgt, height,  0., 50000., altitude|from|surface, m, rainbow
+zhgtrel, zhgtrel, height_relative,  0., 50000., altitude|relative|to|surface, m, rainbow
 zg, zg, geopotential_height, 0., 80000., geopotential|height, m2s-2, rainbow
 WRFgeop, zg, geopotential_height, 0., 80000., geopotential|height, m2s-2, rainbow