1 | ! This subroutine calculates the emissions of BLACK CARBON and ORGANIC |
---|
2 | ! MATTER |
---|
3 | SUBROUTINE 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 |
---|
18 | USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, & |
---|
19 | ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm |
---|
20 | USE 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 |
---|
32 | IMPLICIT 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 | ! |
---|
159 | END SUBROUTINE finemission |
---|