*IF DEF,A05_2C CONV22C.2
C ******************************COPYRIGHT****************************** CONV22C.3
C (c) CROWN COPYRIGHT 1996, METEOROLOGICAL OFFICE, All Rights Reserved. CONV22C.4
C CONV22C.5
C Use, duplication or disclosure of this code is subject to the CONV22C.6
C restrictions as set forth in the contract. CONV22C.7
C CONV22C.8
C Meteorological Office CONV22C.9
C London Road CONV22C.10
C BRACKNELL CONV22C.11
C Berkshire UK CONV22C.12
C RG12 2SZ CONV22C.13
C CONV22C.14
C If no contract has been raised with this copy of the code, the use, CONV22C.15
C duplication or disclosure of it is strictly prohibited. Permission CONV22C.16
C to do so must first be obtained in writing from the Head of Numerical CONV22C.17
C Modelling at the above address. CONV22C.18
C ******************************COPYRIGHT****************************** CONV22C.19
C CONV22C.20
CLL SUBROUTINE CONVEC2------------------------------------------------ CONV22C.21
CLL CONV22C.22
CLL PURPOSE : COMPLETES LIFTING OF THE PARCEL FROM LAYER K TO K+1 CONV22C.23
CLL CONV22C.24
CLL CALL SUBROUTINE PARCEL AND ENVIRON CONV22C.25
CLL CONV22C.26
CLL SUBROUTINE PARCEL CALCULATES AN INITIAL MASS FLUX, CONV22C.27
CLL CARRIES OUT THE DETRAINMENT CALCULATION, TESTS CONV22C.28
CLL TO SEE IF CONVECTION IS TERMINATING AND CALCULATES THE CONV22C.29
CLL PRECIPITATION RATE FROM LAYER K+1 CONV22C.30
CLL CONV22C.31
CLL SUBROUTINE ENVIRON CALCULATES THE EFFECT OF CONVECTION CONV22C.32
CLL UPON THE LARGE-SCALE ATMOSPHERE CONV22C.33
CLL CONV22C.34
CLL SUITABLE FOR SINGLE COLUMN MODEL USE CONV22C.35
CLL CONV22C.36
CLL MODEL MODIFICATION HISTORY: CONV22C.37
CLL VERSION DATE CONV22C.38
CLL 4.2 1/11/96 New deck version based on CONV22A with HADCM2 CONV22C.39
CLL specific modifications: R Jones CONV22C.40
!LL 4.5 19/5/98 Correction of CAPE diagnostic. Julie Gregory AJX2F405.6
!LL 4.5 23/02/98 Call comdecks. D. Robinson ADR1F405.41
CLL CONV22C.41
CLL PROGRAMMING STANDARDS : UNIFIED MODEL DOCUMENTATION PAPER NO. 3 CONV22C.42
CLL CONV22C.43
CLL LOGICAL COMPONENT NUMBER: P27 CONV22C.44
CLL CONV22C.45
CLL SYSTEM TASK : P27 CONV22C.46
CLL CONV22C.47
CLL DOCUMENTATION : UNIFIED MODEL DOCUMENTATION PAPER 27 CONV22C.48
CLL CONV22C.49
CLLEND----------------------------------------------------------------- CONV22C.50
C CONV22C.51
C*L ARGUMENTS--------------------------------------------------------- CONV22C.52
C CONV22C.53
SUBROUTINE CONVEC2 (NPNTS,NLEV,K,THEK,THEKP1,QEK,QEKP1,QSEKP1, 4,6CONV22C.54
* DQSKP1,PSTAR,THPK,QPK,THPKP1,QPKP1,XSQKP1, CONV22C.55
* RBUOY,QSEK,DQSK,THPI,QPI,XPK,FLXK,BWKP1,BGMKP1, CONV22C.56
* BGMK,BLOWST,BLAND,BTERM,DEPTH,PREKP1,DTHEK, CONV22C.57
* DQEK,DTHEKP1,DQEKP1,BINIT,CCA,ICCB,ICCT, CONV22C.58
* TCW,EKP14,EKP34,AMDETK,PK,PKP1, CONV22C.59
* EXK,EXKP1,DELEXKP1,DELPK,DELPKP1, CONV22C.60
* CCLWP,CCW,LCCA,LCBASE,LCTOP,LCCLWP, CONV22C.61
* CAPE) CONV22C.62
C CONV22C.63
IMPLICIT NONE CONV22C.64
C CONV22C.65
C---------------------------------------------------------------------- CONV22C.66
C MODEL CONSTANTS CONV22C.67
C---------------------------------------------------------------------- CONV22C.68
*CALL C_R_CP 
ADR1F405.42
*CALL C_EPSLON ADR1F405.43
C CONV22C.90
C CONV22C.91
C---------------------------------------------------------------------- CONV22C.92
C VECTOR LENGTHS AND LOOP COUNTERS CONV22C.93
C---------------------------------------------------------------------- CONV22C.94
C CONV22C.95
INTEGER NPNTS ! IN VECTOR LENGTH CONV22C.96
C CONV22C.97
INTEGER NLEV ! IN NUMBER OF MODEL LAYERS CONV22C.98
C CONV22C.99
INTEGER I ! LOOP COUNTER CONV22C.100
C CONV22C.101
INTEGER K ! PRESENT MODEL LAYER CONV22C.102
C CONV22C.103
C CONV22C.104
C---------------------------------------------------------------------- CONV22C.105
C VARIABLES WHICH ARE INPUT CONV22C.106
C---------------------------------------------------------------------- CONV22C.107
C CONV22C.108
REAL THEK(NPNTS) ! IN POTENTIAL TEMPERATURE OF CLOUD CONV22C.109
! ENVIRONMENT IN LAYER K (K) CONV22C.110
C CONV22C.111
REAL THEKP1(NPNTS) ! IN POTENTIAL TEMPERATURE OF CLOUD CONV22C.112
! ENVIRONMENT IN LAYER K+1 (K) CONV22C.113
C CONV22C.114
REAL QEK(NPNTS) ! IN MIXING RATIO OF CLOUD CONV22C.115
! ENVIRONMENT IN LAYER K (KG/KG) CONV22C.116
C CONV22C.117
REAL QEKP1(NPNTS) ! IN MIXING RATIO OF CLOUD CONV22C.118
! ENVIRONMENT IN LAYER K+1 (KG/KG) CONV22C.119
C CONV22C.120
REAL QSEKP1(NPNTS) ! IN SATURATION MIXING RATIO OF CLOUD CONV22C.121
! ENVIRONMENT IN LAYER K+1 (KG/KG) CONV22C.122
C CONV22C.123
REAL DQSKP1(NPNTS) ! IN GRADIENT OF SATURATION MIXING RATIO CONV22C.124
! WITH POTENTIAL TEMPERATURE FOR THE CONV22C.125
! CLOUD ENVIRONMENT IN LAYER K+1 CONV22C.126
! (KG/KG) CONV22C.127
C CONV22C.128
REAL PSTAR(NPNTS) ! IN SURFACE PRESSURE (PA) CONV22C.129
C CONV22C.130
REAL THPKP1(NPNTS) ! IN PARCEL POTENTIAL TEMPERATURE CONV22C.131
! IN LAYER K (K) CONV22C.132
C CONV22C.133
REAL QPKP1(NPNTS) ! IN PARCEL MIXING RATIO IN LAYER K+1 CONV22C.134
! (KG/KG) CONV22C.135
C CONV22C.136
REAL XSQKP1(NPNTS) ! IN EXCESS WATER IN PARCEL AFTER LIFTING CONV22C.137
! LAYER K TO K+1 (KG/KG) CONV22C.138
C CONV22C.139
REAL RBUOY(NPNTS) ! IN PARCEL BUOYANCY IN LAYER K+1 (K) CONV22C.140
C CONV22C.141
REAL QSEK(NPNTS) ! IN SATURATION MIXING RATIO OF CLOUD CONV22C.142
! ENVIRONMENT IN LAYER K (KG/KG) CONV22C.143
C CONV22C.144
REAL DQSK(NPNTS) ! IN GRADIENT OF SATURATION MIXING RATIO CONV22C.145
! WITH POTENTIAL TEMPERATURE FOR THE CONV22C.146
! CLOUD ENVIRONMENT OF LAYER K CONV22C.147
! (KG/KG/K) CONV22C.148
C CONV22C.149
REAL THPI(NPNTS) ! IN INITIAL PARCEL POTENTIAL TEMPERATURE CONV22C.150
! (K) CONV22C.151
C CONV22C.152
REAL QPI(NPNTS) ! IN INITIAL PARCEL MIXING RATIO CONV22C.153
! (KG/KG) CONV22C.154
C CONV22C.155
LOGICAL BWKP1(NPNTS) ! IN MASK FOR WHETHER CONDENSATE IS CONV22C.156
! LIQUID IN LAYER K+1 CONV22C.157
C CONV22C.158
LOGICAL BGMKP1(NPNTS) ! IN MASK FOR PARCELS WHICH ARE CONV22C.159
! SATURATED IN LAYER K+1 CONV22C.160
C CONV22C.161
LOGICAL BLAND(NPNTS) ! IN LAND/SEA MASK CONV22C.162
C CONV22C.163
LOGICAL BINIT(NPNTS) ! IN MASK FOR THOSE POINTS AT WHICH CONV22C.164
! CONVECTION IS OCCURING CONV22C.165
C CONV22C.166
LOGICAL BLOWST(NPNTS) ! IN MASK FOR THOSE POINTS AT WHICH CONV22C.167
! STABILITY IS LOW ENOUGH FOR CONV22C.168
! CONVECTION TO OCCUR CONV22C.169
C CONV22C.170
REAL EKP14(NPNTS) ! IN ENTRAINMENT COEFFICIENT AT LEVEL CONV22C.171
! K+1/4 MULTIPLIED BY APPROPRIATE CONV22C.172
! LAYER THICKNESS CONV22C.173
C CONV22C.174
REAL EKP34(NPNTS) ! IN ENTRAINMENT COEFFICIENT AT LEVEL CONV22C.175
! K+1/4 MULTIPLIED BY APPROPRIATE CONV22C.176
! LAYER THICKNESS CONV22C.177
C CONV22C.178
REAL AMDETK(NPNTS) ! IN MIXING DETRAINMENT COEFFICIENT CONV22C.179
! AT LEVEL K MULTIPLIED BY APPROPRIATE CONV22C.180
! LAYER THICKNESS CONV22C.181
C CONV22C.182
REAL DELPKP12(NPNTS) ! IN PRESSURE DIFFERENCE BETWEEN CONV22C.183
! MID-POINTS OF LAYERS K AND K+1 CONV22C.184
! (PA) CONV22C.185
C CONV22C.186
REAL PK(NPNTS) ! IN PRESSURE AT MID-POINT OF LAYER K CONV22C.187
! (PA) CONV22C.188
C CONV22C.189
REAL PKP1(NPNTS) ! IN PRESSURE AT MID-POINT OF LAYER K+1 CONV22C.190
! (PA) CONV22C.191
C CONV22C.192
REAL EXK(NPNTS) ! IN EXNER RATIO AT MID-POINT OF LAYER K CONV22C.193
C CONV22C.194
REAL EXKP1(NPNTS) ! IN EXNER RATIO AT MID-POINT OF LAYER K+1 CONV22C.195
C CONV22C.196
REAL DELEXKP1(NPNTS) ! IN DIFFERENCE IN EXNER RATIO BETWEEN CONV22C.197
! MID-POINTS OF LAYERS K AND K+1 CONV22C.198
C CONV22C.199
REAL DELPK(NPNTS) ! IN DIFFERENCE IN PRESSURE ACROSS LAYER K CONV22C.200
! (PA) CONV22C.201
C CONV22C.202
REAL DELPKP1(NPNTS) ! IN DIFFERENCE IN PRESSURE ACROSS CONV22C.203
! LAYER K+1 (PA) CONV22C.204
C CONV22C.205
C CONV22C.206
C---------------------------------------------------------------------- CONV22C.207
C VARIABLES WHICH ARE INPUT BUT WHICH ARE ALSO UPDATED IN THIS ROUTINE CONV22C.208
C---------------------------------------------------------------------- CONV22C.209
C CONV22C.210
REAL THPK(NPNTS) ! INOUT CONV22C.211
! IN PARCEL POTENTIAL TEMPERATURE CONV22C.212
! IN LAYER K (K) CONV22C.213
! OUT PARCEL POTENTIAL TEMPERATURE CONV22C.214
! IN LAYER K+1 (K) CONV22C.215
C CONV22C.216
REAL QPK(NPNTS) ! INOUT CONV22C.217
! IN PARCEL MIXING RATIO IN LAYER K CONV22C.218
! (KG/KG) CONV22C.219
! OUT PARCEL MIXING RATIO IN LAYER K+1 CONV22C.220
! (KG/KG) CONV22C.221
C CONV22C.222
REAL XPK(NPNTS) ! INOUT CONV22C.223
! IN PARCEL CLOUD WATER IN LAYER K CONV22C.224
! (KG/KG) CONV22C.225
! OUT PARCEL CLOUD WATER IN LAYER K+1 CONV22C.226
! (KG/KG) CONV22C.227
C CONV22C.228
REAL FLXK(NPNTS) ! INOUT CONV22C.229
! IN PARCEL MASSFLUX IN LAYER K (PA/S) CONV22C.230
! OUT PARCEL MASSFLUX IN LAYER K+1 (PA/S) CONV22C.231
C CONV22C.232
LOGICAL BGMK(NPNTS) ! INOUT CONV22C.233
! IN MASK FOR PARCELS WHICH ARE CONV22C.234
! SATURATED IN LAYER K CONV22C.235
! OUT MASK FOR PARCELS WHICH ARE CONV22C.236
! SATURATED IN LAYER K+1 CONV22C.237
C CONV22C.238
REAL DTHEK(NPNTS) ! INOUT CONV22C.239
! IN INCREMENT TO MODEL POTENTIAL CONV22C.240
! TEMPERATURE IN LAYER K DUE TO CONV22C.241
! CONVECTION (MAY BE NONE ZERO CONV22C.242
! DUE TO A PREVIOUS SPLIT FINAL CONV22C.243
! DETRAINEMNT CALCULATION) (K/S) CONV22C.244
! OUT UPDATED INCREMENT TO MODEL POTENTIAL CONV22C.245
! TEMPERATURE IN LAYER K DUE TO CONV22C.246
! CONVECTION (K/S) CONV22C.247
C CONV22C.248
REAL DQEK(NPNTS) ! INOUT CONV22C.249
! IN INCREMENT TO MODEL MIXING RATIO CONV22C.250
! IN LAYER K DUE TO CONVECTION CONV22C.251
! (MAY BE NONE ZERO DUE TO A CONV22C.252
! PREVIOUS SPLIT FINAL DETRAINEMNT CONV22C.253
! CALCULATION) (KG/KG/S) CONV22C.254
! OUT UPDATED INCREMENT TO MODEL MIXING CONV22C.255
! RATIO IN LAYER K DUE TO CONVECTION CONV22C.256
! (KG/KG/S) CONV22C.257
C CONV22C.258
REAL TCW(NPNTS) ! INOUT CONV22C.259
! IN TOTAL CONDENSED WATER SUMMED TO CONV22C.260
! LAYER K (KG/M**2/S) CONV22C.261
! OUT UPDATED TOTAL CONDENSED WATER CONV22C.262
! SUMMED TO LAYER K+1 (KG/M**2/S) CONV22C.263
C CONV22C.264
REAL DEPTH(NPNTS) ! INOUT CONV22C.265
! IN DEPTH OF CONVECTIVE CLOUD TO CONV22C.266
! LAYER K (M) CONV22C.267
! OUT UPDATED DEPTH OF CONVECTIVE CLOUD CONV22C.268
! TO LAYER K+1 (M) CONV22C.269
C CONV22C.270
REAL CCLWP(NPNTS) ! INOUT CONV22C.271
! IN CONDENSED WATER PATH SUMMED TO CONV22C.272
! LAYER K (KG/M**2) CONV22C.273
! OUT UPDATED CONDENSED WATER PATH CONV22C.274
! SUMMED TO LAYER K+1 (KG/M**2) CONV22C.275
C CONV22C.276
REAL CAPE(NPNTS) ! IN CONVECTIVE AVAILABLE POTENTIAL CONV22C.277
! ENERGY UP TO THE CURRENT CONVECTING CONV22C.278
! LAYER CONV22C.279
! OUT CONVECTIVE AVAILABLE POTENTIAL CONV22C.280
! ENERGY INCLUDING ADDITION DUE TO CONV22C.281
! THE CAPE WITHIN THE CURRENT LAYER CONV22C.282
C CONV22C.283
C CONV22C.284
C CONV22C.285
C---------------------------------------------------------------------- CONV22C.286
C VARIABLES WHICH ARE OUTPUT CONV22C.287
C---------------------------------------------------------------------- CONV22C.288
C CONV22C.289
LOGICAL BTERM(NPNTS) ! OUT MASK FOR PARCELS WHICH TERMINATE IN CONV22C.290
! LAYER K+1 CONV22C.291
C CONV22C.292
REAL PREKP1(NPNTS) ! OUT PRECIPITATION FROM PARCEL AS IT CONV22C.293
! RISES FROM LAYER K TO K+1 (KG/M**2/S) CONV22C.294
C CONV22C.295
REAL DTHEKP1(NPNTS) ! OUT INCREMENT TO MODEL POTENTIAL CONV22C.296
! TEMPERATURE IN LAYER K+1 DUE TO CONV22C.297
! CONVECTION (K/S) CONV22C.298
C CONV22C.299
REAL DQEKP1(NPNTS) ! OUT INCREMENT TO MODEL MIXING RATIO CONV22C.300
! IN LAYER K+1 DUE TO CONVECTION CONV22C.301
! (KG/KG/S) CONV22C.302
C CONV22C.303
REAL CCA(NPNTS) ! OUT CONVECTIVE CLOUD AMOUNT (%) CONV22C.304
C CONV22C.305
INTEGER ICCB(NPNTS) ! OUT CONVECTIVE CLOUD BASE LEVEL CONV22C.306
C CONV22C.307
INTEGER ICCT(NPNTS) ! OUT CONVECTIVE CLOUD TOP LEVEL CONV22C.308
C CONV22C.309
REAL CCW(NPNTS) ! OUT CONVECTIVE CLOUD WATER(G/KG) ON CONV22C.310
! MODEL LEVELS CONV22C.311
C CONV22C.312
REAL LCCA(NPNTS) ! OUT LOWEST CONV.CLOUD AMOUNT (%) CONV22C.313
C CONV22C.314
INTEGER LCBASE(NPNTS) ! OUT LOWEST CONV.CLOUD BASE LEVEL CONV22C.315
C CONV22C.316
INTEGER LCTOP(NPNTS) ! OUT LOWEST CONV.CLOUD TOP LEVEL CONV22C.317
C CONV22C.318
REAL LCCLWP(NPNTS) ! OUT LOWEST CONV.CLOUD LIQ.WATER PATH CONV22C.319
C CONV22C.320
C CONV22C.321
C---------------------------------------------------------------------- CONV22C.322
C VARIABLES WHICH ARE DEFINED LOCALLY CONV22C.323
C CONV22C.324
C ON THE IBM ARRAYS ARE ALLOCATED USING A PARAMETER STATEMENT CONV22C.325
C CONV22C.326
C ON THE CRAY ARRAYS ARE DYNAMICALLY ALLOCATED CONV22C.327
C---------------------------------------------------------------------- CONV22C.328
C CONV22C.329
REAL THRK(NPNTS) ! PARCEL DETRAINMENT POTENTIAL CONV22C.330
! TEMPERATURE IN LAYER K (K) CONV22C.331
C CONV22C.332
REAL QRK(NPNTS) ! PARCEL DETRAINMENT MIXING RATIO CONV22C.333
! IN LAYER K (KG/KG) CONV22C.334
C CONV22C.335
REAL XPKP1(NPNTS) ! PARCEL CLOUD WATER IN LAYER K+1 (KG/KG) CONV22C.336
C CONV22C.337
REAL FLXKP1(NPNTS ) ! PARCEL MASSFLUX IN LAYER K+1 (PA/S) CONV22C.338
C CONV22C.339
REAL DELTAK(NPNTS) ! PARCEL FORCED DETRAINMENT RATE CONV22C.340
! IN LAYER K MULTIPLIED BY APPROPRIATE CONV22C.341
! LAYER THICKNESS CONV22C.342
C CONV22C.343
REAL THVP,THVE,RHO ! VIRTUAL TEMPERATURE OF PARCEL, VIRTUAL CONV22C.344
! TEMPERATURE OF ENVIRONMENT AND DENSITY, CONV22C.345
! USED IN CALCULATION OF CAPE. CONV22C.346
C CONV22C.347
LOGICAL L_SHALLOW ! \ DUMMY Variables used in parcel in AJX0F404.392
&, L_CCW ! | versions 3A and 3B. Added here to AJX0F404.393
PARAMETER (L_CCW=.FALSE.) ! / keep argument lists consistent. AJX0F404.394
! AJX0F404.395
C---------------------------------------------------------------------- CONV22C.348
C EXTERNAL ROUTINES CALLED CONV22C.349
C---------------------------------------------------------------------- CONV22C.350
C CONV22C.351
EXTERNAL PARCEL,ENVIRON CONV22C.352
C CONV22C.353
C*--------------------------------------------------------------------- CONV22C.354
CL CONV22C.355
CL---------------------------------------------------------------------- CONV22C.356
CL COMPLETE LIFTING PARCELS TO LAYER K+1 CONV22C.357
CL CONV22C.358
CL SUBROUTINE PARCEL CONV22C.359
CL CONV22C.360
CL UM DOCUMENTATION PAPER P27 CONV22C.361
CL SECTIONS (5),(6),(7),(8),(9) CONV22C.362
CL---------------------------------------------------------------------- CONV22C.363
CL CONV22C.364
CALL PARCEL (K,NPNTS,NLEV,PSTAR,THEKP1,THEK,QEKP1,QEK, CONV22C.365
* QSEK,QSEKP1,DQSK,DQSKP1,BLAND,BWKP1, CONV22C.366
* DELTAK,FLXK,THPK,QPK,THRK,QRK, CONV22C.367
* BTERM,THPKP1,QPKP1,PREKP1,XPK,XPKP1,FLXKP1, CONV22C.368
* XSQKP1,THPI,QPI,BGMK,BGMKP1,BLOWST,RBUOY, CONV22C.369
* CCA,ICCB,ICCT,TCW,DEPTH, CONV22C.370
* EKP14,EKP34,AMDETK,DELPKP1,PK,PKP1, CONV22C.371
* EXK,EXKP1,DELEXKP1,CCLWP,CCW, CONV22C.372
& LCCA,LCBASE,LCTOP,LCCLWP,L_SHALLOW,L_CCW) AJX0F404.396
CL CONV22C.374
CL---------------------------------------------------------------------- CONV22C.375
CL CALCULATE THE EFFECT ON THE ENVIRONMENT (EXCEPT FOR THE CONV22C.376
CL THE EVAPORATION OF PRECIPITATION AND CHANGE OF PHASE) CONV22C.377
CL CONV22C.378
CL SUBROUTINE ENVIRON CONV22C.379
CL CONV22C.380
CL UM DOCUMENTATION PAPER P27 CONV22C.381
CL SECTION (10) CONV22C.382
CL---------------------------------------------------------------------- CONV22C.383
CL CONV22C.384
CALL ENVIRON (NPNTS,DTHEK,DQEK,DTHEKP1,DQEKP1, CONV22C.385
* THEK,QEK,DELTAK,FLXK, CONV22C.386
* THPK,QPK,THRK,QRK,THEKP1,QEKP1, CONV22C.387
* BTERM,THPKP1,QPKP1,XPK,XPKP1,BWKP1,FLXKP1, CONV22C.388
* BLOWST,EKP14,EXK,EXKP1,DELPK,DELPKP1, CONV22C.389
* AMDETK) CONV22C.390
C CONV22C.391
DO 10 I=1,NPNTS CONV22C.392
C CONV22C.393
C----------------------------------------------------------------------- CONV22C.394
C RESET BINIT WHERE CONVECTION HAS TERMINATED CONV22C.395
C----------------------------------------------------------------------- CONV22C.396
C CONV22C.397
BINIT(I) = .NOT.BTERM(I) CONV22C.398
10 CONTINUE CONV22C.399
C CONV22C.400
CL--------------------------------------------------------------------- CONV22C.401
CL CALCULATE CONTRIBUTION TO CAPE AND RATE OF CHANGE OF CAPE DUE TO CONV22C.402
CL THE UPDRAUGHT CONV22C.403
CL--------------------------------------------------------------------- CONV22C.404
C CONV22C.405
DO I=1,NPNTS CONV22C.406
THVP=THPK(I)*(1.0+C_VIRTUAL*QPK(I)) CONV22C.407
THVE=THEK(I)*(1.0+C_VIRTUAL*QEK(I)) CONV22C.408
RHO=PK(I)/(R*THEK(I)*EXK(I)) CONV22C.409
C CONV22C.410
CAPE(I)=CAPE(I)+(THVP-THVE)*DELPK(I)/(RHO*THVE) AJX2F405.7
C CONV22C.412
IF(BTERM(I))THEN CONV22C.413
THVP=THPKP1(I)*(1.0+C_VIRTUAL*QPKP1(I)) CONV22C.414
THVE=THEKP1(I)*(1.0+C_VIRTUAL*QEKP1(I)) CONV22C.415
RHO=PKP1(I)/(R*THEKP1(I)*EXKP1(I)) CONV22C.416
C CONV22C.417
CAPE(I)=CAPE(I)+(THVP-THVE)*DELPKP1(I)/(RHO*THVE) AJX2F405.8
C CONV22C.419
END IF CONV22C.420
C CONV22C.421
END DO CONV22C.422
C CONV22C.423
CL CONV22C.424
CL--------------------------------------------------------------------- CONV22C.425
CL SWAP PARCEL VALUES READY FOR THE NEXT PART OF ASCENT CONV22C.426
CL FROM LAYER K+1 TO K+2 CONV22C.427
CL--------------------------------------------------------------------- CONV22C.428
CL CONV22C.429
IF (K .LT. NLEV-1) THEN CONV22C.430
DO 30 I=1,NPNTS CONV22C.431
IF (BINIT(I)) THEN CONV22C.432
THPK(I) = THPKP1(I) CONV22C.433
QPK(I) = QPKP1(I) CONV22C.434
XPK(I) = XPKP1(I) CONV22C.435
FLXK(I) = FLXKP1(I) CONV22C.436
BGMK(I) = BGMKP1(I) CONV22C.437
END IF CONV22C.438
30 CONTINUE CONV22C.439
END IF CONV22C.440
C CONV22C.441
RETURN CONV22C.442
END CONV22C.443
*ENDIF CONV22C.444