Update mozbc mapping for mozart-mosaic 4bin

[Naser Mohammadzadeh]

Hello,

I want to use use mozbc for chem_opt=201.

I found this mapping:

spc_map = 'o3 -> O3', 'no -> NO',

'no2 -> NO2', 'no3 -> NO3', 'nh3 -> NH3', 'hno3 -> HNO3', 'hno4 -> HO2NO2',

'n2o5 -> N2O5', 'ho -> OH', 'ho2 -> HO2', 'h2o2 -> H2O2', 'ch4 -> CH4', 'co -> CO', 'ch3o2 -> CH3O2', 'ch3ooh -> CH3OOH', 'hcho -> CH2O', 'ch3oh -> CH3OH', 'c2h4 -> C2H4',

'ald -> CH3CHO', 'ch3cooh -> CH3COOH', 'acet -> CH3COCH3', 'mgly -> CH3COCHO', 'paa -> CH3COOOH',

'pan -> PAN', 'mpan -> MPAN', 'macr -> MACR',

'mvk -> MVK', 'c2h6 -> C2H6', 'c3h6 -> C3H6', 'c3h8 -> C3H8', 'c2h5oh -> C2H5OH',

'onit -> ONIT', 'onitr -> ONITR', 'isopr -> ISOP',

'isopn -> ISOPNO3', 'acetol -> HYAC', 'glyald -> GLYALD', 'hydrald -> HYDRALD', 'mek -> MEK',

'bigene -> BIGENE', 'open -> BIGALD', 'bigalk -> BIGALK', 'tol -> TOLUENE',

'cres -> CRESOL', 'dms -> DMS', 'so2 -> SO2', 'sulf -> SO4',

'oc_a01->9.60024e-2*OC1+9.60024e-2*OC2+1.152026*SOA;1.e9' 'oc_a02->.254356*OC1+.254356*OC2+3.05226*SOA;1.e9', 'oc_a03->6.00723e-2*OC1+6.00723e-2*OC2+.720865*SOA;1.e9', 'oc_a04->.00104*OC1+.00104*OC2+.012481*SOA;1.e9', 'bc_a01->6.5325E-01*CB1+6.5325E-01*CB2;0.413793103e9', 'bc_a02->2.5458E-01*CB1+2.5458E-01*CB2;0.413793103e9', 'bc_a03->4.0521E-03*CB1+4.0521E-03*CB2;0.413793103e9', 'bc_a04->1.6794E-06*CB1+1.6794E-06*CB2;0.413793103e9',

'so4_a01->2.4433e-2*SO4;3.31423e9', 'so4_a02->.47111*SO4;3.31423e9', 'so4_a03->.47857*SO4;3.31423e9', 'so4_a04->2.5727e-2*SO4;3.31423e9', 'no3_a01->2.4433e-2*NH4NO3;2.76186e9', 'no3_a02->.47111*NH4NO3;2.76186e9', 'no3_a03->.47857*NH4NO3;2.76186e9', 'no3_a04->2.5727e-2*NH4NO3;2.76186e9', 'nh4_a01->2.4433e-2*NH4;6.22876e8', 'nh4_a02->.47111*NH4;6.22876e8', 'nh4_a03->.47857*NH4;6.22876e8', 'nh4_a04->2.5727e-2*NH4;6.22876e8', 'na_a01->.2*SA1;1.e9',

'na_a02->.4*SA1;1.e9', 'na_a03->.2*SA1+.8*SA2;1.e9', 'na_a04->.8*SA3;1.e9', 'cl_a01->.3*SA1;1.e9', 'cl_a02->.6*SA1;1.e9', 'cl_a03->.3*SA1+1.2*SA2;1.e9', 'cl_a04->1.2*SA3;1.e9', 'oin_a01->0.5*[DUST1];1.e9',

'oin_a02->0.5*[DUST1]+0.24*[DUST2];1.e9', 'oin_a03->0.76*[DUST2]+1.0*[DUST3];1.e9', 'oin_a04->1.0*[DUST4];1.e9',

'num_a01- >1.73e17*OC1+1.73e17*OC2+1.73e17*SOA+5.64e17*CB1+5.64e17*CB2+7.67e16*SO4+6.39 e16*NH4NO3+1.44e16*NH4+3.22e17*SA1+4.83e17*SA1+0.5*3.93e17*[DUST1];1.0',

'num_a02- >1.71e16*OC1+1.71e16*OC2+1.71e16*SOA+9.91e15*CB1+9.91e15*CB2+6.06e16*SO4+5.05 e16*NH4NO3+1.14e16*NH4+2.90e16*SA1+4.36e16*SA1+0.5*1.17e16*[DUST1]+0.24*1.17e16 *[DUST2];1.0',

