-
Notifications
You must be signed in to change notification settings - Fork 0
/
BiomassBurningMapping_FINN.inc
61 lines (58 loc) · 3.23 KB
/
BiomassBurningMapping_FINN.inc
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
!-------------------------------------------------------------------------!
! Mapping of emissions from FINNv1 to EMEP EmChem09 species
! Included by ForestFire_ml
!-------------------------------------------------------------------------!
character(len=*), parameter :: &
BiomassBurningMapping = "FINNv1.5toEmChem09soa"
!FINNv1.5 uses 7 more GEOS-CHEM species and 2 more EMEP:
integer, parameter :: &
!A2018 NBB_DEFS = 27 & ! No mapping lines below
NBB_DEFS = 27-1 & ! No mapping lines below
,NEMEPSPECS = 17 ! No EMEP chemical mech specs used
!A2018 ,NEMEPSPECS = 18 ! No EMEP chemical mech specs used
!-----------------------------------------------------------------------!
! Column-1 gives FINNv1 species ,
! Column-2 gives MW if needed to get to kg/day. Not needed for FINNv1
! (Some FINN emissions are given with mass basis, kg/day, so just
! set to 1.0. Others are given as mole/day, so multiply by MW/1000
! to get kg/day)
! Column-3 gives mass fraction of FINNv1 emissions assign to EMEP species
! Column-4 gives EMEP species
!-----------------------------------------------------------------------!
! BBname unitsfac frac emep:
type(bbtype), dimension(NBB_DEFS) :: FF_defs = (/ &
bbtype("CO ",0.028 , 1.000, CO ) & ! Tracer, will add to CO too
!A2018 tmp ,bbtype("CO ",0.028 , 1.000, FFIRE_CO ) & ! Tracer, will add to CO too
,bbtype("NO ",0.030 , 1.000, NO ) &
,bbtype("NO2 ",0.046 , 1.000, NO2 ) &
,bbtype("SO2 ",0.064 , 1.000, SO2 ) &
,bbtype("NH3 ",0.017 , 1.000, NH3 ) &
,bbtype("ACET",0.058 , 1.000, C2H6 ) & !acetone
,bbtype("ALD2",0.044 , 1.000, CH3CHO ) &
,bbtype("ALK4",0.058 , 1.000, NC4H10 ) &
,bbtype("C2H6",0.030 , 1.000, C2H6 ) &
,bbtype("C3H8",0.044 , 0.700, NC4H10 ) & ! obs
,bbtype("CH2O",0.030 , 1.000, HCHO ) &
,bbtype("MEK ",0.072 , 1.000, MEK ) &
,bbtype("PRPE",0.042 , 1.000, C3H6 ) &
! We read in OC and PM25, but want OM and REMPPM25
,bbtype("PM25",1.0 , 1.000, FFIRE_REMPPM25 ) & ! Will need to subtract OM, BC
,bbtype("OC ",1.7 , 1.000, FFIRE_OM ) & ! Put OM/OC=1.7 in fac
,bbtype("BC ",1.0 , 1.000, FFIRE_BC ) &
! Subtract, assuming OM/OC=1.7. ForestFire_ml will pevent zeros
,bbtype("OC ",-1.7 , 1.000, FFIRE_REMPPM25 ) & ! Will subtract OM
,bbtype("BC ",-1.0 , 1.000, FFIRE_REMPPM25 ) & ! Will subtract BC
! FINN v1.5, GEOS-CHEM 2015 changes: excludes 1 species:
! and needs one less EMEP species, C5H8
! ,bbtype("ISOP",0.068 , 1.000, C5H8 ) &
! FINN v1.5, GEOS-CHEM 2015 changes: added 8 species:
! and needs 3 more EMEP species = C2H4, OXYL, MGLYOX
,bbtype("C2H4",0.028 , 1.000, C2H4 ) & ! v1.5, new EMEP
,bbtype("C3H8",0.044 , 1.000, NC4H10 ) & ! v1.5
,bbtype("GLYC",0.060 , 1.000, CH3CHO ) & ! v1.5 hydroxy-aceteldehyed
,bbtype("HAC", 0.060 , 1.000, C2H6 ) & ! v1.5 hydroxy-acetone. Query surrogate?
,bbtype("BENZ", 0.078 , 1.000, OXYL ) & ! v1.5, new EMEP
,bbtype("TOLU", 0.092 , 1.000, OXYL ) & ! v1.5
,bbtype("XYLE", 0.106 , 1.000, OXYL ) & ! v1.5
,bbtype("MGLY", 0.072 , 1.000, MGLYOX ) & ! v1.5, new EMEP
/)