*IF DEF,A05_2A,OR,DEF,A05_2C,OR,DEF,A05_3B,OR,DEF,A05_3C AJX1F405.147
C ******************************COPYRIGHT****************************** GTS2F400.10117
C (c) CROWN COPYRIGHT 1995, METEOROLOGICAL OFFICE, All Rights Reserved. GTS2F400.10118
C GTS2F400.10119
C Use, duplication or disclosure of this code is subject to the GTS2F400.10120
C restrictions as set forth in the contract. GTS2F400.10121
C GTS2F400.10122
C Meteorological Office GTS2F400.10123
C London Road GTS2F400.10124
C BRACKNELL GTS2F400.10125
C Berkshire UK GTS2F400.10126
C RG12 2SZ GTS2F400.10127
C GTS2F400.10128
C If no contract has been raised with this copy of the code, the use, GTS2F400.10129
C duplication or disclosure of it is strictly prohibited. Permission GTS2F400.10130
C to do so must first be obtained in writing from the Head of Numerical GTS2F400.10131
C Modelling at the above address. GTS2F400.10132
C ******************************COPYRIGHT****************************** GTS2F400.10133
C GTS2F400.10134
CLL SUBROUTINE TERM_CON----------------------------------------------- TERMCO1A.3
CLL TERMCO1A.4
CLL PURPOSE : RETURENS A MASK FOR POINTS AT WHICH CONVECTION TERMCO1A.5
CLL IS TERMINATING TERMCO1A.6
CLL TERMCO1A.7
CLL SUITABLE FOR SINGLE COLUMN MODEL USE TERMCO1A.8
CLL TERMCO1A.9
CLL CODE REWORKED FOR CRAY Y-MP BY D.GREGORY AUTUMN/WINTER 1989/90 TERMCO1A.10
CLL TERMCO1A.11
CLL MODEL MODIFICATION HISTORY FROM MODEL VERSION 3.0: TERMCO1A.12
CLL VERSION DATE TERMCO1A.13
CLL TERMCO1A.14
CLL PROGRAMMING STANDARDS : UNIFIED MODEL DOCUMENTATION PAPER NO. 4 TERMCO1A.15
CLL VERSION NO. 1 TERMCO1A.16
CLL TERMCO1A.17
CLL LOGICAL COMPONENTS COVERED: P27 TERMCO1A.18
CLL TERMCO1A.19
CLL DOCUMENTATION : UNIFIED MODEL DOCUMENTATION PAPER P27 TERMCO1A.20
CLL TERMCO1A.21
CLLEND----------------------------------------------------------------- TERMCO1A.22
C TERMCO1A.23
C*L ARGUMENTS--------------------------------------------------------- TERMCO1A.24
C TERMCO1A.25
SUBROUTINE TERM_CON(NPNTS,NLEV,K,BTERM,BWKP1,FLXKP1,THEKP1,QEKP1, 2TERMCO1A.26
* THPI,QPI,QSEKP1,DELTAK,EXKP1,EKP14,EKP34, TERMCO1A.27
* PSTAR) TERMCO1A.28
C TERMCO1A.29
IMPLICIT NONE TERMCO1A.30
C TERMCO1A.31
C----------------------------------------------------------------------- TERMCO1A.32
C MODEL CONSTANTS TERMCO1A.33
C----------------------------------------------------------------------- TERMCO1A.34
C TERMCO1A.35
*CALL C_R_CP 
TERMCO1A.36
*CALL MPARFL TERMCO1A.37
*CALL XSBMIN TERMCO1A.38
*CALL C_LHEAT TERMCO1A.39
*CALL C_EPSLON TERMCO1A.40
*CALL QSTICE TERMCO1A.41
C TERMCO1A.42
C----------------------------------------------------------------------- TERMCO1A.43
C VECTOR LENGTHS AND LOOP COUNTERS TERMCO1A.44
C----------------------------------------------------------------------- TERMCO1A.45
C TERMCO1A.46
INTEGER NPNTS ! IN VECTOR LENGTH TERMCO1A.47
C TERMCO1A.48
INTEGER NLEV ! IN NUMBER OF MODEL LAYER TERMCO1A.49
C TERMCO1A.50
INTEGER K ! IN PRESENT MODEL LAYER TERMCO1A.51
C TERMCO1A.52
INTEGER I ! LOOP COUNTER TERMCO1A.53
C TERMCO1A.54
C TERMCO1A.55
C----------------------------------------------------------------------- TERMCO1A.56
C VARIABLES THAT ARE INPUT TERMCO1A.57
C----------------------------------------------------------------------- TERMCO1A.58
C TERMCO1A.59
REAL THEKP1(NPNTS) ! IN POTENTIAL TEMPERATURE OF CLOUD TERMCO1A.60
! ENVIRONMENT IN LAYER K+1 (K) TERMCO1A.61
C TERMCO1A.62
REAL QEKP1(NPNTS) ! IN MIXING RATIO OF CLOUD TERMCO1A.63
! ENVIRONMENT IN LAYER K+1 (KG/KG) TERMCO1A.64
C TERMCO1A.65
REAL QSEKP1(NPNTS) ! IN SATURATION MIXING RATIO OF CLOUD TERMCO1A.66
! ENVIRONMENT IN LAYER K+1 (KG/KG) TERMCO1A.67
C TERMCO1A.68
REAL THPI(NPNTS) ! IN INITIAL PARCEL POTENTIAL TEMPERATURE TERMCO1A.69
! (K) TERMCO1A.70
C TERMCO1A.71
REAL QPI(NPNTS) ! IN INITIAL PARCEL MIXING RATIO (KG/KG) TERMCO1A.72
C TERMCO1A.73
REAL FLXKP1(NPNTS) ! IN PARCEL MASSFLUX IN LAYER K+1 (PA/S) TERMCO1A.74
C TERMCO1A.75
LOGICAL BWKP1(NPNTS) ! IN MASK FOR WHETHER CONDENSATE IS TERMCO1A.76
! LIQUID IN LAYER K+1 TERMCO1A.77
C TERMCO1A.78
REAL DELTAK(NPNTS) ! IN FORCED DETRAINMENT IN LAYER K TERMCO1A.79
! MULTIPLIED BY APPROPRIATE TERMCO1A.80
! LAYER THICKNESS TERMCO1A.81
C TERMCO1A.82
REAL EXKP1(NPNTS) ! IN EXNER RATIO FOR LEVEL K+1 TERMCO1A.83
C TERMCO1A.84
REAL EKP14(NPNTS) ! IN ENTRAINMENT RATE FOR LEVEL K+1/4 TERMCO1A.85
! MULTIPLIED BY APPROPRIATE TERMCO1A.86
! LAYER THICKNESS TERMCO1A.87
C TERMCO1A.88
REAL EKP34(NPNTS) ! IN ENTRAINMENT RATE FOR LEVEL K+3/4 TERMCO1A.89
! MULTIPLIED BY APPROPRIATE TERMCO1A.90
! LAYER THICKNESS TERMCO1A.91
C TERMCO1A.92
REAL PSTAR(NPNTS) ! IN SURFACE PRESSURE (PA) TERMCO1A.93
C TERMCO1A.94
C TERMCO1A.95
C----------------------------------------------------------------------- TERMCO1A.96
C VARIABLES THAT ARE OUTPUT TERMCO1A.97
C----------------------------------------------------------------------- TERMCO1A.98
C TERMCO1A.99
LOGICAL BTERM(NPNTS) ! OUT MASK OF THOSE POINTS AT WHICH TERMCO1A.100
! CONVECTION IS ENDING TERMCO1A.101
C TERMCO1A.102
C TERMCO1A.103
C----------------------------------------------------------------------- TERMCO1A.104
C VARIABLES THAT ARE DEFINED LOCALLY TERMCO1A.105
C----------------------------------------------------------------------- TERMCO1A.106
C TERMCO1A.107
REAL EL ! LATENT HEAT OF CONDENSATION OR TERMCO1A.108
! (CONDENSATION + FUSION) (J/KG) TERMCO1A.109
C TERMCO1A.110
REAL FLXMIN ! MINIMUM CONVECTIVE MASSFLUX BELOW TERMCO1A.111
! WHICH TERMINAL DETRAINMENT OCCURS TERMCO1A.112
! (PA/S) TERMCO1A.113
C TERMCO1A.114
REAL THVUNDI ! POTENTIAL TEMPERATURE OF AN TERMCO1A.115
! UNDILUTE PARCEL IN LAYER K+1 TERMCO1A.116
! FROM THE STARTING LAYER OF TERMCO1A.117
! CONVECTION (K) TERMCO1A.118
C TERMCO1A.119
REAL THVEKP1 ! VIRTUAL POTENTIAL TEMPERATURE TERMCO1A.120
! OF ENVIRONMENT IN LAYER K+1 (K) TERMCO1A.121
C TERMCO1A.122
C*--------------------------------------------------------------------- TERMCO1A.123
C TERMCO1A.124
C---------------------------------------------------------------------- TERMCO1A.125
C CALCULATE MINIMUM MASS FLUX BELOW WHICH CONVECTION IS TERMINATED TERMCO1A.126
C---------------------------------------------------------------------- TERMCO1A.127
C TERMCO1A.128
DO 10 I=1,NPNTS TERMCO1A.129
FLXMIN = MPARFL*(1.+EKP14(I))*(1.+EKP34(I))*PSTAR(I) TERMCO1A.130
C TERMCO1A.131
C----------------------------------------------------------------------- TERMCO1A.132
C CREATE A VECTOR OF LATENT HEATS TERMCO1A.133
C----------------------------------------------------------------------- TERMCO1A.134
C TERMCO1A.135
IF (BWKP1(I)) THEN TERMCO1A.136
EL = LC TERMCO1A.137
ELSE TERMCO1A.138
EL = LC + LF TERMCO1A.139
ENDIF TERMCO1A.140
CL TERMCO1A.141
CL---------------------------------------------------------------------- TERMCO1A.142
CL PARCELS ARE ONLY CONSIDERED FOR TERMINATION IF THEY ARE DETRAINING TERMCO1A.143
CL EXCEPT AT THE TOP MODEL LAYER, WHERE ALL CONVECTION TERMINATES TERMCO1A.144
CL TERMCO1A.145
CL IF THE PARCEL HAS A POTENTIAL TEMPETURE GREATER THAN THE TERMCO1A.146
CL POTENTIAL TEMPERATURE OF AN UNDILUTE PARCEL FORM THE STARTING TERMCO1A.147
CL LAYER OF CONVECION IN LAYER K+1 THEN CONVECTION IS TERMINATED TERMCO1A.148
CL TERMCO1A.149
CL UM DOCUMENTATION PAPER P27 TERMCO1A.150
CL SECTION (7), EQUATION (32) TERMCO1A.151
CL TERMCO1A.152
CL CONVECTION IS ALSO TERMINATED IF MASS FLUX IN LAYER K+1 IS LESS TERMCO1A.153
CL IS LESS THAN A MINIMUM VALUE TERMCO1A.154
CL TERMCO1A.155
CL UM DOCUMENTATION PAPER P27 TERMCO1A.156
CL SECTION (7), EQUATION (33) TERMCO1A.157
CL---------------------------------------------------------------------- TERMCO1A.158
CL TERMCO1A.159
THVUNDI=( THPI(I) + (EL/(EXKP1(I)*CP)) *(QPI(I) - QSEKP1(I)) TERMCO1A.160
* +((LC-EL)/(EXKP1(I)*CP))*MAX(0.0,(QPI(I)-QSTICE)) TERMCO1A.161
* )*(1.+C_VIRTUAL*QSEKP1(I)) TERMCO1A.162
C TERMCO1A.163
THVEKP1 = (THEKP1(I)*(1.+C_VIRTUAL*QEKP1(I)) + XSBMIN) TERMCO1A.164
C TERMCO1A.165
BTERM(I) = (((FLXKP1(I) .LT. FLXMIN) .OR. (THVUNDI .LT. THVEKP1)) TERMCO1A.166
* .AND. (DELTAK(I).GT.0.0)) .OR. (K+1) .EQ. NLEV TERMCO1A.167
C TERMCO1A.168
10 CONTINUE TERMCO1A.169
C TERMCO1A.170
RETURN TERMCO1A.171
END TERMCO1A.172
*ENDIF TERMCO1A.173