Changeset 3123

Timestamp:

Dec 13, 2017, 2:16:38 PM (6 years ago)

Author:

oboucher

Message:

Compute z_tropopause for diagnostic purpose from geoid rather than from the surface as per CF convention

Location:

LMDZ6/trunk/libf/phylmd

Files:

• physiq_mod.F90 (modified) (2 diffs)
• stratosphere_mask.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -3537 +3537 @@
        !
        !--WMO criterion to determine tropopause
-       CALL stratosphere_mask(missing_val, t_seri, pplay, latitude_deg)
+       CALL stratosphere_mask(missing_val, pphis, t_seri, pplay, latitude_deg)
        !
        !--STRAT AEROSOL
@@ -3577 +3577 @@
 #ifdef CPP_StratAer
        !--compute stratospheric mask
-       CALL stratosphere_mask(missing_val, t_seri, pplay, latitude_deg)
+       CALL stratosphere_mask(missing_val, pphis, t_seri, pplay, latitude_deg)
        !--interactive strat aerosols
        CALL calcaerosolstrato_rrtm(pplay,t_seri,paprs,debut)

LMDZ6/trunk/libf/phylmd/stratosphere_mask.F90

@@ -2 +2 @@
 ! $Id$
 !
-SUBROUTINE stratosphere_mask(missing_val, t_seri, pplay, xlat)
+SUBROUTINE stratosphere_mask(missing_val, pphis, t_seri, pplay, xlat)
 
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -39 +39 @@
 IMPLICIT NONE
 
+INCLUDE "YOMCST.h"
+
 REAL, INTENT(IN)                       :: missing_val ! missing value, also XIOS
+REAL,DIMENSION(klon),INTENT(IN)        :: pphis   ! Geopotentiel de surface
 REAL,DIMENSION(klon,klev),INTENT(IN)   :: t_seri  ! Temperature
 REAL,DIMENSION(klon,klev),INTENT(IN)   :: pplay   ! pression pour le mileu de chaque couche (en Pa)
@@ -62 +65 @@
   ENDDO
   psrf=pplay(i,1)
-  zsrf=0.0
+  zsrf=pphis(i)/RG           !--altitude de la surface
   call twmo(missing_val, klev, t, p, psrf, zsrf, plimu, pliml, gamma, ptrp, ttrp, ztrp)
   tp(i)=ptrp