Changeset 2100 in lmdz_wrf

Timestamp:

Aug 25, 2018, 10:57:00 PM (7 years ago)

Author:

lfita

Message:

Adding:

• `convini': convective initialization (pr(t) > 0.0001) in time units

Location:

trunk/tools

Files:

• diag_tools.py (modified) (2 diffs)
• diagnostics.inf (modified) (1 diff)
• diagnostics.py (modified) (2 diffs)
• variables_values.dat (modified) (1 diff)

trunk/tools/diag_tools.py

@@ -55 +55 @@
 # var_clt: Function to compute the total cloud fraction following 'newmicro.F90' from
 #   LMDZ using 1D vertical column values
+# var_convini: Function returns convective initialization (pr(t) > 0.0001) in time units
 # var_hur: Function to compute relative humidity following 'August - Roche - Magnus' 
 #   formula
@@ -663 +664 @@
 
     return virtual
+
+def var_convini(pr, time, dimns, dimvns):
+    """ This function returns convective initialization (pr(t) > 0.0001) in time units
+      pr: precipitation fux [kgm-2s-1]
+      time: time in CF coordinates
+    """
+    fname = 'var_convini'
+
+    dt = pr.shape[0]
+    dy = pr.shape[1]
+    dx = pr.shape[2]
+
+    vardims = dimns[:]
+    varvdims = dimvns[:]
+
+    vardims.pop(0)
+    varvdims.pop(0)
+
+    prmin = 0.0001
+    convini = np.ones((dy, dx), dtype=np.float)*gen.fillValueF
+    for it in range(dt):
+# NOT working ?
+#        convini = np.where(convini == gen.fillValueF and pr[it,:,:] >= prmin,        \
+#          time[it], fillValueF)
+        for j in range(dy):
+            for i in range(dx):
+                if convini[j,i] == gen.fillValueF and pr[it,j,i] >= prmin:
+                    convini[j,i] = time[it]
+                    break
+
+    return convini, vardims, varvdims
 
 def Forcompute_zint(var, zinterlev, zweights, dimns, dimvns):

trunk/tools/diagnostics.inf

@@ -11 +11 @@
 clivi, WRFclivi, WRFdens@QICE@QHAIL@QGRAUPEL
 clwvi, WRFclwvi, WRFdens@QCLOUD@QICE@QHAIL@QGRAUPEL
+convini, convini, pr@time
 fog, fog_FRAML50, QVAPOR@WRFt@WRFp
 fog, fog_FRAML50, Q2@T2@PSFC

trunk/tools/diagnostics.py

@@ -84 +84 @@
 ## MAIN
     #######
-availdiags = ['ACRAINTOT', 'accum', 'clt', 'cllmh', 'deaccum', 'fog_K84', 'fog_RUC', \
-  'LMDZrh', 'mslp', 'OMEGAw', 'RAINTOT',                                             \
+availdiags = ['ACRAINTOT', 'accum', 'clt', 'cllmh', 'convini', 'deaccum', 'fog_K84', \
+  'fog_RUC', 'LMDZrh', 'mslp', 'OMEGAw', 'RAINTOT',                                  \
   'rvors', 'td', 'turbulence', 'uavaFROMwswd', 'WRFcape_afwa', 'WRFclivi',           \
   'WRFclwvi', 'WRF_denszint', 'WRFgeop',                                             \
@@ -520 +520 @@
         ncvar.insert_variable(ncobj, 'clt', diagout, diagoutd, diagoutvd, newnc)
 
+# convini (pr, time)
+    elif diagn == 'convini':
+            
+        var0 = ncobj.variables[depvars[0]][:]
+        var1 = ncobj.variables[depvars[1]][:]
+        otime = ncobj.variables[depvars[1]]
+
+        dnamesvar = ncobj.variables[depvars[0]].dimensions
+        dvnamesvar = ncvar.var_dim_dimv(dnamesvar,dnames,dvnames)
+
+        diagout, diagoutd, diagoutvd  = diag.var_convini(var0, var1, dnames, dvnames)
+
+        ncvar.insert_variable(ncobj, 'convini', diagout, diagoutd, diagoutvd, newnc, \
+          fill=gen.fillValueF)
+        # Getting the right units
+        ovar = newnc.variables['convini']
+        if gen.searchInlist(otime.ncattrs(), 'units'): 
+            tunits = otime.getncattr('units')
+            ncvar.set_attribute(ovar, 'units', tunits)
+            newnc.sync()
+
 # deaccum: deacumulation of any variable as (Variable, time [as [tunits] 
 #   from/since ....], newvarname)

trunk/tools/variables_values.dat

@@ -92 +92 @@
 CLW, clw, total_cloud_properties, 0., 3., total|column-integrated|water|vapor|and|cloud|variables, 1, Purples, $clw$, clw
 clw, clw, total_cloud_properties, 0., 3., total|column-integrated|water|vapor|and|cloud|variables, 1, Purples, $clw$, clw
+convini, convini, convection_initialization, 0., 300000., initialization|of|convection|(pr>=0.00001), time, rainbow, $conv^{ini}$, convini
 cpt, cpt, convective_point, 0., 1., convective|point, 1, seismic, $cpt$, cpt
 ptconv, cpt, convective_point, 0., 1., convective|point, 1, seismic, $cpt$, cpt