anchor.diff:  
                  ! Volume gained from cloud collection of aerosols
                  DO ll = in1a, fn2b
  -                  zplusterm(1:8) = zplusterm(1:8) + zccca(ll,cc)*aero(ii,jj,ll)%volc(1:8)
  +                  zplusterm(1:nspec) = zplusterm(1:nspec) + zccca(ll,cc)*aero(ii,jj,ll)%volc(1:nspec)
                  END DO
   
                  ! Volume gained from smaller droplets in a
                  DO ll = ica%cur, cc-1
  -                  zplusterm(1:8) = zplusterm(1:8) + zcccc(ll,cc)*cloud(ii,jj,ll)%volc(1:8)
  +                  zplusterm(1:nspec) = zplusterm(1:nspec) + zcccc(ll,cc)*cloud(ii,jj,ll)%volc(1:nspec)
                  END DO
   
                  ! Volume gained from smaller or equal droplets in b
                  DO ll = icb%cur, kk
  -                  zplusterm(1:8) = zplusterm(1:8) + zcccc(ll,cc)*cloud(ii,jj,ll)%volc(1:8)
  +                  zplusterm(1:nspec) = zplusterm(1:nspec) + zcccc(ll,cc)*cloud(ii,jj,ll)%volc(1:nspec)
                  END DO
   
                  ! Update the hydrometeor volume concentrations
  -               cloud(ii,jj,cc)%volc(1:8) = max(0.,( cloud(ii,jj,cc)%volc(1:8) +  &
  -                                                ptstep*zplusterm(1:8)*cloud(ii,jj,cc)%numc ) /  &
  +               cloud(ii,jj,cc)%volc(1:nspec) = max(0.,( cloud(ii,jj,cc)%volc(1:nspec) +  &
  +                                                ptstep*zplusterm(1:nspec)*cloud(ii,jj,cc)%numc ) /  &
                                                   (1. + ptstep*zminusterm) )
   
                  ! Update the hydrometeor number concentration (Removal by coagulation with larger bins and self)
  
positive.code: 
           END DO ! kbdim
        END DO ! klev
  
     END SUBROUTINE coagulation
  
  
negative.code:                        sum(zcolrateia(1:nice))  +  &       ! and ice particles
                         sum(zcolratesa(1:nsnw))             ! and snow particles
  
              IF(pcocsv(ii,jj) > 1.e-10 .AND. zcs_ocsv > 1.e-30 .AND. ioc > 0) THEN
  
  
                 zcvap_new3 = pcocsv(ii,jj)/(1.+ptstep*zcs_ocsv)   ! new gas phase concentration [#/m3]
  