UPP (develop)
|
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). | |
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.
[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. |
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
Definition in file CALWXT_BOURG.f.
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).
[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().