*IF DEF,A05_2A,OR,DEF,A05_2C,OR,DEF,A05_3B,OR,DEF,A05_3C AJX1F405.181
C ******************************COPYRIGHT****************************** GTS2F400.3043
C (c) CROWN COPYRIGHT 1995, METEOROLOGICAL OFFICE, All Rights Reserved. GTS2F400.3044
C GTS2F400.3045
C Use, duplication or disclosure of this code is subject to the GTS2F400.3046
C restrictions as set forth in the contract. GTS2F400.3047
C GTS2F400.3048
C Meteorological Office GTS2F400.3049
C London Road GTS2F400.3050
C BRACKNELL GTS2F400.3051
C Berkshire UK GTS2F400.3052
C RG12 2SZ GTS2F400.3053
C GTS2F400.3054
C If no contract has been raised with this copy of the code, the use, GTS2F400.3055
C duplication or disclosure of it is strictly prohibited. Permission GTS2F400.3056
C to do so must first be obtained in writing from the Head of Numerical GTS2F400.3057
C Modelling at the above address. GTS2F400.3058
C ******************************COPYRIGHT****************************** GTS2F400.3059
C GTS2F400.3060
CLL SUBROUTINE FLX_INIT----------------------------------------------- FLXINI2A.3
CLL FLXINI2A.4
CLL PURPOSE : CALCULATE INITIAL DOWNDRAUGHT MASSFLUX FLXINI2A.5
CLL FLXINI2A.6
CLL SUITABLE FOR SINGLE COLUMN MODEL USE FLXINI2A.7
CLL FLXINI2A.8
CLL CODE WRITTEN FOR CRAY Y-MP BY S.BETT AND D.GREGORY SUMMER 1992 FLXINI2A.9
CLL FLXINI2A.10
CLL MODEL MODIFICATION HISTORY FROM MODEL VERSION 3.0: FLXINI2A.11
CLL VERSION DATE FLXINI2A.12
CLL 3.3 23/12/93 : DG020893 : DUE TO CHANGE IN WAY CLOUD TOP IS DG020893.53
CLL ESTIMATED BECAUSE OF CHANGES TO THE DG020893.54
CLL CALCULATION OF FORCED DETRAINMENT DG020893.55
CLL DG020893.56
CLL FLXINI2A.13
CLL PROGRAMMING STANDARDS : UNIFIED MODEL DOCUMENTATION PAPER NO. 3 FLXINI2A.14
CLL VERSION NO. 4 DATED 5/2/92 FLXINI2A.15
CLL FLXINI2A.16
CLL SYSTEM TASK : P27 FLXINI2A.17
CLL FLXINI2A.18
CLL DOCUMENTATION : UNIFIED MODEL DOCUMENTATION PAPER P27 FLXINI2A.19
CLL FLXINI2A.20
CLLEND----------------------------------------------------------------- FLXINI2A.21
C FLXINI2A.22
C*L ARGUMENTS--------------------------------------------------------- FLXINI2A.23
C FLXINI2A.24
SUBROUTINE FLX_INIT (NPNTS,KCT,ICCB,ICCT,FLX,FLX_DD_K,BDDI, 4FLXINI2A.25
* FLX_STRT) FLXINI2A.26
C FLXINI2A.27
IMPLICIT NONE FLXINI2A.28
C FLXINI2A.29
C---------------------------------------------------------------------- FLXINI2A.30
C VECTOR LENGTHS AND LOOP COUNTERS FLXINI2A.31
C---------------------------------------------------------------------- FLXINI2A.32
C FLXINI2A.33
INTEGER I ! LOOP COUNTER FLXINI2A.34
C FLXINI2A.35
INTEGER NPNTS ! IN NUMBER OF POINTS FLXINI2A.36
C FLXINI2A.37
INTEGER KCT ! IN CONVECTIVE CLOUD TOP FLXINI2A.38
C FLXINI2A.39
C---------------------------------------------------------------------- FLXINI2A.40
C VARIABLES WHICH ARE INPUT FLXINI2A.41
C---------------------------------------------------------------------- FLXINI2A.42
C FLXINI2A.43
INTEGER ICCB(NPNTS) ! IN CONVECTIVE CLOUD BASE FLXINI2A.44
C FLXINI2A.45
INTEGER ICCT(NPNTS) ! IN CONVECTIVE CLOUD TOP FLXINI2A.46
C FLXINI2A.47
REAL FLX(NPNTS,KCT+1) ! IN CONVECTIVE MASSFLUX (PA/S) FLXINI2A.48
C FLXINI2A.49
LOGICAL BDDI(NPNTS) ! IN MASK FOR THOSE POINTS WHERE FLXINI2A.50
! DOWNDRAUGHT MAY INITIATE FLXINI2A.51
C FLXINI2A.52
C---------------------------------------------------------------------- FLXINI2A.53
C VARIABLES WHICH ARE OUTPUT FLXINI2A.54
C---------------------------------------------------------------------- FLXINI2A.55
C FLXINI2A.56
REAL FLX_DD_K(NPNTS) ! OUT DOWNDRAUGHT MASSFLUX OF LAYER K FLXINI2A.57
! (PA/S) FLXINI2A.58
C FLXINI2A.59
REAL FLX_STRT(NPNTS) ! OUT UPDRAUGHT MASSFLUX AT LEVEL FLXINI2A.60
! DOWNDRAUGHT STARTS (PA/S) FLXINI2A.61
C FLXINI2A.62
C---------------------------------------------------------------------- FLXINI2A.63
C VARIABLES WHICH ARE DEFINED LOCALLY FLXINI2A.64
C---------------------------------------------------------------------- FLXINI2A.65
C FLXINI2A.66
INTEGER KDDREF ! REFERENCE LEVEL FOR DOWNDRAUGHT FLXINI2A.67
! MASSFLUX FLXINI2A.68
C FLXINI2A.69
C---------------------------------------------------------------------- FLXINI2A.70
C CALCULATE DOWNDRAUGHT MASSFLUX BASED ON A REFERENCE LEVEL WHICH IS FLXINI2A.71
C 3/4 CLOUD DEPTH FLXINI2A.72
C---------------------------------------------------------------------- FLXINI2A.73
C FLXINI2A.74
DO I=1,NPNTS FLXINI2A.75
IF (BDDI(I)) THEN FLXINI2A.76
KDDREF = ICCB(I) + 0.75*(ICCT(I) - ICCB(I)) FLXINI2A.77
IF (KDDREF .GE. ICCT(I)-1) KDDREF=ICCT(I)-1 DG020893.57
FLX_STRT(I) = FLX(I,KDDREF) FLXINI2A.79
FLX_DD_K(I) = FLX_STRT(I) * 0.05 FLXINI2A.80
END IF FLXINI2A.81
END DO FLXINI2A.82
C FLXINI2A.83
RETURN FLXINI2A.84
END FLXINI2A.85
C FLXINI2A.86
*ENDIF FLXINI2A.87