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