BEGSIM output unit = 7 file = inca.doc version=CH4_AER COMMENTS LMDz4 INCA version CH4_AER Didier Hauglustaine 01/2003 Michael Schulz, Susanne Bauer, Yves Balkanski, Christiane Textor, Sarah Guibert 07/2003 Anne Cozic 03/2005 ENDPAR SPECIES Solution Rn222 Pb210a -> Pb210 Pb210b -> Pb210 Be7 CO2FF -> CO2 Kr85 SF6 CFC11 -> CCl3F CFC12 -> CCl2F2 CFC113 -> CF2ClCFCl2 HCFC22 -> CHClF2 HFC134a -> CH2FCF3 HCFC141 -> CH3CCl2F HCFC142 -> CH3CClF2 Ha1301 -> CF3Br Ha1211 -> CF2BrCl MCF -> CH3CCl3 CCl4, CH3Cl, CH3Br CH4, N2O O3I -> O3 O3S -> O3 O3, O, O1D -> O H, OH, HO2, H2O2 N, NO, NO2, NO3, HNO3, HNO4, N2O5 CO CH3O2, CH3OOH, CH2O H2 CH3O, CH3OH SO2 H2S DMS ->CH3SCH3 DMSO ->CH3SOCH3 NH3 CO2OC -> CO2 CO2BIM -> CO2 CO2BIH -> CO2 * * nous n'utilisons pas de masse molaire pour cidustm - le travail est fait tout le temps * avec mmr - pour eviter que le preproc se plante nous lui indiquons une masse molaire * de 0. - idem pour CIN et cie. CIDUSTM ->@0. CINO3M ->NO3 CIN ->@0. * * AIBCM is for the hydrophobic part of BC, resp. AIPOMM hydrophobic part of POM * AIBCM -> C AIPOMM -> C AIN ->@0. * CSSSM ->NaCl CSSO4M ->SO4 CSNH4M ->NH4 CSMSAM ->CH3SO3H CSNO3M ->NO3 CSN ->@0. * ASSSM ->NaCl ASBCM ->C ASPOMM ->C ASSO4M ->SO4 ASNH4M ->NH4 ASMSAM ->CH3SO3H ASNO3M ->NO3 ASN ->@0. * SSSSM ->NaCl SSN ->@0. Endlst Fixed M, N2, O2, H2O, CO2 Endlst Col-int O3 = 0. O2 = 0. Endent Groups OX = O3 + O1D + O Endent ENDPAR Solution classes Explicit Rn222, Pb210a, Pb210b, Be7, CO2FF, Kr85 SF6 CFC11, CFC12, CFC113 HCFC22, HFC134a, HCFC141, HCFC142 Ha1301, Ha1211 MCF, CCl4, CH3Cl, CH3Br CH4, N2O CO, H2 CO2OC,CO2BIM,CO2BIH * CIDUSTM CINO3M CIN * AIBCM AIPOMM AIN * CSSSM CSSO4M CSNH4M CSMSAM CSNO3M CSN * ASSSM ASNH4M ASMSAM ASNO3M ASBCM ASPOMM ASN * SSSSM SSN *WARNING O3I and O3S always last explicit species unless * exp_slv is modified O3I, O3S Endlst QSSA Endlst Rodas Endlst EBI Endlst Implicit OX H, OH, HO2, H2O2 N, NO, NO2, NO3, HNO3, HNO4, N2O5 CH3O2, CH3OOH, CH2O CH3O, CH3OH SO2,H2S,DMS,DMSO,ASSO4M Endlst ENDPAR CHEMISTRY Photolysis *WARNING Ordering imposed by the input photorate table [jO2] O2 + hv -> 2 * O [jO3O1D] O3 + hv -> O1D + O2 [jO3O] O3 + hv -> O + O2 [jNO2] NO2 + hv -> NO + O [jNO3NO] NO3 + hv -> NO + O2 [jNO3NO2] NO3 + hv -> NO2 + O [jN2O5] N2O5 + hv -> NO2 + NO3 [jHNO3] HNO3 + hv -> NO2 + OH [jHNO4] HNO4 + hv -> NO2 + HO2 [jCH3OOH] CH3OOH + hv -> CH2O + H + OH [jCH2OCO2H] CH2O + hv -> CO + 2 * H [jCH2OCOH2] CH2O + hv -> CO + H2 [jH2O2] H2O2 + hv -> 2 * OH [jN2O] N2O + hv -> O1D + N2 [jCFC11] CFC11 + hv -> [jCFC12] CFC12 + hv -> [jCFC113] CFC113 + hv -> [jHCFC22] HCFC22 + hv -> [jHCFC141] HCFC141 + hv -> [jHCFC142] HCFC142 + hv -> [jHa1301] Ha1301 + hv -> [jHa1211] Ha1211 + hv -> [jMCF] MCF + hv -> [jCCl4] CCl4 + hv -> [jCH3Cl] CH3Cl + hv -> [jCH3Br] CH3Br + hv -> [jH2O] H2O + hv -> OH + H *WARNING NO is always the last photorate [jNO] NO + hv -> N + O Endent Reactions *Gas-phase reactions [kOO2M] O + O2 + M -> O3 + M [kOOM] O + O + M -> O2 + M [kNO2NO3M] NO2 + NO3 + M -> N2O5 + M [kN2O5M] N2O5 + M -> NO2 + NO3 + M [kNO2OHM] NO2 + OH + M -> HNO3 + M [kNO2HO2M] NO2 + HO2 + M -> HNO4 + M [kHNO4M] HNO4 + M -> HO2 + NO2 + M [kHNO3OH] HNO3 + OH -> NO3 + H2O [kCOOH] CO + OH -> CO2 + H [kHO2HO2] HO2 + HO2 -> H2O2 [kHO2M] H + O2 + M -> HO2 + M * === inorganic species === [kSO2OHM] SO2 + OH + M -> ASSO4M + HO2 + M [kDMSOHDMSO] DMS + OH -> 0.6*SO2 + 0.4*DMSO + 1.2*CH2O *Aerosol reactions [kN2O5] N2O5 -> 2 * HNO3 [kNO2] NO2 -> 0.5*HNO3+0.5*HNO2 [kNO3] NO3 -> HNO3 [kHO2] HO2 -> 0.5*H2O2+0.5*O2 *Aqueous phase reactions *(are in chimiaq and must be hardcoded) *Arrhenius type reactions reaction; A, -Ea [kRn222] Rn222 -> Pb210a + Pb210b ; 2.100e-06 [kBe7] Be7 -> ; 1.514e-07 [kKr85] Kr85 -> ; 2.043e-09 [kSF6O1D] SF6 + O1D -> ; 1.800e-14 [kMCFOH] MCF + OH -> H2O ; 1.800e-12, -1550 [kCH3ClOH] CH3Cl + OH -> H2O + HO2; 4.000e-12, -1400 [kCH3BrOH] CH3Br + OH -> H2O + HO2; 4.000e-12, -1450 [kCH3BrO1D] CH3Br + O1D -> ; 1.800e-10 [kCFC11O1D] CFC11 + O1D -> ; 2.300e-10 [kCFC12O1D] CFC12 + O1D -> ; 1.400e-10 [kCFC113O1D] CFC113 + O1D -> ; 2.000e-10 [kHCFC22OH] HCFC22 + OH -> ; 1.000e-12, -1600 [kHCFC22O1D] HCFC22 + O1D -> H2O; 1.000e-10 [kHCFC141OH] HCFC141 + OH -> ; 1.700e-12, -1700 [kHCFC141O1D] HCFC141 + O1D -> ; 2.600e-10 [kHCFC142OH] HCFC142 + OH -> ; 1.300e-12, -1800 [kHCFC142O1D] HCFC142 + O1D -> ; 2.200e-10 [kHCFC134aOH] HFC134a + OH -> ; 1.500e-12, -1750 [kHFC134aO1D] HFC134a + O1D -> ; 4.900e-11 [kHa1301O1D] Ha1301 + O1D -> ; 1.000e-10 [kHa1211O1D] Ha1211 + O1D -> ; 1.500e-10 [kOO3] O + O3 -> 2 * O2 ; 8e-12, -2060 [kO1DO3] O1D + O3 -> 2 * O2 ; 1.2e-10 [kO1DN2] O1D + N2 -> O + N2 ; 1.8e-11, 110 [kO1DO2] O1D + O2 -> O + O2 ; 3.2e-11, 70 [kO1DH2O] O1D + H2O -> 2 * OH ; 2.2e-10 [kN2OO1D2NO] N2O + O1D -> 2 * NO ; 6.7e-11 [kN2OO1DN2] N2O + O1D -> N2 + O2 ; 4.9e-11 [kNOHO2] NO + HO2 -> NO2 + OH ; 3.5e-12, 250 [kNOO3] NO + O3 -> NO2 + O2 ; 3e-12, -1500 [kNO2O] NO2 + O -> NO + O2 ; 5.6e-12, 180 [kNO2O3] NO2 + O3 -> NO3 + O2 ; 1.2e-13, -2450 [kNO3HO2] NO3 + HO2 -> 0.4 * HNO3 + 0.6 * OH + 0.6 * NO2 ; 2.3e-12, 170. [kHNO4OH] HNO4 + OH -> H2O + NO2 + O2 ; 1.3e-12, 380 [kCH4OH] CH4 + OH -> CH3O2 + H2O ; 2.45e-12, -1775 [kCH4O1DCH3O2] CH4 + O1D -> CH3O2 + OH ; 1.125e-10 [kCH4O1DH2] CH4 + O1D -> H2 + CH2O ; 3.75e-11 [kCH3O2NO] CH3O2 + NO -> CH3O + NO2 ; 3.0e-12, 280 [kCH3OO2] CH3O + O2 -> CH2O + HO2 ; 3.9e-14, -900 [kCH3OHOH] CH3OH + OH -> CH2O + HO2 + H2O ; 6.7e-12, -600 [kCH3O2CH3O2CH3OH] CH3O2 + CH3O2 -> CH3OH + CH2O + O2 ; 1.5e-13, 190 [kCH3O2CH3O22CH2O] CH3O2 + CH3O2 -> 2 * CH3O + O2 ; 1.0e-13, 190 [kCH3O2HO2] CH3O2 + HO2 -> CH3OOH + O2 ; 3.8e-13, 800 [kCH2OOHOHCH2O] CH3OOH + OH -> CH2O + OH + H2O ; 1.14e-12, 200 [kCH3OOHOHCH3O2] CH3OOH + OH -> CH3O2 + H2O ; 2.66e-12, 200 [kCH2OOH] CH2O + OH -> CO + H2O + HO2 ; 1e-11 [kH2O1D] H2 + O1D -> H + OH ; 1.1e-10 [kOOH] O + OH -> H + O2 ; 2.2e-11, 120 [kHO2O] HO2 + O -> OH + O2 ; 3e-11, 200 [kOHO3] OH + O3 -> HO2 + O2 ; 1.5e-12, -880 [kHO2O3] HO2 + O3 -> OH + 2 * O2 ; 2.0e-14, -680 [kH2O2OH] H2O2 + OH -> H2O + HO2 ; 2.9e-12, -160 [kOHHO2] OH + HO2 -> H2O + O2 ; 4.8e-11, 250 [kOHOH] OH + OH -> H2O + O ; 4.2e-12, -240 [kH2OH] H2 + OH -> H2O + H; 5.5e-12, -2000 [kNO2NO] N + O2 -> NO + O ; 1.5e-11, -3600 [kNNO] N + NO -> N2 + O ; 2.1e-11, 100 [kNNO2] N + NO2 -> N2O + O ; 5.8e-12, 220 [kHO3] H + O3 -> OH + O2 ; 1.4e-10, -470 [kHNO2] H + NO2 -> OH + NO ; 4.0e-10, -340 [kHHO22OH] H + HO2 -> 2*OH ; 7.05e-11 [kHHO2H2] H + HO2 -> H2 + O2 ; 7.29e-12 [kHHO2H2O] H + HO2 -> H2O + O ; 1.62e-12 [kCH2OO] CH2O + O -> HO2 + OH + CO ; 3.4e-11, -1600. [kH2SOH] H2S + OH -> SO2 + HO2 ; 6.e-12, -75. [kDMSOHCH2O] DMS + OH -> SO2 + 2 CH2O ; 9.6e-12, -234. [kDMSNO3] DMS + NO3 -> SO2 +NO2 + 2 CH2O ; 1.9e-13,500 [kDMSOOH] DMSO + OH -> 0.6*SO2 + 0.4*ASMSAM + 1.5*CH2O ; 5.8e-11, 0. [kNONO3] NO + NO3 -> 2*NO2 ; 1.7e-11, 133 [kCH3O2NO3] CH3O2 + NO3 -> CH3O + NO2 + O2 ; 3.1e-12 [kCH2O] CH2O + NO3 -> CO + HO2 + HNO3 ; 5.8e-16 [kAIBCM] AIBCM -> ASBCM ; 9.645e-06 [kAIPOMM] AIPOMM -> ASPOMM ; 9.645e-06 [kAIN] AIN -> ASN ; 9.645e-06 Endent * Wet scavenging Heterogeneous HNO3, H2O2, CH2O, CH3OOH, HNO4 SO2,H2S,DMS,DMSO Pb210a, Pb210b Endent Ext Forcing NO Endent ENDPAR TRANSPORT hor-advect OX, Rn222, Pb210a CO2FF MCF CH4, N2O, O3S,H2O2 N, NO, NO2, NO3, HNO3 HNO4, N2O5, CO, CH3OOH CH2O H2, SO2, H2S, DMS, DMSO CO2OC, CO2BIM CO2BIH CIDUSTM, CIN CSSSM, CSSO4M, CSN AIBCM, AIPOMM, AIN ASSSM, ASBCM, ASPOMM, ASSO4M, ASMSAM, ASN SSSSM, SSN * All endlst ver-advect OX, Rn222, Pb210a CO2FF MCF CH4, N2O, O3S,H2O2 N, NO, NO2, NO3, HNO3 HNO4, N2O5, CO, CH3OOH CH2O H2, SO2, H2S, DMS, DMSO CO2OC, CO2BIM CO2BIH CIDUSTM, CIN CSSSM, CSSO4M, CSN AIBCM, AIPOMM, AIN ASSSM, ASBCM, ASPOMM, ASSO4M, ASMSAM, ASN SSSSM, SSN * All endlst convection OX, Rn222, Pb210a CO2FF MCF CH4, N2O, O3S,H2O2 N, NO, NO2, NO3, HNO3 HNO4, N2O5, CO, CH3OOH CH2O H2, SO2, H2S, DMS, DMSO CO2OC, CO2BIM CO2BIH CIDUSTM, CIN CSSSM, CSSO4M, CSN AIBCM, AIPOMM, AIN ASSSM, ASBCM, ASPOMM, ASSO4M, ASMSAM, ASN SSSSM, SSN * All endlst pblayer OX, Rn222, Pb210a CO2FF MCF CH4, N2O, O3S,H2O2 N, NO, NO2, NO3, HNO3 HNO4, N2O5, CO, CH3OOH CH2O H2, SO2, H2S, DMS, DMSO CO2OC, CO2BIM CO2BIH CIDUSTM, CIN CSSSM, CSSO4M, CSN AIBCM, AIPOMM, AIN ASSSM, ASBCM, ASPOMM, ASSO4M, ASMSAM, ASN SSSSM, SSN * All endlst ENDPAR SIMULATION PARAMETERS Spatial Dimensions Longitude points = XNBPOINT Latitude points = 1 Vertical points = XLEV GCM long = XLON GCM lat = XLAT Endent Numerical Control Implicit Iterations = 7 Endent Surface Flux Rn222 SF6, CO2FF CFC11, CFC12, CFC113 HCFC22, HFC134a, HCFC141, HCFC142 Ha1301, Ha1211 MCF, CCl4, CH3Cl, CH3Br CH4, N2O CO, NO, H2 SO2,H2S,DMS CO2OC,CO2BIM,CO2BIH ASSO4M Endent Version Options netcdfdouble = on machine = vpp work dir = $QSUB_WORKDIR/run subfile = lmdz.job multitask = off norms = off fixer = off cpus = 1 relhum = on geoheight = on userhook = on namemod = on transmod = on tavgprnt = on modules = on Endent Execution Options *-- Not activated Qsub file = yes Endent Job Control *-- Not activated Simulation length = 180d Simulation Time Step = 1200 Cray Time Limit = 200h Cray Memory = 120 Case = test Restart = true Crayque = reg Endent Inputs *-- Not activated Start time = 2Y:01Mon:4d Dynamics MSS file = /STACY/ccm2/fo02/hist/h0023 Dynhisttape = short Multi level fields = 15 Single level fields = 6 Ic MSS file = /DIDIER/ctm/mozv1std/hist/h000 Endent Outputs *-- Not activated Write frequency = 1d Start File Number = 001 Density = 4 Retention time = 1y Species Inst * All Endlst Avgr Pb210a HNO3,NO2,NO,NO3,N2O5,H2O2 SO2,DMS CH4,OH,CO CIDUSTM, CIN CSSSM, CSSO4M, CSN AIBCM, AIPOMM, AIN ASSSM, ASBCM, ASPOMM, ASSO4M, ASMSAM, ASN SSSSM, SSN Endlst Endent Group Members Inst * All Endlst Avgr All Endlst Endent Surface Flux Inst * All Endlst Avgr Pb210a CO,NO,NO2,SO2,DMS CIDUSTM, CIN CSSSM, CSSO4M, CSN AIBCM, AIPOMM, AIN ASSSM, ASBCM, ASPOMM, ASSO4M, ASMSAM, ASN SSSSM, SSN Endlst Endent Deposition Velocity Inst * OX, Pb210a Endlst Avgr OX, Pb210a HNO3,NO,NO2,SO2,DMS CIDUSTM, CIN CSSSM, CSSO4M, CSN AIBCM, AIPOMM, AIN ASSSM, ASBCM, ASPOMM, ASSO4M, ASMSAM, ASN SSSSM, SSN Endlst Endent Temperature = avgr, inst Watervapor = avgr, inst Surface Pressure = avgr, inst Photorates Inst * All Endlst Avgr * All Endlst Endent Reaction rates Inst Endlst Avgr Endlst Endent Washout rates Inst * HNO3, H2O2 Endlst Avgr HNO3, H2O2 Endlst Endent External forcing Inst NO Endlst Avgr NO Endlst Endent Print time step = 24h Endent ENDPAR ENDSIM