*IF DEF,A05_2A,OR,DEF,A05_2C AJX1F405.156
C ******************************COPYRIGHT****************************** GTS2F400.1765
C (c) CROWN COPYRIGHT 1995, METEOROLOGICAL OFFICE, All Rights Reserved. GTS2F400.1766
C GTS2F400.1767
C Use, duplication or disclosure of this code is subject to the GTS2F400.1768
C restrictions as set forth in the contract. GTS2F400.1769
C GTS2F400.1770
C Meteorological Office GTS2F400.1771
C London Road GTS2F400.1772
C BRACKNELL GTS2F400.1773
C Berkshire UK GTS2F400.1774
C RG12 2SZ GTS2F400.1775
C GTS2F400.1776
C If no contract has been raised with this copy of the code, the use, GTS2F400.1777
C duplication or disclosure of it is strictly prohibited. Permission GTS2F400.1778
C to do so must first be obtained in writing from the Head of Numerical GTS2F400.1779
C Modelling at the above address. GTS2F400.1780
C ******************************COPYRIGHT****************************** GTS2F400.1781
C GTS2F400.1782
CLL SUBROUTINE DDRAUGHT----------------------------------------------- DDRAUG2A.3
CLL DDRAUG2A.4
CLL PURPOSE : DOWNDRAUGHT ROUTINE DDRAUG2A.5
CLL DDRAUG2A.6
CLL CONVECTIVE DOWNDRAUGHT BASED ON PARCEL THEORY DDRAUG2A.7
CLL DDRAUG2A.8
CLL CARRY OUT DRY DESCENT DDRAUG2A.9
CLL DDRAUG2A.10
CLL CALCULATE SUBSATURATION DDRAUG2A.11
CLL DDRAUG2A.12
CLL CALCULATE EFFECT ON THE ENVIRONMENT DDRAUG2A.13
CLL DDRAUG2A.14
CLL SUITABLE FOR SINGLE COLUMN MODEL USE DDRAUG2A.15
CLL DDRAUG2A.16
CLL CODE WRITTEN FOR CRAY Y-MP BY S.BETT AND D.GREGORY AUTUMN 1991 DDRAUG2A.17
CLL DDRAUG2A.18
CLL MODEL MODIFICATION HISTORY FROM MODEL VERSION 3.0: DDRAUG2A.19
CLL VERSION DATE DDRAUG2A.20
CLL 4.5 Jul. 98 Kill the IBM specific lines (JCThil) AJC1F405.6
CLL DDRAUG2A.21
CLL PROGRAMMING STANDARDS : UNIFIED MODEL DOCUMENTATION PAPER NO. 3 DDRAUG2A.22
CLL VERSION NO. 4 DATED 5/2/92 DDRAUG2A.23
CLL DDRAUG2A.24
CLL LOGICAL COMPONENTS COVERED: DDRAUG2A.25
CLL DDRAUG2A.26
CLL SYSTEM TASK : P27 DDRAUG2A.27
CLL DDRAUG2A.28
CLL DOCUMENTATION : UNIFIED MODEL DOCUMENTATION PAPER P27 DDRAUG2A.29
CLL DDRAUG2A.30
CLLEND----------------------------------------------------------------- DDRAUG2A.31
C DDRAUG2A.32
C*L ARGUMENTS--------------------------------------------------------- DDRAUG2A.33
C DDRAUG2A.34
SUBROUTINE DDRAUGHT (NPNTS,K,KCT,THDD_K,QDD_K,THE_K,THE_KM1,QE_K, 4,21DDRAUG2A.35
& QE_KM1,DTHBYDT_K,DTHBYDT_KM1,DQBYDT_K, DDRAUG2A.36
& DQBYDT_KM1,FLX_DD_K,P_KM1,DELPK, DDRAUG2A.37
& DELPKM1,EXK,EXKM1,DELTD,DELQD,AMDETK,EKM14, DDRAUG2A.38
& EKM34,RAIN,SNOW,BDD_START,BDDWT_K, DDRAUG2A.39
& BDDWT_KM1,BDD_ON,B_DD_END,CCA) DDRAUG2A.40
C DDRAUG2A.41
IMPLICIT NONE DDRAUG2A.42
C DDRAUG2A.43
C----------------------------------------------------------------------- DDRAUG2A.44
C MODEL CONSTANTS DDRAUG2A.45
C----------------------------------------------------------------------- DDRAUG2A.46
C DDRAUG2A.47
*CALL C_0_DG_C 
DDRAUG2A.48
*CALL DDKMDET DDRAUG2A.49
C DDRAUG2A.50
C----------------------------------------------------------------------- DDRAUG2A.51
C VECTOR LENGTHS AND LOOP COUNTERS DDRAUG2A.52
C----------------------------------------------------------------------- DDRAUG2A.53
C DDRAUG2A.54
C DDRAUG2A.58
INTEGER I ! LOOP COUNTER DDRAUG2A.59
C DDRAUG2A.60
INTEGER NPNTS ! IN NUMBER OF POINTS DDRAUG2A.61
C DDRAUG2A.62
INTEGER K ! IN PRESENT MODEL LAYER DDRAUG2A.63
C DDRAUG2A.64
C----------------------------------------------------------------------- DDRAUG2A.65
C VARIABLES WHICH ARE INPUT DDRAUG2A.66
C----------------------------------------------------------------------- DDRAUG2A.67
C DDRAUG2A.68
INTEGER KCT ! IN CONVECTIVE CLOUD TOP DDRAUG2A.69
C DDRAUG2A.70
REAL THE_KM1(NPNTS) ! IN POTENTIAL TEMPERATURE OF DDRAUG2A.71
! ENVIRONMENT IN LAYER K-1 (K) DDRAUG2A.72
C DDRAUG2A.73
REAL QE_KM1(NPNTS) ! IN MIXING RATIO OF ENVIRONMENT IN DDRAUG2A.74
! LAYER K-1 (KG/KG) DDRAUG2A.75
C DDRAUG2A.76
REAL P_KM1(NPNTS) ! IN PRESSURE OF LAYER K-1 (PA) DDRAUG2A.77
C DDRAUG2A.78
REAL DELPK(NPNTS) ! IN CHANGE IN PRESSURE ACROSS DDRAUG2A.79
! LAYER K (PA) DDRAUG2A.80
C DDRAUG2A.81
REAL DELPKM1(NPNTS) ! IN CHANGE IN PRESSURE ACROSS DDRAUG2A.82
! LAYER K-1 (PA) DDRAUG2A.83
C DDRAUG2A.84
REAL EXK(NPNTS) ! IN EXNER RATIO IN LAYER K DDRAUG2A.85
C DDRAUG2A.86
REAL EXKM1(NPNTS) ! IN EXNER RATIO IN LAYER K-1 DDRAUG2A.87
C DDRAUG2A.88
REAL AMDETK(NPNTS) ! IN MIXING DETRAINMENT RATE DDRAUG2A.89
C DDRAUG2A.90
REAL EKM14(NPNTS) ! IN EXNER RATIO AT LAYER K-1/4 DDRAUG2A.91
C DDRAUG2A.92
REAL EKM34(NPNTS) ! IN EXNER RATIO AT LAYER K-3/4 DDRAUG2A.93
C DDRAUG2A.94
REAL DELTD(NPNTS) ! IN COOLING NECESSARY TO ACHIEVE DDRAUG2A.95
! SATURATION (K) DDRAUG2A.96
C DDRAUG2A.97
REAL DELQD(NPNTS) ! IN MOISTENING NECESSARY TO ACHIEVE DDRAUG2A.98
! SATURATION (KG/KG) DDRAUG2A.99
C DDRAUG2A.100
LOGICAL BDDWT_K(NPNTS) ! IN MASK FOR THOSE POINTS IN DDRAUG2A.101
! DOWNDRAUGHT WHERE PRECIPITATION DDRAUG2A.102
! IS LIQUID IN LAYER K DDRAUG2A.103
C DDRAUG2A.104
LOGICAL BDDWT_KM1(NPNTS) ! IN MASK FOR THOSE POINTS IN DDRAUG2A.105
! DOWNDRAUGHT WHERE PRECIPITATION DDRAUG2A.106
! IS LIQUID IN LAYER K-1 DDRAUG2A.107
C DDRAUG2A.108
REAL CCA(NPNTS) ! IN CONVECTIVE CLOUD AMOUNT DDRAUG2A.109
C DDRAUG2A.110
C----------------------------------------------------------------------- DDRAUG2A.111
C VARIABLES WHICH ARE INPUT AND OUTPUT DDRAUG2A.112
C----------------------------------------------------------------------- DDRAUG2A.113
C DDRAUG2A.114
REAL THDD_K(NPNTS) ! INOUT DDRAUG2A.115
! IN POTENTIAL TEMPERATURE OF DDRAUG2A.116
! DOWNDRAUGHT IN LAYER K (K) DDRAUG2A.117
! OUT POTENTIAL TEMPERATURE RESET DDRAUG2A.118
! FOR NEXT LAYER (K) DDRAUG2A.119
C DDRAUG2A.120
REAL QDD_K(NPNTS) ! INOUT DDRAUG2A.121
! IN DOWNDRAUGHT MIXING RATIO OF DDRAUG2A.122
! LAYER K (KG/KG) DDRAUG2A.123
! OUT MIXING RATIO RESET FOR NEXT DDRAUG2A.124
! LAYER (KG/KG) DDRAUG2A.125
C DDRAUG2A.126
REAL THE_K(NPNTS) ! INOUT DDRAUG2A.127
! IN POTENTIAL TEMPERATURE OF DDRAUG2A.128
! ENVIRONMENT IN LAYER K (K) DDRAUG2A.129
! OUT ENVIRONMENT POTENTIAL DDRAUG2A.130
! TEMPERATURE RESET FOR NEXT DDRAUG2A.131
! LAYER (K) DDRAUG2A.132
C DDRAUG2A.133
REAL QE_K(NPNTS) ! INOUT DDRAUG2A.134
! IN MIXING RATIO OF ENVIRONMENT DDRAUG2A.135
! LAYER K (KG/KG) DDRAUG2A.136
! OUT ENVIRONMENT MIXING RATIO DDRAUG2A.137
! RESET FOR NEXT LAYER (KG/KG) DDRAUG2A.138
C DDRAUG2A.139
REAL FLX_DD_K(NPNTS) ! INOUT DDRAUG2A.140
! IN DOWNDRAUGHT MASS FLUX OF DDRAUG2A.141
! LAYER K (PA/S) DDRAUG2A.142
! OUT DOWNDRAUGHT MASS FLUX RESET DDRAUG2A.143
! FOR NEXT LAYER (PA/S) DDRAUG2A.144
C DDRAUG2A.145
REAL RAIN(NPNTS) ! INOUT DDRAUG2A.146
! IN AMOUNT OF RAIN (KG/M**2/S) DDRAUG2A.147
! OUT UPDATED RAINFALL (KG/M**2/S) DDRAUG2A.148
C DDRAUG2A.149
REAL SNOW(NPNTS) ! INOUT DDRAUG2A.150
! IN AMOUNT OF SNOW(KG/M**2/S) DDRAUG2A.151
! OUT UPDATED SNOWFALL (KG/M**2/S) DDRAUG2A.152
C DDRAUG2A.153
REAL DTHBYDT_K(NPNTS) ! INOUT DDRAUG2A.154
! IN INCREMENT TO MODEL POTENTIAL DDRAUG2A.155
! TEMPERATURE OF LAYER K (K/S) DDRAUG2A.156
! OUT UPDATED INCREMENT TO MODEL DDRAUG2A.157
! POTENTIAL TEMPERATURE IN DDRAUG2A.158
! LAYER K (K/S) DDRAUG2A.159
C DDRAUG2A.160
REAL DTHBYDT_KM1(NPNTS) ! INOUT DDRAUG2A.161
! IN INCREMENT TO MODEL POTENTIAL DDRAUG2A.162
! TEMPERATURE IN LAYER K-1 (K/S) DDRAUG2A.163
! OUT UPDATED INCREMENT TO MODEL DDRAUG2A.164
! POTENTIAL TEMPERATURE IN DDRAUG2A.165
! LAYER K-1 (K/S) DDRAUG2A.166
C DDRAUG2A.167
REAL DQBYDT_K(NPNTS) ! INOUT DDRAUG2A.168
! IN INCREMENT TO MODEL MIXING DDRAUG2A.169
! RATIO IN LAYER K (KG/KG/S) DDRAUG2A.170
! OUT UPDATED INCREMENT TO MODEL DDRAUG2A.171
! MIXING RATIO IN LAYER K DDRAUG2A.172
! (KG/KG/S) DDRAUG2A.173
C DDRAUG2A.174
REAL DQBYDT_KM1(NPNTS) ! INOUT DDRAUG2A.175
! IN INCREMENT TO MODEL MIXING DDRAUG2A.176
! RATIO IN LAYER K-1 (KG/KG/S) DDRAUG2A.177
! OUT UPDATED INCREMENT TO MODEL DDRAUG2A.178
! MIXING RATIO IN LAYER K-1 DDRAUG2A.179
! (KG/KG/S) DDRAUG2A.180
C DDRAUG2A.181
LOGICAL BDD_ON(NPNTS) ! INOUT DDRAUG2A.182
! IN MASK FOR THOSE POINTS WHERE DD DDRAUG2A.183
! HAS CONTINUED FROM LAYER K+1 DDRAUG2A.184
! OUT MASK FOR THOSE POINTS WHERE DD DDRAUG2A.185
! CONTINUES TO LAYER K-1 DDRAUG2A.186
C DDRAUG2A.187
C----------------------------------------------------------------------- DDRAUG2A.188
C VARIABLES WHICH ARE OUTPUT DDRAUG2A.189
C----------------------------------------------------------------------- DDRAUG2A.190
C DDRAUG2A.191
LOGICAL BDD_START(NPNTS) ! OUT MASK FOR THOSE POINTS WHERE DDRAUG2A.192
! DOWNDRAUGHT MAY START IN DDRAUG2A.193
! LAYER K-1 DDRAUG2A.194
C DDRAUG2A.195
LOGICAL B_DD_END(NPNTS) ! OUT MASK FOR THOSE POINTS WHERE DDRAUG2A.196
! DOWNDRAUGHT IS ENDING IN DDRAUG2A.197
! LAYER K-1 DDRAUG2A.198
C DDRAUG2A.199
C----------------------------------------------------------------------- DDRAUG2A.200
C VARIABLES WHICH ARE DEFINED LOCALLY DDRAUG2A.201
C----------------------------------------------------------------------- DDRAUG2A.202
C DDRAUG2A.203
C DDRAUG2A.232
REAL THDD_KM1(NPNTS) ! POTENTIAL TEMPERATURE OF DDRAUG2A.233
! DOWNDRAUGHT IN LAYER K-1 (K) DDRAUG2A.234
C DDRAUG2A.235
REAL QDD_KM1(NPNTS) ! DOWNDRAUGHT MIXING RATIO OF DDRAUG2A.236
! LAYER K-1 (KG/KG) DDRAUG2A.237
C DDRAUG2A.238
REAL QSATDD(NPNTS) ! SATURATED DOWNDRAUGHT MIXING DDRAUG2A.239
! RATIO (KG/KG) DDRAUG2A.240
C DDRAUG2A.241
REAL TDD_KM1(NPNTS) ! TEMPERATURE OF DOWNDRAUGHT DDRAUG2A.242
! IN LAYER K-1 (K) DDRAUG2A.243
C DDRAUG2A.244
REAL THDDS(NPNTS) ! POTENTIAL TEMPERATURE OF DDRAUG2A.245
! SATURATED DOWNDRAUGHT (K) DDRAUG2A.246
C DDRAUG2A.247
REAL QDDS(NPNTS) ! SATURATED DOWNDRAUGHT MIXING DDRAUG2A.248
! RATIO (KG/KG) DDRAUG2A.249
C DDRAUG2A.250
REAL FLX_DD_KM1(NPNTS) ! DOWNDRAUGHT MASS FLUX IN DDRAUG2A.251
! LAYER K-1 (PA/S) DDRAUG2A.252
C DDRAUG2A.253
REAL RAIN_TMP(NPNTS) ! LIQUID PRECIPITATION STORE DDRAUG2A.254
C DDRAUG2A.255
REAL SNOW_TMP(NPNTS) ! SNOW STORE DDRAUG2A.256
C DDRAUG2A.257
C DDRAUG2A.259
C----------------------------------------------------------------------- DDRAUG2A.260
C EXTERNAL ROUTINES CALLED DDRAUG2A.261
C----------------------------------------------------------------------- DDRAUG2A.262
C DDRAUG2A.263
EXTERNAL SATCAL, CRS_FRZL, QSAT, DEVAP, TERMDD, DDRAUG2A.264
* DD_ENV, EVP DDRAUG2A.265
C DDRAUG2A.266
C----------------------------------------------------------------------- DDRAUG2A.267
C CALCULATE MASK FOR THOSE POINTS IN DOWNDRAUGHT WHERE PRECIPITATION DDRAUG2A.268
C IS LIQUID DDRAUG2A.269
C DDRAUG2A.270
C STORE PRECIPITATION IN LAYER K IN TEMPORARY VARIABLES DDRAUG2A.271
C----------------------------------------------------------------------- DDRAUG2A.272
C DDRAUG2A.273
DO I=1,NPNTS DDRAUG2A.274
IF (K .EQ. KCT+1 .OR. BDD_START(I)) THEN DDRAUG2A.275
BDDWT_K(I) = THDD_K(I) .GT. TM/EXK(I) DDRAUG2A.276
ELSE DDRAUG2A.277
BDDWT_K(I) = BDDWT_KM1(I) DDRAUG2A.278
END IF DDRAUG2A.279
RAIN_TMP(I) = RAIN(I) DDRAUG2A.280
SNOW_TMP(I) = SNOW(I) DDRAUG2A.281
C DDRAUG2A.282
C----------------------------------------------------------------------- DDRAUG2A.283
C DRY DESCENT FROM LAYER K TO K-1 DDRAUG2A.284
C DDRAUG2A.285
C ENTRAINMENT CALCULATION DDRAUG2A.286
C----------------------------------------------------------------------- DDRAUG2A.287
C DDRAUG2A.288
THDD_KM1(I) = (THDD_K(I)+(EKM14(I)*THE_K(I)) + DDRAUG2A.289
* (1.0+EKM14(I))*EKM34(I)*THE_KM1(I)) / DDRAUG2A.290
* ((1.0+EKM14(I))*(1.0+EKM34(I))) DDRAUG2A.291
QDD_KM1(I) = (QDD_K(I)+(EKM14(I)*QE_K(I)) + DDRAUG2A.292
* (1.0+EKM14(I))*EKM34(I)*QE_KM1(I))/ DDRAUG2A.293
* ((1.0+EKM14(I))*(1.0+EKM34(I))) DDRAUG2A.294
C DDRAUG2A.295
C----------------------------------------------------------------------- DDRAUG2A.296
C UPDATE MASS FLUX AND CALCULATE TEMPERATURE OF LAYER K-1 DDRAUG2A.297
C----------------------------------------------------------------------- DDRAUG2A.298
C DDRAUG2A.299
FLX_DD_KM1(I) = FLX_DD_K(I)*(1.0+EKM34(I))*(1.0+EKM14(I))* DDRAUG2A.300
* (1.0-AMDETK(I)) DDRAUG2A.301
C DDRAUG2A.302
TDD_KM1(I) = THDD_KM1(I)*EXKM1(I) DDRAUG2A.303
END DO DDRAUG2A.304
C DDRAUG2A.305
C----------------------------------------------------------------------- DDRAUG2A.306
C CALCULATE SUBSATURATION DDRAUG2A.307
C CALCULATE TEMPERATURE IF BROUGHT TO SATURATION DDRAUG2A.308
C----------------------------------------------------------------------- DDRAUG2A.309
C DDRAUG2A.310
CALL SATCAL(NPNTS,TDD_KM1,THDD_KM1,P_KM1,QDDS,THDDS, DDRAUG2A.311
& K,EXKM1,QDD_KM1,THE_KM1) DDRAUG2A.312
C DDRAUG2A.313
DO I=1,NPNTS DDRAUG2A.314
BDDWT_KM1(I) = THDDS(I) .GT. TM/EXKM1(I) DDRAUG2A.315
END DO DDRAUG2A.316
C DDRAUG2A.317
C----------------------------------------------------------------------- DDRAUG2A.318
C CALCULATE CHANGE OF PHASE DUE TO DOWNDRAUGHT SATURATION TEMPERATURE DDRAUG2A.319
C----------------------------------------------------------------------- DDRAUG2A.320
C DDRAUG2A.321
CALL CRS_FRZL (NPNTS,RAIN,SNOW,THDD_KM1,EXKM1,FLX_DD_KM1, DDRAUG2A.322
& BDDWT_KM1) DDRAUG2A.323
C DDRAUG2A.324
DO I=1,NPNTS DDRAUG2A.325
TDD_KM1(I) = THDD_KM1(I)*EXKM1(I) DDRAUG2A.326
END DO DDRAUG2A.327
C DDRAUG2A.328
C----------------------------------------------------------------------- DDRAUG2A.329
C RECALCULATE SUBSATURATION TEMPERATURE DDRAUG2A.330
C----------------------------------------------------------------------- DDRAUG2A.331
C DDRAUG2A.332
CALL SATCAL(NPNTS,TDD_KM1,THDD_KM1,P_KM1,QDDS,THDDS, DDRAUG2A.333
& K,EXKM1,QDD_KM1,THE_KM1) DDRAUG2A.334
C DDRAUG2A.335
C----------------------------------------------------------------------- DDRAUG2A.336
C CALCULATE MOISTURE SUBSATURATION DDRAUG2A.337
C----------------------------------------------------------------------- DDRAUG2A.338
C DDRAUG2A.339
CALL QSAT(QSATDD,TDD_KM1,P_KM1,NPNTS) DDRAUG2A.340
C DDRAUG2A.341
C----------------------------------------------------------------------- DDRAUG2A.342
C EVAPORATION CALCULATION AND ADJUSTMENT OF DOWNDRAUGHT TEMPERATURE DDRAUG2A.343
C AND MOISTURE DDRAUG2A.344
C----------------------------------------------------------------------- DDRAUG2A.345
C DDRAUG2A.346
CALL DEVAP (NPNTS,THDD_K,THDD_KM1,QDD_KM1,THDDS,QDDS, DDRAUG2A.347
& FLX_DD_KM1,EXK,EXKM1,QSATDD,RAIN,SNOW, DDRAUG2A.348
& DELPKM1,BDDWT_KM1,CCA,P_KM1) DDRAUG2A.349
C DDRAUG2A.350
C----------------------------------------------------------------------- DDRAUG2A.351
C CHECK IF PARCEL STILL NEGATIVELY BUOYANT SUCH THAT DOWNDRAUGHT CAN DDRAUG2A.352
C CONTINUE TO K-1 DDRAUG2A.353
C----------------------------------------------------------------------- DDRAUG2A.354
C DDRAUG2A.355
CALL TERMDD (NPNTS,BDD_START,THDD_KM1,QDD_KM1,THE_KM1, DDRAUG2A.356
& QE_KM1,K,B_DD_END,BDD_ON) DDRAUG2A.357
C DDRAUG2A.358
C----------------------------------------------------------------------- DDRAUG2A.359
C CALCULATE THE EFFECT ON THE ENVIRONMENT IN LAYER K DDRAUG2A.360
C----------------------------------------------------------------------- DDRAUG2A.361
C DDRAUG2A.362
CALL DD_ENV (NPNTS,THDD_K,THDD_KM1,QDD_K,QDD_KM1,THE_K,THE_KM1, DDRAUG2A.363
& QE_K,QE_KM1,DTHBYDT_K,DTHBYDT_KM1,DQBYDT_K, DDRAUG2A.364
& DQBYDT_KM1,FLX_DD_K,FLX_DD_KM1,DELPK,DELPKM1, DDRAUG2A.365
& DELTD,DELQD,AMDETK,EKM14,B_DD_END, DDRAUG2A.366
& BDD_START,BDD_ON) DDRAUG2A.367
C DDRAUG2A.368
C----------------------------------------------------------------------- DDRAUG2A.369
C RESET DOWNDRAUGHT BIT VECTORS DDRAUG2A.370
C DDRAUG2A.371
C----------------------------------------------------------------------- DDRAUG2A.372
C DDRAUG2A.373
DO I=1,NPNTS DDRAUG2A.374
BDD_START(I) = .FALSE. DDRAUG2A.375
IF (.NOT. BDD_ON(I)) THEN DDRAUG2A.376
RAIN(I) = RAIN_TMP(I) DDRAUG2A.377
SNOW(I) = SNOW_TMP(I) DDRAUG2A.378
END IF DDRAUG2A.379
IF (B_DD_END(I)) BDD_ON(I) = .FALSE. DDRAUG2A.380
END DO DDRAUG2A.381
C DDRAUG2A.382
C----------------------------------------------------------------------- DDRAUG2A.383
C SWITCH POTENTIAL TEMPERATURE, MIXING RATIO AND MASS FLUX FOR DDRAUG2A.384
C CALCULATION AT NEXT MODEL LAYER DDRAUG2A.385
C----------------------------------------------------------------------- DDRAUG2A.386
C DDRAUG2A.387
IF (K.GT.2) THEN DDRAUG2A.388
DO I=1,NPNTS DDRAUG2A.389
IF (BDD_ON(I)) THEN DDRAUG2A.390
THDD_K(I) = THDD_KM1(I) DDRAUG2A.391
QDD_K(I) = QDD_KM1(I) DDRAUG2A.392
FLX_DD_K(I) = FLX_DD_KM1(I) DDRAUG2A.393
END IF DDRAUG2A.394
END DO DDRAUG2A.395
END IF DDRAUG2A.396
DDRAUG2A.397
RETURN DDRAUG2A.398
END DDRAUG2A.399
C DDRAUG2A.400
*ENDIF DDRAUG2A.401