anchor.diff:  !
         IF(fdbk.ge.1.0e-9) THEN
         DO 20 K=1,km
  -      DO 20 j=1,lq
  +      DO j=1,lq
         If(PCTE(j,k).GT.0.0) then
           ZTPP1=pt(j,k)+PTTE(j,k)*ZTMST
           if(ZTPP1.ge.t000) then
  
positive.code: !              CHECK FOR BUOYANCY AND SET FLAGS
  !-------------------------------------------------------
        DO 200 JK=1,KLEV
        DO 200 JL=1,KLON
          ZQOLD(JL,JK)=0.0
    200 CONTINUE
        DO 290 JK=KLEVM1,2,-1
          IS=0
  
negative.code:           ZQSAT(j,k)=ZQSAT(j,k)/(1.-VTMPC1*ZQSAT(j,k))
            PQTE(j,k)=pqvf(j,k)+pqvbl(j,k)
            ZQQ(j,k)=PQTE(j,k)
          end do
     10 CONTINUE
  !
  !-----------------------------------------------------------------------
  