source: LMDZ5/trunk/libf/phylmd/Dust/deposition.F @ 4021

Last change on this file since 4021 was 2630, checked in by fhourdin, 8 years ago

Importation du modèle d'aérosols de Boucher, Escribano et al.

File size: 2.1 KB
Line 
1c Subroutine that estimates the Deposition velocities and the depostion
2C for the different tracers
3      subroutine deposition(vdep_oce,vdep_sic,vdep_ter,vdep_lic,pctsrf,
4     .                      zrho,zdz,pdtphys,RHcl,masse,t_seri,pplay,
5     .                      paprs,lminmax,qmin,qmax,
6     .                        his_ds,source_tr,tr_seri)
7
8      USE dimphy
9      USE infotrac
10      USE indice_sol_mod
11
12      IMPLICIT NONE
13
14#include "dimensions.h"
15#include "chem.h"
16#include "YOMCST.h"
17#include "paramet.h"
18
19c----------------------------- INPUT -----------------------------------
20      LOGICAL lminmax
21      REAL qmin, qmax
22      REAL vdep_oce(nbtr), vdep_sic(nbtr)
23      REAL vdep_ter(nbtr), vdep_lic(nbtr)     
24      REAL pctsrf(klon,nbsrf)
25      REAL zrho(klon,klev)        !Density of air at mid points of Z (kg/m3)
26      REAL zdz(klon,klev)       
27      REAL pdtphys  ! pas d'integration pour la physique (seconde)
28      REAL RHcl(klon,klev)  ! humidite relativen ciel clair
29      REAL t_seri(klon,klev)  ! temperature
30      REAL pplay(klon,klev)  ! pression pour le mileu de chaque couche (en Pa)
31      REAL paprs(klon, klev+1)    !pressure at interface of layers Z (Pa)
32      REAL masse(nbtr)
33                                         
34c----------------------------- OUTPUT ----------------------------------
35      REAL his_ds(klon,nbtr)                                         
36      REAL source_tr(klon,nbtr)
37      REAL tr_seri(klon, klev,nbtr) !conc of tracers
38c--------------------- INTERNAL VARIABLES ------------------------------     
39      INTEGER i, it
40      REAL vdep        !sed. velocity
41
42      DO it=1, nbtr
43      DO i=1, klon
44          vdep=vdep_oce(it)*pctsrf(i,is_oce)+
45     .         vdep_sic(it)*pctsrf(i,is_sic)+
46     .         vdep_ter(it)*pctsrf(i,is_ter)+
47     .         vdep_lic(it)*pctsrf(i,is_lic)
48c--Unit: molec/m2/s for it=1 to nbtr-3, mg/m2/s for it=nbtr-2 to nbtr
49          source_tr(i,it)=source_tr(i,it)
50     .                    -vdep*tr_seri(i,1,it)*zrho(i,1)/1.e2
51          his_ds(i,it)=vdep*tr_seri(i,1,it)*zrho(i,1)/1.e2
52     .                 /RNAVO*masse(it)*1.e3               ! mg/m2/s
53      ENDDO
54      ENDDO
55c
56      END
Note: See TracBrowser for help on using the repository browser.