/[lmdze]/trunk/phylmd/Conflx/flxasc.f90
ViewVC logotype

Diff of /trunk/phylmd/Conflx/flxasc.f90

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 77 by guez, Fri Nov 15 18:45:49 2013 UTC revision 78 by guez, Wed Feb 5 17:51:07 2014 UTC
# Line 24  contains Line 24  contains
24      REAL, intent(in):: pqen(klon, klev)      REAL, intent(in):: pqen(klon, klev)
25      REAL, intent(in):: pqsen(klon, klev)      REAL, intent(in):: pqsen(klon, klev)
26      REAL, intent(in):: pgeo(klon, klev), pgeoh(klon, klev)      REAL, intent(in):: pgeo(klon, klev), pgeoh(klon, klev)
27      REAL pap(klon, klev), paph(klon, klev+1)      REAL, intent(in):: pap(klon, klev), paph(klon, klev+1)
28      REAL pqte(klon, klev)      REAL, intent(in):: pqte(klon, klev)
29      REAL pvervel(klon, klev) ! vitesse verticale en Pa/s      REAL, intent(in):: pvervel(klon, klev) ! vitesse verticale en Pa/s
30      LOGICAL ldland(klon)      LOGICAL, intent(in):: ldland(klon)
31      LOGICAL, intent(inout):: ldcum(klon)      LOGICAL, intent(inout):: ldcum(klon)
32      INTEGER, intent(inout):: ktype(klon)      INTEGER, intent(inout):: ktype(klon)
33      integer klab(klon, klev)      integer klab(klon, klev)
# Line 59  contains Line 59  contains
59    
60      REAL zwmax(klon), zzzmb      REAL zwmax(klon), zzzmb
61      INTEGER klwmin(klon) ! level of maximum vertical velocity      INTEGER klwmin(klon) ! level of maximum vertical velocity
62        real fact
63    
64      !----------------------------------------------------------------------      !----------------------------------------------------------------------
65    
# Line 95  contains Line 96  contains
96            pmful(i, k)=0.            pmful(i, k)=0.
97            plude(i, k)=0.            plude(i, k)=0.
98            pdmfup(i, k)=0.            pdmfup(i, k)=0.
99            IF(.NOT. ldcum(i).OR.ktype(i) == 3) klab(i, k)=0            IF (.NOT. ldcum(i) .OR. ktype(i) == 3) klab(i, k)=0
100            IF(.NOT. ldcum(i).AND.paph(i, k) < 4.E4) kctop0(i)=k            IF (.NOT. ldcum(i) .AND. paph(i, k) < 4e4) kctop0(i) = k
101         ENDDO         ENDDO
102      ENDDO      ENDDO
103    
# Line 142  contains Line 143  contains
143                  pqu(i, k+1) = pqen(i, k)                  pqu(i, k+1) = pqen(i, k)
144                  plu(i, k+1) = 0.0                  plu(i, k+1) = 0.0
145                  zzzmb = MAX(CMFCMIN, -pvervel(i, k)/RG)                  zzzmb = MAX(CMFCMIN, -pvervel(i, k)/RG)
146                  zmfmax = (paph(i, k)-paph(i, k-1))/(RG * pdtime)                  zmfmax = (paph(i, k) - paph(i, k-1)) / (RG * pdtime)
147                  pmfub(i) = MIN(zzzmb, zmfmax)                  pmfub(i) = MIN(zzzmb, zmfmax)
148                  pmfu(i, k+1) = pmfub(i)                  pmfu(i, k+1) = pmfub(i)
149                  pmfus(i, k+1) = pmfub(i) * (RCPD * ptu(i, k+1)+pgeoh(i, k+1))                  pmfus(i, k+1) = pmfub(i) * (RCPD * ptu(i, k+1)+pgeoh(i, k+1))
# Line 171  contains Line 172  contains
172         DO i = 1, klon         DO i = 1, klon
173            pen_u(i, k) = 0.0            pen_u(i, k) = 0.0
174            pde_u(i, k) = 0.0            pde_u(i, k) = 0.0
175            zrho(i)=paph(i, k+1)/(RD * ptenh(i, k+1))            zrho(i) = paph(i, k + 1) / (RD * ptenh(i, k + 1))
176            zpbot(i)=paph(i, kcbot(i))            zpbot(i) = paph(i, kcbot(i))
177            zptop(i)=paph(i, kctop0(i))            zptop(i) = paph(i, kctop0(i))
178         ENDDO         ENDDO
179    
180         DO i = 1, klon         DO i = 1, klon
181            IF(ldcum(i)) THEN            IF (ldcum(i)) THEN
182               zdprho=(paph(i, k+1)-paph(i, k))/(RG * zrho(i))               zdprho = (paph(i, k + 1) - paph(i, k)) / (RG * zrho(i))
183               zentr=pentr(i) * pmfu(i, k+1) * zdprho               zentr=pentr(i) * pmfu(i, k+1) * zdprho
184               llo1=k < kcbot(i)               llo1=k < kcbot(i)
185               IF(llo1) pde_u(i, k)=zentr               IF (llo1) pde_u(i, k)=zentr
186               zpmid=0.5 * (zpbot(i)+zptop(i))               zpmid=0.5 * (zpbot(i)+zptop(i))
187               llo2=llo1.AND.ktype(i) == 2.AND. &               llo2 = llo1 .AND. ktype(i) == 2 &
188                    (zpbot(i)-paph(i, k) < 0.2E5.OR. &                    .AND. (zpbot(i) - paph(i, k) < 0.2E5 .OR. paph(i, k) > zpmid)
189                    paph(i, k) > zpmid)               IF (llo2) pen_u(i, k)=zentr
190               IF(llo2) pen_u(i, k)=zentr               llo2 = llo1 .AND. (ktype(i) == 1 .OR. ktype(i) == 3) .AND. &
191               llo2=llo1.AND.(ktype(i) == 1.OR.ktype(i) == 3).AND. &                    (k >= MAX(klwmin(i), kctop0(i) + 2) .OR. pap(i, k) > zpmid)
192                    (k >= MAX(klwmin(i), kctop0(i)+2).OR.pap(i, k) > zpmid)               IF (llo2) pen_u(i, k)=zentr
193               IF(llo2) pen_u(i, k)=zentr               llo1=pen_u(i, k) > 0. .AND. (ktype(i) == 1 .OR. ktype(i) == 2)
194               llo1=pen_u(i, k) > 0..AND.(ktype(i) == 1.OR.ktype(i) == 2)               IF (llo1) THEN
195               IF(llo1) THEN                  fact = 1. + 3. * (1. - MIN(1., (zpbot(i) - pap(i, k)) / 1.5E4))
196                  zentr=zentr * (1.+3. * (1.-MIN(1., (zpbot(i)-pap(i, k))/1.5E4)))                  zentr = zentr * fact
197                  pen_u(i, k)=pen_u(i, k) * (1.+3. * (1.-MIN(1., &                  pen_u(i, k)=pen_u(i, k) * fact
198                       (zpbot(i)-pap(i, k))/1.5E4)))                  pde_u(i, k)=pde_u(i, k) * fact
                 pde_u(i, k)=pde_u(i, k) * (1.+3. * (1.-MIN(1., &  
                      (zpbot(i)-pap(i, k))/1.5E4)))  
