Merge remote-tracking branch 'origin/release-v4.6.1'
[WRF.git] / chem / module_data_megan2.F
blob6bfc3a0badfac2191ee1dfd3d9825ae830c79940
1 MODULE module_data_megan2
3   !
4   ! This module contains parameters for biogenic emissions model
5   ! MEGAN v2.04
6   !
7   !  Reference: 
8   !
9   !    Estimates of global terrestial isoprene emissions using MEGAN
10   !    (Model of Emissions of Gases and Aerosols from Nature )
11   !    A. Guenther, T. Karl, P. Harley, C. Wiedinmyer, P.I. Palmer, and C. Geron
12   !    Atmospheric Chemistry and Physics, 6, 3181-3210, 2006
13   !
14   !    MEGAN v2.0 Documentation
15   !
16   ! August, 2007
17   !
18   ! Serena H. Chung          Washington State University
19   ! Tan Sakulyanontvittaya   University of Colorado
20   ! Christine Wiedinmyer     National Center for Atmospheric Research
21   !
23   IMPLICIT NONE
25   !=======================================================================
26   ! The following contains 20 clases of species for MEGAN v2.04
27   !  Originally from Tan's SPC_MGN.EXT (12/02/2006)
28   !=======================================================================
30   INTEGER, PARAMETER :: n_mgn_spc  = 20
31   INTEGER, PARAMETER :: n_spc_char = 16
32   CHARACTER (len=n_spc_char), DIMENSION (n_mgn_spc) :: mgn_spc
33   
34   ! Indices for 20 Classes of MEGAN v2.04 Species
35   !    The ordering here follows that of Table 2 of MEGAN v2.0
36   !    Documentation.  For programming, the ordering does not matter, 
37   !    except that imgn_isop should always be 1.
38   INTEGER, PARAMETER ::     &
39        &  imgn_isop  =  1   & ! isoprene
40        & ,imgn_mbo   =  2   & ! MBO (only (2-methyl-3-ene-2-ol)-methylbutenol) 
41        & ,imgn_myrc  =  3   & ! myrcene
42        & ,imgn_sabi  =  4   & ! sabinene
43        & ,imgn_limo  =  5   & ! limonene
44        & ,imgn_3car  =  6   & ! 3 carene
45        & ,imgn_ocim  =  7   & ! trans-beta-ocimene
46        & ,imgn_bpin  =  8   & ! beta-pinene
47        & ,imgn_apin  =  9   & ! alpha-pinene
48        & ,imgn_afarn = 10   & ! alpha-farnescene
49        & ,imgn_bcar  = 11   & ! beta-caryophyllene
50        & ,imgn_meoh  = 12   & ! methanol
51        & ,imgn_acto  = 13   & ! acetone
52        & ,imgn_acta  = 14   & ! acetaldehyde and ethanol
53        & ,imgn_form  = 15   & ! formic acid, formaldehyde, acetic acid
54        & ,imgn_ch4   = 16   & ! methane
55        & ,imgn_no    = 17   & ! nitric oxide (and NO2 and NH3 )
56        & ,imgn_omtp  = 18   & ! other monoterpenes
57        & ,imgn_osqt  = 19   & ! other sesquiterpenes
58        & ,imgn_co    = 20     ! CO and other VOC                                          
60   !
61   ! For reference, below is how Tan's ordereing the 20 classes of species
62   ! in his SPC_MGN.EXT.  They are different from the ordering above.
63   ! The following is provide for convenience when copying/pasting
64   ! Tan's code to here.
65   !
66   ! Tan's Order        Species Class
67   !   1                isop
68   !   2                myrc
69   !   3                sabi
70   !   4                limo
71   !   5                3car
72   !   6                ocim
73   !   7                bpin
74   !   8                apin
75   !   9                omtp
76   !  10                afarn
77   !  11                bcar
78   !  12                osqt
79   !  13                mbo
80   !  14                meoh
81   !  15                acto
82   !  16                ch4
83   !  17                no
84   !  18                acta
85   !  19                form
86   !  20                co        
89   ! Names of the 20 calses of Species
90   ! (These are not used, except for testing stage)
91   DATA  mgn_spc( imgn_isop  )   / 'isoprene        ' /  ! isoprene                                              
92   DATA  mgn_spc( imgn_mbo   )   / 'MBO             ' /  ! MBO (only (2-methyl-3-ene-2-ol)-methylbutenol)               
93   DATA  mgn_spc( imgn_myrc  )   / 'myrcene         ' /  ! myrcene                                               
94   DATA  mgn_spc( imgn_sabi  )   / 'sabinene        ' /  ! sabinene                                              
95   DATA  mgn_spc( imgn_limo  )   / 'limonene        ' /  ! limonene                                              
96   DATA  mgn_spc( imgn_3car  )   / '3carene         ' /  ! 3 carene                                              
97   DATA  mgn_spc( imgn_ocim  )   / 'ocimene         ' /  ! trans-beta-ocimene                                    
98   DATA  mgn_spc( imgn_bpin  )   / 'b-pinene        ' /  ! beta-pinene                                           
99   DATA  mgn_spc( imgn_apin  )   / 'a-pininen       ' /  ! alpha-pinene                                          
100   DATA  mgn_spc( imgn_afarn )   / 'a-farnescene    ' /  ! alpha-farnescene                                      
101   DATA  mgn_spc( imgn_bcar  )   / 'b-caryophyllene ' /  ! beta-caryophyllene                                    
102   DATA  mgn_spc( imgn_meoh  )   / 'methanol        ' /  ! methanol                                              
103   DATA  mgn_spc( imgn_acto  )   / 'actone          ' /  ! acetone                                               
104   DATA  mgn_spc( imgn_acta  )   / 'acta            ' /  ! acetaldehyde and ethanol
105   DATA  mgn_spc( imgn_form  )   / 'form            ' /  ! formic acid, formaldehyde, acetic acid  
106   DATA  mgn_spc( imgn_ch4   )   / 'methane         ' /  ! methane                                               
107   DATA  mgn_spc( imgn_no    )   / 'NO              ' /  ! nitric oxide  (and NO2 and NH3 )
108   DATA  mgn_spc( imgn_omtp  )   / 'other_monoterp  ' /  ! other monoterpenes                                    
109   DATA  mgn_spc( imgn_osqt  )   / 'other_sesquiterp' /  ! other sesquiterpenes                                  
110   DATA  mgn_spc( imgn_co    )   / 'CO              ' /  ! co and other voc                                      
112   !=====================================================================
113   !  The following contains "temperature dependent" parameter for
114   !  non-isoprene species.  This is the parameter "beta" provided in
115   !  Table 2 of MEGAN v2.0 Documentaion.
116   !  Originally from Tan's TEMPD_PRM.EXT (11/30/2006)
117   !  Updated on 08/13/2007 based on Alex Guenther's email.
118   !=====================================================================
120   REAL , DIMENSION(n_mgn_spc) :: tdf_prm
122   DATA  tdf_prm ( imgn_isop  )  / 0.09  / ! isoprene  (but not used for isoprene)
123   DATA  tdf_prm ( imgn_mbo   )  / 0.09  / ! MBO (only (2-methyl-3-ene-2-ol)-methylbutenol)
124   DATA  tdf_prm ( imgn_myrc  )  / 0.10  / ! myrcene                                
125   DATA  tdf_prm ( imgn_sabi  )  / 0.10  / ! sabinene                               
126   DATA  tdf_prm ( imgn_limo  )  / 0.10  / ! limonene                               
127   DATA  tdf_prm ( imgn_3car  )  / 0.10  / ! 3 carene                               
128   DATA  tdf_prm ( imgn_ocim  )  / 0.10  / ! trans-beta-ocimene                     
129   DATA  tdf_prm ( imgn_bpin  )  / 0.10  / ! beta-pinene                            
130   DATA  tdf_prm ( imgn_apin  )  / 0.10  / ! alpha-pinene                           
131   DATA  tdf_prm ( imgn_afarn )  / 0.17  / ! alpha-farnescene                       
132   DATA  tdf_prm ( imgn_bcar  )  / 0.17  / ! beta-caryophyllene                     
133   DATA  tdf_prm ( imgn_meoh  )  / 0.08  / ! methanol                               
134   DATA  tdf_prm ( imgn_acto  )  / 0.11  / ! acetone                                
135   DATA  tdf_prm ( imgn_acta  )  / 0.13  / ! acetaldehyde and ethanol               
136   DATA  tdf_prm ( imgn_form  )  / 0.09  / ! formic acid, formaldehyde, acetic acid
137   DATA  tdf_prm ( imgn_ch4   )  / 0.05  / ! methane                                
138   DATA  tdf_prm ( imgn_no    )  / 0.11  / ! nitric oxide (and NO2 and NH3)                          
139   DATA  tdf_prm ( imgn_omtp  )  / 0.10  / ! other monoterpenes                     
140   DATA  tdf_prm ( imgn_osqt  )  / 0.17  / ! other sesquiterpenes                   
141   DATA  tdf_prm ( imgn_co    )  / 0.09  / ! co and other voc                       
143   !=====================================================================
144   !  The following contains the fraction of total emisssions
145   !  for which light-dependency is applied (via GAMMA_P).  This is the
146   !  parameter "LDF" provided in Table 2 of MEGAN v2.0 Documentaion.
147   !  Originally from Tan's LD_FCT.EXT (11/30/2006)
148   !  Updated 08/13/2007 based on Alex Guenther's email.
149   !=====================================================================
151   REAL , DIMENSION(n_mgn_spc) :: ldf_fct
153   DATA  ldf_fct ( imgn_isop  )  / 0.9999 /  ! isoprene                               
154   DATA  ldf_fct ( imgn_mbo   )  / 0.9999 /  ! MBO (only (2-methyl-3-ene-2-ol)-methylbutenol)
155   DATA  ldf_fct ( imgn_myrc  )  / 0.05   /  ! myrcene                                
156   DATA  ldf_fct ( imgn_sabi  )  / 0.10   /  ! sabinene                               
157   DATA  ldf_fct ( imgn_limo  )  / 0.05   /  ! limonene                               
158   DATA  ldf_fct ( imgn_3car  )  / 0.05   /  ! 3 carene                               
159   DATA  ldf_fct ( imgn_ocim  )  / 0.80   /  ! trans-beta-ocimene                     
160   DATA  ldf_fct ( imgn_bpin  )  / 0.10   /  ! beta-pinene                            
161   DATA  ldf_fct ( imgn_apin  )  / 0.10   /  ! alpha-pinene                           
162   DATA  ldf_fct ( imgn_afarn )  / 0.50   /  ! alpha-farnescene                       
163   DATA  ldf_fct ( imgn_bcar  )  / 0.50   /  ! beta-caryophyllene                     
164   DATA  ldf_fct ( imgn_meoh  )  / 0.75   /  ! methanol                         
165   DATA  ldf_fct ( imgn_acto  )  / 0.25   /  ! acetone                                
166   DATA  ldf_fct ( imgn_acta  )  / 0.50   /  ! acetaldehyde and ethanol ??          
167   DATA  ldf_fct ( imgn_form  )  / 0.50   /  ! formic acid, formaldehyde, acetic acid ??
168   DATA  ldf_fct ( imgn_ch4   )  / 0.75   /  ! methane                                
169   DATA  ldf_fct ( imgn_no    )  / 0.00   /  ! nitric oxide (and NO2 and NH3 )                          
170   DATA  ldf_fct ( imgn_omtp  )  / 0.10   /  ! other monoterpenes                     
171   DATA  ldf_fct ( imgn_osqt  )  / 0.50   /  ! other sesquiterpenes                   
172   DATA  ldf_fct ( imgn_co    )  / 0.50   /  ! CO and other VOC                       
174   !=====================================================================
175   !  The following contains "production and loss within canopy" factors.
176   !  See Eqn (1)  and Section 3.3 of Guenther et al. [2006].
177   !  Originally from Tan's PDT_LOS_CP.EXT (11/30/2006)
178   !=====================================================================
180   REAL , DIMENSION(n_mgn_spc) :: rho_fct
181   
182   ! RHO_FCT are all 1.0 right now
183   DATA  rho_fct ( imgn_isop  )  / 1.0  / ! isoprene                               
184   DATA  rho_fct ( imgn_mbo   )  / 1.0  / ! MBO (only (2-methyl-3-ene-2-ol)-methylbutenol)
185   DATA  rho_fct ( imgn_myrc  )  / 1.0  / ! myrcene                                
186   DATA  rho_fct ( imgn_sabi  )  / 1.0  / ! sabinene                               
187   DATA  rho_fct ( imgn_limo  )  / 1.0  / ! limonene                               
188   DATA  rho_fct ( imgn_3car  )  / 1.0  / ! 3 carene                               
189   DATA  rho_fct ( imgn_ocim  )  / 1.0  / ! trans-beta-ocimene                     
190   DATA  rho_fct ( imgn_bpin  )  / 1.0  / ! beta-pinene                            
191   DATA  rho_fct ( imgn_apin  )  / 1.0  / ! alpha-pinene                           
192   DATA  rho_fct ( imgn_afarn )  / 1.0  / ! alpha-farnescene                       
193   DATA  rho_fct ( imgn_bcar  )  / 1.0  / ! beta-caryophyllene                     
194   DATA  rho_fct ( imgn_meoh  )  / 1.0  / ! methanol                               
195   DATA  rho_fct ( imgn_acto  )  / 1.0  / ! acetone                                
196   DATA  rho_fct ( imgn_acta  )  / 1.0  / ! acetaldehyde and ethanol               
197   DATA  rho_fct ( imgn_form  )  / 1.0  / ! formic acid, formaldehyde, acetic acid 
198   DATA  rho_fct ( imgn_ch4   )  / 1.0  / ! methane                                
199   DATA  rho_fct ( imgn_no    )  / 1.0  / ! nitric oxide (and NO2 and NH3 )                           
200   DATA  rho_fct ( imgn_omtp  )  / 1.0  / ! other monoterpenes                     
201   DATA  rho_fct ( imgn_osqt  )  / 1.0  / ! other sesquiterpenes                   
202   DATA  rho_fct ( imgn_co    )  / 1.0  / ! CO and other VOC                       
205   !=====================================================================
206   !  The following contains relative emission rates assigned to each
207   !  canopy fraction.  See Eqn (16) of Guenther et al. [2006] and
208   !  MEGAN v2.0 Documentation.
209   !  Originally Tan's REL_EM_ACT.EXT (11/30/06)
210   !=====================================================================
212   INTEGER, PARAMETER ::        n_cat=5
213   REAL, DIMENSION(n_cat) ::    anew, agro, amat, aold
215   DATA  anew (1) , agro (1) , amat (1) , aold (1)  /  1.0  , 1.0 , 1.0   , 1.0  /
216   DATA  anew (2) , agro (2) , amat (2) , aold (2)  /  2.0  , 1.8 , 0.95  , 1.0  /
217   DATA  anew (3) , agro (3) , amat (3) , aold (3)  /  0.4  , 0.6 , 1.075 , 1.0  /
218   DATA  anew (4) , agro (4) , amat (4) , aold (4)  /  3.0  , 2.6 , 0.85  , 1.0  /
219   DATA  anew (5) , agro (5) , amat (5) , aold (5)  /  0.05 , 0.6 , 1.125 , 1.0  /
221   !=====================================================================
222   !  The following contains the base emission factors (microgram m-2 hr-1) for
223   !  n_mgn_spc= 20 MEGAN v2.04 classes of compounds for the following
224   !  four plant functional types:
225   !     BT = breadleaf tree
226   !     NT = needle leaf tree
227   !     SB = shrub and brush
228   !     HB = herb/crops/grasses
229   !  Originally Tan's EF_MGN20_Alex.EXT
230   !  Updated 08/13/2007 based on Alex's email.
231   !=====================================================================
232 !shc  real, dimension (n_mgn_spc) :: EF_BT, EF_NT, EF_SB, EF_HB
233   INTEGER, PARAMETER :: n_pft = 4 ! number of plant functional types
234   REAL, DIMENSION (n_mgn_spc,n_pft) :: EF
236   ! Plant functional type indices
237   !  k_bt, k_nt, k_sb, and k_bh are hardwired by the column-order
238   !  of the data entry for EF and EF_frac.  Do not change their values or
239   !  order of the columns in the data entry.
241   INTEGER, PARAMETER :: k_bt = 1, k_nt = 2, k_sb=3, k_hb=4
242   
243   !                                   BT        NT        SB          HB
244   DATA EF( imgn_isop  , 1:n_pft )  /  13000.00 , 2000.00 , 11000.00 ,  400.00  / !not used for isoprene
245   DATA EF( imgn_mbo   , 1:n_pft )  /      0.10 ,  100.00 ,     1.00 ,    0.01  /
246   DATA EF( imgn_myrc  , 1:n_pft )  /     20.00 ,   75.00 ,    22.00 ,    0.30  /
247   DATA EF( imgn_sabi  , 1:n_pft )  /     45.00 ,   70.00 ,    50.00 ,    0.70  /
248   DATA EF( imgn_limo  , 1:n_pft )  /     45.00 ,  100.00 ,    52.00 ,    0.70  /
249   DATA EF( imgn_3car  , 1:n_pft )  /     18.00 ,  160.00 ,    25.00 ,    0.30  /
250   DATA EF( imgn_ocim  , 1:n_pft )  /     90.00 ,   60.00 ,    85.00 ,    1.00  /
251   DATA EF( imgn_bpin  , 1:n_pft )  /     90.00 ,  300.00 ,   100.00 ,    1.50  /
252   DATA EF( imgn_apin  , 1:n_pft )  /    180.00 ,  450.00 ,   200.00 ,    2.00  /
253   DATA EF( imgn_afarn , 1:n_pft )  /     35.00 ,   30.00 ,    30.00 ,    0.50  /
254   DATA EF( imgn_bcar  , 1:n_pft )  /     30.00 ,   60.00 ,    45.00 ,    0.90  /
255   DATA EF( imgn_meoh  , 1:n_pft )  /    800.00 ,  800.00 ,   800.00 ,  800.00  /
256   DATA EF( imgn_acto  , 1:n_pft )  /    240.00 ,  240.00 ,   240.00 ,   80.00  /
257   DATA EF( imgn_acta  , 1:n_pft )  /    240.00 ,  240.00 ,   240.00 ,   80.00  /
258   DATA EF( imgn_form  , 1:n_pft )  /     70.00 ,   70.00 ,    70.00 ,   70.00  /
259   DATA EF( imgn_ch4   , 1:n_pft )  /     30.00 ,   30.00 ,    30.00 ,   30.00  /
260   DATA EF( imgn_no    , 1:n_pft )  /      5.00 ,    6.00 ,    30.00 ,   70.00  /
261   DATA EF( imgn_omtp  , 1:n_pft )  /     90.00 ,  180.00 ,   110.00 ,    4.80  /
262   DATA EF( imgn_osqt  , 1:n_pft )  /     75.00 ,  110.00 ,    85.00 ,    1.40  /
263   DATA EF( imgn_co    , 1:n_pft )  /   1000.00 , 1000.00 ,  1000.00 , 1000.00  /
266   !=======================================================================
267   ! The following contains MEGAN 138 speciated species for MEGAN v2.04.
268   ! See Table 3 of MEGAN v2.0 Documentation.
269   ! Originally from Tan's SPC_SPCAT.EXT (12/02/2006)
270   !=======================================================================
273   INTEGER, PARAMETER :: n_spca_spc = 138                         ! number of speciated species
274   CHARACTER(len=n_spc_char), DIMENSION(n_spca_spc) :: spca_spc   ! speciated species names
275   REAL                     , DIMENSION(n_spca_spc) :: spca_mwt   ! speciated species molecular weight
277   !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
278   ! _a  = alpha     , _b  = beta      , _c  = cis       , _al = allo    ,
279   ! _g  = gamma     , _d  = delta     , _t  = trans     , _m  = methyl  ,
280   ! _p  = para      , _o  = ortho     , _e  = ene       , _ol = ol      ,
281   ! met = methyl    , 2met= dimethyl  , MBO = methylbutenol             ,
282   ! 2s  = disulfide , s   = sulfide   , OXD = oxide     , ACT = acetate ,
283   ! PPPP= propenylpropyl              , DCTT= decatetraene              ,
284   ! CCO = acetaldehyde
285   !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
287   ! Indices for 138 MEGAN speciated species
288   !    The ordering below follows that of Tan's SPC_SPCA.TXT file.
289   !    For programming, the ordering does not matter, except that 
290   !    is_isoprene should always be 1.
291   INTEGER, PARAMETER ::              &
292        &  is_isoprene         =   1  & ! isoprene
293        & ,is_myrcene          =   2  & ! myrcene
294        & ,is_sabinene         =   3  & ! sabiene
295        & ,is_limonene         =   4  & ! limonene
296        & ,is_carene_3         =   5  & ! 3-carene
297        & ,is_ocimene_t_b      =   6  & ! trans-beta-ocimene
298        & ,is_pinene_b         =   7  & ! beta-pinene
299        & ,is_pinene_a         =   8  & ! alpha-pinene
300        & ,is_2met_styrene     =   9  & ! dimethyl styrene
301        & ,is_cymene_p         =  10  & ! para-cymene
302        & ,is_cymene_o         =  11  & ! ortho-cymene
303        & ,is_phellandrene_a   =  12  & ! alpha-phellandrene
304        & ,is_thujene_a        =  13  & ! alpha-thujene
305        & ,is_terpinene_a      =  14  & ! alpha-terpinene
306        & ,is_terpinene_g      =  15  & ! gamma-terpinene
307        & ,is_terpinolene      =  16  & ! terpinolene
308        & ,is_phellandrene_b   =  17  & ! beta-phellandrene
309        & ,is_camphene         =  18  & ! camphene
310        & ,is_bornene          =  19  & ! bornene
311        & ,is_fenchene_a       =  20  & ! alpha-fenchene
312        & ,is_ocimene_al       =  21  & ! allo ocimene
313        & ,is_ocimene_c_b      =  22  & ! cis beta-ocimene
314        & ,is_tricyclene       =  23  & ! tricyclene
315        & ,is_estragole        =  24  & ! estrogole
316        & ,is_camphor          =  25  & ! camphor
317        & ,is_fenchone         =  26  & ! fenchone
318        & ,is_piperitone       =  27  & ! piperitone
319        & ,is_thujone_a        =  28  & ! alpha-thujone
320        & ,is_thujone_b        =  29  & ! beta thujone
321        & ,is_cineole_1_8      =  30  & ! (1,8)-cineole
322        & ,is_borneol          =  31  & ! borneol
323        & ,is_linalool         =  32  & ! linalool
324        & ,is_terpineol_4      =  33  & ! 4-terpineol
325        & ,is_terpineol_a      =  34  & ! alpha-terpineol
326        & ,is_linalool_OXD_c   =  35  & ! cis linalool oxide
327        & ,is_linalool_OXD_t   =  36  & ! trans linalool oxide
328        & ,is_ionone_b         =  37  & ! beta-ionone
329        & ,is_bornyl_ACT       =  38  & ! bornyl acetate
330        & ,is_farnescene_a     =  39  & ! alpha-farnescene
331        & ,is_caryophyllene_b  =  40  & ! beta-caryophyllene
332        & ,is_acoradiene       =  41  & ! acoradiene
333        & ,is_aromadendrene    =  42  & ! aromadendrene
334        & ,is_bergamotene_a    =  43  & ! alpha-bergamotene
335        & ,is_bergamotene_b    =  44  & ! beta-bergamotene
336        & ,is_bisabolene_a     =  45  & ! alpha-bisabolene
337        & ,is_bisabolene_b     =  46  & ! beta-bisabolene
338        & ,is_bourbonene_b     =  47  & ! beta-bourbonene
339        & ,is_cadinene_d       =  48  & ! delta-cadinene
340        & ,is_cadinene_g       =  49  & ! gamma-cadinene
341        & ,is_cedrene_a        =  50  & ! alpha-cedrene
342        & ,is_copaene_a        =  51  & ! alpha-copaene
343        & ,is_cubebene_a       =  52  & ! alpha-cubebene 
344        & ,is_cubebene_b       =  53  & ! beta-cubebene
345        & ,is_elemene_b        =  54  & ! beta-elemene
346        & ,is_farnescene_b     =  55  & ! beta-farnescene
347        & ,is_germacrene_B     =  56  & ! B germacrene
348        & ,is_germacrene_D     =  57  & ! D germacrene
349        & ,is_gurjunene_b      =  58  & ! beta-gurjunene
350        & ,is_humulene_a       =  59  & ! alph-humulene
351        & ,is_humulene_g       =  60  & ! gamma-humulene
352        & ,is_isolongifolene   =  61  & ! isolongifolene
353        & ,is_longifolene      =  62  & ! longifolene
354        & ,is_longipinene      =  63  & ! longipinene
355        & ,is_muurolene_a      =  64  & ! alpha-muurolene
356        & ,is_muurolene_g      =  65  & ! gamma-muurolene
357        & ,is_selinene_b       =  66  & ! beta-selinene
358        & ,is_selinene_d       =  67  & ! delta-selinene
359        & ,is_nerolidol_c      =  68  & ! cis-nerolidol
360        & ,is_nerolidol_t      =  69  & ! trans-nerolidol
361        & ,is_cedrol           =  70  & ! cedrol
362        & ,is_MBO_2m3e2ol      =  71  & !  
363        & ,is_methanol         =  72  & ! methanol
364        & ,is_acetone          =  73  & ! acetone
365        & ,is_methane          =  74  & ! methane
366        & ,is_ammonia          =  75  & ! ammonia
367        & ,is_nitrous_OXD      =  76  & ! nitrous oxide
368        & ,is_nitric_OXD       =  77  & ! nitric oxide
369        & ,is_acetaldehyde     =  78  & ! acetaldehyde
370        & ,is_ethanol          =  79  & ! ethanol
371        & ,is_formic_acid      =  80  & ! formic acid
372        & ,is_formaldehyde     =  81  & ! formaldehye
373        & ,is_acetic_acid      =  82  & ! acetic acid
374        & ,is_MBO_3m2e1ol      =  83  & ! (3-methyl-2-ene-1-ol)-methybutenol 
375        & ,is_MBO_3m3e1ol      =  84  & ! (3-methyl-3-ene-1-ol)-methybutenol 
376        & ,is_benzaldehyde     =  85  & ! benzaldehye
377        & ,is_butanone_2       =  86  & ! 2-butonene
378        & ,is_decanal          =  87  & ! decanal
379        & ,is_dodecene_1       =  88  & ! 1-dodecene
380        & ,is_geranyl_acetone  =  89  & ! geranyl acetone
381        & ,is_heptanal         =  90  & ! heptanal
382        & ,is_heptane          =  91  & ! heptane
383        & ,is_hexane           =  92  & ! hexane
384        & ,is_met_benzoate     =  93  & ! methyl-benzoate
385        & ,is_met_heptenone    =  94  & ! methyl-heptonene
386        & ,is_neryl_acetone    =  95  & ! neryl-acetone
387        & ,is_nonanal          =  96  & ! nonanal
388        & ,is_nonenal          =  97  & ! nonenal
389        & ,is_octanal          =  98  & ! octanal
390        & ,is_octanol          =  99  & ! octanol
391        & ,is_octenol_1e3ol    = 100  & ! (1-ene-3-ol)-octenol
392        & ,is_oxopentanal      = 101  & ! oxopentanal
393        & ,is_pentane          = 102  & ! pentane
394        & ,is_phenyl_CCO       = 103  & ! phenylacetaldehyde
395        & ,is_pyruvic_acid     = 104  & ! pyruvic acid
396        & ,is_terpinyl_ACT_a   = 105  & ! alpha-terpinyl-acetate
397        & ,is_tetradecene_1    = 106  & ! 1-tetradecene
398        & ,is_toluene          = 107  & ! toluene
399        & ,is_carbon_monoxide  = 108  & ! CO
400        & ,is_butene           = 109  & ! butene
401        & ,is_ethane           = 110  & ! ethane
402        & ,is_ethene           = 111  & ! ethene
403        & ,is_hydrogen_cyanide = 112  & ! hydrogen cyanide
404        & ,is_propane          = 113  & ! propane
405        & ,is_propene          = 114  & ! propene
406        & ,is_carbon_2s        = 115  & ! carbon disulfide
407        & ,is_carbonyl_s       = 116  & ! caronyl sulfide
408        & ,is_diallyl_2s       = 117  & ! diallyl disulfide
409        & ,is_2met_2s          = 118  & ! dimethyl disulfide
410        & ,is_2met_s           = 119  & ! dimethyl sulfide
411        & ,is_met_chloride     = 120  & ! methyl chloride
412        & ,is_met_bromide      = 121  & ! methyl bromide
413        & ,is_met_iodide       = 122  & ! methyl iodide
414        & ,is_hydrogen_s       = 123  & ! hydrogen sulfide
415        & ,is_met_mercaptan    = 124  & ! methyl mercaptan
416        & ,is_met_propenyl_2s  = 125  & ! methyl propenyl disulfide
417        & ,is_PPPP_2s          = 126  & ! propenylpropyldisulfdie
418        & ,is_2met_nonatriene  = 127  & ! dimethyl nonatriene
419        & ,is_met_salicylate   = 128  & ! methyl salicylate
420        & ,is_indole           = 129  & ! indole
421        & ,is_jasmone          = 130  & ! jasmone
422        & ,is_met_jasmonate    = 131  & ! methyl jasmonate
423        & ,is_3met_3DCTT       = 132  & ! trimethyl tridecatetraene
424        & ,is_hexanal          = 133  & ! hexanal
425        & ,is_hexanol_1        = 134  & ! 1-hexanol
426        & ,is_hexenal_c3       = 135  & ! cis-3-hexenal
427        & ,is_hexenal_t2       = 136  & ! trans-2-hexenal
428        & ,is_hexenol_c3       = 137  & ! cis-3-hexenol
429        & ,is_hexenyl_ACT_c3   = 138    ! cis-3-hexenyl acetate
431   ! Names and Molecular Weights of 138 MEGAN speciated species  
432   ! (The names are not used in WRF-Chem, except for testing stage)
433   DATA spca_spc(is_isoprene        ), spca_mwt(is_isoprene        ) / 'isoprene        ' ,  68.12 /
434   DATA spca_spc(is_myrcene         ), spca_mwt(is_myrcene         ) / 'myrcene         ' , 136.23 /
435   DATA spca_spc(is_sabinene        ), spca_mwt(is_sabinene        ) / 'sabinene        ' , 136.23 /
436   DATA spca_spc(is_limonene        ), spca_mwt(is_limonene        ) / 'limonene        ' , 136.23 /
437   DATA spca_spc(is_carene_3        ), spca_mwt(is_carene_3        ) / 'carene_3        ' , 136.23 /
438   DATA spca_spc(is_ocimene_t_b     ), spca_mwt(is_ocimene_t_b     ) / 'ocimene_t_b     ' , 136.23 /
439   DATA spca_spc(is_pinene_b        ), spca_mwt(is_pinene_b        ) / 'pinene_b        ' , 136.23 /
440   DATA spca_spc(is_pinene_a        ), spca_mwt(is_pinene_a        ) / 'pinene_a        ' , 136.23 /
441   DATA spca_spc(is_2met_styrene    ), spca_mwt(is_2met_styrene    ) / '2met_styrene    ' , 132.20 /
442   DATA spca_spc(is_cymene_p        ), spca_mwt(is_cymene_p        ) / 'cymene_p        ' , 134.22 /
443   DATA spca_spc(is_cymene_o        ), spca_mwt(is_cymene_o        ) / 'cymene_o        ' , 134.22 /
444   DATA spca_spc(is_phellandrene_a  ), spca_mwt(is_phellandrene_a  ) / 'phellandrene_a  ' , 136.23 /
445   DATA spca_spc(is_thujene_a       ), spca_mwt(is_thujene_a       ) / 'thujene_a       ' , 136.23 /
446   DATA spca_spc(is_terpinene_a     ), spca_mwt(is_terpinene_a     ) / 'terpinene_a     ' , 136.23 /
447   DATA spca_spc(is_terpinene_g     ), spca_mwt(is_terpinene_g     ) / 'terpinene_g     ' , 136.23 /
448   DATA spca_spc(is_terpinolene     ), spca_mwt(is_terpinolene     ) / 'terpinolene     ' , 136.23 /
449   DATA spca_spc(is_phellandrene_b  ), spca_mwt(is_phellandrene_b  ) / 'phellandrene_b  ' , 136.23 /
450   DATA spca_spc(is_camphene        ), spca_mwt(is_camphene        ) / 'camphene        ' , 136.23 /
451   DATA spca_spc(is_bornene         ), spca_mwt(is_bornene         ) / 'bornene         ' , 136.23 /
452   DATA spca_spc(is_fenchene_a      ), spca_mwt(is_fenchene_a      ) / 'fenchene_a      ' , 136.23 /
453   DATA spca_spc(is_ocimene_al      ), spca_mwt(is_ocimene_al      ) / 'ocimene_al      ' , 136.23 /
454   DATA spca_spc(is_ocimene_c_b     ), spca_mwt(is_ocimene_c_b     ) / 'ocimene_c_b     ' , 136.23 /
455   DATA spca_spc(is_tricyclene      ), spca_mwt(is_tricyclene      ) / 'tricyclene      ' , 136.23 /
456   DATA spca_spc(is_estragole       ), spca_mwt(is_estragole       ) / 'estragole       ' , 148.20 /
457   DATA spca_spc(is_camphor         ), spca_mwt(is_camphor         ) / 'camphor         ' , 152.23 /
458   DATA spca_spc(is_fenchone        ), spca_mwt(is_fenchone        ) / 'fenchone        ' , 152.23 /
459   DATA spca_spc(is_piperitone      ), spca_mwt(is_piperitone      ) / 'piperitone      ' , 152.23 /
460   DATA spca_spc(is_thujone_a       ), spca_mwt(is_thujone_a       ) / 'thujone_a       ' , 152.23 /
461   DATA spca_spc(is_thujone_b       ), spca_mwt(is_thujone_b       ) / 'thujone_b       ' , 152.23 /
462   DATA spca_spc(is_cineole_1_8     ), spca_mwt(is_cineole_1_8     ) / 'cineole_1_8     ' , 154.25 /
463   DATA spca_spc(is_borneol         ), spca_mwt(is_borneol         ) / 'borneol         ' , 154.25 /
464   DATA spca_spc(is_linalool        ), spca_mwt(is_linalool        ) / 'linalool        ' , 154.25 /
465   DATA spca_spc(is_terpineol_4     ), spca_mwt(is_terpineol_4     ) / 'terpineol_4     ' , 154.25 /
466   DATA spca_spc(is_terpineol_a     ), spca_mwt(is_terpineol_a     ) / 'terpineol_a     ' , 154.25 /
467   DATA spca_spc(is_linalool_oxd_c  ), spca_mwt(is_linalool_oxd_c  ) / 'linalool_oxd_c  ' , 170.25 /
468   DATA spca_spc(is_linalool_oxd_t  ), spca_mwt(is_linalool_oxd_t  ) / 'linalool_oxd_t  ' , 170.25 /
469   DATA spca_spc(is_ionone_b        ), spca_mwt(is_ionone_b        ) / 'ionone_b        ' , 192.30 /
470   DATA spca_spc(is_bornyl_act      ), spca_mwt(is_bornyl_act      ) / 'bornyl_act      ' , 196.29 /
471   DATA spca_spc(is_farnescene_a    ), spca_mwt(is_farnescene_a    ) / 'farnescene_a    ' , 204.35 /
472   DATA spca_spc(is_caryophyllene_b ), spca_mwt(is_caryophyllene_b ) / 'caryophyllene_b ' , 204.35 /
473   DATA spca_spc(is_acoradiene      ), spca_mwt(is_acoradiene      ) / 'acoradiene      ' , 204.35 /
474   DATA spca_spc(is_aromadendrene   ), spca_mwt(is_aromadendrene   ) / 'aromadendrene   ' , 204.35 /
475   DATA spca_spc(is_bergamotene_a   ), spca_mwt(is_bergamotene_a   ) / 'bergamotene_a   ' , 204.35 /
476   DATA spca_spc(is_bergamotene_b   ), spca_mwt(is_bergamotene_b   ) / 'bergamotene_b   ' , 204.35 /
477   DATA spca_spc(is_bisabolene_a    ), spca_mwt(is_bisabolene_a    ) / 'bisabolene_a    ' , 204.35 /
478   DATA spca_spc(is_bisabolene_b    ), spca_mwt(is_bisabolene_b    ) / 'bisabolene_b    ' , 204.35 /
479   DATA spca_spc(is_bourbonene_b    ), spca_mwt(is_bourbonene_b    ) / 'bourbonene_b    ' , 204.35 /
480   DATA spca_spc(is_cadinene_d      ), spca_mwt(is_cadinene_d      ) / 'cadinene_d      ' , 204.35 /
481   DATA spca_spc(is_cadinene_g      ), spca_mwt(is_cadinene_g      ) / 'cadinene_g      ' , 204.35 /
482   DATA spca_spc(is_cedrene_a       ), spca_mwt(is_cedrene_a       ) / 'cedrene_a       ' , 204.35 /
483   DATA spca_spc(is_copaene_a       ), spca_mwt(is_copaene_a       ) / 'copaene_a       ' , 204.35 /
484   DATA spca_spc(is_cubebene_a      ), spca_mwt(is_cubebene_a      ) / 'cubebene_a      ' , 204.35 /
485   DATA spca_spc(is_cubebene_b      ), spca_mwt(is_cubebene_b      ) / 'cubebene_b      ' , 204.35 /
486   DATA spca_spc(is_elemene_b       ), spca_mwt(is_elemene_b       ) / 'elemene_b       ' , 204.35 /
487   DATA spca_spc(is_farnescene_b    ), spca_mwt(is_farnescene_b    ) / 'farnescene_b    ' , 204.35 /
488   DATA spca_spc(is_germacrene_b    ), spca_mwt(is_germacrene_B    ) / 'germacrene_b    ' , 204.35 /
489   DATA spca_spc(is_germacrene_d    ), spca_mwt(is_germacrene_D    ) / 'germacrene_d    ' , 204.35 /
490   DATA spca_spc(is_gurjunene_b     ), spca_mwt(is_gurjunene_b     ) / 'gurjunene_b     ' , 204.35 /
491   DATA spca_spc(is_humulene_a      ), spca_mwt(is_humulene_a      ) / 'humulene_a      ' , 204.35 /
492   DATA spca_spc(is_humulene_g      ), spca_mwt(is_humulene_g      ) / 'humulene_g      ' , 204.35 /
493   DATA spca_spc(is_isolongifolene  ), spca_mwt(is_isolongifolene  ) / 'isolongifolene  ' , 204.35 /
494   DATA spca_spc(is_longifolene     ), spca_mwt(is_longifolene     ) / 'longifolene     ' , 204.35 /
495   DATA spca_spc(is_longipinene     ), spca_mwt(is_longipinene     ) / 'longipinene     ' , 204.35 /
496   DATA spca_spc(is_muurolene_a     ), spca_mwt(is_muurolene_a     ) / 'muurolene_a     ' , 204.35 /
497   DATA spca_spc(is_muurolene_g     ), spca_mwt(is_muurolene_g     ) / 'muurolene_g     ' , 204.35 /
498   DATA spca_spc(is_selinene_b      ), spca_mwt(is_selinene_b      ) / 'selinene_b      ' , 204.35 /
499   DATA spca_spc(is_selinene_d      ), spca_mwt(is_selinene_d      ) / 'selinene_d      ' , 204.35 /
500   DATA spca_spc(is_nerolidol_c     ), spca_mwt(is_nerolidol_c     ) / 'nerolidol_c     ' , 222.37 /
501   DATA spca_spc(is_nerolidol_t     ), spca_mwt(is_nerolidol_t     ) / 'nerolidol_t     ' , 222.37 /
502   DATA spca_spc(is_cedrol          ), spca_mwt(is_cedrol          ) / 'cedrol          ' , 222.37 /
503   DATA spca_spc(is_mbo_2m3e2ol     ), spca_mwt(is_mbo_2m3e2ol     ) / 'mbo_2m3e2ol     ' , 86.13  /
504   DATA spca_spc(is_methanol        ), spca_mwt(is_methanol        ) / 'methanol        ' , 32.04  /
505   DATA spca_spc(is_acetone         ), spca_mwt(is_acetone         ) / 'acetone         ' , 58.08  /
506   DATA spca_spc(is_methane         ), spca_mwt(is_methane         ) / 'methane         ' , 16.04  /
507   DATA spca_spc(is_ammonia         ), spca_mwt(is_ammonia         ) / 'ammonia         ' , 17.03  /
508   DATA spca_spc(is_nitrous_oxd     ), spca_mwt(is_nitrous_oxd     ) / 'nitrous_oxd     ' , 44.01  /
509   DATA spca_spc(is_nitric_oxd      ), spca_mwt(is_nitric_oxd      ) / 'nitric_oxd      ' , 30.01  /
510   DATA spca_spc(is_acetaldehyde    ), spca_mwt(is_acetaldehyde    ) / 'acetaldehyde    ' , 44.05  /
511   DATA spca_spc(is_ethanol         ), spca_mwt(is_ethanol         ) / 'ethanol         ' , 46.07  /
512   DATA spca_spc(is_formic_acid     ), spca_mwt(is_formic_acid     ) / 'formic_acid     ' , 46.03  /
513   DATA spca_spc(is_formaldehyde    ), spca_mwt(is_formaldehyde    ) / 'formaldehyde    ' , 30.03  /
514   DATA spca_spc(is_acetic_acid     ), spca_mwt(is_acetic_acid     ) / 'acetic_acid     ' , 60.05  /
515   DATA spca_spc(is_mbo_3m2e1ol     ), spca_mwt(is_mbo_3m2e1ol     ) / 'mbo_3m2e1ol     ' , 86.13  /
516   DATA spca_spc(is_mbo_3m3e1ol     ), spca_mwt(is_mbo_3m3e1ol     ) / 'mbo_3m3e1ol     ' , 86.13  /
517   DATA spca_spc(is_benzaldehyde    ), spca_mwt(is_benzaldehyde    ) / 'benzaldehyde    ' , 106.12 /
518   DATA spca_spc(is_butanone_2      ), spca_mwt(is_butanone_2      ) / 'butanone_2      ' , 72.11  /
519   DATA spca_spc(is_decanal         ), spca_mwt(is_decanal         ) / 'decanal         ' , 156.27 /
520   DATA spca_spc(is_dodecene_1      ), spca_mwt(is_dodecene_1      ) / 'dodecene_1      ' , 168.32 /
521   DATA spca_spc(is_geranyl_acetone ), spca_mwt(is_geranyl_acetone ) / 'geranyl_acetone ' , 194.31 /
522   DATA spca_spc(is_heptanal        ), spca_mwt(is_heptanal        ) / 'heptanal        ' , 114.19 /
523   DATA spca_spc(is_heptane         ), spca_mwt(is_heptane         ) / 'heptane         ' , 100.20 /
524   DATA spca_spc(is_hexane          ), spca_mwt(is_hexane          ) / 'hexane          ' , 86.18  /
525   DATA spca_spc(is_met_benzoate    ), spca_mwt(is_met_benzoate    ) / 'met_benzoate    ' , 136.15 /
526   DATA spca_spc(is_met_heptenone   ), spca_mwt(is_met_heptenone   ) / 'met_heptenone   ' , 126.20 /
527   DATA spca_spc(is_neryl_acetone   ), spca_mwt(is_neryl_acetone   ) / 'neryl_acetone   ' , 194.31 /
528   DATA spca_spc(is_nonanal         ), spca_mwt(is_nonanal         ) / 'nonanal         ' , 142.24 /
529   DATA spca_spc(is_nonenal         ), spca_mwt(is_nonenal         ) / 'nonenal         ' , 140.22 /
530   DATA spca_spc(is_octanal         ), spca_mwt(is_octanal         ) / 'octanal         ' , 128.21 /
531   DATA spca_spc(is_octanol         ), spca_mwt(is_octanol         ) / 'octanol         ' , 130.23 /
532   DATA spca_spc(is_octenol_1e3ol   ), spca_mwt(is_octenol_1e3ol   ) / 'octenol_1e3ol   ' , 128.21 /
533   DATA spca_spc(is_oxopentanal     ), spca_mwt(is_oxopentanal     ) / 'oxopentanal     ' , 100.12 /
534   DATA spca_spc(is_pentane         ), spca_mwt(is_pentane         ) / 'pentane         ' , 72.15  /
535   DATA spca_spc(is_phenyl_cco      ), spca_mwt(is_phenyl_cco      ) / 'phenyl_cco      ' , 120.15 /
536   DATA spca_spc(is_pyruvic_acid    ), spca_mwt(is_pyruvic_acid    ) / 'pyruvic_acid    ' , 88.06  /
537   DATA spca_spc(is_terpinyl_act_a  ), spca_mwt(is_terpinyl_act_a  ) / 'terpinyl_act_a  ' , 196.29 /
538   DATA spca_spc(is_tetradecene_1   ), spca_mwt(is_tetradecene_1   ) / 'tetradecene_1   ' , 196.37 /
539   DATA spca_spc(is_toluene         ), spca_mwt(is_toluene         ) / 'toluene         ' , 92.14  /
540   DATA spca_spc(is_carbon_monoxide ), spca_mwt(is_carbon_monoxide ) / 'carbon_monoxide ' , 28.01  /
541   DATA spca_spc(is_butene          ), spca_mwt(is_butene          ) / 'butene          ' , 56.11  /
542   DATA spca_spc(is_ethane          ), spca_mwt(is_ethane          ) / 'ethane          ' , 30.07  /
543   DATA spca_spc(is_ethene          ), spca_mwt(is_ethene          ) / 'ethene          ' , 28.05  /
544   DATA spca_spc(is_hydrogen_cyanide), spca_mwt(is_hydrogen_cyanide) / 'hydrogen_cyanide' , 27.03  /
545   DATA spca_spc(is_propane         ), spca_mwt(is_propane         ) / 'propane         ' , 44.10  /
546   DATA spca_spc(is_propene         ), spca_mwt(is_propene         ) / 'propene         ' , 42.08  /
547   DATA spca_spc(is_carbon_2s       ), spca_mwt(is_carbon_2s       ) / 'carbon_2s       ' , 76.14  /
548   DATA spca_spc(is_carbonyl_s      ), spca_mwt(is_carbonyl_s      ) / 'carbonyl_s      ' , 60.08  /
549   DATA spca_spc(is_diallyl_2s      ), spca_mwt(is_diallyl_2s      ) / 'diallyl_2s      ' , 146.28 /
550   DATA spca_spc(is_2met_2s         ), spca_mwt(is_2met_2s         ) / '2met_2s         ' , 94.20  /
551   DATA spca_spc(is_2met_s          ), spca_mwt(is_2met_s          ) / '2met_s          ' , 62.14  /
552   DATA spca_spc(is_met_chloride    ), spca_mwt(is_met_chloride    ) / 'met_chloride    ' , 50.49  /
553   DATA spca_spc(is_met_bromide     ), spca_mwt(is_met_bromide     ) / 'met_bromide     ' , 94.94  /
554   DATA spca_spc(is_met_iodide      ), spca_mwt(is_met_iodide      ) / 'met_iodide      ' , 141.94 /
555   DATA spca_spc(is_hydrogen_s      ), spca_mwt(is_hydrogen_s      ) / 'hydrogen_s      ' , 34.08  /
556   DATA spca_spc(is_met_mercaptan   ), spca_mwt(is_met_mercaptan   ) / 'met_mercaptan   ' , 48.11  /
557   DATA spca_spc(is_met_propenyl_2s ), spca_mwt(is_met_propenyl_2s ) / 'met_propenyl_2s ' , 120.24 /
558   DATA spca_spc(is_pppp_2s         ), spca_mwt(is_pppp_2s         ) / 'pppp_2s         ' , 148.29 /
559   DATA spca_spc(is_2met_nonatriene ), spca_mwt(is_2met_nonatriene ) / '2met_nonatriene ' , 150.26 /
560   DATA spca_spc(is_met_salicylate  ), spca_mwt(is_met_salicylate  ) / 'met_salicylate  ' , 152.15 /
561   DATA spca_spc(is_indole          ), spca_mwt(is_indole          ) / 'indole          ' , 117.15 /
562   DATA spca_spc(is_jasmone         ), spca_mwt(is_jasmone         ) / 'jasmone         ' , 164.24 /
563   DATA spca_spc(is_met_jasmonate   ), spca_mwt(is_met_jasmonate   ) / 'met_jasmonate   ' , 224.30 /
564   DATA spca_spc(is_3met_3dctt      ), spca_mwt(is_3met_3dctt      ) / '3met_3dctt      ' , 218.38 /
565   DATA spca_spc(is_hexanal         ), spca_mwt(is_hexanal         ) / 'hexanal         ' , 100.16 /
566   DATA spca_spc(is_hexanol_1       ), spca_mwt(is_hexanol_1       ) / 'hexanol_1       ' , 102.17 /
567   DATA spca_spc(is_hexenal_c3      ), spca_mwt(is_hexenal_c3      ) / 'hexenal_c3      ' , 98.14  /
568   DATA spca_spc(is_hexenal_t2      ), spca_mwt(is_hexenal_t2      ) / 'hexenal_t2      ' , 98.14  /
569   DATA spca_spc(is_hexenol_c3      ), spca_mwt(is_hexenol_c3      ) / 'hexenol_c3      ' , 100.16 /
570   DATA spca_spc(is_hexenyl_act_c3  ), spca_mwt(is_hexenyl_act_c3  ) / 'hexenyl_act_c3  ' , 142.20 /
573   !=====================================================================
574   ! The following is a "table" for mapping nmegan=n_spca_spc=138 
575   ! MEGAN v2.04 species to n_mgn_spc=20 classes of species.  
576   ! See Table 3 of MEGAN v2.0 Documentation.
577   ! Originally from Tan's MAP_MGN20T138.EXT.EXT (12/02/2006)
578   !=====================================================================
579                                                                                                      
580   INTEGER, DIMENSION(n_spca_spc) :: mg20_map
581                                                                                    
582   !---------------------------------------------------------------------
583   ! Some naming convention
584   ! _a  = alpha     , _b  = beta      , _c  = cis       , _al = allo    ,
585   ! _g  = gamma     , _d  = delta     , _t  = trans     , _m  = methyl  ,
586   ! _p  = para      , _o  = ortho     , _e  = ene       , _ol = ol      ,
587   ! met = methyl    , 2met= dimethyl  , MBO = methylbutenol             ,
588   ! 2s  = disulfide , s   = sulfide   , OXD = oxide     , ACT = acetate ,
589   ! PPPP= propenylpropyl              , DCTT= decatetraene              ,
590   ! CCO = acetaldehyde
591   !--------------------------------------------------------------------
593   DATA    mg20_map ( is_isoprene        )        /  imgn_isop     /
594   DATA    mg20_map ( is_myrcene         )        /  imgn_myrc     /
595   DATA    mg20_map ( is_sabinene        )        /  imgn_sabi     /
596   DATA    mg20_map ( is_limonene        )        /  imgn_limo     /
597   DATA    mg20_map ( is_carene_3        )        /  imgn_3car     /
598   DATA    mg20_map ( is_ocimene_t_b     )        /  imgn_ocim     /
599   DATA    mg20_map ( is_pinene_b        )        /  imgn_bpin     /
600   DATA    mg20_map ( is_pinene_a        )        /  imgn_apin     /
601   DATA    mg20_map ( is_2met_styrene    )        /  imgn_omtp     /
602   DATA    mg20_map ( is_cymene_p        )        /  imgn_omtp     /
603   DATA    mg20_map ( is_cymene_o        )        /  imgn_omtp     /
604   DATA    mg20_map ( is_phellandrene_a  )        /  imgn_omtp     /
605   DATA    mg20_map ( is_thujene_a       )        /  imgn_omtp     /
606   DATA    mg20_map ( is_terpinene_a     )        /  imgn_omtp     /
607   DATA    mg20_map ( is_terpinene_g     )        /  imgn_omtp     /
608   DATA    mg20_map ( is_terpinolene     )        /  imgn_omtp     /
609   DATA    mg20_map ( is_phellandrene_b  )        /  imgn_omtp     /
610   DATA    mg20_map ( is_camphene        )        /  imgn_omtp     /
611   DATA    mg20_map ( is_bornene         )        /  imgn_omtp     /
612   DATA    mg20_map ( is_fenchene_a      )        /  imgn_omtp     /
613   DATA    mg20_map ( is_ocimene_al      )        /  imgn_omtp     /
614   DATA    mg20_map ( is_ocimene_c_b     )        /  imgn_omtp     /
615   DATA    mg20_map ( is_tricyclene      )        /  imgn_omtp     /
616   DATA    mg20_map ( is_estragole       )        /  imgn_omtp     /
617   DATA    mg20_map ( is_camphor         )        /  imgn_omtp     /
618   DATA    mg20_map ( is_fenchone        )        /  imgn_omtp     /
619   DATA    mg20_map ( is_piperitone      )        /  imgn_omtp     /
620   DATA    mg20_map ( is_thujone_a       )        /  imgn_omtp     /
621   DATA    mg20_map ( is_thujone_b       )        /  imgn_omtp     /
622   DATA    mg20_map ( is_cineole_1_8     )        /  imgn_omtp     /
623   DATA    mg20_map ( is_borneol         )        /  imgn_omtp     /
624   DATA    mg20_map ( is_linalool        )        /  imgn_omtp     /
625   DATA    mg20_map ( is_terpineol_4     )        /  imgn_omtp     /
626   DATA    mg20_map ( is_terpineol_a     )        /  imgn_omtp     /
627   DATA    mg20_map ( is_linalool_oxd_c  )        /  imgn_omtp     /
628   DATA    mg20_map ( is_linalool_oxd_t  )        /  imgn_omtp     /
629   DATA    mg20_map ( is_ionone_b        )        /  imgn_omtp     /
630   DATA    mg20_map ( is_bornyl_act      )        /  imgn_omtp     /
631   DATA    mg20_map ( is_farnescene_a    )        /  imgn_afarn    /
632   DATA    mg20_map ( is_caryophyllene_b )        /  imgn_bcar     /
633   DATA    mg20_map ( is_acoradiene      )        /  imgn_osqt     /
634   DATA    mg20_map ( is_aromadendrene   )        /  imgn_osqt     /
635   DATA    mg20_map ( is_bergamotene_a   )        /  imgn_osqt     /
636   DATA    mg20_map ( is_bergamotene_b   )        /  imgn_osqt     /
637   DATA    mg20_map ( is_bisabolene_a    )        /  imgn_osqt     /
638   DATA    mg20_map ( is_bisabolene_b    )        /  imgn_osqt     /
639   DATA    mg20_map ( is_bourbonene_b    )        /  imgn_osqt     /
640   DATA    mg20_map ( is_cadinene_d      )        /  imgn_osqt     /
641   DATA    mg20_map ( is_cadinene_g      )        /  imgn_osqt     /
642   DATA    mg20_map ( is_cedrene_a       )        /  imgn_osqt     /
643   DATA    mg20_map ( is_copaene_a       )        /  imgn_osqt     /
644   DATA    mg20_map ( is_cubebene_a      )        /  imgn_osqt     /
645   DATA    mg20_map ( is_cubebene_b      )        /  imgn_osqt     /
646   DATA    mg20_map ( is_elemene_b       )        /  imgn_osqt     /
647   DATA    mg20_map ( is_farnescene_b    )        /  imgn_osqt     /
648   DATA    mg20_map ( is_germacrene_b    )        /  imgn_osqt     /
649   DATA    mg20_map ( is_germacrene_d    )        /  imgn_osqt     /
650   DATA    mg20_map ( is_gurjunene_b     )        /  imgn_osqt     /
651   DATA    mg20_map ( is_humulene_a      )        /  imgn_osqt     /
652   DATA    mg20_map ( is_humulene_g      )        /  imgn_osqt     /
653   DATA    mg20_map ( is_isolongifolene  )        /  imgn_osqt     /
654   DATA    mg20_map ( is_longifolene     )        /  imgn_osqt     /
655   DATA    mg20_map ( is_longipinene     )        /  imgn_osqt     /
656   DATA    mg20_map ( is_muurolene_a     )        /  imgn_osqt     /
657   DATA    mg20_map ( is_muurolene_g     )        /  imgn_osqt     /
658   DATA    mg20_map ( is_selinene_b      )        /  imgn_osqt     /
659   DATA    mg20_map ( is_selinene_d      )        /  imgn_osqt     /
660   DATA    mg20_map ( is_nerolidol_c     )        /  imgn_osqt     /
661   DATA    mg20_map ( is_nerolidol_t     )        /  imgn_osqt     /
662   DATA    mg20_map ( is_cedrol          )        /  imgn_osqt     /
663   DATA    mg20_map ( is_mbo_2m3e2ol     )        /  imgn_mbo      /
664   DATA    mg20_map ( is_methanol        )        /  imgn_meoh     /
665   DATA    mg20_map ( is_acetone         )        /  imgn_acto     /
666   DATA    mg20_map ( is_methane         )        /  imgn_ch4      /
667   DATA    mg20_map ( is_ammonia         )        /  imgn_no       /
668   DATA    mg20_map ( is_nitrous_oxd     )        /  imgn_no       /
669   DATA    mg20_map ( is_nitric_oxd      )        /  imgn_no       /
670   DATA    mg20_map ( is_acetaldehyde    )        /  imgn_acta     /
671   DATA    mg20_map ( is_ethanol         )        /  imgn_acta     /
672   DATA    mg20_map ( is_formic_acid     )        /  imgn_form     /
673   DATA    mg20_map ( is_formaldehyde    )        /  imgn_form     /
674   DATA    mg20_map ( is_acetic_acid     )        /  imgn_form     /
675   DATA    mg20_map ( is_mbo_3m2e1ol     )        /  imgn_co       /
676   DATA    mg20_map ( is_mbo_3m3e1ol     )        /  imgn_co       /
677   DATA    mg20_map ( is_benzaldehyde    )        /  imgn_co       /
678   DATA    mg20_map ( is_butanone_2      )        /  imgn_co       /
679   DATA    mg20_map ( is_decanal         )        /  imgn_co       /
680   DATA    mg20_map ( is_dodecene_1      )        /  imgn_co       /
681   DATA    mg20_map ( is_geranyl_acetone )        /  imgn_co       /
682   DATA    mg20_map ( is_heptanal        )        /  imgn_co       /
683   DATA    mg20_map ( is_heptane         )        /  imgn_co       /
684   DATA    mg20_map ( is_hexane          )        /  imgn_co       /
685   DATA    mg20_map ( is_met_benzoate    )        /  imgn_co       /
686   DATA    mg20_map ( is_met_heptenone   )        /  imgn_co       /
687   DATA    mg20_map ( is_neryl_acetone   )        /  imgn_co       /
688   DATA    mg20_map ( is_nonanal         )        /  imgn_co       /
689   DATA    mg20_map ( is_nonenal         )        /  imgn_co       /
690   DATA    mg20_map ( is_octanal         )        /  imgn_co       /
691   DATA    mg20_map ( is_octanol         )        /  imgn_co       /
692   DATA    mg20_map ( is_octenol_1e3ol   )        /  imgn_co       /
693   DATA    mg20_map ( is_oxopentanal     )        /  imgn_co       /
694   DATA    mg20_map ( is_pentane         )        /  imgn_co       /
695   DATA    mg20_map ( is_phenyl_cco      )        /  imgn_co       /
696   DATA    mg20_map ( is_pyruvic_acid    )        /  imgn_co       /
697   DATA    mg20_map ( is_terpinyl_act_a  )        /  imgn_co       /
698   DATA    mg20_map ( is_tetradecene_1   )        /  imgn_co       /
699   DATA    mg20_map ( is_toluene         )        /  imgn_co       /
700   DATA    mg20_map ( is_carbon_monoxide )        /  imgn_co       /
701   DATA    mg20_map ( is_butene          )        /  imgn_co       /
702   DATA    mg20_map ( is_ethane          )        /  imgn_co       /
703   DATA    mg20_map ( is_ethene          )        /  imgn_co       /
704   DATA    mg20_map ( is_hydrogen_cyanide)        /  imgn_co       /
705   DATA    mg20_map ( is_propane         )        /  imgn_co       /
706   DATA    mg20_map ( is_propene         )        /  imgn_co       /
707   DATA    mg20_map ( is_carbon_2s       )        /  imgn_co       /
708   DATA    mg20_map ( is_carbonyl_s      )        /  imgn_co       /
709   DATA    mg20_map ( is_diallyl_2s      )        /  imgn_co       /
710   DATA    mg20_map ( is_2met_2s         )        /  imgn_co       /
711   DATA    mg20_map ( is_2met_s          )        /  imgn_co       /
712   DATA    mg20_map ( is_met_chloride    )        /  imgn_co       /
713   DATA    mg20_map ( is_met_bromide     )        /  imgn_co       /
714   DATA    mg20_map ( is_met_iodide      )        /  imgn_co       /
715   DATA    mg20_map ( is_hydrogen_s      )        /  imgn_co       /
716   DATA    mg20_map ( is_met_mercaptan   )        /  imgn_co       /
717   DATA    mg20_map ( is_met_propenyl_2s )        /  imgn_co       /
718   DATA    mg20_map ( is_pppp_2s         )        /  imgn_co       /
719   DATA    mg20_map ( is_2met_nonatriene )        /  imgn_co       /
720   DATA    mg20_map ( is_met_salicylate  )        /  imgn_co       /
721   DATA    mg20_map ( is_indole          )        /  imgn_co       /
722   DATA    mg20_map ( is_jasmone         )        /  imgn_co       /
723   DATA    mg20_map ( is_met_jasmonate   )        /  imgn_co       /
724   DATA    mg20_map ( is_3met_3dctt      )        /  imgn_co       /
725   DATA    mg20_map ( is_hexanal         )        /  imgn_co       /
726   DATA    mg20_map ( is_hexanol_1       )        /  imgn_co       /
727   DATA    mg20_map ( is_hexenal_c3      )        /  imgn_co       /
728   DATA    mg20_map ( is_hexenal_t2      )        /  imgn_co       /
729   DATA    mg20_map ( is_hexenol_c3      )        /  imgn_co       /
730   DATA    mg20_map ( is_hexenyl_act_c3  )        /  imgn_co       /
732   !====================================================================
733   ! The following is a table of for the speciation fractions 
734   ! (dimensionless) of the base emission factors EF of n_mgn_spc=20 
735   ! classes of species into n_spca_spc=138 species.  
736   ! See Table 3 of MEGAN v2.0 Documentation.
737   ! Originally from Tan's EFFS_MGN20T138.EXT (12/02/2006)
738   ! Updated on 08/13/2007 based on Alex Guenther's email.
739   !====================================================================
740   REAL, DIMENSION (n_spca_spc, n_pft) :: EF_frac
742   ! The row order of the data statements below can be changed; however
743   ! the column order is hardwreid to k_bt=1, k_nt=2, k_sb=3, and
744   ! k_hb=4 (declared above).  Do not change the column order
745   ! of the data entry below.
747   ! so many decimal places?
748   !                                                   BT            NT            SB            HB
749   DATA  EF_frac( is_isoprene        , k_bt:k_hb )   / 1.000000    , 1.000000    , 1.000000    , 1.000000     / !not used for isoprene
750   DATA  EF_frac( is_myrcene         , k_bt:k_hb )   / 1.000000    , 1.000000    , 1.000000    , 1.000000     /
751   DATA  EF_frac( is_sabinene        , k_bt:k_hb )   / 1.000000    , 1.000000    , 1.000000    , 1.000000     /
752   DATA  EF_frac( is_limonene        , k_bt:k_hb )   / 1.000000    , 1.000000    , 1.000000    , 1.000000     /
753   DATA  EF_frac( is_carene_3        , k_bt:k_hb )   / 1.000000    , 1.000000    , 1.000000    , 1.000000     /
754   DATA  EF_frac( is_ocimene_t_b     , k_bt:k_hb )   / 1.000000    , 1.000000    , 1.000000    , 1.000000     /
755   DATA  EF_frac( is_pinene_b        , k_bt:k_hb )   / 1.000000    , 1.000000    , 1.000000    , 1.000000     /
756   DATA  EF_frac( is_pinene_a        , k_bt:k_hb )   / 1.000000    , 1.000000    , 1.000000    , 1.000000     /
757   DATA  EF_frac( is_2met_styrene    , k_bt:k_hb )   / 0.011470    , 0.005540    , 0.009240    , 0.010400     /
758   DATA  EF_frac( is_cymene_p        , k_bt:k_hb )   / 0.057340    , 0.055430    , 0.046210    , 0.041580     /
759   DATA  EF_frac( is_cymene_o        , k_bt:k_hb )   / 0.034400    , 0.016630    , 0.027730    , 0.031190     /
760   DATA  EF_frac( is_phellandrene_a  , k_bt:k_hb )   / 0.045870    , 0.055430    , 0.046210    , 0.041580     /
761   DATA  EF_frac( is_thujene_a       , k_bt:k_hb )   / 0.011470    , 0.033260    , 0.036970    , 0.041580     /
762   DATA  EF_frac( is_terpinene_a     , k_bt:k_hb )   / 0.057340    , 0.055430    , 0.046210    , 0.041580     /
763   DATA  EF_frac( is_terpinene_g     , k_bt:k_hb )   / 0.057340    , 0.055430    , 0.046210    , 0.041580     /
764   DATA  EF_frac( is_terpinolene     , k_bt:k_hb )   / 0.057340    , 0.066520    , 0.055450    , 0.062370     /
765   DATA  EF_frac( is_phellandrene_b  , k_bt:k_hb )   / 0.057340    , 0.166300    , 0.092420    , 0.103950     /
766   DATA  EF_frac( is_camphene        , k_bt:k_hb )   / 0.172020    , 0.249450    , 0.184840    , 0.145530     /
767   DATA  EF_frac( is_bornene         , k_bt:k_hb )   / 0.010320    , 0.004990    , 0.008320    , 0.008320     /
768   DATA  EF_frac( is_fenchene_a      , k_bt:k_hb )   / 0.003440    , 0.001660    , 0.002770    , 0.004160     /
769   DATA  EF_frac( is_ocimene_al      , k_bt:k_hb )   / 0.011470    , 0.005540    , 0.009240    , 0.010400     /
770   DATA  EF_frac( is_ocimene_c_b     , k_bt:k_hb )   / 0.045870    , 0.022170    , 0.036970    , 0.041580     /
771   DATA  EF_frac( is_tricyclene      , k_bt:k_hb )   / 0.011470    , 0.005540    , 0.009240    , 0.010400     /
772   DATA  EF_frac( is_estragole       , k_bt:k_hb )   / 0.003440    , 0.001660    , 0.002770    , 0.004160     /
773   DATA  EF_frac( is_camphor         , k_bt:k_hb )   / 0.034400    , 0.033260    , 0.046210    , 0.041580     /
774   DATA  EF_frac( is_fenchone        , k_bt:k_hb )   / 0.011470    , 0.005540    , 0.009240    , 0.010400     /
775   DATA  EF_frac( is_piperitone      , k_bt:k_hb )   / 0.003440    , 0.001660    , 0.002770    , 0.004160     /
776   DATA  EF_frac( is_thujone_a       , k_bt:k_hb )   / 0.011470    , 0.027720    , 0.046210    , 0.041580     /
777   DATA  EF_frac( is_thujone_b       , k_bt:k_hb )   / 0.002290    , 0.005540    , 0.009240    , 0.010400     /
778   DATA  EF_frac( is_cineole_1_8     , k_bt:k_hb )   / 0.057340    , 0.011090    , 0.036970    , 0.041580     /
779   DATA  EF_frac( is_borneol         , k_bt:k_hb )   / 0.008030    , 0.003880    , 0.006470    , 0.006240     /
780   DATA  EF_frac( is_linalool        , k_bt:k_hb )   / 0.137610    , 0.066520    , 0.110910    , 0.124740     /
781   DATA  EF_frac( is_terpineol_4     , k_bt:k_hb )   / 0.006880    , 0.003330    , 0.005550    , 0.006240     /
782   DATA  EF_frac( is_terpineol_a     , k_bt:k_hb )   / 0.034400    , 0.016630    , 0.027730    , 0.031190     /
783   DATA  EF_frac( is_linalool_oxd_c  , k_bt:k_hb )   / 0.006880    , 0.003330    , 0.005550    , 0.006240     /
784   DATA  EF_frac( is_linalool_oxd_t  , k_bt:k_hb )   / 0.034400    , 0.016630    , 0.027730    , 0.031190     /
785   DATA  EF_frac( is_ionone_b        , k_bt:k_hb )   / 0.002290    , 0.001110    , 0.001850    , 0.002080     /
786   DATA  EF_frac( is_bornyl_act      , k_bt:k_hb )   / 0.001150    , 0.002770    , 0.002770    , 0.002080     /
787   DATA  EF_frac( is_farnescene_a    , k_bt:k_hb )   / 1.000000    , 1.000000    , 1.000000    , 1.000000     /
788   DATA  EF_frac( is_caryophyllene_b , k_bt:k_hb )   / 1.000000    , 1.000000    , 1.000000    , 1.000000     /
789   DATA  EF_frac( is_acoradiene      , k_bt:k_hb )   / 0.018570    , 0.015950    , 0.019160    , 0.021860     /
790   DATA  EF_frac( is_aromadendrene   , k_bt:k_hb )   / 0.007430    , 0.006380    , 0.007660    , 0.010930     /
791   DATA  EF_frac( is_bergamotene_a   , k_bt:k_hb )   / 0.083570    , 0.143540    , 0.095790    , 0.098360     /
792   DATA  EF_frac( is_bergamotene_b   , k_bt:k_hb )   / 0.001860    , 0.001590    , 0.001920    , 0.005460     /
793   DATA  EF_frac( is_bisabolene_a    , k_bt:k_hb )   / 0.001860    , 0.001590    , 0.001920    , 0.005460     /
794   DATA  EF_frac( is_bisabolene_b    , k_bt:k_hb )   / 0.055710    , 0.119620    , 0.067050    , 0.076500     /
795   DATA  EF_frac( is_bourbonene_b    , k_bt:k_hb )   / 0.027860    , 0.023920    , 0.028740    , 0.032790     /
796   DATA  EF_frac( is_cadinene_d      , k_bt:k_hb )   / 0.013930    , 0.011960    , 0.014370    , 0.016390     /
797   DATA  EF_frac( is_cadinene_g      , k_bt:k_hb )   / 0.009290    , 0.007970    , 0.009580    , 0.010930     /
798   DATA  EF_frac( is_cedrene_a       , k_bt:k_hb )   / 0.005570    , 0.004780    , 0.005750    , 0.005460     /
799   DATA  EF_frac( is_copaene_a       , k_bt:k_hb )   / 0.009290    , 0.007970    , 0.009580    , 0.010930     /
800   DATA  EF_frac( is_cubebene_a      , k_bt:k_hb )   / 0.013930    , 0.011960    , 0.014370    , 0.016390     /
801   DATA  EF_frac( is_cubebene_b      , k_bt:k_hb )   / 0.009290    , 0.007970    , 0.009580    , 0.010930     /
802   DATA  EF_frac( is_elemene_b       , k_bt:k_hb )   / 0.018570    , 0.015950    , 0.019160    , 0.021860     /
803   DATA  EF_frac( is_farnescene_b    , k_bt:k_hb )   / 0.278550    , 0.239230    , 0.287360    , 0.218580     /
804   DATA  EF_frac( is_germacrene_b    , k_bt:k_hb )   / 0.009290    , 0.007970    , 0.009580    , 0.010930     /
805   DATA  EF_frac( is_germacrene_d    , k_bt:k_hb )   / 0.027860    , 0.023920    , 0.028740    , 0.032790     /
806   DATA  EF_frac( is_gurjunene_b     , k_bt:k_hb )   / 0.004640    , 0.003990    , 0.004790    , 0.005460     /
807   DATA  EF_frac( is_humulene_a      , k_bt:k_hb )   / 0.139280    , 0.199360    , 0.172410    , 0.163930     /
808   DATA  EF_frac( is_humulene_g      , k_bt:k_hb )   / 0.001860    , 0.001590    , 0.001920    , 0.005460     /
809   DATA  EF_frac( is_isolongifolene  , k_bt:k_hb )   / 0.001860    , 0.001590    , 0.001920    , 0.005460     /
810   DATA  EF_frac( is_longifolene     , k_bt:k_hb )   / 0.001860    , 0.001590    , 0.001920    , 0.005460     /
811   DATA  EF_frac( is_longipinene     , k_bt:k_hb )   / 0.001860    , 0.001590    , 0.001920    , 0.005460     /
812   DATA  EF_frac( is_muurolene_a     , k_bt:k_hb )   / 0.013930    , 0.011960    , 0.014370    , 0.016390     /
813   DATA  EF_frac( is_muurolene_g     , k_bt:k_hb )   / 0.046430    , 0.039870    , 0.047890    , 0.054640     /
814   DATA  EF_frac( is_selinene_b      , k_bt:k_hb )   / 0.185700    , 0.079740    , 0.114940    , 0.109290     /
815   DATA  EF_frac( is_selinene_d      , k_bt:k_hb )   / 0.001860    , 0.001590    , 0.001920    , 0.005460     /
816   DATA  EF_frac( is_nerolidol_c     , k_bt:k_hb )   / 0.001860    , 0.001590    , 0.001920    , 0.005460     /
817   DATA  EF_frac( is_nerolidol_t     , k_bt:k_hb )   / 0.004640    , 0.001590    , 0.001920    , 0.005460     /
818   DATA  EF_frac( is_cedrol          , k_bt:k_hb )   / 0.001860    , 0.001590    , 0.001920    , 0.005460     /
819   DATA  EF_frac( is_mbo_2m3e2ol     , k_bt:k_hb )   / 1.000000    , 1.000000    , 1.000000    , 1.000000     /
820   DATA  EF_frac( is_methanol        , k_bt:k_hb )   / 1.000000    , 1.000000    , 1.000000    , 1.000000     /
821   DATA  EF_frac( is_acetone         , k_bt:k_hb )   / 1.000000    , 1.000000    , 1.000000    , 1.000000     /
822   DATA  EF_frac( is_methane         , k_bt:k_hb )   / 1.000000    , 1.000000    , 1.000000    , 1.000000     /
823   DATA  EF_frac( is_ammonia         , k_bt:k_hb )   / 0.000000    , 0.000000    , 0.000000    , 0.000000     /
824   DATA  EF_frac( is_nitrous_oxd     , k_bt:k_hb )   / 0.000000    , 0.000000    , 0.000000    , 0.000000     /
825   DATA  EF_frac( is_nitric_oxd      , k_bt:k_hb )   / 1.000000    , 1.000000    , 1.000000    , 1.000000     /
826   DATA  EF_frac( is_acetaldehyde    , k_bt:k_hb )   / 0.500000    , 0.500000    , 0.500000    , 0.500000     /
827   DATA  EF_frac( is_ethanol         , k_bt:k_hb )   / 0.500000    , 0.500000    , 0.500000    , 0.500000     /
828   DATA  EF_frac( is_formic_acid     , k_bt:k_hb )   / 0.285710    , 0.285710    , 0.285710    , 0.285710     /
829   DATA  EF_frac( is_formaldehyde    , k_bt:k_hb )   / 0.428570    , 0.428570    , 0.428570    , 0.428570     /
830   DATA  EF_frac( is_acetic_acid     , k_bt:k_hb )   / 0.285710    , 0.285710    , 0.285710    , 0.285710     /
831   DATA  EF_frac( is_mbo_3m2e1ol     , k_bt:k_hb )   / 0.000520    , 0.000520    , 0.000520    , 0.000520     /
832   DATA  EF_frac( is_mbo_3m3e1ol     , k_bt:k_hb )   / 0.000520    , 0.000520    , 0.000520    , 0.000520     /
833   DATA  EF_frac( is_benzaldehyde    , k_bt:k_hb )   / 0.000210    , 0.000210    , 0.000210    , 0.000210     /
834   DATA  EF_frac( is_butanone_2      , k_bt:k_hb )   / 0.001030    , 0.001030    , 0.001030    , 0.001030     /
835   DATA  EF_frac( is_decanal         , k_bt:k_hb )   / 0.000210    , 0.000210    , 0.000210    , 0.000210     /
836   DATA  EF_frac( is_dodecene_1      , k_bt:k_hb )   / 0.000210    , 0.000210    , 0.000210    , 0.000210     /
837   DATA  EF_frac( is_geranyl_acetone , k_bt:k_hb )   / 0.003100    , 0.003100    , 0.003100    , 0.003100     /
838   DATA  EF_frac( is_heptanal        , k_bt:k_hb )   / 0.000210    , 0.000210    , 0.000210    , 0.000210     /
839   DATA  EF_frac( is_heptane         , k_bt:k_hb )   / 0.000210    , 0.000210    , 0.000210    , 0.000210     /
840   DATA  EF_frac( is_hexane          , k_bt:k_hb )   / 0.000100    , 0.000100    , 0.000100    , 0.000100     /
841   DATA  EF_frac( is_met_benzoate    , k_bt:k_hb )   / 0.000100    , 0.000100    , 0.000100    , 0.000100     /
842   DATA  EF_frac( is_met_heptenone   , k_bt:k_hb )   / 0.000210    , 0.000210    , 0.000210    , 0.000210     /
843   DATA  EF_frac( is_neryl_acetone   , k_bt:k_hb )   / 0.000100    , 0.000100    , 0.000100    , 0.000100     /
844   DATA  EF_frac( is_nonanal         , k_bt:k_hb )   / 0.000210    , 0.000210    , 0.000210    , 0.000210     /
845   DATA  EF_frac( is_nonenal         , k_bt:k_hb )   / 0.000100    , 0.000100    , 0.000100    , 0.000100     /
846   DATA  EF_frac( is_octanal         , k_bt:k_hb )   / 0.000210    , 0.000210    , 0.000210    , 0.000210     /
847   DATA  EF_frac( is_octanol         , k_bt:k_hb )   / 0.000210    , 0.000210    , 0.000210    , 0.000210     /
848   DATA  EF_frac( is_octenol_1e3ol   , k_bt:k_hb )   / 0.000210    , 0.000210    , 0.000210    , 0.000210     /
849   DATA  EF_frac( is_oxopentanal     , k_bt:k_hb )   / 0.000210    , 0.000210    , 0.000210    , 0.000210     /
850   DATA  EF_frac( is_pentane         , k_bt:k_hb )   / 0.000100    , 0.000100    , 0.000100    , 0.000100     /
851   DATA  EF_frac( is_phenyl_cco      , k_bt:k_hb )   / 0.000100    , 0.000100    , 0.000100    , 0.000100     /
852   DATA  EF_frac( is_pyruvic_acid    , k_bt:k_hb )   / 0.002060    , 0.002060    , 0.002060    , 0.002060     /
853   DATA  EF_frac( is_terpinyl_act_a  , k_bt:k_hb )   / 0.000100    , 0.000100    , 0.000100    , 0.000100     /
854   DATA  EF_frac( is_tetradecene_1   , k_bt:k_hb )   / 0.000100    , 0.000100    , 0.000100    , 0.000100     /
855   DATA  EF_frac( is_toluene         , k_bt:k_hb )   / 0.002060    , 0.002060    , 0.002060    , 0.002060     /
856   DATA  EF_frac( is_carbon_monoxide , k_bt:k_hb )   / 0.619070    , 0.619070    , 0.619070    , 0.619070     /
857   DATA  EF_frac( is_butene          , k_bt:k_hb )   / 0.036110    , 0.036110    , 0.036110    , 0.036110     /
858   DATA  EF_frac( is_ethane          , k_bt:k_hb )   / 0.002060    , 0.002060    , 0.002060    , 0.002060     /
859   DATA  EF_frac( is_ethene          , k_bt:k_hb )   / 0.134130    , 0.134130    , 0.134130    , 0.134130     /
860   DATA  EF_frac( is_hydrogen_cyanide, k_bt:k_hb )   / 0.004130    , 0.004130    , 0.004130    , 0.004130     /
861   DATA  EF_frac( is_propane         , k_bt:k_hb )   / 0.001030    , 0.001030    , 0.001030    , 0.001030     /
862   DATA  EF_frac( is_propene         , k_bt:k_hb )   / 0.082540    , 0.082540    , 0.082540    , 0.082540     /
863   DATA  EF_frac( is_carbon_2s       , k_bt:k_hb )   / 0.000310    , 0.000310    , 0.000310    , 0.000310     /
864   DATA  EF_frac( is_carbonyl_s      , k_bt:k_hb )   / 0.000620    , 0.000620    , 0.000620    , 0.000620     /
865   DATA  EF_frac( is_diallyl_2s      , k_bt:k_hb )   / 0.000100    , 0.000100    , 0.000100    , 0.000100     /
866   DATA  EF_frac( is_2met_2s         , k_bt:k_hb )   / 0.000310    , 0.000310    , 0.000310    , 0.000310     /
867   DATA  EF_frac( is_2met_s          , k_bt:k_hb )   / 0.001240    , 0.001240    , 0.001240    , 0.001240     /
868   DATA  EF_frac( is_met_chloride    , k_bt:k_hb )   / 0.001030    , 0.001030    , 0.001030    , 0.001030     /
869   DATA  EF_frac( is_met_bromide     , k_bt:k_hb )   / 0.000210    , 0.000210    , 0.000210    , 0.000210     /
870   DATA  EF_frac( is_met_iodide      , k_bt:k_hb )   / 0.000100    , 0.000100    , 0.000100    , 0.000100     /
871   DATA  EF_frac( is_hydrogen_s      , k_bt:k_hb )   / 0.000520    , 0.000520    , 0.000520    , 0.000520     /
872   DATA  EF_frac( is_met_mercaptan   , k_bt:k_hb )   / 0.000100    , 0.000100    , 0.000100    , 0.000100     /
873   DATA  EF_frac( is_met_propenyl_2s , k_bt:k_hb )   / 0.000100    , 0.000100    , 0.000100    , 0.000100     /
874   DATA  EF_frac( is_pppp_2s         , k_bt:k_hb )   / 0.000100    , 0.000100    , 0.000100    , 0.000100     /
875   DATA  EF_frac( is_2met_nonatriene , k_bt:k_hb )   / 0.020640    , 0.020640    , 0.020640    , 0.020640     /
876   DATA  EF_frac( is_met_salicylate  , k_bt:k_hb )   / 0.002060    , 0.002060    , 0.002060    , 0.002060     /
877   DATA  EF_frac( is_indole          , k_bt:k_hb )   / 0.000210    , 0.000210    , 0.000210    , 0.000210     /
878   DATA  EF_frac( is_jasmone         , k_bt:k_hb )   / 0.000210    , 0.000210    , 0.000210    , 0.000210     /
879   DATA  EF_frac( is_met_jasmonate   , k_bt:k_hb )   / 0.002060    , 0.002060    , 0.002060    , 0.002060     /
880   DATA  EF_frac( is_3met_3dctt      , k_bt:k_hb )   / 0.000210    , 0.000210    , 0.000210    , 0.000210     /
881   DATA  EF_frac( is_hexanal         , k_bt:k_hb )   / 0.003100    , 0.003100    , 0.003100    , 0.003100     /
882   DATA  EF_frac( is_hexanol_1       , k_bt:k_hb )   / 0.003100    , 0.003100    , 0.003100    , 0.003100     /
883   DATA  EF_frac( is_hexenal_c3      , k_bt:k_hb )   / 0.025790    , 0.025790    , 0.025790    , 0.025790     /
884   DATA  EF_frac( is_hexenal_t2      , k_bt:k_hb )   / 0.015480    , 0.015480    , 0.015480    , 0.015480     /
885   DATA  EF_frac( is_hexenol_c3      , k_bt:k_hb )   / 0.015480    , 0.015480    , 0.015480    , 0.015480     /
886   DATA  EF_frac( is_hexenyl_act_c3  , k_bt:k_hb )   / 0.015480    , 0.015480    , 0.015480    , 0.015480     / 
889 END MODULE module_data_megan2