anchor.diff:      !***
       !***  CLOUD MUST BE AT LEAST TWO LAYERS THICK
       !***
  -    !        IF(LBOT-LTOP<2)LTOP=LBOT-2  (eliminate this criterion)
  +        !    IF(LBOT-LTOP<2)LTOP=LBOT-2  (eliminate this criterion)
       !
       !-- End: Buoyancy check (24 Aug 2006)
       !
  
positive.code:             APEKXY=APEK(LBOT)
      !
              LMID=.5*(LBOT+LTOP)
      !
              DO L=LBOT,LTOP,-1
              TREFKX=((PKXXXY-PKXXXX)*SMIX                                    &
              &          +TREFKX*APEKXX)/APEKXY
  
negative.code:             IF(L<=LMID) TREFK(L)=MAX(TREFK(L), TK(L)+DTSHAL)
              APEKXX=APEKXY
              PKXXXX=PKXXXY
              APEKXY=APEK(L-1)  ! error in ICAR; index below lower bound of 1
              PKXXXY=PK(L-1)
              ENDDO
      !
  