Changeset 2669 in lmdz_wrf

Timestamp:

Jul 12, 2019, 6:07:25 PM (5 years ago)

Author:

lfita

Message:

Advancing on the implementation of the frontogenesis diagnostic

File:

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

trunk/tools/module_ForDiagnosticsVars.f90

@@ -2056 +2056 @@
     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(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
@@ -2065 +2065 @@
     ! 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
 
 !!!!!!! Variables
@@ -2079 +2082 @@
     fname = 'var_Frontogenesis'
 
+    ! Computing components separately by time-step
+    DO it=1, dt
+      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 deformation3D(dx, dy, dz, thetagrad, uagrad, vagrad, thetadef)
+      CALL tilting3D(dx, dy, dz, thetagrad, wagrad, thetatilt)
+      xdef(:,:,:,it) = thetadef(:,:,:,1)
+      ydef(:,:,:,it) = thetadef(:,:,:,2)
+      zdef(:,:,:,it) = thetadef(:,:,:,3)
+      xtilt(:,:,:,it) = thetatilt(:,:,:,1)
+      ytilt(:,:,:,it) = thetatilt(:,:,:,2)
+      zdiv(:,:,:,it) = thetatilt(:,:,:,3)
+
+    END DO
+
 
   END SUBROUTINE var_Frontogenesis