CALWXT_BOURG.f File Reference

Subroutine that calculate precipitation type (Bourgouin). More...

Go to the source code of this file.

Functions/Subroutines
subroutine	calwxt_bourg_post (im, jm, jsta_2l, jend_2u, jsta, jend, lm, lp1, iseed, g, pthresh, t, q, pmid, pint, lmh, prec, zint, ptype, me)

Detailed Description

Subroutine that calculate precipitation type (Bourgouin).

This routine computes precipitation type. using a decision tree approach that uses the so-called "energy method" of Bourgouin of AES (Canada) 1992.

Parameters

[in]	im	integer i dimension.
[in]	jm	integer j dimension.
[in]	jsta_2l	integer j dimension start point (including haloes).
[in]	jend_2u	integer j dimension end point (including haloes).
[in]	jsta	integer j dimension start point (excluding haloes).
[in]	jend	integer j dimension end point (excluding haloes).
[in]	lm	integer k dimension.
[in]	lp1	integer k dimension plus 1.
[in]	iseed	integer random number seed.
[in]	g	real gravity (m/s**2).
[in]	pthresh	real precipitation threshold (m).
[in]	t	real(im,jsta_2l:jend_2u,lm) mid layer temp (K).
[in]	q	real(im,jsta_2l:jend_2u,lm) specific humidity (kg/kg).
[in]	pmid	real(im,jsta_2l:jend_2u,lm) mid layer pressure (Pa).
[in]	pint	real(im,jsta_2l:jend_2u,lp1) interface pressure (Pa).
[in]	lmh	real(im,jsta_2l:jend_2u) max number of layers.
[in]	prec	real(im,jsta_2l:jend_2u) precipitation (m).
[in]	zint	real(im,jsta_2l:jend_2u,lp1) interface height (m).
[out]	ptype	integer(im,jm) instantaneous weather type () acts like a 4 bit binary 1111 = rain/freezing rain/ice pellets/snow. where the one's digit is for snow the two's digit is for ice pellets the four's digit is for freezing rain and the eight's digit is for rain in other words... ptype=1 snow ptype=2 ice pellets/mix with ice pellets ptype=4 freezing rain/mix with freezing rain ptype=8 rain

Program history log:

Date	Programmer	Comments
1999-07-06	M Baldwin	Initial
1999-09-20	M Baldwin	make more consistent with bourgouin (1992)
2005-08-24	G Manikin	added to wrf post
2007-06-19	M Iredell	mersenne twister, best practices
2015-??-??	S Moorthi	changed random number call and optimization and cleanup

Remarks: vertical order of arrays must be layer 1 = top and layer lmh = bottom

Author

M Baldwin np22

Date

1999-07-06

Definition in file CALWXT_BOURG.f.