! This subroutine calculates the emissions of BLACK CARBON and ORGANIC ! MATTER SUBROUTINE finemission(zdz,pdtphys,zalt,kminbc,kmaxbc, & scale_param_bb,scale_param_ff, & iregion_ind,iregion_bb, & nbreg_ind,nbreg_bb, & lmt_bcff,lmt_bcnff,lmt_bcbb_l,lmt_bcbb_h, & lmt_bcba,lmt_omff,lmt_omnff,lmt_ombb_l, & lmt_ombb_h,lmt_omnat,lmt_omba,id_fine, & flux_sparam_bb,flux_sparam_ff, & source_tr,flux_tr,tr_seri) USE dimphy USE indice_sol_mod USE infotrac ! USE phytracr_spl_mod, ONLY : nbreg_dust, nbreg_ind, nbreg_bb USE dimensions_mod, ONLY: iim, jjm, llm, ndm USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, & ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm USE yomcst_mod_h, ONLY: RPI, RCLUM, RHPLA, RKBOL, RNAVO & , RDAY, REA, REPSM, RSIYEA, RSIDAY, ROMEGA & , R_ecc, R_peri, R_incl & , RA, RG, R1SA & , RSIGMA & , R, RMD, RMV, RD, RV, RCPD & , RMO3, RMCO2, RMC, RMCH4, RMN2O, RMCFC11, RMCFC12 & , RCPV, RCVD, RCVV, RKAPPA, RETV, eps_w & , RCW, RCS & , RLVTT, RLSTT, RLMLT, RTT, RATM & , RESTT, RALPW, RBETW, RGAMW, RALPS, RBETS, RGAMS & , RALPD, RBETD, RGAMD IMPLICIT NONE INCLUDE "chem.h" INTEGER :: i, k, kminbc, kmaxbc !============================= INPUT =================================== REAL :: pdtphys ! pas d'integration pour la physique (seconde) REAL :: zalt(klon,klev) REAL :: zdz(klon,klev) ! !------------------------- Scaling Parameters -------------------------- ! INTEGER :: nbreg_ind,nbreg_bb INTEGER :: iregion_ind(klon) !Defines regions for SO2, BC & OM INTEGER :: iregion_bb(klon) !Defines regions for SO2, BC & OM REAL :: scale_param_bb(nbreg_bb) !Scaling parameter for biomas burning REAL :: scale_param_ff(nbreg_ind) !Scaling parameter for industrial emissions (fossil fuel) INTEGER :: id_fine !============================= OUTPUT ================================== REAL :: source_tr(klon,nbtr) REAL :: flux_tr(klon,nbtr) REAL :: tr_seri(klon,klev,nbtr) ! traceur REAL :: flux_sparam_bb(klon), flux_sparam_ff(klon) !========================= LOCAL VARIABLES ============================= REAL :: zzdz !------------------------- BLACK CARBON emissions ---------------------- REAL :: lmt_bcff(klon) ! emissions de BC fossil fuels REAL :: lmt_bcnff(klon) ! emissions de BC non-fossil fuels REAL :: lmt_bcbb_l(klon) ! emissions de BC biomass basses REAL :: lmt_bcbb_h(klon) ! emissions de BC biomass hautes REAL :: lmt_bcba(klon) ! emissions de BC bateau !------------------------ ORGANIC MATTER emissions --------------------- REAL :: lmt_omff(klon) ! emissions de OM fossil fuels REAL :: lmt_omnff(klon) ! emissions de OM non-fossil fuels REAL :: lmt_ombb_l(klon) ! emissions de OM biomass basses REAL :: lmt_ombb_h(klon) ! emissions de OM biomass hautes REAL :: lmt_omnat(klon) ! emissions de OM Natural REAL :: lmt_omba(klon) ! emissions de OM bateau EXTERNAL condsurfc !======================================================================== ! LOW LEVEL EMISSIONS !======================================================================== ! corresponds to bc_source.EQ.3 DO i=1,klon IF (iregion_ind(i).GT.0) THEN IF(id_fine>0) source_tr(i,id_fine)=source_tr(i,id_fine)+ & (scale_param_ff(iregion_ind(i))*lmt_bcff(i)+ & !g/m2/s scale_param_ff(iregion_ind(i))*lmt_omff(i) & ) * 1.e4 !g/m2/s ! IF(id_fine>0) flux_tr(i,id_fine)=flux_tr(i,id_fine)+ & (scale_param_ff(iregion_ind(i))*lmt_bcff(i)+ & !mg/m2/s scale_param_ff(iregion_ind(i))*lmt_omff(i) & ) * 1.e4 *1.e3 !mg/m2/s ! flux_sparam_ff(i)= flux_sparam_ff(i) + & scale_param_ff(iregion_ind(i))* & ( lmt_bcff(i)+lmt_omff(i)) & *1.e4*1.e3 ENDIF IF (iregion_bb(i).GT.0) THEN IF(id_fine>0) source_tr(i,id_fine)=source_tr(i,id_fine)+ & (scale_param_bb(iregion_bb(i))*lmt_bcbb_l(i)+ & !g/m2/s scale_param_bb(iregion_bb(i))*lmt_ombb_l(i) & !g/m2/s ) * 1.e4 !g/m2/s ! IF(id_fine>0) flux_tr(i,id_fine)=flux_tr(i,id_fine)+ & (scale_param_bb(iregion_bb(i))*lmt_bcbb_l(i)+ & !mg/m2/s scale_param_bb(iregion_bb(i))*lmt_ombb_l(i)+ & !mg/m2/s scale_param_bb(iregion_bb(i))*lmt_bcbb_h(i)+ & !mg/m2/s scale_param_bb(iregion_bb(i))*lmt_ombb_h(i) & !mg/m2/s ) * 1.e4 *1.e3 !mg/m2/s ! flux_sparam_bb(i)=flux_sparam_bb(i) + & scale_param_bb(iregion_bb(i))*(lmt_bcbb_l(i) + & lmt_bcbb_h(i) + lmt_ombb_l(i) + lmt_ombb_h(i)) & *1.e4*1.e3 ENDIF IF(id_fine>0) source_tr(i,id_fine)=source_tr(i,id_fine)+ & (lmt_bcnff(i)+lmt_bcba(i)+lmt_omnff(i)+ & lmt_omnat(i)+lmt_omba(i)) * 1.e4 !g/m2/s ! IF(id_fine>0) flux_tr(i,id_fine)=flux_tr(i,id_fine)+ & (lmt_bcnff(i)+lmt_omnff(i)+lmt_omnat(i)+ & lmt_omba(i)+lmt_bcba(i)) * 1.e4 *1.e3 !mg/m2/s ! flux_sparam_ff(i)= flux_sparam_ff(i) + & (lmt_omba(i)+lmt_bcba(i))*1.e4*1.e3 ENDDO !======================================================================== ! HIGH LEVEL EMISSIONS !======================================================================== ! Sources hautes de BC/OM ! ! HIGH LEVEL EMISSIONS OF SO2 ARE IN PRECUREMISSION.F ! k=2 !introducing emissions in level 2 !nhl DO i = 1, klon ! !nhl tr_seri(i,k,id_fine)=tr_seri(i,k,id_fine)+scale_param_ff(iregion_ind(i))* !nhl . (lmt_bcff_h(i)+lmt_omff_h(i))/zdz(i,k)/100.*pdtphys ! !nhl ENDDO DO k=kminbc, kmaxbc DO i = 1, klon zzdz=zalt(i,kmaxbc+1)-zalt(i,kminbc) ! IF (iregion_bb(i) .GT.0) THEN IF(id_fine>0) tr_seri(i,k,id_fine)=tr_seri(i,k,id_fine)+ & (scale_param_bb(iregion_bb(i))*lmt_bcbb_h(i)+ & scale_param_bb(iregion_bb(i))*lmt_ombb_h(i)) & /zzdz/100.*pdtphys ENDIF ! ENDDO ENDDO ! END SUBROUTINE finemission