Argument Definitions (+ indicates altered content)
Description SUBROUTINE CONVECT ** VERSION 4.3c ** 20 May, 2002 ** Kerry Emanuel ** ******************************************************************* C. Forster, November 2003 - May 2004: The subroutine has been downloaded from Kerry Emanuel's homepage, where further infos on the convection scheme can be found http://www-paoc.mit.edu/~emanuel/home.html The following changes have been made to integrate this subroutine into FLEXPART Putting most of the variables in a new common block renaming eps to eps0 because there is some eps already in includepa removing the arrays U,V,TRA and related arrays renaming the original arrays T,Q,QS,P,PH to TCONV,QCONV,QSCONV,PCONV_HPA,PHCONV_HPA Initialization of variables has been put into parameter statements of assignment of values at each call, in order to save computation *** On input: *** T: Array of absolute temperature (K) of dimension ND, with first index corresponding to lowest model level. Note that this ar will be altered by the subroutine if dry convective adjustme occurs and if IPBL is not equal to 0. Q: Array of specific humidity (gm/gm) of dimension ND, with firs index corresponding to lowest model level. Must be defined at same grid levels as T. Note that this array will be alte if dry convective adjustment occurs and if IPBL is not equa QS: Array of saturation specific humidity of dimension ND, with f index corresponding to lowest model level. Must be defined at same grid levels as T. Note that this array will be alte if dry convective adjustment occurs and if IPBL is not equa U: Array of zonal wind velocity (m/s) of dimension ND, witth fir index corresponding with the lowest model level. Defined at same levels as T. Note that this array will be altered if dry convective adjustment occurs and if IPBL is not equal t V: Same as U but for meridional velocity. TRA: Array of passive tracer mixing ratio, of dimensions (ND,NTRA) where NTRA is the number of different tracers. If no convective tracer transport is needed, define a dummy input array of dimension (ND,1). Tracers are defined at same vertical levels as T. Note that this array will be alt if dry convective adjustment occurs and if IPBL is not equa P: Array of pressure (mb) of dimension ND, with first index corresponding to lowest model level. Must be defined at same grid levels as T. PH: Array of pressure (mb) of dimension ND+1, with first index corresponding to lowest level. These pressures are defined levels intermediate between those of P, T, Q and QS. The fi value of PH should be greater than (i.e. at a lower level t the first value of the array P. ND: The dimension of the arrays T,Q,QS,P,PH,FT and FQ NL: The maximum number of levels to which convection can penetrate, plus 1. NL MUST be less than or equal to ND-1. NTRA:The number of different tracers. If no tracer transport is needed, set this equal to 1. (On most compilers, setting NTRA to 0 will bypass tracer calculation, saving some CPU.) DELT: The model time step (sec) between calls to CONVECT *** On Output: *** IFLAG: An output integer whose value denotes the following: VALUE INTERPRETATION ----- -------------- 0 No moist convection; atmosphere is no unstable, or surface temperature is l than 250 K or surface specific humidi is non-positive. 1 Moist convection occurs. 2 No moist convection: lifted condensat level is above the 200 mb level. 3 No moist convection: cloud base is hi then the level NL-1. 4 Moist convection occurs, but a CFL co on the subsidence warming is violated does not cause the scheme to terminat FT: Array of temperature tendency (K/s) of dimension ND, defined grid levels as T, Q, QS and P. FQ: Array of specific humidity tendencies ((gm/gm)/s) of dimensi defined at same grid levels as T, Q, QS and P. FU: Array of forcing of zonal velocity (m/s^2) of dimension ND, defined at same grid levels as T. FV: Same as FU, but for forcing of meridional velocity. FTRA: Array of forcing of tracer content, in tracer mixing ratio p second, defined at same levels as T. Dimensioned (ND,NTRA) PRECIP: Scalar convective precipitation rate (mm/day). WD: A convective downdraft velocity scale. For use in surface flux parameterizations. See convect.ps file for details. TPRIME: A convective downdraft temperature perturbation scale (K). For use in surface flux parameterizations. See convect.ps file for details. QPRIME: A convective downdraft specific humidity perturbation scale (gm/gm). For use in surface flux parameterizations. See convect.ps file for details. CBMF: The cloud base mass flux ((kg/m**2)/s). THIS SCALAR VALUE BE STORED BY THE CALLING PROGRAM AND RETURNED TO CONVECT ITS NEXT CALL. That is, the value of CBMF must be "rememb by the calling program between calls to CONVECT. *** THE PARAMETER NA SHOULD IN GENERAL BE GREATER THAN *** *** OR EQUAL TO ND + 1 ***
Source file: convect43c.f
Intrinsic Functions Called
External Functions and Subroutines Called
Parameter Variables Used
Local Variables (+ indicates altered content)
Referenced Common Block Variables (+ indicates altered content)
Argument Definitions (+ indicates altered content)
Description
Source file: convect43c.f
Intrinsic Functions Called
Parameter Variables Used
Local Variables (+ indicates altered content)
Referenced Common Block Variables (+ indicates altered content)