anchor.diff:     ENDDO
      lastdt = dt
   
  -     DO J = jts,jte  
  +     DO J = jts,jte
        DO I= its,ite
           CU_ACT_FLAG(i,j) = .true.
        ENDDO
        ENDDO
   
  -     DO J = jts,jte  
  +     DO J = jts,jte
           DO I=its,ite
   ! if (i.eq. 110 .and. j .eq. 59 ) then
   !   write(0,*) 'nca = ',nca(i,j),' CU_ACT_FLAG = ',CU_ACT_FLAG(i,j)
  
positive.code:             QLFXOUT,QIPA,QIFXIN,QIFXOUT,QRPA,              &
              QRFXIN,QRFXOUT,QSPA,QSFXIN,QSFXOUT,            &
              QL0,QLG,QI0,QIG,QR0,QRG,QS0,QSG
                                           
        REAL, DIMENSION( kts:kte+1 ) :: OMG
        REAL, DIMENSION( kts:kte ) :: RAINFB,SNOWFB
  
  
negative.code:                                                    PRATEC, &
                                                        NCA
  !
  !...DEFINE LOCAL VARIABLES...
  !
        REAL, DIMENSION( kts:kte ) ::                        &
              Q0,Z0,TV0,TU,TVU,QU,TZ,TVD,                    &
              QD,QES,THTES,TG,TVG,QG,WU,WD,W0,EMS,EMSD,      &
  