*IF DEF,A05_2A,OR,DEF,A05_2C,OR,DEF,A05_3B AJX1F405.130
C ******************************COPYRIGHT****************************** GTS2F400.5221
C (c) CROWN COPYRIGHT 1995, METEOROLOGICAL OFFICE, All Rights Reserved. GTS2F400.5222
C GTS2F400.5223
C Use, duplication or disclosure of this code is subject to the GTS2F400.5224
C restrictions as set forth in the contract. GTS2F400.5225
C GTS2F400.5226
C Meteorological Office GTS2F400.5227
C London Road GTS2F400.5228
C BRACKNELL GTS2F400.5229
C Berkshire UK GTS2F400.5230
C RG12 2SZ GTS2F400.5231
C GTS2F400.5232
C If no contract has been raised with this copy of the code, the use, GTS2F400.5233
C duplication or disclosure of it is strictly prohibited. Permission GTS2F400.5234
C to do so must first be obtained in writing from the Head of Numerical GTS2F400.5235
C Modelling at the above address. GTS2F400.5236
C ******************************COPYRIGHT****************************** GTS2F400.5237
C GTS2F400.5238
CLL SUBROUTINE LATENT_H----------------------------------------------- LATENT1A.3
CLL LATENT1A.4
CLL PURPOSE : CALCULATES A NEW PARCEL TEMPERATURE AFTER LATENT1A.5
CLL CONDENSATION HAS OCCURRED LATENT1A.6
CLL SUITABLE FOR SINGLE COLUMN MODEL USE LATENT1A.7
CLL LATENT1A.8
CLL CODE REWORKED FOR CRAY Y-MP BY D.GREGORY AUTUMN/WINTER 1989/90 LATENT1A.9
CLL LATENT1A.10
CLL MODEL MODIFICATION HISTORY FROM MODEL VERSION 3.0: LATENT1A.11
CLL VERSION DATE LATENT1A.12
CLL LATENT1A.13
CLL PROGRAMMING STANDARDS : UNIFIED MODEL DOCUMENTATION PAPER NO. 4 LATENT1A.14
CLL VERSION NO. 1 LATENT1A.15
CLL LATENT1A.16
CLL LOGICAL COMPONENTS COVERED: P27 LATENT1A.17
CLL LATENT1A.18
CLL DOCUMENTATION : UNIFIED MODEL DOCUMENTATION PAPER P27 LATENT1A.19
CLL LATENT1A.20
CLLEND----------------------------------------------------------------- LATENT1A.21
C LATENT1A.22
C*L ARGUMENTS--------------------------------------------------------- LATENT1A.23
C LATENT1A.24
SUBROUTINE LATENT_H (NPNTS,THPKP1,QPKP1,THEKP1, 2LATENT1A.25
* QSEKP1,DQSKP1,BGMKP1,BWKP1,EXKP1) LATENT1A.26
C LATENT1A.27
C----------------------------------------------------------------------- LATENT1A.28
C ADJUSTS PARCEL POTENTIAL TEMPERATURES TO ACCOUNT FOR THE LATENT LATENT1A.29
C HEAT RELEASE FROM SUPERSATURATED PARCELS CONDENSING/DEPOSITING OUT LATENT1A.30
C----------------------------------------------- AUTHOR: M FISHER 1987 - LATENT1A.31
C LATENT1A.32
IMPLICIT NONE LATENT1A.33
C LATENT1A.34
C---------------------------------------------------------------------- LATENT1A.35
C MODEL CONSTANTS LATENT1A.36
C---------------------------------------------------------------------- LATENT1A.37
C LATENT1A.38
*CALL C_R_CP 
LATENT1A.39
*CALL C_LHEAT LATENT1A.40
C LATENT1A.41
C---------------------------------------------------------------------- LATENT1A.42
C VECTOR LENGTHS AND LOOP COUNTERS LATENT1A.43
C---------------------------------------------------------------------- LATENT1A.44
C LATENT1A.45
INTEGER NPNTS ! VECTOR LENGTHS LATENT1A.46
C LATENT1A.47
INTEGER I ! LOOP COUNTER LATENT1A.48
C LATENT1A.49
C LATENT1A.50
C---------------------------------------------------------------------- LATENT1A.51
C VARIABLES WHICH ARE INPUT LATENT1A.52
C---------------------------------------------------------------------- LATENT1A.53
C LATENT1A.54
REAL QPKP1(NPNTS) ! IN PARCEL MIXING RATIO IN LAYER K+1 LATENT1A.55
! (KG/KG/S) LATENT1A.56
C LATENT1A.57
REAL THEKP1(NPNTS) ! IN ENVIRONMENT POTENTIAL FOR LAYER K+1 LATENT1A.58
! (K/S) LATENT1A.59
C LATENT1A.60
REAL QSEKP1(NPNTS) ! IN SATURATION MIXING RATIO OF THE LATENT1A.61
! ENVIRONMENT IN LAYER K+1 (KG/KG/S) LATENT1A.62
C LATENT1A.63
REAL DQSKP1(NPNTS) ! IN GRADIENT OF SATURATION MIXING RATIO LATENT1A.64
! WITH POTENTIAL TEMPERATURE FOR THE LATENT1A.65
! ENVIRONMENT IN LAYER K+1 (KG/KG/K) LATENT1A.66
C LATENT1A.67
LOGICAL BGMKP1(NPNTS) ! IN MASK FOR PARCELS WHICH ARE LATENT1A.68
! SATURATED LAYER K+1 LATENT1A.69
C LATENT1A.70
LOGICAL BWKP1(NPNTS) ! IN MASK FOR POINTS IN WHICH LATENT1A.71
! CONDENSATE IS LIQUID IN LAYER K+1 LATENT1A.72
C LATENT1A.73
REAL EXKP1(NPNTS) ! IN EXNER RATIO AT MID-POINT OF LATENT1A.74
! LAYER K+1 LATENT1A.75
C LATENT1A.76
C LATENT1A.77
C---------------------------------------------------------------------- LATENT1A.78
C VARIABLES WHICH ARE INPUT AND OUTPUT LATENT1A.79
C---------------------------------------------------------------------- LATENT1A.80
C LATENT1A.81
REAL THPKP1(NPNTS) ! INOUT LATENT1A.82
! IN INITIAL ESTIMATE OF PARCEL LATENT1A.83
! POTENTIAL TEMPERATURE IN LATENT1A.84
! LAYER K+1 (K/S) LATENT1A.85
! OUT PARCEL POTENTIAL TEMPERATURE LATENT1A.86
! IN LAYER K+1 AFTER LATENT LATENT1A.87
! HEATING (K/S) LATENT1A.88
C LATENT1A.89
C LATENT1A.90
C---------------------------------------------------------------------- LATENT1A.91
C VARIABLES WHICH ARE DEFINED LOCALLY LATENT1A.92
C---------------------------------------------------------------------- LATENT1A.93
C LATENT1A.94
REAL EL ! LATENT HEAT OF CONDENSATION OR LATENT1A.95
! (CONDENSATION + FUSION) (J/KG) LATENT1A.96
C LATENT1A.97
C*---------------------------------------------------------------------- LATENT1A.98
CL LATENT1A.99
CL---------------------------------------------------------------------- LATENT1A.100
CL ADJUST PARCEL POTENTIAL TEMPERATURES TO ACCOUNT FOR LATENT HEATING LATENT1A.101
CL LATENT1A.102
CL UM DOCUMENTATION PAPER P27 LATENT1A.103
CL SECTION (4), EQUATION(21) LATENT1A.104
CL---------------------------------------------------------------------- LATENT1A.105
CL LATENT1A.106
DO 10 I=1,NPNTS LATENT1A.107
IF (BWKP1(I)) THEN LATENT1A.108
EL = LC LATENT1A.109
ELSE LATENT1A.110
EL = LC + LF LATENT1A.111
ENDIF LATENT1A.112
C LATENT1A.113
IF (BGMKP1(I)) THEN LATENT1A.114
THPKP1(I) = ( THPKP1(I) + LATENT1A.115
* (EL/(EXKP1(I)*CP))*(QPKP1(I)-QSEKP1(I)+THEKP1(I)*DQSKP1(I)) LATENT1A.116
* ) / (1.+(EL/(EXKP1(I)*CP))*DQSKP1(I)) LATENT1A.117
ENDIF LATENT1A.118
10 CONTINUE LATENT1A.119
C LATENT1A.120
RETURN LATENT1A.121
END LATENT1A.122
*ENDIF LATENT1A.123