anchor.diff:                 ! [m3(SO4)/m3(air)] by condensation
   
                  !-- Change of volume concentration of sulphate in aerosol [fxm]
  -               aero(ii,jj,in1a:fn2b)%volc(1) = aero(ii,jj,in1a:fn2b)%volc(1) + zdvolsa
  +               aero(ii,jj,in1a:fn2b)%volc(iso4) = aero(ii,jj,in1a:fn2b)%volc(iso4) + zdvolsa
   
                  !-- Clouds
  -               cloud(ii,jj,1:ncld)%volc(1) = cloud(ii,jj,1:ncld)%volc(1)  +  &
  +               cloud(ii,jj,1:ncld)%volc(iso4) = cloud(ii,jj,1:ncld)%volc(iso4)  +  &
                                                 zcolrateca(1:ncld)/zcs_su*mvsu*zdvap1
   
                  ! Rain drops
  -               precp(ii,jj,1:nprc)%volc(1) = precp(ii,jj,1:nprc)%volc(1)  +  &
  +               precp(ii,jj,1:nprc)%volc(iso4) = precp(ii,jj,1:nprc)%volc(iso4)  +  &
                                                 zcolratepa(1:nprc)/zcs_su*mvsu*zdvap1
   
                  !-- Ice clouds
  -               ice(ii,jj,1:nice)%volc(1) = ice(ii,jj,1:nice)%volc(1)  +  &
  +               ice(ii,jj,1:nice)%volc(iso4) = ice(ii,jj,1:nice)%volc(iso4)  +  &
                                               zcolrateia(1:nice)/zcs_su*mvsu*zdvap1
   
                  ! Snow particles
  -               snow(ii,jj,1:nsnw)%volc(1) = snow(ii,jj,1:nsnw)%volc(1)  +  &
  +               snow(ii,jj,1:nsnw)%volc(iso4) = snow(ii,jj,1:nsnw)%volc(iso4)  +  &
                                                zcolratesa(1:nsnw)/zcs_su*mvsu*zdvap1
   
                  !-- Change of number concentration in the smallest bin caused by nucleation
  
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: 
           END DO ! kbdim
        END DO ! klev
  
        DEALLOCATE(zplusterm)
  
     END SUBROUTINE coagulation
  
  