CALWXT_BOURG.f File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	calwxt_bourg_post (im, ista_2l, iend_2u, ista, iend, jm, jsta_2l, jend_2u, jsta, jend, lm, lp1, iseed, g, pthresh, t, q, pmid, pint, lmh, prec, zint, ptype, me)
	calwxt_bourg_post Subroutine that calculates precipitation type (Bourgouin).

Detailed Description

Subroutine that calculates 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
[in]	me	integer Identifier for the processor used in the current instance.

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
2021-10-31	J Meng	2D DECOMPOSITION

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.

Function/Subroutine Documentation

calwxt_bourg_post()

subroutine calwxt_bourg_post	(	integer, intent(in)	im,
		integer, intent(in)	ista_2l,
		integer, intent(in)	iend_2u,
		integer, intent(in)	ista,
		integer, intent(in)	iend,
		integer, intent(in)	jm,
		integer, intent(in)	jsta_2l,
		integer, intent(in)	jend_2u,
		integer, intent(in)	jsta,
		integer, intent(in)	jend,
		integer, intent(in)	lm,
		integer, intent(in)	lp1,
		integer, intent(in)	iseed,
		real, intent(in)	g,
		real, intent(in)	pthresh,
		real, dimension(ista_2l:iend_2u,jsta_2l:jend_2u,lm), intent(in)	t,
		real, dimension(ista_2l:iend_2u,jsta_2l:jend_2u,lm), intent(in)	q,
		real, dimension(ista_2l:iend_2u,jsta_2l:jend_2u,lm), intent(in)	pmid,
		real, dimension(ista_2l:iend_2u,jsta_2l:jend_2u,lp1), intent(in)	pint,
		real, dimension(ista_2l:iend_2u,jsta_2l:jend_2u), intent(in)	lmh,
		real, dimension(ista_2l:iend_2u,jsta_2l:jend_2u), intent(in)	prec,
		real, dimension(ista_2l:iend_2u,jsta_2l:jend_2u,lp1), intent(in)	zint,
		integer, dimension(ista:iend,jsta:jend), intent(out)	ptype,
		integer, intent(in)	me
	)

calwxt_bourg_post Subroutine that calculates precipitation type (Bourgouin).

Parameters

[in]	im	integer i dimension.
[in]	ista_2l	integer i dimension start point (including haloes).
[in]	iend_2u	integer i dimension end point (including haloes).
[in]	ista	integer i dimension start point (excluding haloes).
[in]	iend	integer i dimension end point (excluding haloes).
[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.
[in]	me	integer Identifier for the processor used in the current instance.

Definition at line 82 of file CALWXT_BOURG.f.

Referenced by surfce().