source: LMDZ6/trunk/libf/phylmd/Dust/incloud_scav.f90 @ 5354

Last change on this file since 5354 was 5337, checked in by Laurent Fairhead, 3 weeks ago

Getting rid of dependance to dynamics

File size: 3.2 KB
Line 
1! Subroutine that calculates the effect of precipitation in scavenging
2! WITHIN the cloud, for large scale as well as convective precipitation
3SUBROUTINE incloud_scav(lminmax,qmin,qmax,masse,henry,kk,prfl, &
4        psfl,pmflxr,pmflxs,zrho,zdz,t_seri,pdtphys, &
5        his_dhlsc,his_dhcon1,tr_seri)
6
7USE chem_mod_h
8    USE dimphy
9  USE infotrac_phy, ONLY: nbtr
10  USE indice_sol_mod
11
12!!USE paramet_mod_h
13USE yomcst_mod_h
14IMPLICIT NONE
15
16
17
18
19
20  !============================= INPUT ===================================
21  REAL :: qmin, qmax
22  REAL :: masse(nbtr)
23  REAL :: henry(nbtr)         !--cste de Henry  mol/l/atm
24  REAL :: kk(nbtr)            !--coefficient de var avec T (K)
25  REAL :: prfl(klon,klev+1),   psfl(klon,klev+1)     !--large-scale
26   ! REAL prfl(klon,klev),   psfl(klon,klev)     !--large-scale
27  REAL :: pmflxr(klon,klev+1), pmflxs(klon,klev+1)   !--convection
28   ! REAL pmflxr(klon,klev), pmflxs(klon,klev)   !--convection
29  REAL :: zrho(klon,klev), zdz(klon,klev)
30  REAL :: t_seri(klon,klev)
31  LOGICAL :: lminmax
32  REAL :: pdtphys
33   ! REAL prfl(klon,klev+1),   psfl(klon,klev+1)     !--large-scale  ! Titane
34   ! REAL pmflxr(klon,klev+1), pmflxs(klon,klev+1)   !--convection   ! Titane
35  !============================= OUTPUT ==================================
36  REAL :: tr_seri(klon,klev,nbtr) ! traceur
37  REAL :: aux_var1(klon,klev) ! traceur
38  REAL :: aux_var2(klon) ! traceur
39  REAL :: aux_var3(klon) ! traceur
40  REAL :: his_dhlsc(klon,nbtr)        ! in-cloud scavenging lsc
41  REAL :: his_dhcon1(klon,nbtr)       ! in-cloud scavenging con
42  !========================= LOCAL VARIABLES =============================
43  INTEGER :: it, i, j
44
45  EXTERNAL minmaxqfi, inscav_spl
46
47  DO it=1, nbtr
48  !
49  DO i=1,klon
50    aux_var2(i)=his_dhlsc(i,it)
51    aux_var3(i)=his_dhcon1(i,it)
52  ENDDO
53  DO j=1,klev
54  DO i=1,klon
55    aux_var1(i,j)=tr_seri(i,j,it)
56  ENDDO
57  ENDDO
58  !
59  IF (lminmax) THEN
60    CALL minmaxqfi(aux_var1,qmin,qmax,'avt inscav')
61  !nhl      CALL minmaxqfi(tr_seri(1,1,it),qmin,qmax,'avt inscav')
62  ENDIF
63  !
64  !nhl      CALL inscav_spl(pdtphys,it,masse(it),henry(it),kk(it),0.5e-3,
65  !nhl     .             prfl,psfl,zrho,zdz,t_seri,tr_seri(1,1,it),
66  !nhl     .             his_dhlsc(1,it))
67  CALL inscav_spl(pdtphys,it,masse(it),henry(it),kk(it),0.5e-3, &
68        prfl,psfl,zrho,zdz,t_seri,aux_var1,aux_var2)
69  !
70  IF (lminmax) THEN
71    CALL minmaxqfi(aux_var1,qmin,qmax,'depot humide lsc')
72  !nhl      CALL minmaxqfi(tr_seri(1,1,it),qmin,qmax,'depot humide lsc')
73  ENDIF
74  !
75  !
76  !-scheme for convective in-cloud scavenging
77  !
78  !nhl      CALL inscav_spl(pdtphys,it,masse(it),henry(it),kk(it),1.e-3,
79  !nhl     .             pmflxr,pmflxs,zrho,zdz,t_seri,tr_seri(1,1,it),
80  !nhl     .             his_dhcon1(1,it))
81  CALL inscav_spl(pdtphys,it,masse(it),henry(it),kk(it),1.e-3, &
82        pmflxr,pmflxs,zrho,zdz,t_seri,aux_var1,aux_var3)
83  !
84  IF (lminmax) THEN
85    CALL minmaxqfi(aux_var1,qmin,qmax,'depot humide con')
86  !nhl      CALL minmaxqfi(tr_seri(1,1,it),qmin,qmax,'depot humide con')
87  ENDIF
88  !
89  DO j=1,klev
90  DO i=1,klon
91    tr_seri(i,j,it)=aux_var1(i,j)
92  ENDDO
93  ENDDO
94  DO i=1,klon
95    his_dhlsc(i,it)=aux_var2(i)
96    his_dhcon1(i,it)=aux_var3(i)
97  ENDDO
98
99  !
100  ENDDO !--boucle sur it
101
102END SUBROUTINE incloud_scav
Note: See TracBrowser for help on using the repository browser.