anchor.diff:         zplusterm(:,:,:) = 0.
          
          ! corresponding index for regime b cloud droplets
  -       cc = MAX(kk-fca%cur+ncld,icb%cur) ! Regime a has more bins than b:
  -                                         ! Set this at minimum to beginnign of b.
  +       index_b = MAX(kk-fca%cur+ncld,icb%cur) ! If regime a has more bins than b:
  +                                              ! Set this at minimum to beginnign of b.
  +
  +       ! Droplets lost to larger aerosol bins
  +       IF (lscgca) THEN
  +          ! Aerosol bins a, parallel aerosol bin plus 1
  +          index_aero = ica%par + kk
  +          CALL accumulateSinkReverse(kbdim,klev,zccca,kk,index_aero,fn2a,aero,ncld,nbins,zminusterm)
  +          ! Aerosol bins b, aerosol bin plus 1
  +          index_aero = icb%par + (index_b-fca%cur)
  +          CALL accumulateSinkReverse(kbdim,klev,zccca,kk,index_aero,fn2b,aero,ncld,nbins,zminusterm)
  +       END IF
   
          ! Droplets lost by those with larger nucleus in regime a
          IF (lscgcc) &
  
positive.code:        
         ! Droplets lost by those with larger nucleus in regime a
         IF (lscgcc) &
              CALL accumulateSink(kbdim,klev,zcccc,kk,cc+1,fca%cur,cloud,ncld,ncld,zminusterm)
         
         ! Droplets lost by collection by rain drops
         IF (lscgpc) &
  
negative.code:        
         ! Droplets lost by those with larger nucleus in regime b
         IF (lscgcc) &
              CALL accumulateSink(kbdim,klev,zcccc,kk,index_b+1,fcb%cur,cloud,ncld,ncld,zminusterm)
         
         ! Droplets lost by collection by rain drops
         IF (lscgpc) &
  