anchor.diff:        REAL :: temppi,pressi
   
         LOGICAL :: any_cloud, any_precp, any_ice, any_snow
  +
  +      INTEGER :: nspec, ndry, iwa  ! Total number of compounds, number of dry compounds (shoudl be total-1) and index of water
  +
  +      nspec = spec%getNSpec(type="wet")
  +      ndry = spec%getNSpec(type="dry")
  +      iwa = spec%getIndex("H2O")
  +
  +      ALLOCATE(zplusterm(nspec))
  +      zplusterm = 0.
   
         !-----------------------------------------------------------------------------
         !-- 1) Coagulation to coarse mode calculated in a simplified way: ------------
  
positive.code:                zdvoloc(in2a:fn2b) = zdvoloc(in2a:fn2b) + &                     ! volume change of particles
                                      zcolrate(in2a:fn2b)/zcs_ocsv*mvoc*zdvap3   !  [m3(OC)/m3(air)]
  
                 aero(ii,jj,in1a:fn2b)%volc(2) = &                   !-- change of volume due
                    aero(ii,jj,in1a:fn2b)%volc(2) + zdvoloc        !   due to condensation in 1a-2b
  
                 ! Condensation on hydromets
                 cloud(ii,jj,1:ncld)%volc(2) = cloud(ii,jj,1:ncld)%volc(2)  +  &
                                                zcolrateca(1:ncld)/zcs_ocsv*mvoc*zdvap3
  
                 ! Condensation on rain drops
                 precp(ii,jj,1:nprc)%volc(2) = precp(ii,jj,1:nprc)%volc(2)  +  &
                                                zcolratepa(1:nprc)/zcs_ocsv*mvoc*zdvap3
  
                 ! Condensation on ice particles
                 ice(ii,jj,1:nice)%volc(2) = ice(ii,jj,1:nice)%volc(2)  +  &
                                              zcolrateia(1:nice)/zcs_ocsv*mvoc*zdvap3
  
                 ! Condensation on snow particles
                 snow(ii,jj,1:nsnw)%volc(2) = snow(ii,jj,1:nsnw)%volc(2)  +  &
                                               zcolratesa(1:nprc)/zcs_ocsv*mvoc*zdvap3
  
              END IF
  
  
              ! ---- 5.2) Sulphate -------------------------------------------
              IF(pcsa(ii,jj) > 1.e-10 .AND. zcs_tot > 1.e-30 .AND. spec%isUsed('SO4')) THEN
  
                 !-- Ratio of mass transfer between nucleation and condensation
  
  
negative.code:             !--- 5.1) Organic vapours ------------------------
  
              !---- 5.1.1) Non-volatile organic compound: condenses onto all bins
              IF(pcocnv(ii,jj) > 1.e-10 .AND. zcs_tot > 1.e-30 .AND. ioc > 0) THEN
  
                 zn_vs_c = 0.
  
  