*IF DEF,A05_2A,OR,DEF,A05_2C,OR,DEF,A05_3B AJX1F405.177
C ******************************COPYRIGHT****************************** GTS2F400.2287
C (c) CROWN COPYRIGHT 1995, METEOROLOGICAL OFFICE, All Rights Reserved. GTS2F400.2288
C GTS2F400.2289
C Use, duplication or disclosure of this code is subject to the GTS2F400.2290
C restrictions as set forth in the contract. GTS2F400.2291
C GTS2F400.2292
C Meteorological Office GTS2F400.2293
C London Road GTS2F400.2294
C BRACKNELL GTS2F400.2295
C Berkshire UK GTS2F400.2296
C RG12 2SZ GTS2F400.2297
C GTS2F400.2298
C If no contract has been raised with this copy of the code, the use, GTS2F400.2299
C duplication or disclosure of it is strictly prohibited. Permission GTS2F400.2300
C to do so must first be obtained in writing from the Head of Numerical GTS2F400.2301
C Modelling at the above address. GTS2F400.2302
C ******************************COPYRIGHT****************************** GTS2F400.2303
C GTS2F400.2304
CLL SUBROUTINE DQSDTH------------------------------------------------- DQSDTH1A.3
CLL DQSDTH1A.4
CLL PURPOSE : CALCULATES GARDIENT OF SATURATION MIXING RATIO DQSDTH1A.5
CLL WITH POTENTIAL TEMPERATURE FORM THE DQSDTH1A.6
CLL CLAUSIUS-CLAPEYRON EQUATION DQSDTH1A.7
CLL DQSDTH1A.8
CLL SUITABLE FOR SINGLE COLUMN MODEL USE DQSDTH1A.9
CLL DQSDTH1A.10
CLL CODE REWORKED FOR CRAY Y-MP BY D.GREGORY AUTUMN/WINTER 1989/90 DQSDTH1A.11
CLL DQSDTH1A.12
CLL MODEL MODIFICATION HISTORY FROM MODEL VERSION 3.0: DQSDTH1A.13
CLL VERSION DATE DQSDTH1A.14
CLL DQSDTH1A.15
CLL PROGRAMMING STANDARDS : UNIFIED MODEL DOCUMENTATION PAPER NO. 4 DQSDTH1A.16
CLL VERSION NO. 1 DQSDTH1A.17
CLL DQSDTH1A.18
CLL LOGICAL COMPONENTS COVERED: P27 DQSDTH1A.19
CLL DQSDTH1A.20
CLL SYSTEM TASK : DQSDTH1A.21
CLL DQSDTH1A.22
CLL DOCUMENTATION : UNIFIED MODEL DOCUMENTATION PAPER P27 DQSDTH1A.23
CLL SECTION(4), EQUATION (20) DQSDTH1A.24
CLL DQSDTH1A.25
CLLEND----------------------------------------------------------------- DQSDTH1A.26
C DQSDTH1A.27
C*L ARGUMENTS--------------------------------------------------------- DQSDTH1A.28
C DQSDTH1A.29
SUBROUTINE DQS_DTH (DQS,K,THEK,QSEK,EXK,NPNTS) 12DQSDTH1A.30
IMPLICIT NONE DQSDTH1A.31
C DQSDTH1A.32
C---------------------------------------------------------------------- DQSDTH1A.33
C MODEL CONSTANTS DQSDTH1A.34
C---------------------------------------------------------------------- DQSDTH1A.35
C DQSDTH1A.36
*CALL RV 
DQSDTH1A.37
*CALL C_LHEAT DQSDTH1A.38
C DQSDTH1A.39
C---------------------------------------------------------------------- DQSDTH1A.40
C VECTOR LENGTHS AND LOOP COUNTERS DQSDTH1A.41
C---------------------------------------------------------------------- DQSDTH1A.42
C DQSDTH1A.43
INTEGER NPNTS ! IN VECTOR LENGTH DQSDTH1A.44
C DQSDTH1A.45
INTEGER K ! IN PRESENT MODEL LAYER DQSDTH1A.46
C DQSDTH1A.47
INTEGER I ! LOOP COUNTER DQSDTH1A.48
C DQSDTH1A.49
C DQSDTH1A.50
C---------------------------------------------------------------------- DQSDTH1A.51
C VARIABLES WHICH ARE INPUT DQSDTH1A.52
C---------------------------------------------------------------------- DQSDTH1A.53
C DQSDTH1A.54
REAL THEK(NPNTS) ! IN POTENTIAL TEMPERATURE FOR LAYER K (K) DQSDTH1A.55
C DQSDTH1A.56
REAL QSEK(NPNTS) ! IN SATURATION MIXING RATIO FOR LAYER K (K) DQSDTH1A.57
C DQSDTH1A.58
REAL EXK(NPNTS) ! IN EXNER RATIO FOR LAYER K DQSDTH1A.59
C DQSDTH1A.60
C DQSDTH1A.61
C---------------------------------------------------------------------- DQSDTH1A.62
C VARIABLES WHICH ARE OUTPUT DQSDTH1A.63
C---------------------------------------------------------------------- DQSDTH1A.64
C DQSDTH1A.65
REAL DQS(NPNTS) ! OUT GRADIENT OF SATURATION MIXING RATIO DQSDTH1A.66
! WITH POTENTIAL TEMPERATURE FOR LAYER K DQSDTH1A.67
! (KG/KG/S) DQSDTH1A.68
C DQSDTH1A.69
C DQSDTH1A.70
C---------------------------------------------------------------------- DQSDTH1A.71
C VARIABLES WHICH ARE DEFINED LOCALLY DQSDTH1A.72
C---------------------------------------------------------------------- DQSDTH1A.73
C DQSDTH1A.74
REAL EL ! LATENT HEATING OF CONDENSATION OR DQSDTH1A.75
! (CONDENSATION PLUS HEATING) (J/KG) DQSDTH1A.76
C DQSDTH1A.77
C*--------------------------------------------------------------------- DQSDTH1A.78
CL DQSDTH1A.79
CL--------------------------------------------------------------------- DQSDTH1A.80
CL NO SIGNIFICANT STRUCTURE DQSDTH1A.81
CL--------------------------------------------------------------------- DQSDTH1A.82
CL DQSDTH1A.83
DO 10 I=1,NPNTS DQSDTH1A.84
C DQSDTH1A.85
C----------------------------------------------------------------------- DQSDTH1A.86
C CREATE A VECTOR OF LATENT HEATS ACCORDING TO WHETHER QSAT IS WRT DQSDTH1A.87
C ICE OR WATER DQSDTH1A.88
C----------------------------------------------------------------------- DQSDTH1A.89
C DQSDTH1A.90
IF (THEK(I) .GT. 273./EXK(I)) THEN DQSDTH1A.91
EL = LC DQSDTH1A.92
ELSE DQSDTH1A.93
EL = LC + LF DQSDTH1A.94
ENDIF DQSDTH1A.95
C DQSDTH1A.96
C----------------------------------------------------------------------- DQSDTH1A.97
C CALCULATE D(QSAT)/D(THETA) DQSDTH1A.98
C----------------------------------------------------------------------- DQSDTH1A.99
C DQSDTH1A.100
DQS(I) = EL*QSEK(I)/(EXK(I)*RV*THEK(I)*THEK(I)) DQSDTH1A.101
10 CONTINUE DQSDTH1A.102
C DQSDTH1A.103
RETURN DQSDTH1A.104
END DQSDTH1A.105
*ENDIF DQSDTH1A.106