
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
