Changeset 2668 in lmdz_wrf

Timestamp:

Jul 12, 2019, 5:58:40 PM (5 years ago)

Author:

lfita

Message:

Adding:

• `deltat3D': Subroutine to compute the temporal derivative of a 3D field

File:

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

trunk/tools/module_scientific.f90

@@ -107 +107 @@
 ! curl2D_c1o: Subroutine to compute the first order centered curl of a 2D vectorial field
 ! deformation3D: Subroutine to compute the deformation of a 3D field
+! deltat3D: Subroutine to compute the temporal derivative of a 3D field
 ! divergence2D_1o: Subroutine to compute the first order divergence of a 2D vectorial field
 ! divergence2D_c1o: Subroutine to compute the first order centered divergence of a 2D vectorial field
@@ -7967 +7968 @@
   SUBROUTINE tilting3D(dx, dy, dz, vargrad, wagrad, tilt)
   ! Subroutine to compute the tilting of a 3D field
-  !   FROM: From: https://en.wikipedia.org/wiki/Finite_difference_coefficient
 
     IMPLICIT NONE
@@ -7994 +7994 @@
   END SUBROUTINE tilting3D
 
+  SUBROUTINE deltat3D(dx, dy, dz, dt, var, ddt, dtkind, vardt)
+  ! Subroutine to compute the temporal derivative of a 3D field
+
+    IMPLICIT NONE
+
+    INTEGER, INTENT(in)                                  :: dx, dy, dz, dt
+    REAL(r_k), DIMENSION(dx,dy,dz,dt), INTENT(in)        :: vargrad
+    CHARACTER(LEN=*), INTENT(in)                         :: dtkind
+    REAL(r_k), DIMENSION(dx,dy,dz,dt), INTENT(out)       :: vardt
+
+! Local 
+    INTEGER                                              :: i, j, k, it
+
+!!!!!!! Variables
+! dx, dy, dz: shape of the 3D field
+! dt: number of time-steps
+! var: variable [X]
+! ddt: delta t [s]
+! dtkind: kind of delta t to compute
+!   'backward': [X(it) - X(it-1)] / ddt
+!   'centered': [X(it+1) - 2*X(it) + X(it-1)] / 2*ddt
+!   'forward': [X(it+1) - X(it)] / ddt
+! vardt: temporal derivative of var [X s-1]
+
+    fname = 'deltat3D'
+
+    vardt = fillval64
+
+    IF (TRIM(dtkind) == 'backward') THEN
+      DO it=2, dt
+        vardt(:,:,:,it) = (var(:,:,:,it) - var(:,:,:,it-1))/ddt
+      END DO
+    ELSE IF (TRIM(dtkind) == 'centered') THEN
+      DO it=2, dt-1
+        vardt(:,:,:,it) = (var(:,:,:,it+1) - 2*var(:,:,:,it) + var(:,:,:,it-1))/(2.*ddt)
+      END DO
+    ELSE IF (TRIM(dtkind) == 'forward') THEN
+      DO it=1, dt-1
+        vardt(:,:,:,it) = (var(:,:,:,it+1) - var(:,:,:,it))/ddt
+      END DO
+    ELSE
+      msg = "kind of delta t '" // TRIM(dtkind) // "' not ready !!"
+      CALL ErrMsg(msg, fname, -1)
+    END IF    
+
+  END SUBROUTINE tilting3D
+
 END MODULE module_scientific
+