phylmd/Conflx/flxdlfs.f90

      SUBROUTINE flxdlfs(ptenh, pqenh, pgeoh, paph, ptu, pqu, &
           ldcum, kcbot, kctop, pmfub, prfl, ptd, pqd, &
           pmfd, pmfds, pmfdq, pdmfdp, kdtop, lddraf)
      use dimens_m
      use dimphy
      use SUPHEC_M
      use yoethf_m
            use yoecumf
      IMPLICIT none
!
!----------------------------------------------------------------------
! THIS ROUTINE CALCULATES LEVEL OF FREE SINKING FOR
! CUMULUS DOWNDRAFTS AND SPECIFIES T,Q,U AND V VALUES
!
! TO PRODUCE LFS-VALUES FOR CUMULUS DOWNDRAFTS
! FOR MASSFLUX CUMULUS PARAMETERIZATION
!
! INPUT ARE ENVIRONMENTAL VALUES OF T,Q,U,V,P,PHI
! AND UPDRAFT VALUES T,Q,U AND V AND ALSO
! CLOUD BASE MASSFLUX AND CU-PRECIPITATION RATE.
! IT RETURNS T,Q,U AND V VALUES AND MASSFLUX AT LFS.
!
! CHECK FOR NEGATIVE BUOYANCY OF AIR OF EQUAL PARTS OF
! MOIST ENVIRONMENTAL AIR AND CLOUD AIR.
!----------------------------------------------------------------------
!
      REAL ptenh(klon,klev)
      REAL pqenh(klon,klev)
      REAL pgeoh(klon,klev), paph(klon,klev+1)
      REAL ptu(klon,klev), pqu(klon,klev)
      REAL pmfub(klon)
      REAL prfl(klon)
!
      REAL ptd(klon,klev), pqd(klon,klev)
      REAL pmfd(klon,klev), pmfds(klon,klev), pmfdq(klon,klev)
      REAL pdmfdp(klon,klev)
      INTEGER  kcbot(klon), kctop(klon), kdtop(klon)
      LOGICAL  ldcum(klon), lddraf(klon)
!
      REAL ztenwb(klon,klev), zqenwb(klon,klev), zcond(klon)
      REAL zttest, zqtest, zbuo, zmftop
      LOGICAL  llo2(klon)
      INTEGER i, k, is, icall
!----------------------------------------------------------------------
      DO i= 1, klon
         lddraf(i)=.FALSE.
         kdtop(i)=klev+1
      ENDDO
!
!----------------------------------------------------------------------
! DETERMINE LEVEL OF FREE SINKING BY
! DOING A SCAN FROM TOP TO BASE OF CUMULUS CLOUDS
!
! FOR EVERY POINT AND PROCEED AS FOLLOWS:
!     (1) DETEMINE WET BULB ENVIRONMENTAL T AND Q
!     (2) DO MIXING WITH CUMULUS CLOUD AIR
!     (3) CHECK FOR NEGATIVE BUOYANCY
!
! THE ASSUMPTION IS THAT AIR OF DOWNDRAFTS IS MIXTURE
! OF 50% CLOUD AIR + 50% ENVIRONMENTAL AIR AT WET BULB
! TEMPERATURE (I.E. WHICH BECAME SATURATED DUE TO
! EVAPORATION OF RAIN AND CLOUD WATER)
!----------------------------------------------------------------------
!
      DO 290 k = 3, klev-3
!
      is=0
      DO 212 i= 1, klon
         ztenwb(i,k)=ptenh(i,k)
         zqenwb(i,k)=pqenh(i,k)
         llo2(i) = ldcum(i).AND.prfl(i).GT.0. &
                   .AND..NOT.lddraf(i) &
                   .AND.(k.LT.kcbot(i).AND.k.GT.kctop(i))
         IF ( llo2(i) ) is = is + 1
  212 CONTINUE
      IF(is.EQ.0) GO TO 290
!
      icall=2
      CALL flxadjtq(paph(1,k), ztenwb(1,k), zqenwb(1,k), llo2, icall)
!
!----------------------------------------------------------------------
! DO MIXING OF CUMULUS AND ENVIRONMENTAL AIR
! AND CHECK FOR NEGATIVE BUOYANCY.
! THEN SET VALUES FOR DOWNDRAFT AT LFS.
!----------------------------------------------------------------------
      DO 222 i= 1, klon
      IF (llo2(i)) THEN
         zttest=0.5*(ptu(i,k)+ztenwb(i,k))
         zqtest=0.5*(pqu(i,k)+zqenwb(i,k))
         zbuo=zttest*(1.+RETV*zqtest)- &
              ptenh(i,k)*(1.+RETV  *pqenh(i,k))
         zcond(i)=pqenh(i,k)-zqenwb(i,k)
         zmftop=-CMFDEPS*pmfub(i)
         IF (zbuo.LT.0..AND.prfl(i).GT.10.*zmftop*zcond(i)) THEN
            kdtop(i)=k
            lddraf(i)=.TRUE.
            ptd(i,k)=zttest
            pqd(i,k)=zqtest
            pmfd(i,k)=zmftop
            pmfds(i,k)=pmfd(i,k)*(RCPD*ptd(i,k)+pgeoh(i,k))
            pmfdq(i,k)=pmfd(i,k)*pqd(i,k)
            pdmfdp(i,k-1)=-0.5*pmfd(i,k)*zcond(i)
            prfl(i)=prfl(i)+pdmfdp(i,k-1)
         ENDIF
      ENDIF
  222 CONTINUE
!
  290 CONTINUE
!
      RETURN
      END
