anchor.diff:  
   !----------------------------------------------------------------------
   
  -!--- CALL CUMULUS PARAMETERIZATION                                        
  -!                                                                        
  -!...TST IS THE NUMBER OF TIME STEPS IN 10 MINUTES...W0AVG IS CLOSE TO A    
  -!...RUNNING MEAN VERTICAL VELOCITY...NOTE THAT IF YOU CHANGE TST, IT WIL  
  -!...CHANGE THE FREQUENCY OF THE CONVECTIVE INTITIATION CHECK (SEE BELOW) 
  +!--- CALL CUMULUS PARAMETERIZATION
  +!
  +!...TST IS THE NUMBER OF TIME STEPS IN 10 MINUTES...W0AVG IS CLOSE TO A
  +!...RUNNING MEAN VERTICAL VELOCITY...NOTE THAT IF YOU CHANGE TST, IT WIL
  +!...CHANGE THE FREQUENCY OF THE CONVECTIVE INTITIATION CHECK (SEE BELOW)
   !...NOTE THAT THE ORDERING OF VERTICAL LAYERS MUST BE REVERSED FOR W0AVG
  -!...BECAUSE THE ORDERING IS REVERSED IN KFPARA...                         
  -!                                                                           
  +!...BECAUSE THE ORDERING IS REVERSED IN KFPARA...
  +!
      DXSQ=DX*DX
      qi_flag = .FALSE.
      qs_flag = .FALSE.
  
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:                  DQDT,DQIDT,DQCDT,DQRDT,DQSDT,DTDT, &
                   RAINCV,PRATEC,NCA,                 &
                   warm_rain,qi_flag,qs_flag,         &
                   ids,ide, jds,jde, kds,kde,         &
                   ims,ime, jms,jme, kms,kme,         &
                   its,ite, jts,jte, kts,kte          )
  
              IF ( PRESENT( RTHCUTEN ) .AND. PRESENT( RQVCUTEN ) ) THEN
                DO K=kts,kte
                   RTHCUTEN(I,K,J)=DTDT(K)/pi(I,K,J)
  