'num_a03- >1.13e14*OC1+1.13e14*OC2+1.13e14*SOA+1.80e12*CB1+1.80e12*CB2+2.44e15*SO4+2.03 e15*NH4NO3+4.69e14*NH4+4.12e14*SA1+1.64e15*SA2+0.76*9.55e13*[DUST2]+9.55e13*[D UST3];1.0',

'num_a04- >4.09e10*OC1+4.09e10*OC2+4.09e10*SOA+3.42e07*CB1+3.42e07*CB2+3.25e12*SO4+2.71 e12*NH4NO3+6.12e11*NH4+7.60e13*SA3+1.49e12*[DUST4];1.0',

But I faced an error to not exist specific variables. 

Do you have any updated version of this schemes?

Thank you 

Naser

[Héctor Jorquera González]

Hi Naser

Try this mapping, it worked for my run of mech 201:

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%   beginning  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

spc_map   = 'o3 -> O3', 'n2o -> N2O', 'no -> NO',
            'no2 -> NO2', 'nh3 -> NH3', 'hno3 -> HNO3', 'hno4 -> HO2NO2',
            'n2o5 -> N2O5', 'h2o2 -> H2O2',
            'ch4 -> CH4', 'co -> CO', 'ch3ooh -> CH3OOH',

            'hcho -> CH2O', 'ch3oh -> CH3OH', 'c2h4 -> C2H4',

            'ald -> CH3CHO', 'acet -> CH3COCH3', 'mgly -> CH3COCHO',

            'pan -> PAN', 'mpan -> MPAN', 'macr -> MACR',
            'mvk -> MVK', 'c2h6 -> C2H6', 'c3h6 -> C3H6', 'c3h8 -> C3H8',
            'c2h5oh -> C2H5OH',

            'isopr -> ISOP','acetol -> HYAC', 'mek -> MEK',
            'bigene -> BIGENE', 'bigalk -> BIGALK',
            'tol -> TOLUENE', 'benzene -> BENZENE',

            'cres -> CRESOL', 'dms -> DMS', 'so2 -> SO2',

           'asoa1_a01=0.0',
           'asoa1_a02=0.0',
           'asoa1_a03=0.0',
           'asoa1_a04=0.0',
           'asoa2_a01=0.0',
           'asoa2_a02=0.0',
           'asoa2_a03=0.0',
           'asoa2_a04=0.0',
           'asoa3_a01=0.0',
           'asoa3_a02=0.0',
           'asoa3_a03=0.0',
           'asoa3_a04=0.0',
           'asoa4_a01=0.0',
           'asoa4_a02=0.0',
           'asoa4_a03=0.0',
           'asoa4_a04=0.0',
           'asoaX_a01=0.0',
           'asoaX_a02=0.0',
           'asoaX_a03=0.0',
           'asoaX_a04=0.0',
           'bsoa1_a01=0.0',
           'bsoa1_a02=0.0',
           'bsoa1_a03=0.0',
           'bsoa1_a04=0.0',
           'bsoa2_a01=0.0',
           'bsoa2_a02=0.0',
           'bsoa2_a03=0.0',
           'bsoa2_a04=0.0',
           'bsoa3_a01=0.0',
           'bsoa3_a02=0.0',
           'bsoa3_a03=0.0',
           'bsoa3_a04=0.0',
           'bsoa4_a01=0.0',
           'bsoa4_a02=0.0',
           'bsoa4_a03=0.0',
           'bsoa4_a04=0.0',
           'bsoaX_a01=0.0',
           'bsoaX_a02=0.0',
           'bsoaX_a03=0.0',
           'bsoaX_a04=0.0',
           'cvasoa1=0.0',
           'cvasoa2=0.0',
           'cvasoa3=0.0',
           'cvasoa4=0.0',
           'cvasoaX=0.0',
           'cvbsoa1=0.0',
           'cvbsoa2=0.0',
           'cvbsoa3=0.0',
           'cvbsoa4=0.0',
           'cvbsoaX=0.0',
          'sulf -> so4_a1+so4_a2+so4_a3;0.30173e6',
          'oc_a01->0.1216*pom_a1+0.9886*soa1_a2+0.1216*soa1_a1+0.9886*soa2_a2+
           0.1216*soa2_a1+0.9886*soa3_a2+0.1216*soa3_a1+0.9886*soa4_a2+0.1216*soa4_a1+
           0.9886*soa5_a2+0.1216*soa5_a1;1.e9          ',
          'oc_a02->0.7618*pom_a1+0.0114*soa1_a2+0.7618*soa1_a1+0.0134*soa2_a2+
           0.7618*soa2_a1+0.0114*soa3_a2+0.7621*soa3_a1+0.0114*soa4_a2+0.7621*soa4_a1+
           0.0114*soa5_a2+ 0.7621*soa5_a1;1.e9          ',
          'oc_a03->0.1164*pom_a1+0.0000*soa1_a2+0.1164*soa1_a1+0.0000*soa2_a2+
           0.1164*soa2_a1+0.0000*soa3_a2+0.1164*soa3_a1+0.0000*soa4_a2+0.1164*soa4_a1+
           0.0000*soa5_a2+ 0.1164*soa5_a1;1.e9          ',
          'oc_a04->0.0002*pom_a1+0.0000*soa1_a2+0.0002*soa1_a1+0.0000*soa2_a2+
           0.0002*soa2_a1+0.0000*soa3_a2+0.0002*soa3_a1+0.0000*soa4_a2+0.0002*soa4_a1+
           0.0000*soa5_a2+ 0.0002*soa5_a1;1.e9          ',
          'bc_a01->0.1216*bc_a1+0.1216*bc_a4;1.e9',
          'bc_a02->0.7618*bc_a1+0.7618*bc_a4;1.e9',
          'bc_a03->0.1164*bc_a1+0.1164*bc_a4;1.e9',
          'bc_a04->0.0002*bc_a1+0.0002*bc_a4;1.e9',
          'so4_a01->0.9886*so4_a2+0.1216*so4_a1+0.0000*so4_a3;1.e9',
          'so4_a02->0.0114*so4_a2+0.7618*so4_a1+0.0002*so4_a3;1.e9',
          'so4_a03->0.0000*so4_a2+0.1164*so4_a1+0.0995*so4_a3;1.e9',
          'so4_a04->0.0000*so4_a2+0.0002*so4_a1+0.9003*so4_a3;1.e9',
          'nh4_a01->0.1856*so4_a2+0.0050*so4_a1+0.0000*so4_a3;1.e9',
          'nh4_a02->0.0021*so4_a2+0.0930*so4_a1+0.0000*so4_a3;1.e9',
          'nh4_a03->0.0000*so4_a2+0.0203*so4_a1+0.0186*so4_a3;1.e9',
          'nh4_a04->0.0000*so4_a2+0.0000*so4_a1+0.1690*so4_a3;1.e9',
          'no3_a01->0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3;1.e9',
          'no3_a02->0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3;1.e9',
          'no3_a03->0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3;1.e9',
          'no3_a04->0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3;1.e9',
          'na_a01->0.3889*ncl_a2+0.0479*ncl_a1+0.0000*ncl_a3;1.e9',
          'na_a02->0.0045*ncl_a2+0.2997*ncl_a1+0.0000*ncl_a3;1.e9',
          'na_a03->0.0000*ncl_a2+0.0458*ncl_a1+0.0391*ncl_a3;1.e9',
          'na_a04->0.0000*ncl_a2+0.0000*ncl_a1+0.3542*ncl_a3;1.e9',
          'cl_a01->0.5996*ncl_a2+0.0737*ncl_a1+0.0000*ncl_a3;1.e9',
          'cl_a02->0.0068*ncl_a2+0.4621*ncl_a1+0.0000*ncl_a3;1.e9',
          'cl_a03->0.0000*ncl_a2+0.0709*ncl_a1+0.0604*ncl_a3;1.e9',
          'cl_a04->0.0000*ncl_a2+0.0001*ncl_a1+0.5462*ncl_a3;1.e9',
          'oin_a01->0.9886*dst_a2+0.1216*dst_a1+0.0000*dst_a3;1.e9',
          'oin_a02->0.0114*dst_a2+0.7618*dst_a1+0.0002*dst_a3;1.e9',
          'oin_a03->0.0000*dst_a2+0.1164*dst_a1+0.0995*dst_a3;1.e9',
          'oin_a04->0.0000*dst_a2+0.0002*dst_a1+0.9003*dst_a3;1.e9',
          'num_a01->0.9996*num_a2+0.7135*num_a1+0.0000*num_a3;1.0',
          'num_a02->0.0004*num_a2+0.2847*num_a1+0.0239*num_a3;1.0',
          'num_a03->0.0000*num_a2+0.0016*num_a1+0.6258*num_a3;1.0',
          'num_a04->0.0000*num_a2+0.0377*num_a1+0.3501*num_a3;1.0',
/

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%     end   %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

I did not have a clear idea about several secondary organic aerosols and vapor precursors, so I zeroed them throughout

I am using camchem output files to do the processing of bdry and init files.

All the best

Hector