Ignore:
Timestamp:
Jan 30, 2020, 10:36:35 AM (5 years ago)
Author:
aboissinot
Message:

The thermal plume model is able to manage several plumes in the same column and
work without the convective adjustment.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.GENERIC/libf/phystd/convadj.F

    r2127 r2232  
    33     &                   pu,pv,ph,pq,
    44     &                   pdufi,pdvfi,pdhfi,pdqfi,
    5      &                   pduadj,pdvadj,pdhadj,
    6      &                   pdqadj,lmax)
     5     &                   pduadj,pdvadj,pdhadj,pdqadj)
    76
    87      USE tracer_h
     
    4140      REAL pu(ngrid,nlay),pdufi(ngrid,nlay),pduadj(ngrid,nlay)
    4241      REAL pv(ngrid,nlay),pdvfi(ngrid,nlay),pdvadj(ngrid,nlay)
    43       INTEGER lmax(ngrid)
    4442
    4543!     Tracers
     
    157155      DO l=2,nlay
    158156        DO ig=1,ngrid
    159           IF (zhc(ig,l).LT.zhc(ig,l-1).and.(l.GT.lmax(ig))) THEN
     157          IF (zhc(ig,l).LT.zhc(ig,l-1)) THEN
    160158            vtest(ig)=.true.
    161159          ENDIF
     
    197195            l2 = l2 + 1
    198196            IF (l2 .GT. nlay) EXIT
    199             IF ((zhc(i, l2).LT.zhc(i, l2-1)).and.(l2.GT.lmax(i))) THEN
     197            IF (zhc(i, l2).LT.zhc(i, l2-1)) THEN
    200198 
    201199!     l2 is the highest level of the unstable column
     
    226224                down = .false.
    227225                IF (l1 .ne. 1) then    !--  and then
    228                   IF ((zhmc.LT.zhc(i, l1-1)).and.(l1.GT.lmax(i))) then
     226                  IF (zhmc.LT.zhc(i, l1-1)) then
    229227                    down = .true.
    230228                  END IF
Note: See TracChangeset for help on using the changeset viewer.