Changeset 1762 in lmdz_wrf for trunk/tools/module_ForDiagnostics.f90

Timestamp:

Jan 24, 2018, 3:54:22 PM (7 years ago)

Author:

lfita

Message:

Adding:

`mrso': total soil moisture
`slw': total liquid water

File:

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

trunk/tools/module_ForDiagnostics.f90

@@ -533 +533 @@
 
   END SUBROUTINE compute_psl_ptarget4d2
+
+  SUBROUTINE compute_massvertint1D(var, mutot, dz, deta, integral)
+    ! Subroutine to vertically integrate a 1D variable in eta vertical coordinates
+
+    IMPLICIT NONE
+
+    INTEGER, INTENT(in)                                  :: dz
+    REAL(r_k), INTENT(in)                                :: mutot
+    REAL(r_k), DIMENSION(dz), INTENT(in)                 :: var, deta
+    REAL(r_k), INTENT(out)                               :: integral
+
+! Local
+    INTEGER                                              :: k
+
+!!!!!!! Variables
+! var: vertical variable to integrate (assuming kgkg-1)
+! mutot: total dry-air mass in column
+! dz: vertical dimension
+! deta: eta-levels difference between full eta-layers
+
+    fname = 'compute_massvertint1D'
+
+!    integral=0.
+!    DO k=1,dz
+!      integral = integral + var(k)*deta(k)
+!    END DO
+     integral = SUM(var*deta)
+
+    integral=integral*mutot/g
+
+    RETURN
+
+  END SUBROUTINE compute_massvertint1D
+
+  SUBROUTINE compute_zint4D(var4D, dlev, zweight, d1, d2, d3, d4, int3D)
+    ! Subroutine to vertically integrate a 4D variable in any vertical coordinates
+
+    IMPLICIT NONE
+
+    INTEGER, INTENT(in)                                  :: d1,d2,d3,d4
+    REAL(r_k), DIMENSION(d1,d2,d3,d4), INTENT(in)        :: var4D, dlev, zweight
+    REAL(r_k), DIMENSION(d1,d2,d4), INTENT(out)          :: int3D
+
+! Local
+    INTEGER                                              :: i,j,l
+
+!!!!!!! Variables
+! var4D: vertical variable to integrate
+! dlev: height of layers
+! zweight: weight for each level to be applied (=1. for no effect)
+
+    fname = 'compute_zint4D'
+
+    DO i=1,d1
+      DO j=1,d2
+        DO l=1,d4
+          CALL compute_vertint1D(var4D(i,j,:,l),d3, dlev(i,j,:,l), zweight(i,j,:,l), &
+            int3D(i,j,l))
+        END DO
+      END DO
+    END DO
+
+    RETURN
+
+  END SUBROUTINE compute_zint4D
+
+  SUBROUTINE compute_vertint1D(var, dz, deta, zweight, integral)
+    ! Subroutine to vertically integrate a 1D variable in any vertical coordinates
+
+    IMPLICIT NONE
+
+    INTEGER, INTENT(in)                                  :: dz
+    REAL(r_k), DIMENSION(dz), INTENT(in)                 :: var, deta, zweight
+    REAL(r_k), INTENT(out)                               :: integral
+
+! Local
+    INTEGER                                              :: k
+
+!!!!!!! Variables
+! var: vertical variable to integrate
+! dz: vertical dimension
+! deta: eta-levels difference between layers
+! zweight: weight for each level to be applied (=1. for no effect)
+
+    fname = 'compute_vertint1D'
+
+!    integral=0.
+!    DO k=1,dz
+!      integral = integral + var(k)*deta(k)
+!    END DO
+    integral = SUM(var*deta*zweight)
+
+    RETURN
+
+  END SUBROUTINE compute_vertint1D
 
   SUBROUTINE compute_tv4d(ta,qv,tv,d1,d2,d3,d4)