Changeset 3141

Timestamp:

Dec 19, 2017, 7:49:05 PM (6 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.

Location:

LMDZ6/trunk/libf/phylmd

Files:

• regr_pr_time_av_m.F90 (modified) (1 diff)
• tropopause_m.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/regr_pr_time_av_m.F90

@@ -312 +312 @@
       Sig_to0 = MIN(SigT_in,SigT_ou)               !--- highest tropopause
       beta    = (Sig_bo0/Sig_to0)**gamm            !--- stretching exponent
-      Sig_bot = MIN(Sig_bo0*beta,0.1*(9.+Sig_bot)) !--- must be <1
+      Sig_bot = MIN(Sig_bo0*beta,0.1*(9.+Sig_bo0)) !--- must be <1
       ibot = locate(Sig_ou(:),Sig_bot)             !--- layer index
       IF(ibot-iout<2) THEN                         !--- at least one layer thick

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