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

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.

File:

• LMDZ6/trunk/libf/phylmd/regr_pr_time_av_m.F90 (modified) (1 diff)

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