! subroutine ratecoef ! ! Chemical reaction rate coefficients (formerly in ratec.F). ! use params,only: nlev use cons,only: cgm use flds_ratecoef use flds_atmos,only: tvib,tn use flds_ionatm,only: te,ti implicit none ! ! Local: integer :: k real,parameter :: phic = 0.75 real :: tr,t1,t2,a1,a2,a3,a4,b1,b2,aa ! ! adl,a1d,asg are in module flds_ratecoef: adl = 2.58e-4 a1d = 6.8e-3 asg = 0.0758 do k=1,nlev tr=(ti(k)+tn(k))*0.5 t1=0.667*ti(k)+0.333*tn(k) t2=0.636*ti(k)+0.364*tn(k) a1=t1/300. a2=a1*a1 a3=a2*a1 a4=a2*a2 b1=t2/300. b2=b1*b1 rk1nv(k)=2.82e-11-7.74e-12*a1+1.073e-12*a2-5.17e-14* | a3+9.65e-16*a4 rk1(k)=1.7e-11*(300./t1)**0.77+8.54e-11*exp(-3467./ | t1) rk2(k)=1.533e-12-5.92e-13*b1+8.6e-14*b2 rk2nv(k)=rk2(k) tvib(k) = tn(k) if (tr >= 1500.) then rk3(k)=5.2e-11*(tr/300.)**0.2 else rk3(k)=1.4e-10*(300./tr)**0.44 endif rk4(k)=1.e-10 rk5(k)=4.4e-10 rk6(k)=4.e-10 rk7(k)=2.e-10 rk8(k)=1.e-12 rk9(k)=6.e-11 rk10(k)=5.0e-11 rk12(k)=6.e-10 rk11(k)=rk12(k)*8.*sqrt((ti(k)+tn(k)/16.)/(tn(k)+ti(k)/16.))/9. rk13(k)=9.4e-10 rk14(k)=2.0e-9 rk15(k)=2.4e-9 rk16(k)=3.4e-10 rk17(k)=1.0e-10 rk18(k)=4.0e-10 rk19(k)=4.0e-8*(300./te(k))**0.5 rk20(k)=1.5e-7*(300./te(k))**0.5 rk21(k)=0.047 rk22(k)=0.172 rk23(k)=4.8e-10 rk24(k)=8.0e-10 rk25(k)=5.0e-12 rk26(k)=6.6e-8*(300./te(k))**0.5 rk27(k)=7.e-10 rk28(k)=7.7e-5 alp1(k)=4.2e-7*(300./te(k))**0.85 if (te(k) >= 1200.) then alp2(k)=1.6e-7*(300./te(k))**0.55 else alp2(k)=2.7e-7*(300./te(k))**0.7 endif alp3(k)=1.8e-7*(300./te(k))**0.39 ! ! Odd nitrogen chemistry: ! ! OLD JPL VALUE ! rb1(k)=4.4e-12*exp(-3200./tn(k)) ! NEW JPL VALUE rb1(k)=1.5e-11*exp(-3600./tn(k)) rb2(k)=5.0e-12 ! ! 3.4E-11 IS THE JPL OLD RATE ! rb3(k)=3.4e-11 ! 1.6E-10*EXP(-460./T(K)) THE SISKIND AND RUSCH VALUE ! rb3(k)=1.6e-10*exp(-460./tn(k)) ! THIS ONE IS FROM JANE FOX FOR VIKING DATA ! rb3(k)=2.5e-10*sqrt(tn(k)/300.)*exp(-600./tn(k)) rb3(k)=3.4e-11*sqrt(tn(k)/300.) ! rb4(k)=7.0e-13 rb5(k)=3.6e-10*sqrt(te(k)/300.) rb6(k)=7.e-11 rb7(k)=1.06e-5 rb8(k)=5.0e-11 rb9(k)=1.8e-12*exp(-1370./tn(k)) rb10(k)=3.5e-12*exp(250./tn(k)) rb11(k)=9.3e-12 ! rb12(k)=1.2e-13*exp(-2450./tn(k)) rb12(k)=0. rb13(k)=3.0e-12 rb14(k)=6.7e-11 ! ! OXYGEN - HYDROGEN CHEMISTRY ! METASTABLES ! rkm1(k)=1.8e-11*exp(110./tn(k)) rkm2a(k)=3.2e-11*exp(70./tn(k))*phic rkm2b(k)=3.2e-11*exp(70./tn(k))*(1.-phic) rkm3(k)=2.2e-10 rkm4(k)=1.0e-10 rkm5a(k)=1.4e-10 rkm5b(k)=1.4e-11 rkm6(k)=7.4e-11*exp(120./tn(k)) rkm7a(k)=1.2e-10 rkm7b(k)=1.2e-10 rkm8a(k)=8.e-12 rkm9(k)=3.6e-18*exp(-220./tn(k)) rkm10(k)=1.0e-20 rkm11(k)=1.3e-16 rkm13(k)=2.1e-15 rkm14(k)=4.2e-13 rkm15(k)=2.2e-11 rkm16(k)=8.0e-14 rkm17(k)=3.9e-17 ! ! HYDROGEN-OXYGEN REACTIONS ! rkm20(k)=9.59e-34*exp(480./tn(k)) aa=6.0e-34*(300./tn(k))*2.3 rkm21(k)=aa rkm22(k)=aa rkm23(k)=aa rkm24(k)=8.0e-12*exp(-2060./tn(k)) rkm25(k)=2.2e-11*exp(120./tn(k)) ! JPL RATE COEFFICIENT ! rkm26(k)=3.0e-11*exp(200./tn(k)) ! CHANGE RKM26 = RKM26*0.5 AS RECOMMENDED BY OH STUDIES rkm26(k)=1.5e-11*exp(200./tn(k)) rkm27(k)=1.4e-12*exp(-2000./tn(k)) rkm28(k)=1.6e-11*exp(-4570./tn(k)) rkm29(k)=1.6e-12*exp(-940./tn(k)) rkm30(k)=4.2e-12*exp(-240./tn(k)) rkm31(k)=4.8e-11*exp(250./tn(k)) rkm32(k)=2.9e-12*exp(-160./tn(k)) rkm33(k)=5.5e-12*exp(-2000./tn(k)) rkm34(k)=1.1e-14*exp(-500./tn(k)) rkm35(k)=2.3e-13*exp(600./tn(k)) rkm36(k)=5.7e-32*(300./tn(k))**1.6 rkm37(k)=1.4e-10*exp(-470./tn(k)) rkm38(k)=8.1e-12 rkm39(k)=4.2e-10*exp(-950./tn(k)) rkm40(k)=8.3e-11*exp(-500./tn(k)) rkm41(k)=5.7e-32*(300./tn(k))**1.6 rkm42(k)=1.35e-13*exp(-100./tn(k)) rkm43(k)=2.4e-12*exp(-1710./tn(k)) rkm44(k)=3.5e-11*exp(-4550./tn(k)) rkm45(k)=1.4e-10 ! ! CARBON REACTIONS ! gam1(k)=2.3e-12*exp(-1700./tn(k)) gam2a(k)=1.4e-10 gam2b(k)=1.4e-11 gam3(k)=4.5e-31*(300./tn(k))**3 gam4(k)=1.0e-10 gam5(k)=4.2e-12*exp(180./tn(k)) gam6(k)=3.3e-13*exp(800./tn(k)) gam7(k)=2.2e-13*exp(220./tn(k)) gam8(k)=3.8e-12*exp(200./tn(k)) gam9(k)=3.9e-14*exp(-900./tn(k)) gam10(k)=1.0e-11 gam11(k)=3.4e-11*exp(-1600./tn(k)) gam12(k)=3.5e-12*exp(140./tn(k)) gam13(k)=3.5e-13*sqrt(tn(k))*exp(-3650./tn(k)) ! TOM SLANGER PROVIDED ME WITH THIS RATE COEFFICIENT ! gam14(k)=3.7e-33*exp(-2190./tn(k)) ! SWAMI PROVIDED ME WITH THIS RATE COEFFICIENT gam14(k)=1.7e-33*exp(-1510./tn(k)) ! STROBEL'S SPECULATIVE RATE ! gam15(k)=1.0e-12 gam15(k)=1.0e-13 ! gam15(k)=0. ! ! CHLORINE REACTIONS ! del1(k)=2.9e-11*exp(-260./tn(k)) del2(k)=3.0e-11*exp(+70./tn(k)) del3(k)=1.1e-11*exp(-1400./tn(k)) ! ! SODIUM REACTIONS ! rkn1(k) = 1.1e-9*exp(-116./tn(k)) rkn2(k) = 2.2e-10*sqrt(tn(k)/200.) rkn3a(k) = 1.1e-9*exp(-568./tn(k)) rkn3b(k) = 3.2e-10*exp(-550./tn(k)) rkn4(k) = 4.4e-10*exp(-507./tn(k)) rkn5(k) = 5.3e-30*(200./tn(k)) rkn6(k) = 1.3e-27*(200./tn(k)) rkn7(k) = 5.0e-30*(200./tn(k))**1.22 rkn8(k) = 5.0e-10*exp(-940./tn(k)) rkn9(k) = 1.e-9*exp(-1000./tn(k)) xjn10(k) = 4.e-3*cgm rkn11(k) = 4.e-11*exp(-550./tn(k)) xjn12(k) = 1.e-3*cgm rkn13(k) = 1.9e-28*(200./tn(k)) xjn14(k) = 2.e-5*cgm rkn15(k) = 1.4e-9 rkn16(k) = 1.0e-9 rkn17(k) = 2.6e-30*(200./tn(k))**2.19 rkn18(k) = 1.e-6*(200./tn(k))**0.5 rkn19(k) = 2.5e-10*sqrt(tn(k)/200.) rkn20(k) = 5.e-10*exp(-1200/tn(k)) rkn21(k) = 1.e-9*exp(-1400./tn(k)) rkn22(k) = 1.e-11*exp(-910./tn(k)) xjn23(k) = 0. rkn24a(k)=1.1e-9*exp(-1100./tn(k)) rkn24b(k)=1.1e-9*exp(-1400./tn(k)) rkn25(k) =3.e-10*exp(-668./tn(k)) xjn26(k) =1.e-4*cgm rkn27a(k)=8.e-10 rkn27b(k)=8.e-10 rkn28(k) =6.e-10 rkn29(k) =8.e-10 rkn30(k) =1.e-6*(200./tn(k))**0.5 enddo ! k=1,nlev end subroutine ratecoef