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

Timestamp:

Mar 12, 2018, 3:37:19 PM (7 years ago)

Author:

lfita

Message:

Adding:

`psl_ecmwf': sea-level pressure computation following ECMWF

File:

• trunk/tools/module_ForDiagnostics.f90 (modified) (2 diffs)

trunk/tools/module_ForDiagnostics.f90

@@ -24 +24 @@
 ! compute_clt3D3: Computation of total cloudiness from a 3D CLDFRA being 3rd dimension the z-dim
 ! compute_clt: Computation of total cloudiness
+! compute_psl_ecmwf: Compute sea level pressure using ECMWF method following Mats Hamrud and Philippe Courtier [Pa]
 ! compute_massvertint1D: Subroutine to vertically integrate a 1D variable in eta vertical coordinates
 ! compute_vertint1D: Subroutine to vertically integrate a 1D variable in any vertical coordinates
@@ -602 +603 @@
   END SUBROUTINE compute_cape_afwa4D
 
+  SUBROUTINE compute_psl_ecmwf(ps, hgt, T, press, unpress, psl, d1, d2, d4)
+! Subroutine to compute sea level pressure using ECMWF method following Mats Hamrud and Philippe Courtier [Pa]
+
+    IMPLICIT NONE
+
+    INTEGER, INTENT(in)                                  :: d1, d2, d4
+    REAL(r_k), DIMENSION(d1,d2,d4), INTENT(in)           :: ps, T, press, unpress
+    REAL(r_k), DIMENSION(d1,d2), INTENT(in)              :: hgt
+    REAL(r_k), DIMENSION(d1,d2,d4), INTENT(out)          :: psl
+ 
+! Local
+    INTEGER                                              :: i, j, it
+
+!!!!!!! Variables
+! ps: surface pressure [Pa]
+! hgt: terrain height [m]
+! T: temperature at first half-mass level [K]
+! press: pressure at first full levels [Pa]
+! unpress: pressure at first mass (half) levels [Pa]
+! psl: sea-level pressure [Pa]
+
+    fname = 'compute_psl_ecmwf'
+
+    DO i=1, d1
+      DO j=1, d2
+        DO it=1, d4
+          CALL var_psl_ecmwf(ps(i,j,it), hgt(i,j), T(i,j,it), unpress(i,j,it), press(i,j,it),         &
+            psl(i,j,it))
+        END DO
+      END DO
+    END DO
+
+    RETURN
+
+  END SUBROUTINE compute_psl_ecmwf
+
   SUBROUTINE compute_zmla_generic4D(tpot, qratio, z, hgt, zmla3D, d1, d2, d3, d4)
 ! Subroutine to compute pbl-height following a generic method