Changeset 2670 in lmdz_wrf

Timestamp:

Jul 12, 2019, 6:14:01 PM (5 years ago)

Author:

lfita

Message:

Adding pressure !

File:

• trunk/tools/module_ForDiagnosticsVars.f90 (modified) (5 diffs)

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 var_Frontogenesis(dx, dy, dz, dt, theta, ua, va, wa, press, ddx, ddy, ddz, ddt, 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
@@ -2055 +2055 @@
 
     INTEGER, INTENT(in)                                  :: dx, dy, dz, dt
+    REAL(r_k), INTENT(in)                                :: ddt
     REAL(r_k), DIMENSION(dx,dy), INTENT(in)              :: ddx,ddy
     REAL(r_k), DIMENSION(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(in)        :: theta, ua, va, wa, press
     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
@@ -2065 +2066 @@
     ! Local
     INTEGER                                              :: i,j,k,l,it
-    REAL(r_k), DIMENSION(dx,dy,dz)                       :: modthetagrad
-    REAL(r_k), DIMENSION(dx,dy,dz,3)                     :: thetagrad, uagrad, vagrad, wagrad,        &
-      thetadef, thetatilt
+    REAL(r_k), DIMENSION(dx,dy,dz)                       :: modthetagrad, thetadt
+    REAL(r_k), DIMENSION(dx,dy,dz,3)                     :: thetagrad, thetadtgrad, uagrad, vagrad,   &
+      wagrad, pressgrad, thetadef, thetatilt
 
 !!!!!!! Variables
@@ -2073 +2074 @@
 ! theta: potential temperature [K]
 ! ua, va, wa: x/y/z wind direction [ms-1]
+! press: pressure [Pa]
+! ddt: temporal delta [s]
 ! ddx, ddy, ddz: x/y/z distances among grid points [m]
 ! xdiab, ydiab, zdiab: x/y/z diabatic term [Ks-1m-1]
@@ -2082 +2085 @@
     fname = 'var_Frontogenesis'
 
+    CALL deltat3D(dx, dy, dz, dt, theta, ddt, 'forward', thetadt)
+
     ! Computing components separately by time-step
     DO it=1, dt
+      CALL gradient3D_1o(dx, dy, dz, thetadt(:,:,:,it), ddx, ddy, ddz, thetadtgrad)
       CALL gradient3D_1o(dx, dy, dz, theta(:,:,:,it), ddx, ddy, ddz, thetagrad)
       CALL gradient3D_1o(dx, dy, dz, ua(:,:,:,it), ddx, ddy, ddz, uagrad)
       CALL gradient3D_1o(dx, dy, dz, va(:,:,:,it), ddx, ddy, ddz, vagrad)
       CALL gradient3D_1o(dx, dy, dz, wa(:,:,:,it), ddx, ddy, ddz, wagrad)
+      CALL gradient3D_1o(dx, dy, dz, press(:,:,:,it), ddx, ddy, ddz, pressgrad)
       CALL deformation3D(dx, dy, dz, thetagrad, uagrad, vagrad, thetadef)
       CALL tilting3D(dx, dy, dz, thetagrad, wagrad, thetatilt)
+      xdiab(:,:,:,it) = thetadtgrad(:,:,:,1)
+
       xdef(:,:,:,it) = thetadef(:,:,:,1)
       ydef(:,:,:,it) = thetadef(:,:,:,2)