*IF DEF,A05_2A AJX1F405.122
C ******************************COPYRIGHT****************************** GTS2F400.1189
C (c) CROWN COPYRIGHT 1995, METEOROLOGICAL OFFICE, All Rights Reserved. GTS2F400.1190
C GTS2F400.1191
C Use, duplication or disclosure of this code is subject to the GTS2F400.1192
C restrictions as set forth in the contract. GTS2F400.1193
C GTS2F400.1194
C Meteorological Office GTS2F400.1195
C London Road GTS2F400.1196
C BRACKNELL GTS2F400.1197
C Berkshire UK GTS2F400.1198
C RG12 2SZ GTS2F400.1199
C GTS2F400.1200
C If no contract has been raised with this copy of the code, the use, GTS2F400.1201
C duplication or disclosure of it is strictly prohibited. Permission GTS2F400.1202
C to do so must first be obtained in writing from the Head of Numerical GTS2F400.1203
C Modelling at the above address. GTS2F400.1204
C ******************************COPYRIGHT****************************** GTS2F400.1205
C GTS2F400.1206
CLL SUBROUTINE CONVEC2------------------------------------------------ CONV21A.3
CLL CONV21A.4
CLL PURPOSE : COMPLETES LIFTING OF THE PARCEL FROM LAYER K TO K+1 CONV21A.5
CLL CONV21A.6
CLL CALL SUBROUTINE PARCEL AND ENVIRON CONV21A.7
CLL CONV21A.8
CLL SUBROUTINE PARCEL CALCULATES AN INITIAL MASS FLUX, CONV21A.9
CLL CARRIES OUT THE DETRAINMENT CALCULATION, TESTS CONV21A.10
CLL TO SEE IF CONVECTION IS TERMINATING AND CALCULATES THE CONV21A.11
CLL PRECIPITATION RATE FROM LAYER K+1 CONV21A.12
CLL CONV21A.13
CLL SUBROUTINE ENVIRON CALCULATES THE EFFECT OF CONVECTION CONV21A.14
CLL UPON THE LARGE-SCALE ATMOSPHERE CONV21A.15
CLL CONV21A.16
CLL SUITABLE FOR SINGLE COLUMN MODEL USE CONV21A.17
CLL CONV21A.18
CLL CODE REWORKED FOR CRAY Y-MP BY D.GREGORY AUTUMN/WINTER 1989/90 CONV21A.19
CLL CONV21A.20
CLL MODEL MODIFICATION HISTORY FROM MODEL VERSION 3.0: CONV21A.21
CLL VERSION DATE CONV21A.22
CLL 3.2 8/07/93 : added convective cloud condensed water diagnostic PI080793.42
CLL : P Inness PI080793.43
CLL CONV21A.23
CLL 3.3 23/12/93 : DG060893 : CORRECTION TO REDUCE OVER PREDICTION DG060893.1
CLL OF CONVECTIVE SNOW; ALLOW MASS FLUX DG060893.2
CLL IN LAYER K-1 WHEN CLOUD UNDERGOING DG060893.3
CLL TERMINAL DETRAINMENT TO BE PASSED DG060893.4
CLL BACK TO CONVECT CORRECTLY. DG060893.5
CLL 3.4 21/03/94 Add lowest conv.cloud diagnostics. R.T.H.Barnes. ARN2F304.76
CLL DG060893.6
CLL 4.0 5/05/95 Added calculation of CAPE for diagnostic output. API2F400.17
CLL Pete Inness. API2F400.18
!LL 4.5 19/5/98 Correction of CAPE diagnostic. Julie Gregory AJX2F405.9
CLL 4.5 22/7/98 Kill the IBM specific lines (JCThil) AJC1F405.2
CLL PROGRAMMING STANDARDS : UNIFIED MODEL DOCUMENTATION PAPER NO. 4 CONV21A.24
CLL VERSION NO. 1 CONV21A.25
CLL CONV21A.26
CLL LOGICAL COMPONENT NUMBER: P27 CONV21A.27
CLL CONV21A.28
CLL SYSTEM TASK : P27 CONV21A.29
CLL CONV21A.30
CLL DOCUMENTATION : UNIFIED MODEL DOCUMENTATION PAPER P27 CONV21A.31
CLL CONV21A.32
CLLEND----------------------------------------------------------------- CONV21A.33
C CONV21A.34
C*L ARGUMENTS--------------------------------------------------------- CONV21A.35
C CONV21A.36
SUBROUTINE CONVEC2 (NPNTS,NLEV,K,THEK,THEKP1,QEK,QEKP1,QSEKP1, 4,6CONV21A.37
* DQSKP1,PSTAR,THPK,QPK,THPKP1,QPKP1,XSQKP1, CONV21A.38
* RBUOY,QSEK,DQSK,THPI,QPI,XPK,FLXK,BWKP1,BGMKP1, CONV21A.39
* BGMK,BLOWST,BLAND,BTERM,DEPTH,PREKP1,DTHEK, CONV21A.40
* DQEK,DTHEKP1,DQEKP1,BINIT,CCA,ICCB,ICCT, CONV21A.41
* TCW,EKP14,EKP34,AMDETK,PK,PKP1, CONV21A.42
* EXK,EXKP1,DELEXKP1,DELPK,DELPKP1, CONV21A.43
* CCLWP,CCW,LCCA,LCBASE,LCTOP,LCCLWP, API2F400.19
* CAPE) API2F400.20
C CONV21A.45
IMPLICIT NONE CONV21A.46
C API2F400.21
C---------------------------------------------------------------------- API2F400.22
C MODEL CONSTANTS API2F400.23
C---------------------------------------------------------------------- API2F400.24
C API2F400.25
*CALL C_R_CP 
API2F400.26
*CALL C_EPSLON API2F400.27
C API2F400.28
C CONV21A.47
C---------------------------------------------------------------------- CONV21A.48
C VECTOR LENGTHS AND LOOP COUNTERS CONV21A.49
C---------------------------------------------------------------------- CONV21A.50
C CONV21A.51
INTEGER NPNTS ! IN VECTOR LENGTH CONV21A.55
C CONV21A.56
INTEGER NLEV ! IN NUMBER OF MODEL LAYERS CONV21A.57
C CONV21A.58
INTEGER I ! LOOP COUNTER CONV21A.59
C CONV21A.60
INTEGER K ! PRESENT MODEL LAYER CONV21A.61
C CONV21A.62
C CONV21A.63
C---------------------------------------------------------------------- CONV21A.64
C VARIABLES WHICH ARE INPUT CONV21A.65
C---------------------------------------------------------------------- CONV21A.66
C CONV21A.67
REAL THEK(NPNTS) ! IN POTENTIAL TEMPERATURE OF CLOUD CONV21A.68
! ENVIRONMENT IN LAYER K (K) CONV21A.69
C CONV21A.70
REAL THEKP1(NPNTS) ! IN POTENTIAL TEMPERATURE OF CLOUD CONV21A.71
! ENVIRONMENT IN LAYER K+1 (K) CONV21A.72
C CONV21A.73
REAL QEK(NPNTS) ! IN MIXING RATIO OF CLOUD CONV21A.74
! ENVIRONMENT IN LAYER K (KG/KG) CONV21A.75
C CONV21A.76
REAL QEKP1(NPNTS) ! IN MIXING RATIO OF CLOUD CONV21A.77
! ENVIRONMENT IN LAYER K+1 (KG/KG) CONV21A.78
C CONV21A.79
REAL QSEKP1(NPNTS) ! IN SATURATION MIXING RATIO OF CLOUD CONV21A.80
! ENVIRONMENT IN LAYER K+1 (KG/KG) CONV21A.81
C CONV21A.82
REAL DQSKP1(NPNTS) ! IN GRADIENT OF SATURATION MIXING RATIO CONV21A.83
! WITH POTENTIAL TEMPERATURE FOR THE CONV21A.84
! CLOUD ENVIRONMENT IN LAYER K+1 CONV21A.85
! (KG/KG) CONV21A.86
C CONV21A.87
REAL PSTAR(NPNTS) ! IN SURFACE PRESSURE (PA) CONV21A.88
C CONV21A.89
REAL THPKP1(NPNTS) ! IN PARCEL POTENTIAL TEMPERATURE CONV21A.90
! IN LAYER K (K) CONV21A.91
C CONV21A.92
REAL QPKP1(NPNTS) ! IN PARCEL MIXING RATIO IN LAYER K+1 CONV21A.93
! (KG/KG) CONV21A.94
C CONV21A.95
REAL XSQKP1(NPNTS) ! IN EXCESS WATER IN PARCEL AFTER LIFTING CONV21A.96
! LAYER K TO K+1 (KG/KG) CONV21A.97
C CONV21A.98
REAL RBUOY(NPNTS) ! IN PARCEL BUOYANCY IN LAYER K+1 (K) CONV21A.99
C CONV21A.100
REAL QSEK(NPNTS) ! IN SATURATION MIXING RATIO OF CLOUD CONV21A.101
! ENVIRONMENT IN LAYER K (KG/KG) CONV21A.102
C CONV21A.103
REAL DQSK(NPNTS) ! IN GRADIENT OF SATURATION MIXING RATIO CONV21A.104
! WITH POTENTIAL TEMPERATURE FOR THE CONV21A.105
! CLOUD ENVIRONMENT OF LAYER K CONV21A.106
! (KG/KG/K) CONV21A.107
C CONV21A.108
REAL THPI(NPNTS) ! IN INITIAL PARCEL POTENTIAL TEMPERATURE CONV21A.109
! (K) CONV21A.110
C CONV21A.111
REAL QPI(NPNTS) ! IN INITIAL PARCEL MIXING RATIO CONV21A.112
! (KG/KG) CONV21A.113
C CONV21A.114
LOGICAL BWKP1(NPNTS) ! IN MASK FOR WHETHER CONDENSATE IS CONV21A.115
! LIQUID IN LAYER K+1 CONV21A.116
C CONV21A.117
LOGICAL BGMKP1(NPNTS) ! IN MASK FOR PARCELS WHICH ARE CONV21A.118
! SATURATED IN LAYER K+1 CONV21A.119
C CONV21A.120
LOGICAL BLAND(NPNTS) ! IN LAND/SEA MASK CONV21A.121
C CONV21A.122
LOGICAL BINIT(NPNTS) ! IN MASK FOR THOSE POINTS AT WHICH CONV21A.123
! CONVECTION IS OCCURING CONV21A.124
C CONV21A.125
LOGICAL BLOWST(NPNTS) ! IN MASK FOR THOSE POINTS AT WHICH CONV21A.126
! STABILITY IS LOW ENOUGH FOR CONV21A.127
! CONVECTION TO OCCUR CONV21A.128
C CONV21A.129
REAL EKP14(NPNTS) ! IN ENTRAINMENT COEFFICIENT AT LEVEL CONV21A.130
! K+1/4 MULTIPLIED BY APPROPRIATE CONV21A.131
! LAYER THICKNESS CONV21A.132
C CONV21A.133
REAL EKP34(NPNTS) ! IN ENTRAINMENT COEFFICIENT AT LEVEL CONV21A.134
! K+1/4 MULTIPLIED BY APPROPRIATE CONV21A.135
! LAYER THICKNESS CONV21A.136
C CONV21A.137
REAL AMDETK(NPNTS) ! IN MIXING DETRAINMENT COEFFICIENT CONV21A.138
! AT LEVEL K MULTIPLIED BY APPROPRIATE CONV21A.139
! LAYER THICKNESS CONV21A.140
C CONV21A.141
REAL DELPKP12(NPNTS) ! IN PRESSURE DIFFERENCE BETWEEN CONV21A.142
! MID-POINTS OF LAYERS K AND K+1 CONV21A.143
! (PA) CONV21A.144
C CONV21A.145
REAL PK(NPNTS) ! IN PRESSURE AT MID-POINT OF LAYER K CONV21A.146
! (PA) CONV21A.147
C CONV21A.148
REAL PKP1(NPNTS) ! IN PRESSURE AT MID-POINT OF LAYER K+1 CONV21A.149
! (PA) CONV21A.150
C CONV21A.151
REAL EXK(NPNTS) ! IN EXNER RATIO AT MID-POINT OF LAYER K CONV21A.152
C CONV21A.153
REAL EXKP1(NPNTS) ! IN EXNER RATIO AT MID-POINT OF LAYER K+1 CONV21A.154
C CONV21A.155
REAL DELEXKP1(NPNTS) ! IN DIFFERENCE IN EXNER RATIO BETWEEN CONV21A.156
! MID-POINTS OF LAYERS K AND K+1 CONV21A.157
C CONV21A.158
REAL DELPK(NPNTS) ! IN DIFFERENCE IN PRESSURE ACROSS LAYER K CONV21A.159
! (PA) CONV21A.160
C CONV21A.161
REAL DELPKP1(NPNTS) ! IN DIFFERENCE IN PRESSURE ACROSS CONV21A.162
! LAYER K+1 (PA) CONV21A.163
C CONV21A.164
C CONV21A.165
C---------------------------------------------------------------------- CONV21A.166
C VARIABLES WHICH ARE INPUT BUT WHICH ARE ALSO UPDATED IN THIS ROUTINE CONV21A.167
C---------------------------------------------------------------------- CONV21A.168
C CONV21A.169
REAL THPK(NPNTS) ! INOUT CONV21A.170
! IN PARCEL POTENTIAL TEMPERATURE CONV21A.171
! IN LAYER K (K) CONV21A.172
! OUT PARCEL POTENTIAL TEMPERATURE CONV21A.173
! IN LAYER K+1 (K) CONV21A.174
C CONV21A.175
REAL QPK(NPNTS) ! INOUT CONV21A.176
! IN PARCEL MIXING RATIO IN LAYER K CONV21A.177
! (KG/KG) CONV21A.178
! OUT PARCEL MIXING RATIO IN LAYER K+1 CONV21A.179
! (KG/KG) CONV21A.180
C CONV21A.181
REAL XPK(NPNTS) ! INOUT CONV21A.182
! IN PARCEL CLOUD WATER IN LAYER K CONV21A.183
! (KG/KG) CONV21A.184
! OUT PARCEL CLOUD WATER IN LAYER K+1 CONV21A.185
! (KG/KG) CONV21A.186
C CONV21A.187
REAL FLXK(NPNTS) ! INOUT CONV21A.188
! IN PARCEL MASSFLUX IN LAYER K (PA/S) CONV21A.189
! OUT PARCEL MASSFLUX IN LAYER K+1 (PA/S) CONV21A.190
C CONV21A.191
LOGICAL BGMK(NPNTS) ! INOUT CONV21A.192
! IN MASK FOR PARCELS WHICH ARE CONV21A.193
! SATURATED IN LAYER K CONV21A.194
! OUT MASK FOR PARCELS WHICH ARE CONV21A.195
! SATURATED IN LAYER K+1 CONV21A.196
C CONV21A.197
REAL DTHEK(NPNTS) ! INOUT CONV21A.198
! IN INCREMENT TO MODEL POTENTIAL CONV21A.199
! TEMPERATURE IN LAYER K DUE TO CONV21A.200
! CONVECTION (MAY BE NONE ZERO CONV21A.201
! DUE TO A PREVIOUS SPLIT FINAL CONV21A.202
! DETRAINEMNT CALCULATION) (K/S) CONV21A.203
! OUT UPDATED INCREMENT TO MODEL POTENTIAL CONV21A.204
! TEMPERATURE IN LAYER K DUE TO CONV21A.205
! CONVECTION (K/S) CONV21A.206
C CONV21A.207
REAL DQEK(NPNTS) ! INOUT CONV21A.208
! IN INCREMENT TO MODEL MIXING RATIO CONV21A.209
! IN LAYER K DUE TO CONVECTION CONV21A.210
! (MAY BE NONE ZERO DUE TO A CONV21A.211
! PREVIOUS SPLIT FINAL DETRAINEMNT CONV21A.212
! CALCULATION) (KG/KG/S) CONV21A.213
! OUT UPDATED INCREMENT TO MODEL MIXING CONV21A.214
! RATIO IN LAYER K DUE TO CONVECTION CONV21A.215
! (KG/KG/S) CONV21A.216
C CONV21A.217
REAL TCW(NPNTS) ! INOUT CONV21A.218
! IN TOTAL CONDENSED WATER SUMMED TO CONV21A.219
! LAYER K (KG/M**2/S) CONV21A.220
! OUT UPDATED TOTAL CONDENSED WATER CONV21A.221
! SUMMED TO LAYER K+1 (KG/M**2/S) CONV21A.222
C CONV21A.223
REAL DEPTH(NPNTS) ! INOUT CONV21A.224
! IN DEPTH OF CONVECTIVE CLOUD TO CONV21A.225
! LAYER K (M) CONV21A.226
! OUT UPDATED DEPTH OF CONVECTIVE CLOUD CONV21A.227
! TO LAYER K+1 (M) CONV21A.228
C CONV21A.229
REAL CCLWP(NPNTS) ! INOUT CONV21A.230
! IN CONDENSED WATER PATH SUMMED TO CONV21A.231
! LAYER K (KG/M**2) CONV21A.232
! OUT UPDATED CONDENSED WATER PATH CONV21A.233
! SUMMED TO LAYER K+1 (KG/M**2) CONV21A.234
C API2F400.29
REAL CAPE(NPNTS) ! IN CONVECTIVE AVAILABLE POTENTIAL API2F400.30
! ENERGY UP TO THE CURRENT CONVECTING API2F400.31
! LAYER API2F400.32
! OUT CONVECTIVE AVAILABLE POTENTIAL API2F400.33
! ENERGY INCLUDING ADDITION DUE TO API2F400.34
! THE CAPE WITHIN THE CURRENT LAYER API2F400.35
C API2F400.36
C CONV21A.235
C CONV21A.236
C---------------------------------------------------------------------- CONV21A.237
C VARIABLES WHICH ARE OUTPUT CONV21A.238
C---------------------------------------------------------------------- CONV21A.239
C CONV21A.240
LOGICAL BTERM(NPNTS) ! OUT MASK FOR PARCELS WHICH TERMINATE IN CONV21A.241
! LAYER K+1 CONV21A.242
C CONV21A.243
REAL PREKP1(NPNTS) ! OUT PRECIPITATION FROM PARCEL AS IT CONV21A.244
! RISES FROM LAYER K TO K+1 (KG/M**2/S) CONV21A.245
C CONV21A.246
REAL DTHEKP1(NPNTS) ! OUT INCREMENT TO MODEL POTENTIAL CONV21A.247
! TEMPERATURE IN LAYER K+1 DUE TO CONV21A.248
! CONVECTION (K/S) CONV21A.249
C CONV21A.250
REAL DQEKP1(NPNTS) ! OUT INCREMENT TO MODEL MIXING RATIO CONV21A.251
! IN LAYER K+1 DUE TO CONVECTION CONV21A.252
! (KG/KG/S) CONV21A.253
C CONV21A.254
REAL CCA(NPNTS) ! OUT CONVECTIVE CLOUD AMOUNT (%) CONV21A.255
C CONV21A.256
INTEGER ICCB(NPNTS) ! OUT CONVECTIVE CLOUD BASE LEVEL CONV21A.257
C CONV21A.258
INTEGER ICCT(NPNTS) ! OUT CONVECTIVE CLOUD TOP LEVEL CONV21A.259
C CONV21A.260
REAL CCW(NPNTS) ! OUT CONVECTIVE CLOUD WATER(G/KG) ON PI080793.45
! MODEL LEVELS PI080793.46
C CONV21A.261
REAL LCCA(NPNTS) ! OUT LOWEST CONV.CLOUD AMOUNT (%) ARN2F304.78
C PI080793.47
INTEGER LCBASE(NPNTS) ! OUT LOWEST CONV.CLOUD BASE LEVEL ARN2F304.79
C ARN2F304.80
INTEGER LCTOP(NPNTS) ! OUT LOWEST CONV.CLOUD TOP LEVEL ARN2F304.81
C ARN2F304.82
REAL LCCLWP(NPNTS) ! OUT LOWEST CONV.CLOUD LIQ.WATER PATH ARN2F304.83
C ARN2F304.84
C ARN2F304.85
C---------------------------------------------------------------------- CONV21A.262
C VARIABLES WHICH ARE DEFINED LOCALLY CONV21A.263
C CONV21A.264
REAL THRK(NPNTS) ! PARCEL DETRAINMENT POTENTIAL CONV21A.286
! TEMPERATURE IN LAYER K (K) CONV21A.287
C CONV21A.288
REAL QRK(NPNTS) ! PARCEL DETRAINMENT MIXING RATIO CONV21A.289
! IN LAYER K (KG/KG) CONV21A.290
C CONV21A.291
REAL XPKP1(NPNTS) ! PARCEL CLOUD WATER IN LAYER K+1 (KG/KG) CONV21A.292
C CONV21A.293
REAL FLXKP1(NPNTS ) ! PARCEL MASSFLUX IN LAYER K+1 (PA/S) CONV21A.294
C CONV21A.295
REAL DELTAK(NPNTS) ! PARCEL FORCED DETRAINMENT RATE CONV21A.296
! IN LAYER K MULTIPLIED BY APPROPRIATE CONV21A.297
! LAYER THICKNESS CONV21A.298
C CONV21A.299
REAL THVP,THVE,RHO ! VIRTUAL TEMPERATURE OF PARCEL, VIRTUAL API2F400.37
! TEMPERATURE OF ENVIRONMENT AND DENSITY, API2F400.38
! USED IN CALCULATION OF CAPE. API2F400.39
C CONV21A.301
LOGICAL L_SHALLOW ! \ DUMMY Variables used in parcel in AJX0F404.382
&, L_CCW ! | versions 3A and 3B. Added here to AJX0F404.383
PARAMETER (L_CCW=.FALSE.) ! / keep argument lists consistent. AJX0F404.384
! AJX0F404.385
C---------------------------------------------------------------------- CONV21A.302
C EXTERNAL ROUTINES CALLED CONV21A.303
C---------------------------------------------------------------------- CONV21A.304
C CONV21A.305
EXTERNAL PARCEL,ENVIRON CONV21A.306
C CONV21A.307
C*--------------------------------------------------------------------- CONV21A.308
CL CONV21A.309
CL---------------------------------------------------------------------- CONV21A.310
CL COMPLETE LIFTING PARCELS TO LAYER K+1 CONV21A.311
CL CONV21A.312
CL SUBROUTINE PARCEL CONV21A.313
CL CONV21A.314
CL UM DOCUMENTATION PAPER P27 CONV21A.315
CL SECTIONS (5),(6),(7),(8),(9) CONV21A.316
CL---------------------------------------------------------------------- CONV21A.317
CL CONV21A.318
CALL PARCEL (K,NPNTS,NLEV,PSTAR,THEKP1,THEK,QEKP1,QEK, CONV21A.319
* QSEK,QSEKP1,DQSK,DQSKP1,BLAND,BWKP1, CONV21A.320
* DELTAK,FLXK,THPK,QPK,THRK,QRK, CONV21A.321
* BTERM,THPKP1,QPKP1,PREKP1,XPK,XPKP1,FLXKP1, CONV21A.322
* XSQKP1,THPI,QPI,BGMK,BGMKP1,BLOWST,RBUOY, CONV21A.323
* CCA,ICCB,ICCT,TCW,DEPTH, CONV21A.324
* EKP14,EKP34,AMDETK,DELPKP1,PK,PKP1, CONV21A.325
* EXK,EXKP1,DELEXKP1,CCLWP,CCW, ARN2F304.86
& LCCA,LCBASE,LCTOP,LCCLWP,L_SHALLOW,L_CCW) AJX0F404.386
CL CONV21A.327
CL---------------------------------------------------------------------- CONV21A.328
CL CALCULATE THE EFFECT ON THE ENVIRONMENT (EXCEPT FOR THE CONV21A.329
CL THE EVAPORATION OF PRECIPITATION AND CHANGE OF PHASE) CONV21A.330
CL CONV21A.331
CL SUBROUTINE ENVIRON CONV21A.332
CL CONV21A.333
CL UM DOCUMENTATION PAPER P27 CONV21A.334
CL SECTION (10) CONV21A.335
CL---------------------------------------------------------------------- CONV21A.336
CL CONV21A.337
CALL ENVIRON (NPNTS,DTHEK,DQEK,DTHEKP1,DQEKP1, CONV21A.338
* THEK,QEK,DELTAK,FLXK, CONV21A.339
* THPK,QPK,THRK,QRK,THEKP1,QEKP1, CONV21A.340
* BTERM,THPKP1,QPKP1,XPK,XPKP1,BWKP1,FLXKP1, CONV21A.341
* BLOWST,EKP14,EXK,EXKP1,DELPK,DELPKP1, CONV21A.342
* AMDETK) CONV21A.343
C CONV21A.344
DO 10 I=1,NPNTS CONV21A.345
C CONV21A.346
C----------------------------------------------------------------------- CONV21A.347
C RESET BINIT WHERE CONVECTION HAS TERMINATED CONV21A.348
C----------------------------------------------------------------------- CONV21A.349
C CONV21A.350
BINIT(I) = .NOT.BTERM(I) CONV21A.351
10 CONTINUE CONV21A.352
C API2F400.40
CL--------------------------------------------------------------------- API2F400.41
CL CALCULATE CONTRIBUTION TO CAPE AND RATE OF CHANGE OF CAPE DUE TO API2F400.42
CL THE UPDRAUGHT API2F400.43
CL--------------------------------------------------------------------- API2F400.44
C API2F400.45
DO I=1,NPNTS API2F400.46
THVP=THPK(I)*(1.0+C_VIRTUAL*QPK(I)) API2F400.47
THVE=THEK(I)*(1.0+C_VIRTUAL*QEK(I)) API2F400.48
RHO=PK(I)/(R*THEK(I)*EXK(I)) API2F400.49
C API2F400.50
CAPE(I)=CAPE(I)+(THVP-THVE)*DELPK(I)/(RHO*THVE) AJX2F405.10
C API2F400.52
IF(BTERM(I))THEN API2F400.53
THVP=THPKP1(I)*(1.0+C_VIRTUAL*QPKP1(I)) API2F400.54
THVE=THEKP1(I)*(1.0+C_VIRTUAL*QEKP1(I)) API2F400.55
RHO=PKP1(I)/(R*THEKP1(I)*EXKP1(I)) API2F400.56
C API2F400.57
CAPE(I)=CAPE(I)+(THVP-THVE)*DELPKP1(I)/(RHO*THVE) AJX2F405.11
C API2F400.59
END IF API2F400.60
C API2F400.61
END DO API2F400.62
C API2F400.63
CL CONV21A.353
CL--------------------------------------------------------------------- CONV21A.354
CL SWAP PARCEL VALUES READY FOR THE NEXT PART OF ASCENT CONV21A.355
CL FROM LAYER K+1 TO K+2 CONV21A.356
CL--------------------------------------------------------------------- CONV21A.357
CL CONV21A.358
DO 30 I=1,NPNTS DG060893.7
THPK(I) = THPKP1(I) CONV21A.362
QPK(I) = QPKP1(I) CONV21A.363
XPK(I) = XPKP1(I) CONV21A.364
FLXK(I) = FLXKP1(I) CONV21A.365
BGMK(I) = BGMKP1(I) CONV21A.366
30 CONTINUE DG060893.8
C CONV21A.370
RETURN CONV21A.371
END CONV21A.372
*ENDIF CONV21A.373