Changeset 3120

Timestamp:

Dec 9, 2017, 6:49:01 PM (7 years ago)

Author:

Laurent Fairhead

Message:

Inclusion of bugfix 3119 from trunk
OB

File:

• LMDZ6/branches/IPSLCM6.0.14/libf/phylmd/stratosphere_mask.F90 (modified) (4 diffs)

LMDZ6/branches/IPSLCM6.0.14/libf/phylmd/stratosphere_mask.F90

@@ -8 +8 @@
 ! determination of tropopause height and temperature from gridded temperature data
 !
-! reference: Reichler, T., M. Dameris, and R. Sausen (2003):
+! reference: Reichler, T., M. Dameris, and R. Sausen (GRL, 10.1029/2003GL018240, 2003) 
 ! modified: 6/28/06 tjr
 ! adapted to LMDZ by C. Kleinschmitt (2016-02-15)
-!
-! input:    temp(nlon,nlat,nlev)    3D-temperature field 
-!       ps(nlon,nlat)       2D-surface pressure field
-!       zs(nlon,nlat)       2D-surface height
-!       nlon            grid points in x
-!       nlat            grid points in y
-!       pfull(nlon,nlat,nlev)   full pressure levels in Pa
-!       plimu           upper limit for tropopause pressure
-!       pliml           lower limit for tropopause pressure 
-!       gamma           tropopause criterion, e.g. -0.002 K/m
-!
-! output:   p_tropopause(klon)     tropopause pressure in Pa with missing values
-!           t_tropopause(klon)     tropopause temperature in K with missing values
-!           z_tropopause(klon)     tropopause height in m with missing values
-!           stratomask             stratospheric mask withtout missing values
-!           ifil                  # of undetermined values
+! committed to LMDz by O. Boucher (2016) with a mistake  
+! mistake corrected by O. Boucher (2017-12-11) 
+! 
+! input:  temp(nlon,nlat,nlev)  3D-temperature field  
+!         ps(nlon,nlat)         2D-surface pressure field 
+!         zs(nlon,nlat)         2D-surface height 
+!         nlon                  grid points in x 
+!         nlat                  grid points in y 
+!         pfull(nlon,nlat,nlev) full pressure levels in Pa 
+!         plimu                 upper limit for tropopause pressure 
+!         pliml                 lower limit for tropopause pressure  
+!         gamma                 tropopause criterion, e.g. -0.002 K/m 
+! 
+! output: p_tropopause(klon)    tropopause pressure in Pa with missing values 
+!         t_tropopause(klon)    tropopause temperature in K with missing values 
+!         z_tropopause(klon)    tropopause height in m with missing values 
+!         stratomask            stratospheric mask withtout missing values 
+!         ifil                  # of undetermined values
 !
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -104 +106 @@
 subroutine twmo(missing_val, level, t, p, ps, zs, plimu, pliml, gamma, ptrp, ttrp, ztrp)
 
+! reference: Reichler, T., M. Dameris, and R. Sausen (GRL, 10.1029/2003GL018240, 2003)
+
 implicit none
 
@@ -128 +132 @@
 ztrp=missing_val
 
-faktor = -RG/R
+faktor = -RG/RD
 
 do j=level,2,-1
@@ -218 +222 @@
 TDLNP = TDLNP + TM*DLNP
 
-ZTRP = ZS + TDLNP*R/RG
-
-if (ZTRP .lt. 0) then
-  print*,ZTRP
-  print*,PS
-  print*,P
-  print*,T
-  print*,ZS
-  stop
-endif
+ZTRP = ZS + TDLNP*RD/RG
+
+!!if (ZTRP .lt. 0) then
+!!  print*,ZTRP
+!!  print*,PS
+!!  print*,P
+!!  print*,T
+!!  print*,ZS
+!!  stop
+!!endif
 
 return