anchor.diff:  
          remcld(:,:,:) = cmdep(:,:,:)
   
  -       ! Account for changes in in liquid water pot temperature
  +       ! Account for changes in liquid water pot temperature
          nc = GetIndex(prtcl,'H2O')
          istr = (nc-1)*ncld+1
          iend = nc*ncld
          DO j = 3,n3-2
             DO i = 3,n2-2
  -             DO k = 1,n1
  +             DO k = 2,n1
                   tlt(k,i,j) = tlt(k,i,j) + SUM(cmdiv(k,i,j,istr:iend))*(alvl/cp)*th(k,i,j)/tk(k,i,j)
                END DO
  +             tlt(2,i,j) = tlt(2,i,j) + SUM(cmdep(i,j,istr:iend))*(alvl/cp)*th(2,i,j)/tk(2,i,j)
             END DO
          END DO
   
  
positive.code: 
         remaer(:,:,:) = amdep(:,:,:)
  
         ! Account for changes in in liquid water pot temperature
         nc = GetIndex(prtcl,'H2O')
         istr = (nc-1)*nbins+1
         iend = nc*nbins
         DO j = 3,n3-2
            DO i = 3,n2-2
               DO k = 1,n1
                  tlt(k,i,j) = tlt(k,i,j) + SUM(amdiv(k,i,j,istr:iend))*(alvl/cp)*th(k,i,j)/tk(k,i,j)
               END DO
            END DO
         END DO
  
  
negative.code: 
         remice(:,:,:) = imdep(:,:,:)
  
         ! Account for changes in liquid water pot temperature
         nc = GetIndex(prtcl,'H2O')
         istr = (nc-1)*nice+1
         iend = nc*nice
         DO j = 3,n3-2
            DO i = 3,n2-2
               DO k = 2,n1
                  tlt(k,i,j) = tlt(k,i,j) + SUM(imdiv(k,i,j,istr:iend))*(alvi/cp)*th(k,i,j)/tk(k,i,j)
               END DO
               tlt(2,i,j) = tlt(2,i,j) + SUM(imdep(i,j,istr:iend))*(alvi/cp)*th(2,i,j)/tk(2,i,j)
            END DO
         END DO
  
  