source: LMDZ6/trunk/libf/phylmd/Dust/finemission.f90 @ 5277

Last change on this file since 5277 was 5274, checked in by abarral, 9 hours ago

Replace yomcst.h by existing module
File size: 6.5 KB
Line 
1! This subroutine calculates the emissions of BLACK CARBON and ORGANIC
2! MATTER
3SUBROUTINE finemission(zdz,pdtphys,zalt,kminbc,kmaxbc, &
4        scale_param_bb,scale_param_ff, &
5        iregion_ind,iregion_bb, &
6        nbreg_ind,nbreg_bb, &
7        lmt_bcff,lmt_bcnff,lmt_bcbb_l,lmt_bcbb_h, &
8        lmt_bcba,lmt_omff,lmt_omnff,lmt_ombb_l, &
9        lmt_ombb_h,lmt_omnat,lmt_omba,id_fine, &
10        flux_sparam_bb,flux_sparam_ff, &
11        source_tr,flux_tr,tr_seri)
12
13  USE dimphy
14  USE indice_sol_mod
15  USE infotrac
16   ! USE phytracr_spl_mod, ONLY : nbreg_dust, nbreg_ind, nbreg_bb
17  USE dimensions_mod, ONLY: iim, jjm, llm, ndm
18USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
19          ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
20USE yomcst_mod_h, ONLY: RPI, RCLUM, RHPLA, RKBOL, RNAVO                   &
21          , RDAY, REA, REPSM, RSIYEA, RSIDAY, ROMEGA                  &
22          , R_ecc, R_peri, R_incl                                      &
23          , RA, RG, R1SA                                         &
24          , RSIGMA                                                     &
25          , R, RMD, RMV, RD, RV, RCPD                    &
26          , RMO3, RMCO2, RMC, RMCH4, RMN2O, RMCFC11, RMCFC12        &
27          , RCPV, RCVD, RCVV, RKAPPA, RETV, eps_w                    &
28          , RCW, RCS                                                 &
29          , RLVTT, RLSTT, RLMLT, RTT, RATM                           &
30          , RESTT, RALPW, RBETW, RGAMW, RALPS, RBETS, RGAMS            &
31          , RALPD, RBETD, RGAMD
32IMPLICIT NONE
33
34
35  INCLUDE "chem.h"
36
37
38
39  INTEGER :: i, k, kminbc, kmaxbc
40  !============================= INPUT ===================================
41  REAL :: pdtphys  ! pas d'integration pour la physique (seconde)
42  REAL :: zalt(klon,klev)
43  REAL :: zdz(klon,klev)
44  !
45  !------------------------- Scaling Parameters --------------------------
46  !
47  INTEGER :: nbreg_ind,nbreg_bb
48  INTEGER :: iregion_ind(klon)  !Defines regions for SO2, BC & OM
49  INTEGER :: iregion_bb(klon)   !Defines regions for SO2, BC & OM
50  REAL :: scale_param_bb(nbreg_bb) !Scaling parameter for biomas burning
51  REAL :: scale_param_ff(nbreg_ind) !Scaling parameter for industrial emissions (fossil fuel)
52  INTEGER :: id_fine
53  !============================= OUTPUT ==================================
54  REAL :: source_tr(klon,nbtr)
55  REAL :: flux_tr(klon,nbtr)
56  REAL :: tr_seri(klon,klev,nbtr) ! traceur
57  REAL :: flux_sparam_bb(klon), flux_sparam_ff(klon)
58  !========================= LOCAL VARIABLES =============================
59  REAL :: zzdz
60  !------------------------- BLACK CARBON emissions ----------------------
61  REAL :: lmt_bcff(klon)       ! emissions de BC fossil fuels
62  REAL :: lmt_bcnff(klon)      ! emissions de BC non-fossil fuels
63  REAL :: lmt_bcbb_l(klon)     ! emissions de BC biomass basses
64  REAL :: lmt_bcbb_h(klon)     ! emissions de BC biomass hautes
65  REAL :: lmt_bcba(klon)       ! emissions de BC bateau
66  !------------------------ ORGANIC MATTER emissions ---------------------
67  REAL :: lmt_omff(klon)     ! emissions de OM fossil fuels
68  REAL :: lmt_omnff(klon)    ! emissions de OM non-fossil fuels
69  REAL :: lmt_ombb_l(klon)   ! emissions de OM biomass basses
70  REAL :: lmt_ombb_h(klon)   ! emissions de OM biomass hautes
71  REAL :: lmt_omnat(klon)    ! emissions de OM Natural
72  REAL :: lmt_omba(klon)     ! emissions de OM bateau
73
74  EXTERNAL condsurfc
75  !========================================================================
76                     ! LOW LEVEL EMISSIONS
77  !========================================================================
78
79  ! corresponds to bc_source.EQ.3
80
81  DO i=1,klon
82     IF (iregion_ind(i).GT.0) THEN
83   IF(id_fine>0)    source_tr(i,id_fine)=source_tr(i,id_fine)+ &
84         (scale_param_ff(iregion_ind(i))*lmt_bcff(i)+ & !g/m2/s
85         scale_param_ff(iregion_ind(i))*lmt_omff(i) &
86         )     * 1.e4                        !g/m2/s
87  !
88  IF(id_fine>0)     flux_tr(i,id_fine)=flux_tr(i,id_fine)+ &
89        (scale_param_ff(iregion_ind(i))*lmt_bcff(i)+ & !mg/m2/s
90        scale_param_ff(iregion_ind(i))*lmt_omff(i) &
91        )     * 1.e4 *1.e3                  !mg/m2/s
92  !
93       flux_sparam_ff(i)= flux_sparam_ff(i) + &
94             scale_param_ff(iregion_ind(i))* &
95             ( lmt_bcff(i)+lmt_omff(i)) &
96             *1.e4*1.e3
97     ENDIF
98     IF (iregion_bb(i).GT.0) THEN
99   IF(id_fine>0)    source_tr(i,id_fine)=source_tr(i,id_fine)+ &
100         (scale_param_bb(iregion_bb(i))*lmt_bcbb_l(i)+ & !g/m2/s
101         scale_param_bb(iregion_bb(i))*lmt_ombb_l(i) & !g/m2/s
102         )     * 1.e4                        !g/m2/s
103  !
104   IF(id_fine>0)    flux_tr(i,id_fine)=flux_tr(i,id_fine)+ &
105         (scale_param_bb(iregion_bb(i))*lmt_bcbb_l(i)+ & !mg/m2/s
106         scale_param_bb(iregion_bb(i))*lmt_ombb_l(i)+ & !mg/m2/s
107         scale_param_bb(iregion_bb(i))*lmt_bcbb_h(i)+ & !mg/m2/s
108         scale_param_bb(iregion_bb(i))*lmt_ombb_h(i) & !mg/m2/s
109         )     * 1.e4 *1.e3                  !mg/m2/s
110  !
111       flux_sparam_bb(i)=flux_sparam_bb(i) + &
112             scale_param_bb(iregion_bb(i))*(lmt_bcbb_l(i) + &
113             lmt_bcbb_h(i) + lmt_ombb_l(i) + lmt_ombb_h(i)) &
114             *1.e4*1.e3
115     ENDIF
116   IF(id_fine>0)  source_tr(i,id_fine)=source_tr(i,id_fine)+ &
117         (lmt_bcnff(i)+lmt_bcba(i)+lmt_omnff(i)+ &
118         lmt_omnat(i)+lmt_omba(i))     * 1.e4           !g/m2/s
119  !
120   IF(id_fine>0)  flux_tr(i,id_fine)=flux_tr(i,id_fine)+ &
121         (lmt_bcnff(i)+lmt_omnff(i)+lmt_omnat(i)+ &
122         lmt_omba(i)+lmt_bcba(i))     * 1.e4 *1.e3      !mg/m2/s
123  !
124     flux_sparam_ff(i)= flux_sparam_ff(i) + &
125           (lmt_omba(i)+lmt_bcba(i))*1.e4*1.e3
126  ENDDO
127
128  !========================================================================
129                     ! HIGH LEVEL EMISSIONS
130  !========================================================================
131
132  !  Sources hautes de BC/OM
133
134  !
135  ! HIGH LEVEL EMISSIONS OF SO2 ARE IN PRECUREMISSION.F
136  !
137  k=2                             !introducing emissions in level 2
138  !nhl      DO i = 1, klon
139  !
140  !nhl      tr_seri(i,k,id_fine)=tr_seri(i,k,id_fine)+scale_param_ff(iregion_ind(i))*
141  !nhl     .               (lmt_bcff_h(i)+lmt_omff_h(i))/zdz(i,k)/100.*pdtphys
142  !
143  !nhl      ENDDO
144
145  DO k=kminbc, kmaxbc
146  DO i = 1, klon
147      zzdz=zalt(i,kmaxbc+1)-zalt(i,kminbc)
148  !
149     IF (iregion_bb(i) .GT.0) THEN
150    IF(id_fine>0)   tr_seri(i,k,id_fine)=tr_seri(i,k,id_fine)+ &
151          (scale_param_bb(iregion_bb(i))*lmt_bcbb_h(i)+ &
152          scale_param_bb(iregion_bb(i))*lmt_ombb_h(i)) &
153          /zzdz/100.*pdtphys
154     ENDIF
155  !
156  ENDDO
157  ENDDO
158  !
159END SUBROUTINE finemission
Note: See TracBrowser for help on using the repository browser.