Changeset 2643 in lmdz_wrf

Timestamp:

Jun 28, 2019, 9:42:25 PM (5 years ago)

Author:

lfita

Message:

Adding:

• Computing of front following Rodrigues et al 2004

Location:

trunk/tools

Files:

• diag_tools.py (modified) (2 diffs)
• diagnostics.py (modified) (2 diffs)

trunk/tools/diag_tools.py

@@ -41 +41 @@
 # Forcompute_fog_RUC: Computation of fog and visibility following RUC method Smirnova, (2000)
 # Forcompute_fog_FRAML50: fog and visibility following Gultepe and Milbrandt, (2010)
+# Forcompute_front_R04: Function to compute front following Rodrigues et al.(2004)
 # Forcompute_potevap_orPM: Function to compute potential evapotranspiration following
 #   Penman-Monteith formulation implemented in ORCHIDEE
@@ -1531 +1532 @@
     return hgtmax, pthgtmax, dhgt, peaks, valleys, origfaces, filtfaces, vardims,    \
       varvdims, ranges, rangeshgtmax, ptrangeshgtmax
+
+def Forcompute_front_R04(tas, uas, vas, dxs, dys, dimns, dimvns):
+    """ Function to compute front following Rodrigues et al.(2004), Rev. Bras. 
+          Geofis. 22, 135-151
+    Forcompute_front_R04(tas, vas, dimns, dimvns)
+      [tas]= surface air-temperature values (assuming [[t],y,x]) [K]
+      [uas]= latitudinal wind (assuming [[t],y,x]) [ms-1]
+      [vas]= meridional wind (assuming [[t],y,x]) [ms-1]
+      [dxs]= grid spacing between grid points along x-axis (assuming [y,x]) [m]
+      [dys]= grid spacing between grid points along y-axis (assuming [y,x]) [m]
+      [dimns]= list of the name of the dimensions of [pa]
+      [dimvns]= list of the name of the variables with the values of the 
+        dimensions of [pa]
+    """
+    fname = 'Forcompute_front_R04'
+
+    psldims = dimns[:]
+    pslvdims = dimvns[:]
+
+    if len(pa.shape) == 3:
+        front = np.zeros((tas.shape[0],tas.shape[1],tas.shape[2]), dtype=np.float)
+
+        dx = tas.shape[2]
+        dy = tas.shape[1]
+        dt = tas.shape[0]
+
+        pfront=fdin.module_fordiagnostics.compute_front_R04d3(tas=tas[:].transpose(),\
+          uas=uas[:].transpose(), vas=vas[:].transpose(), dxs=dxs[:].transpose(),    \
+          dys=dys[:].transpose(), d1=dx, d2=dy, d3=dt)
+        front = pfront.transpose()
+    else:
+        print errormsg
+        print '  ' + fname + ': rank', len(pa.shape), 'not ready !!'
+        print '  it only computes 3D [t,y,x] rank values'
+        quit(-1)
+
+    return front
 
 ####### ###### ##### #### ### ## # END Fortran diagnostics

trunk/tools/diagnostics.py

@@ -90 +90 @@
     #######
 availdiags = ['ACRAINTOT', 'accum', 'clt', 'cllmh', 'convini', 'deaccum', 'fog_K84', \
-  'fog_RUC', 'rhs_tas_tds', 'LMDZrh', 'mslp', 'OMEGAw', 'RAINTOT', 'range_faces',    \
+  'fog_RUC', 'front_R04', 'rhs_tas_tds', 'LMDZrh', 'mslp', 'OMEGAw', 'RAINTOT',      \
+  'range_faces',                                                                     \
   'rvors', 'td', 'timemax', 'timeoverthres', 'turbulence', 'tws', 'uavaFROMwswd',    \
   'WRFbnds', 'WRFcape_afwa', 'WRFclivi', 'WRFclwvi', 'WRF_denszint', 'WRFgeop',      \
@@ -773 +774 @@
         ncvar.insert_variable(ncobj, 'fogvisblty', diag2, diagoutd, diagoutvd, newnc)
 
+# front_R04: Computation of the presence of a front as defined by Rodrigues et al. 
+#     (2004) (vas, tas, dxs, dys) 
+#   Rodrigues et al. 2004: Climatologia de frentes frias no litoral de Santa Catarina,
+#     Rev. Bras. Geof. v.22 n.2 Sao Paulo maio/ago. DOI: 10.1590/S0102-261X2004000200004  
+    elif diagn == 'front_R04':
+
+        var0 = ncobj.variables[depvars[0]]
+        var1 = ncobj.variables[depvars[1]]
+        var2 = ncobj.variables[depvars[2]]
+        var3 = ncobj.variables[depvars[3]]
+
+        dnamesvar = list(var0.dimensions)
+        dvnamesvar = ncvar.var_dim_dimv(dnamesvar,dnames,dvnames)
+
+        diag1, diagoutd, diagoutvd = diag.Forcompute_front_R04(var0, var1,           \
+          var2, var3, 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, 'front', diag1, diagoutd, diagoutvd, newnc)
+
 # LMDZrh (pres, t, r)
     elif diagn == 'LMDZrh':