anchor.diff:              END IF
   
               ! Self collection
  -            CALL accumulateSink(kbdim,klev,nbins,nbins,kk,kk,kk,zccaa,aero,zminusterm,multp=0.5)
  +            CALL accumulateSink(kbdim,klev,nbins,nbins,kk,kk,kk,zccaa,aero,zminus_self,multp=0.5)
            END IF
   
            ! Collection by larger or equal cloud droplet bins
  
positive.code:          END IF
  
           !-- Volume and number concentrations after coagulation update [fxm]
           CALL applyCoag(kbdim,klev,ncld,nspec,kk,ptstep,cloud,zplusterm,zminusterm)
  
        END DO
  
  
negative.code: 
           !-- 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,zminus_self)
  
        END DO
  
  