TRPAUS_NAM.f File Reference

trpaus() computes tropopause data. More...

Go to the source code of this file.

Functions/Subroutines
subroutine	trpaus (ptrop, ttrop, ztrop, utrop, vtrop, shtrop)
	This routine computes tropopause data.

Detailed Description

trpaus() computes tropopause data.

Author

Russ Treadon W/NP2

Date

1992-12-22

Definition in file TRPAUS_NAM.f.

Function/Subroutine Documentation

trpaus()

subroutine trpaus	(	real, dimension(ista:iend,jsta:jend)	ptrop,
		real, dimension(ista:iend,jsta:jend)	ttrop,
		real, dimension(i real vtrop(ista:iend,jsta:jend),shtrop(ista:iend,jsta:jend) real tlapse(lm),dz2(lm),delt2(lm),tlapse2(lm) ! integer i,j real pm,delt,dz,rsqdif ! !***************************************************************************** ! start trpaus here. ! ! loop over the horizontal grid. ! do j=jsta,jend do i=ista,iend ptrop(i,j) = spval ttrop(i,j) = spval ztrop(i,j) = spval utrop(i,j) = spval vtrop(i,j) = spval shtrop(i,j) = spval enddo enddo ! !$omp parallel do !$omp& private(delt,delt2,dz,dz2,ie,iw,l,llmh,pm,rsqdif, !$omp& tlapse,tlapse2,u0,u0l,uh,uh0,ul, !$omp& v0,v0l,vh,vh0) do j=jsta,jend loopi:do i=ista,iend ! ! compute the temperature lapse rate (-dt/dz) between eta ! layers moving up from the ground. the first eta layer ! above pressure "pstart" in which the lapse rate is less ! than the critcal lapse rate is labeled the tropopause. ! llmh=nint(lmh(i,j)) ! loopl: do l=llmh-1,2,-1 pm = pint(i,j,l) delt = t(i,j,l-1)-t(i,j,l) dz = d50*(zint(i,j,l-1)-zint(i,j,l+1)) tlapse(l) = -delt/dz ! if ((tlapse(l)<crtlap).and.(pm<pstart))	ztrop,
			utrop,
			vtrop,
			shtrop
	)

This routine computes tropopause data.

At each mass point a surface up search is made for the first occurrence of a three layer mean lapse rate less than or equal to a critical lapse rate. This critcal lapse rate is 2deg/km. This is in accord with the WMO definition of a tropopause. A maximum tropopause pressure of 500mb is enforced. Onc the tropopause is located in a column, pressure, temperature, u and v winds, and vertical wind shear are computed.

Parameters

[out]	PTROP	Tropopause pressure.
[out]	TTROP	Tropopause temperature.
[out]	ZTROP	Tropopause height.
[out]	UTROP	Tropopause u wind component.
[out]	VTROP	Tropopause v wind component.
[out]	SHTROP	Vertical wind shear at tropopause.

Program history log:

Date	Programmer	Comments
1992-12-22	Russ Treadon	Initial
1997-03-06	Geoff Manikin	Changed criteria for determining the tropopause and added height
1998-06-15	T Black	Conversion from 1-D TO 2-D
2000-01-04	Jim Tuccillo	MPI Version
2002-04-23	Mike Baldwin	WRF Version
2019-10-30	Bo Cui	ReMOVE "GOTO" STATEMENT
2021-09-13	JESSE MENG	2D DECOMPOSITION

Author

Russ Treadon W/NP2

Date

1992-12-22

Definition at line 35 of file TRPAUS_NAM.f.