! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ! ! Initialization File ! ! Generated by KPP-2.2.3 symbolic chemistry Kinetics PreProcessor ! (http://www.cs.vt.edu/~asandu/Software/KPP) ! KPP is distributed under GPL, the general public licence ! (http://www.gnu.org/copyleft/gpl.html) ! (C) 1995-1997, V. Damian & A. Sandu, CGRER, Univ. Iowa ! (C) 1997-2005, A. Sandu, Michigan Tech, Virginia Tech ! With important contributions from: ! M. Damian, Villanova University, USA ! R. Sander, Max-Planck Institute for Chemistry, Mainz, Germany ! ! File : cbm_Initialize.f90 ! Time : Fri Mar 15 14:04:58 2013 ! Working directory : /home/sandu/kpp-2.2.3/examples/Cbm_fortran ! Equation file : cbm.kpp ! Output root filename : cbm ! ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ MODULE cbm_Initialize USE cbm_Parameters, ONLY: dp, NVAR, NFIX IMPLICIT NONE CONTAINS ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ! ! Initialize - function to initialize concentrations ! Arguments : ! ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ SUBROUTINE Initialize ( ) USE cbm_Global INTEGER :: i REAL(kind=dp) :: x CFACTOR = 2.550000e+10_dp x = (1.0E-8)*CFACTOR DO i = 1, NVAR VAR(i) = x END DO x = (1.0E-8)*CFACTOR DO i = 1, NFIX FIX(i) = x END DO VAR(3) = (1.0)*CFACTOR VAR(5) = (10.0)*CFACTOR VAR(7) = (10.0)*CFACTOR VAR(9) = (1.0)*CFACTOR VAR(16) = (300.0)*CFACTOR VAR(17) = (10.0)*CFACTOR VAR(20) = (50.0)*CFACTOR VAR(21) = (10.0)*CFACTOR VAR(22) = (10.0)*CFACTOR VAR(23) = (10.0)*CFACTOR VAR(24) = (10)*CFACTOR VAR(25) = (100.0)*CFACTOR VAR(26) = (20.0)*CFACTOR VAR(31) = (50.0)*CFACTOR FIX(1) = (1.25E+8)*CFACTOR ! constant rate coefficients RCONST(4) = 9.3e-12 RCONST(11) = 2.2e-10 RCONST(18) = 1.3e-21 RCONST(21) = 4.39999e-40 RCONST(24) = 6.6e-12 RCONST(25) = 1e-20 RCONST(36) = 2.2e-13 RCONST(37) = 1e-11 RCONST(41) = 6.3e-16 RCONST(44) = 2.5e-15 RCONST(49) = 2e-12 RCONST(50) = 6.5e-12 RCONST(52) = 8.1e-13 RCONST(54) = 1600 RCONST(55) = 1.5e-11 RCONST(59) = 7.7e-15 RCONST(64) = 8.1e-12 RCONST(65) = 4.2 RCONST(66) = 4.1e-11 RCONST(67) = 2.2e-11 RCONST(68) = 1.4e-11 RCONST(70) = 3e-11 RCONST(73) = 1.7e-11 RCONST(75) = 1.8e-11 RCONST(76) = 9.6e-11 RCONST(77) = 1.2e-17 RCONST(78) = 3.2e-13 RCONST(79) = 8.1e-12 RCONST(81) = 6.8e-13 ! END constant rate coefficients ! INLINED initializations TSTART = 12.D0*3600.D0 TEND = TSTART + 24.D0*3600.D0 * 5 DT = 3600.D0 TEMP = 288.15 ! End INLINED initializations END SUBROUTINE Initialize ! End of Initialize function ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ END MODULE cbm_Initialize