Changeset 3141 for LMDZ6/trunk/libf/phylmd/tropopause_m.F90

Timestamp:

Dec 19, 2017, 7:49:05 PM (7 years ago)

Author:

dcugnet

Message:

• Fix: Sig_bot was used before being initialized in "regr_pr_time_av".
• Add an optional output argument to the "dyn_tropopause" function: the table of the last tropospheric layer in each column.
• A filter for irrealistic tropopause pressure (lower than 80hPa or greater than 400hPa) is included in "dyn_tropopause". This is safer if this function is not used together with "regr_pr_time_av" (which contains already this filter), which is the case for the LMDZ version with coupled REPROBUS and INCA chemistry.

File:

• LMDZ6/trunk/libf/phylmd/tropopause_m.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/tropopause_m.F90

@@ -9 +9 @@
 !-------------------------------------------------------------------------------
 !
-FUNCTION dyn_tropopause(t, ts, paprs, pplay, rot, thet0, pvor0)
+FUNCTION dyn_tropopause(t, ts, paprs, pplay, rot, itrop, thet0, pvor0)
 !
 !-------------------------------------------------------------------------------
@@ -26 +26 @@
   REAL, INTENT(IN)  ::  pplay(:,:) !--- Cells-centers pressure
   REAL, INTENT(IN)  ::    rot(:,:) !--- Cells-centers relative vorticity
-  REAL, INTENT(IN), OPTIONAL :: thet0, pvor0
+  INTEGER, INTENT(OUT), OPTIONAL :: itrop(klon) !--- Last tropospheric layer idx
+  REAL,    INTENT(IN),  OPTIONAL :: thet0, pvor0
 !-------------------------------------------------------------------------------
 ! Local variables:
   include "YOMCST.h"
+  REAL, PARAMETER :: DynPTrMin =8.E+3 !--- Thresholds for minimum and maximum
+  REAL, PARAMETER :: DynPTrMax =4.E+4 !    dynamical tropopause pressure (Pa).
   CHARACTER(LEN=80)  :: sub
   INTEGER :: i, k, kb, kt, kp, ib, ie, nw
@@ -100 +103 @@
     dyn_tropopause(i) = MAX( pplay(i,kp+1)*(pplay(i,kp)/pplay(i,kp+1))**al,    &
                             dyn_tropopause(i) )
+    !--- UNREALISTIC VALUES DETECTION
+    IF(dyn_tropopause(i)<DynPTrMin.OR.dyn_tropopause(i)>DynPTrMax) THEN
+      dyn_tropopause(i)=MIN(MAX(dyn_tropopause(i),DynPTrMax),DynPTrMin)
+      DO kt=1,klev-1; IF(pplay(i,kt+1)>dyn_tropopause(i)) EXIT; END DO; kp=kt
+    END IF
+    !--- LAST TROPOSPHERIC LAYER INDEX NEEDED
+    IF(PRESENT(itrop)) itrop(i)=MAX(kt,kp)
   END DO