API tools faq
paste
Advertisement
Guest User

geostuv

a guest
Jun 22nd, 2018
92
0
Never
Add comment
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Fortran 10.45 KB | None | 0 0
raw download clone embed print report
  1. SELECT CASE(jl)
  2.  
  3.  
  4.     CASE(  1) !O2 -> O + O
  5.        j(  1) =  seval(szabin,theta,tmp,tmp2,b,c,d)!O2 + hv = 2.000O
  6.  
  7.     CASE(  3) !O3 -> O2 + O(3P)
  8.        j(  2) =  seval(szabin,theta,tmp,tmp2,b,c,d)!O3 + hv = O + O2
  9.  
  10.     CASE(  2) !O3 -> O2 + O(1D)
  11.        j(  3) =  seval(szabin,theta,tmp,tmp2,b,c,d)!O3 + hv = O1D + O2
  12.  
  13.     CASE(357) !NO -> N + O(3P)
  14.        j(  6) =  seval(szabin,theta,tmp,tmp2,b,c,d)!NO + hv = O + N
  15.  
  16.     CASE( 23) !CH2O -> H2 + CO
  17.        j(  7) =  seval(szabin,theta,tmp,tmp2,b,c,d)!CH2O + hv = HO2 + H + CO
  18.        j(  8) =  seval(szabin,theta,tmp,tmp2,b,c,d)!CH2O + hv = H2 + CO
  19.  
  20.     CASE(  5) !H2O2 -> 2 OH
  21.        j(  9) =  seval(szabin,theta,tmp,tmp2,b,c,d)!H2O2 + hv = OH + OH
  22.  
  23. CASE(33) !unknown
  24. j( 10) =  seval(szabin,theta,tmp,tmp2,b,c,d)!MP + hv = CH2O + HO2 + OH
  25.  
  26.     CASE(  6) !NO2 -> NO + O(3P)
  27.        j( 11) =  seval(szabin,theta,tmp,tmp2,b,c,d)!NO2 + hv = NO + O
  28.  
  29.     CASE(  8) !NO3 -> NO2 + O(3P)
  30.        j( 12) =  seval(szabin,theta,tmp,tmp2,b,c,d)!NO3 + hv = NO2 + O
  31.  
  32.     CASE(  7) !NO3 -> NO + O2
  33.        j( 13) =  seval(szabin,theta,tmp,tmp2,b,c,d)!NO3 + hv = NO + O2
  34.  
  35.     CASE( 11) !N2O5 -> NO3 + NO2
  36.        j( 14) =  seval(szabin,theta,tmp,tmp2,b,c,d)!N2O5 + hv = NO3 + NO2
  37.  
  38.     CASE( 12) !HNO2 -> OH + NO
  39.        j( 15) =  seval(szabin,theta,tmp,tmp2,b,c,d)!HNO2 + hv = OH + NO
  40.  
  41.     CASE( 13) !HNO3 -> OH + NO2
  42.        j( 16) =  seval(szabin,theta,tmp,tmp2,b,c,d)!HNO3 + hv = OH + NO2
  43.  
  44.     CASE( 19) !HNO4 -> OH + NO3
  45.        j( 17) =  seval(szabin,theta,tmp,tmp2,b,c,d)!HNO4 + hv = OH + NO3
  46.  
  47.     CASE( 18) !HNO4 -> HO2 + NO2
  48.        j( 18) =  seval(szabin,theta,tmp,tmp2,b,c,d)!HNO4 + hv = HO2 + NO2
  49.  
  50.     CASE( 97) !ClONO2 -> Cl + NO3
  51.        j( 19) =  seval(szabin,theta,tmp,tmp2,b,c,d)!ClNO3 + hv = Cl + NO3
  52.  
  53.     CASE( 98) !ClONO2 -> ClO + NO2
  54.        j( 20) =  seval(szabin,theta,tmp,tmp2,b,c,d)!ClNO3 + hv = ClO + NO2
  55.  
  56.     CASE( 95) !ClNO2 -> Cl + NO2
  57.        j( 21) =  seval(szabin,theta,tmp,tmp2,b,c,d)!ClNO2 + hv = Cl + NO2
  58.  
  59.     CASE( 86) !Cl2 -> Cl + Cl
  60.        j( 22) =  seval(szabin,theta,tmp,tmp2,b,c,d)!Cl2 + hv = 2.000Cl
  61.  
  62.     CASE(118) !Br2 -> Br + Br
  63.        j( 23) =  seval(szabin,theta,tmp,tmp2,b,c,d)!Br2 + hv = 2.000Br
  64.  
  65.     CASE( 93) !HOCl -> OH + Cl
  66.        j( 24) =  seval(szabin,theta,tmp,tmp2,b,c,d)!HOCl + hv = Cl + OH
  67.  
  68.     CASE( 90) !OClO -> Products
  69.        j( 25) =  seval(szabin,theta,tmp,tmp2,b,c,d)!OClO + hv = ClO + O
  70.  
  71.     CASE( 91) !ClOOCl -> Cl + ClOO
  72.        j( 26) =  seval(szabin,theta,tmp,tmp2,b,c,d)!Cl2O2 + hv = Cl + ClOO
  73.  
  74.     CASE( 88) !ClO -> Cl + O(3P)
  75.        j( 27) =  seval(szabin,theta,tmp,tmp2,b,c,d)!ClO + hv = Cl + O
  76.  
  77.     CASE(119) !BrO -> Br + O
  78.        j( 28) =  seval(szabin,theta,tmp,tmp2,b,c,d)!BrO + hv = Br + O
  79.  
  80.     CASE(126) !BrONO2 -> Br + NO3
  81.        j( 29) =  seval(szabin,theta,tmp,tmp2,b,c,d)!BrNO3 + hv = Br + NO3
  82.  
  83.     CASE(125) !BrONO2 -> BrO + NO2
  84.        j( 30) =  seval(szabin,theta,tmp,tmp2,b,c,d)!BrNO3 + hv = BrO + NO2
  85.  
  86.     CASE(124) !BrNO2 -> Br + NO2
  87.        j( 31) =  seval(szabin,theta,tmp,tmp2,b,c,d)!BrNO2 + hv = Br + NO2
  88.  
  89.     CASE(120) !HOBr -> OH + Br
  90.        j( 32) =  seval(szabin,theta,tmp,tmp2,b,c,d)!HOBr + hv = Br + OH
  91.  
  92.     CASE(127) !BrCl -> Br + Cl
  93.        j( 33) =  seval(szabin,theta,tmp,tmp2,b,c,d)!BrCl + hv = Br + Cl
  94.  
  95.     CASE(358) !OCS -> CO + SO2
  96.        j( 34) =  seval(szabin,theta,tmp,tmp2,b,c,d)!OCS + hv = SO2 + CO
  97.  
  98.     CASE(  9) !N2O -> N2 + O(1D)
  99.        j( 36) =  seval(szabin,theta,tmp,tmp2,b,c,d)!N2O + hv = N2 + O1D
  100.  
  101. CASE(109) !Unknown
  102. j( 37) =  seval(szabin,theta,tmp,tmp2,b,c,d)!CFC11 + hv = 3.000Cl
  103.  
  104. CASE(110) !unknown
  105. j( 38) =  seval(szabin,theta,tmp,tmp2,b,c,d)!CFC12 + hv = 2.000Cl
  106.  
  107. CASE(111) !unknown
  108. j( 39) =  seval(szabin,theta,tmp,tmp2,b,c,d)!CFC113 + hv = 3.000Cl
  109.  
  110. CASE(112) !unknown
  111. j( 40) =  seval(szabin,theta,tmp,tmp2,b,c,d)!CFC114 + hv = 2.000Cl
  112.  
  113. CASE(113) !unknown
  114. j( 41) =  seval(szabin,theta,tmp,tmp2,b,c,d)!CFC115 + hv = Cl
  115.  
  116.     CASE( 99) !CCl4 -> Products
  117.        j( 42) =  seval(szabin,theta,tmp,tmp2,b,c,d)!CCl4 + hv = 4.000Cl
  118.  
  119.     CASE(102) !CH3Cl -> Products
  120.        j( 43) =  seval(szabin,theta,tmp,tmp2,b,c,d)!CH3Cl + hv = MO2 + Cl
  121.  
  122.     CASE(103) !CH3CCl3 -> Products
  123.        j( 44) =  seval(szabin,theta,tmp,tmp2,b,c,d)!CH3CCl3 + hv = 3.000Cl
  124.  
  125. !CASE() !unknown
  126. !j( 45) =  seval(szabin,theta,tmp,tmp2,b,c,d)!CH2Cl2 + hv = 2.000Cl
  127.  
  128. CASE(114) !unknown
  129. j( 46) =  seval(szabin,theta,tmp,tmp2,b,c,d)!HCFC22 + hv = 2.000Cl
  130.  
  131. CASE(115) !unknown
  132. j( 47) =  seval(szabin,theta,tmp,tmp2,b,c,d)!HCFC123 + hv = 2.000Cl
  133.  
  134. CASE(117) !unknown
  135. j( 48) =  seval(szabin,theta,tmp,tmp2,b,c,d)!HCFC141b + hv = 2.000Cl
  136.  
  137. CASE(118) !unknown
  138. j( 49) =  seval(szabin,theta,tmp,tmp2,b,c,d)!HCFC142b + hv = 2.000Cl
  139.  
  140.     CASE(128) !CH3Br -> Products
  141.        j( 50) =  seval(szabin,theta,tmp,tmp2,b,c,d)!CH3Br + hv = MO2 + Br
  142.  
  143. CASE(134) !unknown
  144. j( 51) =  seval(szabin,theta,tmp,tmp2,b,c,d)!H1211 + hv = Cl + Br
  145.  
  146. CASE(135) !unknown
  147. j( 53) =  seval(szabin,theta,tmp,tmp2,b,c,d)!H1301 + hv = Br
  148.  
  149. CASE(136) !unknown
  150. j( 54) =  seval(szabin,theta,tmp,tmp2,b,c,d)!H2402 + hv = 2.000Br
  151.  
  152. CASE(362) !unknown
  153. j( 55) =  seval(szabin,theta,tmp,tmp2,b,c,d)!CH2Br2 + hv = 2.000Br
  154.  
  155.     CASE(129) !CHBr3 -> Products
  156.        j( 56) =  seval(szabin,theta,tmp,tmp2,b,c,d)!CHBr3 + hv = 3.000Br
  157.  
  158.     CASE( 25) !CH3CHO -> CH4 + CO
  159.        j( 62) =  seval(szabin,theta,tmp,tmp2,b,c,d)!ALD2 + hv = CH4 + CO
  160.  
  161.     CASE( 63) !CH3COCH=CH2 -> C3H6 + CO
  162.        j( 63) =  seval(szabin,theta,tmp,tmp2,b,c,d)!MVK + hv = PRPE + CO
  163.  
  164. !CASE() !unknown
  165. !j( 64) =  seval(szabin,theta,tmp,tmp2,b,c,d)!MVK + hv = MCO3 + CH2O + CO + HO2
  166.  
  167. !CASE() !unknown
  168. !j( 65) =  seval(szabin,theta,tmp,tmp2,b,c,d)!MVK + hv = MO2 + RCO3
  169.  
  170.     CASE(360) !MACR -> products
  171.        j( 66) =  seval(szabin,theta,tmp,tmp2,b,c,d)!MACR + hv = CO + HO2 + CH2O + MCO3
  172.  
  173. !CASE() !unknown
  174. !j( 70) =  seval(szabin,theta,tmp,tmp2,b,c,d)!RCHO + hv = ETO2 + HO2 + CO
  175.  
  176. !CASE() !unknown
  177. !j( 71) =  seval(szabin,theta,tmp,tmp2,b,c,d)!MGLY + hv = MCO3 + CO + HO2
  178.  
  179.     CASE( 70) !CHOCHO -> 2 HO2 + 2 CO
  180.        j( 72) =  seval(szabin,theta,tmp,tmp2,b,c,d)!GLYX + hv = 2.000HO2 + 2.000CO
  181.  
  182.     CASE( 71) !CHOCHO -> H2 + 2 CO
  183.        j( 73) =  seval(szabin,theta,tmp,tmp2,b,c,d)!GLYX + hv = H2 + 2.000CO
  184.  
  185.     CASE( 72) !CHOCHO -> CH2O + CO
  186.        j( 74) =  seval(szabin,theta,tmp,tmp2,b,c,d)!GLYX + hv = CH2O + CO
  187.  
  188. !CASE() !unknown
  189. !j( 75) =  seval(szabin,theta,tmp,tmp2,b,c,d)!HAC + hv = MCO3 + CH2O + HO2
  190.  
  191. !CASE() !unknown
  192. !j( 76) =  seval(szabin,theta,tmp,tmp2,b,c,d)!ACET + hv = MCO3 + MO2
  193.  
  194. !CASE() !unknown
  195. !j( 77) =  seval(szabin,theta,tmp,tmp2,b,c,d)!ACET + hv = 2.000MO2 + CO
  196.  
  197. !CASE() !unknown
  198. !j( 78) =  seval(szabin,theta,tmp,tmp2,b,c,d)!INPN + hv = OH + HO2 + RCHO + NO2
  199.  
  200. !CASE() !unknown
  201. !j( 79) =  seval(szabin,theta,tmp,tmp2,b,c,d)!PRPN + hv = OH + HO2 + RCHO + NO2
  202.  
  203. !CASE() !unknown
  204. !j( 80) =  seval(szabin,theta,tmp,tmp2,b,c,d)!ETP + hv = OH + HO2 + ALD2
  205.  
  206. !CASE() !unknown
  207. !j( 81) =  seval(szabin,theta,tmp,tmp2,b,c,d)!RA3P + hv = OH + HO2 + RCHO
  208.  
  209. !CASE() !unknown
  210. !j( 82) =  seval(szabin,theta,tmp,tmp2,b,c,d)!RB3P + hv = OH + HO2 + ACET
  211.  
  212. !CASE() !unknown
  213. !j( 83) =  seval(szabin,theta,tmp,tmp2,b,c,d)!R4P + hv = OH + HO2 + RCHO
  214.  
  215. !CASE() !unknown
  216. !j( 84) =  seval(szabin,theta,tmp,tmp2,b,c,d)!PP + hv = OH + HO2 + ALD2 + CH2O
  217.  
  218. !CASE() !unknown
  219. !j( 85) =  seval(szabin,theta,tmp,tmp2,b,c,d)!RP + hv = OH + HO2 + ALD2
  220.  
  221. !CASE() !unknown
  222. !j( 88) =  seval(szabin,theta,tmp,tmp2,b,c,d)!ISNP + hv = OH + HO2 + RCHO + NO2
  223.  
  224. !CASE() !unknown
  225. !j( 91) =  seval(szabin,theta,tmp,tmp2,b,c,d)!MAOP + hv = OH + CH2O + MCO3
  226.  
  227. !CASE() !unknown
  228. !j( 93) =  seval(szabin,theta,tmp,tmp2,b,c,d)!MVKN + hv = GLYC + NO2 + MCO3
  229.  
  230. !CASE() !unknown
  231. !j( 94) =  seval(szabin,theta,tmp,tmp2,b,c,d)!ISOPNB + hv = HC5 + NO2 + HO2
  232.  
  233. !CASE() !unknown
  234. !j( 95) =  seval(szabin,theta,tmp,tmp2,b,c,d)!ISOPND + hv = HC5 + NO2 + HO2
  235.  
  236. !CASE() !unknown
  237. !j( 96) =  seval(szabin,theta,tmp,tmp2,b,c,d)!PROPNN + hv = CH2O + NO2 + CO + MO2
  238.  
  239. !CASE() !unknown
  240. !j( 97) =  seval(szabin,theta,tmp,tmp2,b,c,d)!ATOOH + hv = OH + CH2O + MCO3
  241.  
  242. !CASE() !unknown
  243. !j( 99) =  seval(szabin,theta,tmp,tmp2,b,c,d)!MAP + hv = OH + MO2
  244.  
  245. !CASE() !unknown
  246. !j(100) =  seval(szabin,theta,tmp,tmp2,b,c,d)!SO4 + hv = SO2 + 2.000OH
  247.  
  248.     CASE( 89) !ClOO -> Products
  249.        j(102) =  seval(szabin,theta,tmp,tmp2,b,c,d)!ClOO + hv = Cl + O2
  250.  
  251. !CASE() !unknown
  252. !j(104) =  seval(szabin,theta,tmp,tmp2,b,c,d)!MPN + hv = CH2O + NO3 + HO2
  253.  
  254. !CASE() !unknown
  255. !j(105) =  seval(szabin,theta,tmp,tmp2,b,c,d)!MPN + hv = MO2 + NO2
  256.  
  257. !CASE() !unknown
  258. !j(106) =  seval(szabin,theta,tmp,tmp2,b,c,d)!PIP + hv = RCHO + OH + HO2
  259.  
  260. !CASE() !unknown
  261. !j(108) =  seval(szabin,theta,tmp,tmp2,b,c,d)!ETHLN + hv = NO2 + CH2O + CO + HO2
  262.  
  263. !CASE() !unknown
  264. !j(109) =  seval(szabin,theta,tmp,tmp2,b,c,d)!DHDC + hv = MGLY + GLYX + 2.000OH
  265.  
  266. !CASE() !unknown
  267. !j(112) =  seval(szabin,theta,tmp,tmp2,b,c,d)!MONITS + hv = MEK + NO2
  268.  
  269. !CASE() !unknown
  270. !j(113) =  seval(szabin,theta,tmp,tmp2,b,c,d)!MONITU + hv = RCHO + NO2
  271.  
  272. !CASE() !unknown
  273. !j(114) =  seval(szabin,theta,tmp,tmp2,b,c,d)!HONIT + hv = HAC + NO2
  274.  
  275. !CASE() !unknown
  276. !j(115) =  seval(szabin,theta,tmp,tmp2,b,c,d)!I2 + hv = 2.000I
  277.  
  278. !CASE() !unknown
  279. !j(116) =  seval(szabin,theta,tmp,tmp2,b,c,d)!HOI + hv = I + OH
  280.  
  281. !CASE() !unknown
  282. !j(117) =  seval(szabin,theta,tmp,tmp2,b,c,d)!IO + hv = I + O
  283.  
  284. !CASE() !unknown
  285. !j(118) =  seval(szabin,theta,tmp,tmp2,b,c,d)!OIO + hv = I + O2
  286.  
  287. !CASE() !unknown
  288. !j(119) =  seval(szabin,theta,tmp,tmp2,b,c,d)!INO + hv = I + NO
  289.  
  290. !CASE() !unknown
  291. !j(120) =  seval(szabin,theta,tmp,tmp2,b,c,d)!IONO + hv = I + NO2
  292.  
  293. !CASE() !unknown
  294. !j(121) =  seval(szabin,theta,tmp,tmp2,b,c,d)!IONO2 + hv = I + NO3
  295.  
  296. !CASE() !unknown
  297. !j(122) =  seval(szabin,theta,tmp,tmp2,b,c,d)!I2O2 + hv = I + OIO
  298.  
  299. !CASE() !unknown
  300. !j(123) =  seval(szabin,theta,tmp,tmp2,b,c,d)!CH3I + hv = I
  301.  
  302. !CASE() !unknown
  303. !j(124) =  seval(szabin,theta,tmp,tmp2,b,c,d)!CH2I2 + hv = 2.000I
  304.  
  305. !CASE() !unknown
  306. !j(125) =  seval(szabin,theta,tmp,tmp2,b,c,d)!CH2ICl + hv = I + Cl
  307.  
  308. !CASE() !unknown
  309. !j(126) =  seval(szabin,theta,tmp,tmp2,b,c,d)!CH2IBr + hv = I + Br
  310.  
  311. !CASE() !unknown
  312. !j(127) =  seval(szabin,theta,tmp,tmp2,b,c,d)!I2O4 + hv = 2.000OIO
  313.  
  314. !CASE() !unknown
  315. !j(128) =  seval(szabin,theta,tmp,tmp2,b,c,d)!I2O3 + hv = OIO + IO
  316.  
  317. !CASE() !unknown
  318. !j(129) =  seval(szabin,theta,tmp,tmp2,b,c,d)!IBr + hv = I + Br
  319.  
  320. !CASE() !unknown
  321. !j(130) =  seval(szabin,theta,tmp,tmp2,b,c,d)!ICl + hv = I + Cl
  322.  
  323. END SELECT
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement
Public Pastes
Advertisement
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy.  OK, I Understand
Not a member of Pastebin yet?
Sign Up, it unlocks many cool features!