! Subroutine that estimates the Deposition velocities and the depostion
! for the different tracers
SUBROUTINE deposition(vdep_oce, vdep_sic, vdep_ter, vdep_lic, pctsrf, &
        zrho, zdz, pdtphys, RHcl, masse, t_seri, pplay, &
        paprs, lminmax, qmin, qmax, &
        his_ds, source_tr, tr_seri)

  USE dimphy
  USE lmdz_infotrac
  USE indice_sol_mod
  USE lmdz_yomcst

USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
  USE lmdz_paramet
USE lmdz_chem, ONLY: idms, iso2, iso4, ih2s, idmso, imsa, ih2o2, &
          n_avogadro, masse_s, masse_so4, rho_water, rho_ice
  IMPLICIT NONE





  !----------------------------- INPUT -----------------------------------
  LOGICAL :: lminmax
  REAL :: qmin, qmax
  REAL :: vdep_oce(nbtr), vdep_sic(nbtr)
  REAL :: vdep_ter(nbtr), vdep_lic(nbtr)
  REAL :: pctsrf(klon, nbsrf)
  REAL :: zrho(klon, klev)        !Density of air at mid points of Z (kg/m3)
  REAL :: zdz(klon, klev)
  REAL :: pdtphys  ! pas d'integration pour la physique (seconde)
  REAL :: RHcl(klon, klev)  ! humidite relativen ciel clair
  REAL :: t_seri(klon, klev)  ! temperature
  REAL :: pplay(klon, klev)  ! pression pour le mileu de chaque couche (en Pa)
  REAL :: paprs(klon, klev + 1)    !pressure at interface of layers Z (Pa)
  REAL :: masse(nbtr)

  !----------------------------- OUTPUT ----------------------------------
  REAL :: his_ds(klon, nbtr)
  REAL :: source_tr(klon, nbtr)
  REAL :: tr_seri(klon, klev, nbtr) !conc of tracers
  !--------------------- INTERNAL VARIABLES ------------------------------
  INTEGER :: i, it
  REAL :: vdep        !sed. velocity

  DO it = 1, nbtr
    DO i = 1, klon
      vdep = vdep_oce(it) * pctsrf(i, is_oce) + &
              vdep_sic(it) * pctsrf(i, is_sic) + &
              vdep_ter(it) * pctsrf(i, is_ter) + &
              vdep_lic(it) * pctsrf(i, is_lic)
      !--Unit: molec/m2/s for it=1 to nbtr-3, mg/m2/s for it=nbtr-2 to nbtr
      source_tr(i, it) = source_tr(i, it) &
              - vdep * tr_seri(i, 1, it) * zrho(i, 1) / 1.e2
      his_ds(i, it) = vdep * tr_seri(i, 1, it) * zrho(i, 1) / 1.e2 &
              / RNAVO * masse(it) * 1.e3               ! mg/m2/s
    ENDDO
  ENDDO

END SUBROUTINE deposition