Changeset 2666 in lmdz_wrf

Timestamp:

Jul 12, 2019, 5:32:52 PM (5 years ago)

Author:

lfita

Message:

Starting to add copmutation of frontogenesis equation

File:

• trunk/tools/module_ForDiagnosticsVars.f90 (modified) (1 diff)

trunk/tools/module_ForDiagnosticsVars.f90

@@ -2047 +2047 @@
   END SUBROUTINE var_front_R04
 
+  SUBROUTINE var_Frontogenesis(dx, dy, dz, dt, theta, ua, va, wa, ddx, ddy, ddz, xdiab, ydiab, zdiab, \
+    xdef, ydef, zdef, xtilt, ytilt, zdiv, f)
+  ! Subroutine to compute the terms of the equation of frontogenesis
+  ! AFTER: https://en.wikipedia.org/wiki/Frontogenesis, http://glossary.ametsoc.org/wiki/Frontogenetical_function
+
+    IMPLICIT NONE
+
+    INTEGER, INTENT(in)                                  :: dx, dy, dz, dt
+    REAL(r_k), DIMENSION(dx,dy), INTENT(in)              :: ddx,ddy
+    REAL(r_k), DIMENSION(dx,dy,dz), INTENT(in)           :: ddz
+    REAL(r_k), DIMENSION(dx,dy,dz,dt), INTENT(in)        :: theta, ua, va, wa
+    REAL(r_k), DIMENSION(dx,dy,dz,dt), INTENT(out)       :: xdiab, ydiab, zdiab
+    REAL(r_k), DIMENSION(dx,dy,dz,dt), INTENT(out)       :: xdef, ydef, zdef
+    REAL(r_k), DIMENSION(dx,dy,dz,dt), INTENT(out)       :: xtilt, ytilt, zdiv
+    REAL(r_k), DIMENSION(dx,dy,dz,dt), INTENT(out)       :: f
+
+    ! Local
+    INTEGER                                              :: i,j,k,l,it
+
+!!!!!!! Variables
+! dx, dy, dz, dt: dimensions of the variables
+! theta: potential temperature [K]
+! ua, va, wa: x/y/z wind direction [ms-1]
+! ddx, ddy, ddz: x/y/z distances among grid points [m]
+! xdiab, ydiab, zdiab: x/y/z diabatic term [Ks-1m-1]
+! xdef, ydef, zdef: x/y/z deformation term [Ks-1m-1]
+! xtilt, ytilt: x/y tilting term [Ks-1m-1]
+! zdiv: vertical divergence term [Ks-1m-1]
+! f: frontogenetical function
+
+    fname = 'var_Frontogenesis'
+
+
+  END SUBROUTINE var_Frontogenesis
+
 END MODULE module_ForDiagnosticsVars