Index | Rundeck: P1SoM40 | Created: Thu May 9 03:30:17 EDT 2024 |
tropwmo | Module: | File: ATM_UTILS.f |
Author : D. Nodorp/T. Reichler/C. Land GISS Modifications by Jean Lerner/Gavin Schmidt
Version:
Algorithm:
WMO Tropopause Definition From A Temperature Lapse Rate Definition of the Tropopause Based on Ozone, J. M. Roe and W. H. Jasperson, 1981 In the following discussion the lapse rate is defined as -dT/dz. The main features of the WMO tropopause definition are as follows: * The first tropopause (i.e., the conventional tropopause) is defined as the lowest level at which the lapse rate decreases to 2 K/km or less, and the average lapse rate from this level to any level within the next higher 2 km does not exceed 2 K/km. * If above the first tropopause the average lapse rate between any level and all higher levels within 1 km exceed 3 K/km, then a second tropopause is defined by the same criterion as under the statement above. This tropopause may be either within or above the 1 km layer. * A level otherwise satisfying the definition of tropopause, but occuring at an altitude below that of the 500 mb level will not be designated a tropopause unless it is the only level satisfying the definition and the average lapse rate fails to exceed 3 K/km over at least 1 km in any higher layer. * (GISS failsafe) Some cases occur when the lapse rate never falls below 2 K/km. In such cases the failsafe level is that where the lapse rate first falls below 3 K/km. If this still doesn't work (ever?), the level is set to the pressure level below 30mb.
Declaration:
subroutine tropwmo(ptm1, papm1, pk, ptropo, ltropp, ierr)
- klev : used from resolution
- Initial Value
= >lm
- zkappa : used from constant
- Initial Value
= >kapa
- zzkap : used from constant
- Initial Value
= >bykapa
- grav : used from constant
- rgas : used from constant
- ptm1 :
real*8, intent(in), dimension(klev)
- papm1 :
real*8, intent(in), dimension(klev)
- pk :
real*8, intent(in), dimension(klev)
- ptropo :
real*8, intent(out)
- ltropp :
integer, intent(out)
- ierr :
integer, intent(out)
- zgwmo :
real*8, parameter
- min lapse rate (* -1) needed for trop. defn. (-K/km)
Initial Value= -2d-3
- zgwmo2 :
real*8, parameter
- GISS failsafe minimum lapse rate (* -1) (-K/km)
Initial Value= -3d-3
- zdeltaz :
real*8, parameter
- distance to check for lapse rate changes (km)
Initial Value= 2000.0
- zfaktor :
real*8, parameter
- factor for caluclating height from pressure (-rgas/grav)
Initial Value= -GRAV/RGAS
- zplimb :
real*8, parameter
- min pressure at which to define tropopause (mb)
Initial Value= 500.
- zpmk :
real*8, dimension(klev)
- zpm :
real*8, dimension(klev)
- za :
real*8, dimension(klev)
- zb :
real*8, dimension(klev)
- ztm :
real*8, dimension(klev)
- zdtdz :
real*8, dimension(klev)
- zptph :
real*8
- zp2km :
real*8
- zag :
real*8
- zbg :
real*8
- zasum :
real*8
- zaquer :
real*8
- zptf :
real*8
- iplimb :
integer
- iplimt :
integer
- jk :
integer
- jj :
integer
- kcount :
integer
- ltset :
integer
- l :
integer
- ldtdz :
logical