anchor.diff:      REAL, INTENT(in) :: numc(n1,n2,n3,nn)    ! Particle number concentration
       REAL, INTENT(in) :: mass(n1,n2,n3,nn*n4) ! Particle mass mixing ratio
       REAL, INTENT(in) :: dzt(n1)              ! Inverse of grid level thickness
  +    REAL, INTENT(in) :: dt                   ! timestep
       REAL, INTENT(IN) :: clim                ! Concentration limit
       INTEGER, INTENT(IN) :: flag         ! An option for identifying aerosol, cloud, precipitation, ice and snow
       REAL, INTENT(OUT) :: flxdivm(n1,n2,n3,nn*n4), flxdivn(n1,n2,n3,nn) ! Mass and number divergency
  
positive.code: 
  
  
    SUBROUTINE DepositionSlow(n1,n2,n3,n4,nn,tk,adn,pdn,ustar,numc,mass,dzt,clim,flxdivn,flxdivm,depflxn,depflxm,flag)
      IMPLICIT NONE
  
      INTEGER, INTENT(in) :: n1,n2,n3,n4       ! Grid numbers, number of chemical species
  
negative.code:       !-------------------------------------------------------
        IF (sed_aero) THEN
  
         CALL DepositionSlow(n1,n2,n3,n4,nbins,tk,a_dn,1500.,ustar,naerop,maerop, &
                             dzt,tstep,nlim,andiv,amdiv,andep,remaer,1            )
  
         naerot = naerot - andiv
         maerot = maerot - amdiv
  