anchor.diff:              CALL part(ii,jj,itrgt)%updateRhomean()
   
               ! Apply the sink term for number concentration
  -            part(ii,jj,itrgt)%numc = part(ii,jj,itrgt)%numc / ( 1. + ptstep*sink(ii,jj) )
  +            part(ii,jj,itrgt)%numc = part(ii,jj,itrgt)%numc / ( 1. + ptstep*(sink(ii,jj) + sink_self(ii,jj)) )
   
            END DO
         END DO
  
positive.code:                           zccpa(kbdim,klev,nbins,nprc), zccia(kbdim,klev,nbins,nice),     &
                            zccsa(kbdim,klev,nbins,nsnw)
        INTEGER :: kk
        REAL :: zplusterm(nspec,kbdim,klev), zminusterm(kbdim,klev)
  
        INTEGER :: index_cld_a, index_cld_b
        INTEGER :: index_a, index_b
  
negative.code:       INTEGER :: iwa
        
        iwa = spec%getIndex("H2O")
  
        DO jj = 1,klev
           DO ii = 1,kbdim
  