source: LMDZ5/branches/IPSLCM6.0.12/libf/phylmd/Dust/blcloud_scav.F @ 4400

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

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

File size: 3.8 KB
Line 
1c Subroutine that calculates the effect of precipitation in scavenging
2c BELOW the cloud, for large scale as well as convective precipitation
3      SUBROUTINE blcloud_scav(lminmax,qmin,qmax,pdtphys,prfl,psfl,
4     .                        pmflxr,pmflxs,zdz,alpha_r,alpha_s,masse,
5     .                                  his_dhbclsc,his_dhbccon,tr_seri)
6
7      USE dimphy
8      USE indice_sol_mod
9      USE infotrac
10      IMPLICIT NONE
11
12#include "dimensions.h"
13#include "chem.h"
14#include "YOMCST.h"
15#include "paramet.h"
16
17c============================= INPUT ===================================
18      REAL qmin,qmax
19      REAL pdtphys  ! pas d'integration pour la physique (seconde)
20!      REAL prfl(klon,klev),   psfl(klon,klev)     !--large-scale
21!      REAL pmflxr(klon,klev), pmflxs(klon,klev)   !--convection
22      REAL alpha_r(nbtr)!--coefficient d'impaction pour la pluie
23      REAL alpha_s(nbtr)!--coefficient d'impaction pour la neige     
24      REAL masse(nbtr)
25      LOGICAL lminmax
26      REAL zdz(klon,klev)
27      REAL prfl(klon,klev+1),   psfl(klon,klev+1)     !--large-scale  ! Titane
28      REAL pmflxr(klon,klev+1), pmflxs(klon,klev+1)   !--convection   ! Titane
29c============================= OUTPUT ==================================
30      REAL tr_seri(klon,klev,nbtr) ! traceur
31      REAL aux_var1(klon,klev) ! traceur
32      REAL aux_var2(klon,klev) ! traceur
33      REAL his_dhbclsc(klon,nbtr), his_dhbccon(klon,nbtr)
34c========================= LOCAL VARIABLES =============================     
35      INTEGER it, k, i, j
36      REAL d_tr(klon,klev,nbtr)
37
38      EXTERNAL minmaxqfi, bcscav_spl
39     
40      DO it=1, nbtr
41c
42      DO j=1,klev
43      DO i=1,klon
44        aux_var1(i,j)=tr_seri(i,j,it)
45        aux_var2(i,j)=d_tr(i,j,it)
46      ENDDO
47      ENDDO
48c
49cnhl      CALL bcscav_spl(pdtphys,prfl,psfl,alpha_r(it),alpha_s(it),
50cnhl     .                tr_seri(1,1,it),d_tr(1,1,it))
51      CALL bcscav_spl(pdtphys,prfl,psfl,alpha_r(it),alpha_s(it),
52     .                aux_var1,aux_var2)
53c
54      DO j=1,klev
55      DO i=1,klon
56        tr_seri(i,j,it)=aux_var1(i,j)
57        d_tr(i,j,it)=aux_var2(i,j)
58      ENDDO
59      ENDDO
60      DO k = 1, klev
61      DO i = 1, klon
62         tr_seri(i,k,it) = tr_seri(i,k,it) + d_tr(i,k,it)
63         his_dhbclsc(i,it)=his_dhbclsc(i,it)-d_tr(i,k,it)/RNAVO*
64     .                masse(it)*1.e3*1.e6*zdz(i,k)/pdtphys  !--mgS/m2/s
65             
66      ENDDO
67      ENDDO
68c
69      DO i=1,klon
70      DO j=1,klev
71        aux_var1(i,j)=tr_seri(i,j,it)
72        aux_var2(i,j)=d_tr(i,j,it)
73      ENDDO
74      ENDDO
75c
76      IF (lminmax) THEN
77        CALL minmaxqfi(aux_var1,qmin,qmax,'depot humide bc lsc')
78cnhl      CALL minmaxqfi(tr_seri(1,1,it),qmin,qmax,'depot humide bc lsc')
79      ENDIF
80c
81c-scheme for convective scavenging
82c
83cnhl      CALL bcscav_spl(pdtphys,pmflxr,pmflxs,alpha_r(it),alpha_s(it),
84cnhl     .                tr_seri(1,1,it),d_tr(1,1,it))
85
86
87      CALL bcscav_spl(pdtphys,pmflxr,pmflxs,alpha_r(it),alpha_s(it),
88     .                aux_var1,aux_var2)
89
90
91c
92      DO i=1,klon
93      DO j=1,klev
94        tr_seri(i,j,it)=aux_var1(i,j)
95        d_tr(i,j,it)=aux_var2(i,j)
96      ENDDO
97      ENDDO
98c
99      DO k = 1, klev
100      DO i = 1, klon
101         tr_seri(i,k,it) = tr_seri(i,k,it) + d_tr(i,k,it)
102         his_dhbccon(i,it)=his_dhbccon(i,it)-d_tr(i,k,it)/RNAVO*
103     .                masse(it)*1.e3*1.e6*zdz(i,k)/pdtphys    !--mgS/m2/s
104      ENDDO
105      ENDDO
106c
107      IF (lminmax) THEN
108        DO j=1,klev
109        DO i=1,klon
110          aux_var1(i,j)=tr_seri(i,j,it)
111        ENDDO
112        ENDDO
113        CALL minmaxqfi(aux_var1,qmin,qmax,'depot humide bc con')
114cnhl      CALL minmaxqfi(tr_seri(1,1,it),qmin,qmax,'depot humide bc con')
115        DO j=1,klev
116        DO i=1,klon
117          tr_seri(i,j,it)=aux_var1(i,j)
118        ENDDO
119        ENDDO
120      ENDIF
121c
122c
123      ENDDO !--boucle sur it
124c
125      END
Note: See TracBrowser for help on using the repository browser.