Changeset 1809 in lmdz_wrf

Timestamp:

Mar 16, 2018, 6:31:51 PM (7 years ago)

Author:

lfita

Message:

Adding:

• `WRF_denszint': integration of a sum of variables with vertical dependency on density
• 'NONcheckdepvars': list of variables without checking their dependencies

Location:

trunk/tools

Files:

• diagnostics.inf (modified) (1 diff)
• diagnostics.py (modified) (5 diffs)
• documentation/other/diagnostics.html (modified) (1 diff)

trunk/tools/diagnostics.inf

@@ -50 +50 @@
 zg, WRFght, PH@PHB
 zmla, WRFzmlagen, T@QVAPOR@WRFgeop@HGT
+
+# Specials (variables which depend on the selection)
+varn WRF_denszint WRFdens@varn=clivi@QCLOUD@QICE@QGRAUPEL

trunk/tools/diagnostics.py

@@ -78 +78 @@
 availdiags = ['ACRAINTOT', 'accum', 'clt', 'cllmh', 'deaccum', 'LMDZrh', 'mslp',     \
   'OMEGAw', 'RAINTOT',                                                               \
-  'rvors', 'td', 'turbulence', 'WRFcape_afwa', 'WRFclivi', 'WRFclwvi', 'WRFgeop',    \
+  'rvors', 'td', 'turbulence', 'WRFcape_afwa', 'WRFclivi', 'WRFclwvi',               \
+  'WRF_denszint', 'WRFgeop',                                                         \
   'WRFmrso', 'WRFpotevap_orPM', 'WRFp', 'WRFpsl_ecmwf',                              \
   'WRFpsl_ptarget', 'WRFrvors', 'WRFslw', 'ws', 'wds', 'wss', 'WRFheight',           \
@@ -89 +90 @@
   'WRFclivi', 'WRFclwvi', 'WRFdens', 'WRFgeop',                                      \
   'WRFp', 'WRFtd',                                                                   \
-  'WRFpos', 'WRFpotevap_orPM', 'WRFprc', 'WRFprls', 'WRFrh', 'LMDZrh', 'LMDZrhs',    \
+  'WRFpos', 'WRFprc', 'WRFprls', 'WRFrh', 'LMDZrh', 'LMDZrhs',                       \
   'WRFrhs', 'WRFrvors',                                                              \
   'WRFt', 'WRFtime', 'WRFua', 'WRFva', 'WRFwds', 'WRFwss', 'WRFheight', 'WRFz']
+
+# diagnostics not to check their dependeny
+NONcheckdepvars = ['WRF_denszint', 'WRFzwind_log', 'WRFzwind', 'WRFzwindMO']
 
 NONchkvardims = ['WRFtime']
@@ -360 +364 @@
     diagn = diags[idiag].split('|')[0]
     depvars = diags[idiag].split('|')[1].split('@')
-    if diags[idiag].split('|')[1].find('@') != -1:
-        depvars = diags[idiag].split('|')[1].split('@')
-        if depvars[0] == 'deaccum': diagn='deaccum'
-        if depvars[0] == 'accum': diagn='accum'
-        for depv in depvars:
-            if not ncobj.variables.has_key(depv) and not                             \
-              gen.searchInlist(NONcheckingvars, depv) and                            \
-              not gen.searchInlist(methods, depv) and not depvars[0] == 'deaccum'    \
-              and not depvars[0] == 'accum' and not depv[0:2] == 'z=':
+    if not gen.searchInlist(NONcheckdepvars, diagn):
+        if diags[idiag].split('|')[1].find('@') != -1:
+            depvars = diags[idiag].split('|')[1].split('@')
+            if depvars[0] == 'deaccum': diagn='deaccum'
+            if depvars[0] == 'accum': diagn='accum'
+            for depv in depvars:
+                if not ncobj.variables.has_key(depv) and not                         \
+                  gen.searchInlist(NONcheckingvars, depv) and                        \
+                  not gen.searchInlist(methods, depv) and not depvars[0] == 'deaccum'\
+                  and not depvars[0] == 'accum' and not depv[0:2] == 'z=':
+                    print errormsg
+                    print '  ' + main + ": file '" + opts.ncfile +                   \
+                      "' does not have variable '" + depv + "' !!"
+                    quit(-1)
+        else:
+            depvars = diags[idiag].split('|')[1]
+            if not ncobj.variables.has_key(depvars) and not                          \
+              gen.searchInlist(NONcheckingvars, depvars) and                         \
+              not gen.searchInlist(methods, depvars):
                 print errormsg
                 print '  ' + main + ": file '" + opts.ncfile +                       \
-                  "' does not have variable '" + depv + "' !!"
+                  "' does not have variable '" + depvars + "' !!"
                 quit(-1)
-    else:
-        depvars = diags[idiag].split('|')[1]
-        if not ncobj.variables.has_key(depvars) and not                              \
-          gen.searchInlist(NONcheckingvars, depvars) and                             \
-          not gen.searchInlist(methods, depvars):
-            print errormsg
-            print '  ' + main + ": file '" + opts.ncfile +                           \
-              "' does not have variable '" + depvars + "' !!"
-            quit(-1)
 
     print "\n    Computing '" + diagn + "' from: ", depvars, '...'
@@ -745 +750 @@
 
         newvarn = depvars + 'turb'
-        print main + '; Lluis newvarn:', newvarn
         ncvar.insert_variable(ncobj, newvarn, diagout, diagoutd, 
           diagoutvd, newnc)
-        print main + '; Lluis variables:', newnc.variables.keys()
         varobj = newnc.variables[newvarn]
         attrv = varobj.long_name
@@ -847 +850 @@
             varsadd.append(nonvd)
         ncvar.insert_variable(ncobj, 'clwvi', diagout, diagoutd, diagoutvd, newnc)
+
+# WRF_denszint WRF vertical integration as WRFdens, sum(Q[water species1], ..., Q[water speciesN]), varn=[varN]
+    elif diagn == 'WRF_denszint':
+            
+        var0 = WRFdens
+        varn = depvars[1].split('=')[1]
+        qtot = ncobj.variables[depvars[2]]
+        qtotv = ncobj.variables[depvars[2]]
+        Nspecies = len(depvars) - 2
+        for iv in range(Nspecies):
+            if ncobj.variables.has_key(depvars[iv+2]):
+                var1 = ncobj.variables[depvars[iv+2]]
+                qtotv = qtotv + var1[:]
+
+        dnamesvar = list(qtot.dimensions)
+        dvnamesvar = ncvar.var_dim_dimv(dnamesvar,dnames,dvnames)
+
+        diagout, diagoutd, diagoutvd = diag.compute_clwvl(var0, qtotv, dnamesvar,dvnamesvar)
+
+        # 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, varn, diagout, diagoutd, diagoutvd, newnc)
 
 # WRFgeop geopotential from WRF as PH + PHB

trunk/tools/documentation/other/diagnostics.html

@@ -113 +113 @@
        </TBODY>
       </TABLE>
+      Special diagnostics (variables which depend on the selection):<BR>
+      <TABLE ID="t01">
+       <THEAD>
+         <TR><TH>CF Name</TH><TH>Name in script</TH><TH>variables</TH><TH>description</TH></TR>
+       </THEAD>
+       <TFOOT>
+       </TFOOT>
+       <TBODY>
+         <TR><TD>[varn]</TD><TD>WRF_denszint</TD><TD>WRFdens@varn=clivi@QCLOUD@QICE@QGRAUPEL</TD><TD> Vertical mass integration (vertical density dependence) of a series of variables from WRF with name [varn]</TD></TR>
+       </TBODY>
     </DIV>
     <DIV CLASS="valins">