BEGSIM
output unit = 7
file = inca.doc
version = CH4

COMMENTS

     LMDz.4
     INCA version CH4

Didier Hauglustaine, 02/2002
LMDz4 update 02/2004

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
      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
      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
      Endlst

   ENDPAR


   CHEMISTRY

      Photolysis
 [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
 [jNO]       NO + hv -> N + O
      Endent

      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
 [kN2O5] 		N2O5 -> 2 * HNO3
 [kNO2] 		NO2 -> 0.5*HNO3+0.5*HNO2
 [kNO3] 		NO3 -> HNO3
 [kHO2] 		HO2 -> 0.5*H2O2+0.5*O2
 [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.
   Endent

   Heterogeneous
   HNO3, H2O2, CH2O, CH3OOH, HNO4
   Pb210a, Pb210b
   Endent

   Ext Forcing
   NO
   Endent

   ENDPAR


   TRANSPORT

     hor-advect

	OX, Rn222, Pb210a
	CO2FF, Kr85, SF6, CFC11
 	MCF
 	CH4, N2O, O3I, O3S
 	H2O2 
 	NO, NO2, NO3, HNO3, HNO4, N2O5
 	CO
 	CH3OOH, CH2O 
 	H2

*      	All
     endlst

     ver-advect

	OX, Rn222, Pb210a
	CO2FF, Kr85, SF6, CFC11
        MCF
        CH4, N2O, O3I, O3S
        H2O2
        NO, NO2, NO3, HNO3, HNO4, N2O5
        CO
        CH3OOH, CH2O
        H2

*     	All
     endlst

     convection

	OX, Rn222, Pb210a
	CO2FF, Kr85, SF6, CFC11
        MCF
        CH4, N2O, O3I, O3S
        H2O2
        NO, NO2, NO3, HNO3, HNO4, N2O5
        CO
        CH3OOH, CH2O
        H2

*     	All
     endlst

     pblayer

	OX, Rn222, Pb210a
	CO2FF, Kr85, SF6, CFC11
        MCF
        CH4, N2O, O3I, O3S
        H2O2
        NO, NO2, NO3, HNO3, HNO4, N2O5
        CO
        CH3OOH, CH2O
        H2

*     	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
     Endent


     Version Options
        netcdfdouble = on
        machine = sx5
	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

        Retention time = 1y
        Write frequency = 1d
        Start File Number = 001
        Density = 4

	Species
	  Inst
            All
          Endlst
	  Avgr
            All
          Endlst
	Endent

        Group Members
	  Inst
            All
          Endlst
          Avgr
            All
          Endlst
	Endent

	Surface Flux
	  Inst
            All
          Endlst
          Avgr
            All
          Endlst
	Endent

	Deposition Velocity
	  Inst
            Pb210a, Pb210b
          Endlst
	  Avgr
	    Pb210a, Pb210b
	  Endlst
	Endent

	Temperature = avgr, inst
	Watervapor = avgr, inst
	Surface Pressure = avgr, inst

	Photorates
	  Inst
	    All
	  Endlst
	  Avgr
	    All
	  Endlst
	Endent

	Reaction rates
	  Inst
	    76,77,78,79
	  Endlst
	  Avgr
  	    76,77,78,79
	  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 = 12h

     Endent


   ENDPAR


ENDSIM
