anchor.diff:  
            !-- Volume and number concentrations after coagulation update. Again the Ice variant. 
            CALL applyCoagIce(kbdim,klev,nice,nspec,ndry,kk,ptstep,ice,     &
  -                           rhoic,rhorime,zplusterm,zplus_rime,zminusterm )
  +                           rhoic,rhorime,zplusterm,zplus_rime,zminusterm,zminus_self)
   
         END DO
   
  
positive.code:                CALL accumulateSink(kbdim,klev,ncld,ncld,kk,kk+1,fcb%cur,zcccc,cloud,zminusterm)
                 CALL accumulateSink(kbdim,klev,ncld,ncld,kk,index_a+1,fca%cur,zcccc,cloud,zminusterm)
              END IF
              CALL accumulateSink(kbdim,klev,ncld,ncld,kk,kk,kk,zcccc,cloud,zminusterm,multp=0.5)
           END IF
           
           ! Collection by precipitation
  
negative.code:          END IF
  
           !-- Volume and number concentrations after coagulation update [fxm]
           CALL applyCoag(kbdim,klev,nbins,nspec,kk,ptstep,aero,zplusterm,zminusterm,zminus_self)
  
        END DO !kk
  
  