199               ENDIF               ENDIF
200               IF(llo2.AND.pqenh(i, k+1) > 1.E-5) &               IF (llo2 .AND. pqenh(i, k+1) > 1e-5) &
201                    pen_u(i, k)=zentr+MAX(pqte(i, k), 0.)/pqenh(i, k+1) * &                    pen_u(i, k)=zentr+MAX(pqte(i, k), 0.)/pqenh(i, k+1) * &
202                    zrho(i) * zdprho                    zrho(i) * zdprho
203            ENDIF            ENDIF
# Line 210  contains Line 209  contains
209            IF (llflag(i)) THEN            IF (llflag(i)) THEN
210               IF (k < kcbot(i)) THEN               IF (k < kcbot(i)) THEN
211                  zmftest = pmfu(i, k+1)+pen_u(i, k)-pde_u(i, k)                  zmftest = pmfu(i, k+1)+pen_u(i, k)-pde_u(i, k)
212                  zmfmax = MIN(zmftest, (paph(i, k)-paph(i, k-1))/(RG * pdtime))                  zmfmax = MIN(zmftest, &
213                         (paph(i, k) - paph(i, k - 1)) / (RG * pdtime))
214                  pen_u(i, k)=MAX(pen_u(i, k)-MAX(0.0, zmftest-zmfmax), 0.0)                  pen_u(i, k)=MAX(pen_u(i, k)-MAX(0.0, zmftest-zmfmax), 0.0)
215               ENDIF               ENDIF
216               pde_u(i, k)=MIN(pde_u(i, k), 0.75 * pmfu(i, k+1))               pde_u(i, k)=MIN(pde_u(i, k), 0.75 * pmfu(i, k+1))
# Line 244  contains Line 244  contains
244         CALL flxadjtq(paph(1, k), ptu(1, k), pqu(1, k), llflag, icall)         CALL flxadjtq(paph(1, k), ptu(1, k), pqu(1, k), llflag, icall)
245    
246         DO i = 1, klon         DO i = 1, klon
247            IF(llflag(i).AND.pqu(i, k).NE.zqold(i)) THEN            IF (llflag(i) .AND. pqu(i, k).NE.zqold(i)) THEN
248               klab(i, k) = 2               klab(i, k) = 2
249               plu(i, k) = plu(i, k)+zqold(i)-pqu(i, k)               plu(i, k) = plu(i, k)+zqold(i)-pqu(i, k)
250               zbuo = ptu(i, k) * (1.+RETV * pqu(i, k))- &               zbuo = ptu(i, k) * (1.+RETV * pqu(i, k))- &
# Line 256  contains Line 256  contains
256                  zdnoprc = 1.5E4                  zdnoprc = 1.5E4
257                  IF (ldland(i)) zdnoprc = zdland(i)                  IF (ldland(i)) zdnoprc = zdland(i)
258                  zprcon = CPRCON                  zprcon = CPRCON
259                  IF ((zpbot(i)-paph(i, k)) < zdnoprc) zprcon = 0.0                  IF ((zpbot(i) - paph(i, k)) < zdnoprc) zprcon = 0.
260                  zlnew=plu(i, k)/(1.+zprcon * (pgeoh(i, k)-pgeoh(i, k+1)))                  zlnew=plu(i, k)/(1.+zprcon * (pgeoh(i, k)-pgeoh(i, k+1)))
261                  pdmfup(i, k)=MAX(0., (plu(i, k)-zlnew) * pmfu(i, k))                  pdmfup(i, k)=MAX(0., (plu(i, k)-zlnew) * pmfu(i, k))
262                  plu(i, k)=zlnew                  plu(i, k)=zlnew
# Line 273  contains Line 273  contains
273               pmfuq(i, k)=pqu(i, k) * pmfu(i, k)               pmfuq(i, k)=pqu(i, k) * pmfu(i, k)
274            ENDIF            ENDIF
275         end DO         end DO
   
276      end DO      end DO
277    
278      ! Determine convective fluxes above non-buoyancy level (note:      ! Determine convective fluxes above non-buoyancy level (note:
Legend:
Removed from v.77  
changed lines
  Added in v.78
  ViewVC Help
Powered by ViewVC 1.1.21