anchor.diff:           IF ( ALL(cloud(:,:,kk)%numc < cloud(:,:,kk)%nlim) ) CYCLE
   
            zminusterm(:,:) = 0.
  +         zminus_self(:,:) = 0.
            zplusterm(:,:,:) = 0.
   
            ! Corresponding index in the regime b droplets
  
positive.code:          !-- Volume and number concentrations after coagulation update. Use the Ice variant, which 
           !   calculates the volume mean densities etc correctly for the pristine and rimed ice contributions.
           CALL applyCoagIce(kbdim,klev,nice,nspec,ndry,kk,ptstep,ice,      &
                             rhoic,rhorime,zplusterm,zplus_rime,zminusterm  )
  
        END DO
  
  
negative.code:          IF ( ALL(cloud(:,:,kk)%numc < cloud(:,:,kk)%nlim) ) CYCLE
  
           zminusterm(:,:) = 0.
           zminus_self(:,:) = 0.
           zplusterm(:,:,:) = 0.
  
           ! Corresponding index in the regime b droplets
  