1	SUBROUTINE flxdlfs(ptenh, pqenh, pgeoh, paph, ptu, pqu, &
2	ldcum, kcbot, kctop, pmfub, prfl, ptd, pqd, &
3	pmfd, pmfds, pmfdq, pdmfdp, kdtop, lddraf)
4	use dimens_m
5	use dimphy
6	use SUPHEC_M
7	use yoethf_m
8	use yoecumf
9	IMPLICIT none
10	!
11	!----------------------------------------------------------------------
12	! THIS ROUTINE CALCULATES LEVEL OF FREE SINKING FOR
13	! CUMULUS DOWNDRAFTS AND SPECIFIES T,Q,U AND V VALUES
14	!
15	! TO PRODUCE LFS-VALUES FOR CUMULUS DOWNDRAFTS
16	! FOR MASSFLUX CUMULUS PARAMETERIZATION
17	!
18	! INPUT ARE ENVIRONMENTAL VALUES OF T,Q,U,V,P,PHI
19	! AND UPDRAFT VALUES T,Q,U AND V AND ALSO
20	! CLOUD BASE MASSFLUX AND CU-PRECIPITATION RATE.
21	! IT RETURNS T,Q,U AND V VALUES AND MASSFLUX AT LFS.
22	!
23	! CHECK FOR NEGATIVE BUOYANCY OF AIR OF EQUAL PARTS OF
24	! MOIST ENVIRONMENTAL AIR AND CLOUD AIR.
25	!----------------------------------------------------------------------
26	!
27	REAL ptenh(klon,klev)
28	REAL pqenh(klon,klev)
29	REAL pgeoh(klon,klev), paph(klon,klev+1)
30	REAL ptu(klon,klev), pqu(klon,klev)
31	REAL pmfub(klon)
32	REAL prfl(klon)
33	!
34	REAL ptd(klon,klev), pqd(klon,klev)
35	REAL pmfd(klon,klev), pmfds(klon,klev), pmfdq(klon,klev)
36	REAL pdmfdp(klon,klev)
37	INTEGER kcbot(klon), kctop(klon), kdtop(klon)
38	LOGICAL ldcum(klon), lddraf(klon)
39	!
40	REAL ztenwb(klon,klev), zqenwb(klon,klev), zcond(klon)
41	REAL zttest, zqtest, zbuo, zmftop
42	LOGICAL llo2(klon)
43	INTEGER i, k, is, icall
44	!----------------------------------------------------------------------
45	DO i= 1, klon
46	lddraf(i)=.FALSE.
47	kdtop(i)=klev+1
48	ENDDO
49	!
50	!----------------------------------------------------------------------
51	! DETERMINE LEVEL OF FREE SINKING BY
52	! DOING A SCAN FROM TOP TO BASE OF CUMULUS CLOUDS
53	!
54	! FOR EVERY POINT AND PROCEED AS FOLLOWS:
55	! (1) DETEMINE WET BULB ENVIRONMENTAL T AND Q
56	! (2) DO MIXING WITH CUMULUS CLOUD AIR
57	! (3) CHECK FOR NEGATIVE BUOYANCY
58	!
59	! THE ASSUMPTION IS THAT AIR OF DOWNDRAFTS IS MIXTURE
60	! OF 50% CLOUD AIR + 50% ENVIRONMENTAL AIR AT WET BULB
61	! TEMPERATURE (I.E. WHICH BECAME SATURATED DUE TO
62	! EVAPORATION OF RAIN AND CLOUD WATER)
63	!----------------------------------------------------------------------
64	!
65	DO 290 k = 3, klev-3
66	!
67	is=0
68	DO 212 i= 1, klon
69	ztenwb(i,k)=ptenh(i,k)
70	zqenwb(i,k)=pqenh(i,k)
71	llo2(i) = ldcum(i).AND.prfl(i).GT.0. &
72	.AND..NOT.lddraf(i) &
73	.AND.(k.LT.kcbot(i).AND.k.GT.kctop(i))
74	IF ( llo2(i) ) is = is + 1
75	212 CONTINUE
76	IF(is.EQ.0) GO TO 290
77	!
78	icall=2
79	CALL flxadjtq(paph(1,k), ztenwb(1,k), zqenwb(1,k), llo2, icall)
80	!
81	!----------------------------------------------------------------------
82	! DO MIXING OF CUMULUS AND ENVIRONMENTAL AIR
83	! AND CHECK FOR NEGATIVE BUOYANCY.
84	! THEN SET VALUES FOR DOWNDRAFT AT LFS.
85	!----------------------------------------------------------------------
86	DO 222 i= 1, klon
87	IF (llo2(i)) THEN
88	zttest=0.5*(ptu(i,k)+ztenwb(i,k))
89	zqtest=0.5*(pqu(i,k)+zqenwb(i,k))
90	zbuo=zttest(1.+RETVzqtest)- &
91	ptenh(i,k)(1.+RETV pqenh(i,k))
92	zcond(i)=pqenh(i,k)-zqenwb(i,k)
93	zmftop=-CMFDEPS*pmfub(i)
94	IF (zbuo.LT.0..AND.prfl(i).GT.10.zmftopzcond(i)) THEN
95	kdtop(i)=k
96	lddraf(i)=.TRUE.
97	ptd(i,k)=zttest
98	pqd(i,k)=zqtest
99	pmfd(i,k)=zmftop
100	pmfds(i,k)=pmfd(i,k)(RCPDptd(i,k)+pgeoh(i,k))
101	pmfdq(i,k)=pmfd(i,k)*pqd(i,k)
102	pdmfdp(i,k-1)=-0.5pmfd(i,k)zcond(i)
103	prfl(i)=prfl(i)+pdmfdp(i,k-1)
104	ENDIF
105	ENDIF
106	222 CONTINUE
107	!
108	290 CONTINUE
109	!
110	RETURN
111	END