BASIC.f90

    &TPK%K020(:)*TPK%O1D(:)*TPK%N2(:)+TPK%K021(:)*TPK%O1D(:)*TPK%O2(:)+&
    &0.00000*TPK%K079(:)*PCONC(:,JP_ALKE)*PCONC(:,JP_O3)+&
    &0.09000*TPK%K080(:)*PCONC(:,JP_BIO)*PCONC(:,JP_O3))/&
   &(TPK%K018(:)*TPK%O2(:)+TPK%K019(:)*PCONC(:,JP_O3)+TPK%K029(:)*PCONC(:,JP_NO)+&
   &TPK%K030(:)*PCONC(:,JP_NO2)+TPK%K031(:)*PCONC(:,JP_NO2)+TPK%K054(:)*PCONC(:,JP_BIO)+&
   &TPK%K055(:)*PCONC(:,JP_CARBO))
! /END_CODE/
CALL SUB0
CALL SUB1
CALL SUB2
CALL SUB3
CALL SUB4
CALL SUB5
CALL SUB6
CALL SUB7
CALL SUB8
CALL SUB9

CONTAINS

SUBROUTINE SUB0
!
!Indices 1 a 10
!
!
!PPROD(O3) = +K018*<O3P>*<O2>+0.17307*K0102*<CARBOP>*<HO2>+KTC21*<WC_O3>+KTR21*
!<WR_O3>
 PPROD(:,1) = +TPK%K018(:)*TPK%O3P(:)*TPK%O2(:)+0.17307*TPK%K0102(:)*PCONC(:,39&
&)*PCONC(:,15)+TPK%KTC21(:)*PCONC(:,42)+TPK%KTR21(:)*PCONC(:,67)
!PLOSS(O3) = +K002+K003+K019*<O3P>+K023*<OH>+K024*<HO2>+K042*<NO>+K043*<NO2>+K0
!79*<ALKE>+K080*<BIO>+K081*<CARBO>+K082*<PAN>+K087*<ADD>+KTC1+KTR1
 PLOSS(:,1) = +TPK%K002(:)+TPK%K003(:)+TPK%K019(:)*TPK%O3P(:)+TPK%K023(:)*PCONC&
&(:,14)+TPK%K024(:)*PCONC(:,15)+TPK%K042(:)*PCONC(:,3)+TPK%K043(:)*PCONC(:,4)+T&
&PK%K079(:)*PCONC(:,19)+TPK%K080(:)*PCONC(:,20)+TPK%K081(:)*PCONC(:,25)+TPK%K08&
&2(:)*PCONC(:,27)+TPK%K087(:)*PCONC(:,37)+TPK%KTC1(:)+TPK%KTR1(:)
!
!PPROD(H2O2) = +K027*<HO2>*<HO2>+K028*<HO2>*<HO2>*<H2O>+0.01833*K079*<ALKE>*<O3
!>+0.00100*K080*<BIO>*<O3>+KTC22*<WC_H2O2>+KTR22*<WR_H2O2>
 PPROD(:,2) = +TPK%K027(:)*PCONC(:,15)*PCONC(:,15)+TPK%K028(:)*PCONC(:,15)*PCON&
&C(:,15)*TPK%H2O(:)+0.01833*TPK%K079(:)*PCONC(:,19)*PCONC(:,1)+0.00100*TPK%K080&
&(:)*PCONC(:,20)*PCONC(:,1)+TPK%KTC22(:)*PCONC(:,43)+TPK%KTR22(:)*PCONC(:,68)
!PLOSS(H2O2) = +K009+K026*<OH>+KTC2+KTR2
 PLOSS(:,2) = +TPK%K009(:)+TPK%K026(:)*PCONC(:,14)+TPK%KTC2(:)+TPK%KTR2(:)
!
!PPROD(NO) = +K001*<NO2>+K004*<HONO>+K007*<NO3>+K030*<O3P>*<NO2>+K046*<NO3>*<NO
!2>+KTC23*<WC_NO>+KTR23*<WR_NO>
 PPROD(:,3) = +TPK%K001(:)*PCONC(:,4)+TPK%K004(:)*PCONC(:,7)+TPK%K007(:)*PCONC(&
&:,5)+TPK%K030(:)*TPK%O3P(:)*PCONC(:,4)+TPK%K046(:)*PCONC(:,5)*PCONC(:,4)+TPK%K&
&TC23(:)*PCONC(:,44)+TPK%KTR23(:)*PCONC(:,69)
!PLOSS(NO) = +K029*<O3P>+K032*<OH>+K035*<HO2>+K042*<O3>+K044*<NO>*<O2>+K044*<NO
!>*<O2>+K045*<NO3>+K090*<MO2>+K091*<ALKAP>+K092*<ALKEP>+K093*<BIOP>+K094*<AROP>
!+K095*<CARBOP>+K096*<OLN>+K130*<XO2>+KTC3+KTR3
 PLOSS(:,3) = +TPK%K029(:)*TPK%O3P(:)+TPK%K032(:)*PCONC(:,14)+TPK%K035(:)*PCONC&
&(:,15)+TPK%K042(:)*PCONC(:,1)+TPK%K044(:)*PCONC(:,3)*TPK%O2(:)+TPK%K044(:)*PCO&
&NC(:,3)*TPK%O2(:)+TPK%K045(:)*PCONC(:,5)+TPK%K090(:)*PCONC(:,32)+TPK%K091(:)*P&
&CONC(:,33)+TPK%K092(:)*PCONC(:,34)+TPK%K093(:)*PCONC(:,35)+TPK%K094(:)*PCONC(:&
&,38)+TPK%K095(:)*PCONC(:,39)+TPK%K096(:)*PCONC(:,40)+TPK%K130(:)*PCONC(:,41)+T&
&PK%KTC3(:)+TPK%KTR3(:)
!
!PPROD(NO2) = +K005*<HNO3>+0.65*K006*<HNO4>+K008*<NO3>+K017*<ONIT>+K029*<O3P>*<
!NO>+K034*<OH>*<NO3>+K035*<HO2>*<NO>+K037*<HNO4>+0.7*K038*<HO2>*<NO3>+K039*<OH>
!*<HONO>+K041*<OH>*<HNO4>+K042*<O3>*<NO>+K044*<NO>*<NO>*<O2>+K044*<NO>*<NO>*<O2
!>+K045*<NO3>*<NO>+K045*<NO3>*<NO>+K046*<NO3>*<NO2>+K048*<N2O5>+K049*<NO3>*<NO3
!>+K049*<NO3>*<NO3>+K071*<ONIT>*<OH>+0.10530*K074*<CARBO>*<NO3>+0.40*K078*<PAN>
!*<NO3>+0.70*K082*<PAN>*<O3>+K089*<PAN>+K090*<MO2>*<NO>+0.91541*K091*<ALKAP>*<N
!O>+K092*<ALKEP>*<NO>+0.84700*K093*<BIOP>*<NO>+0.95115*K094*<AROP>*<NO>+K095*<C
!ARBOP>*<NO>+1.81599*K096*<OLN>*<NO>+0.32440*K110*<OLN>*<MO2>+0.00000*K116*<OLN
!>*<CARBOP>+0.00000*K118*<OLN>*<OLN>+K119*<MO2>*<NO3>+K120*<ALKAP>*<NO3>+K121*<
!ALKEP>*<NO3>+K122*<BIOP>*<NO3>+K123*<AROP>*<NO3>+K124*<CARBOP>*<NO3>+1.74072*K
!125*<OLN>*<NO3>+K130*<XO2>*<NO>+K131*<XO2>*<NO3>+KTC24*<WC_NO2>+KTR24*<WR_NO2>
 PPROD(:,4) = +TPK%K005(:)*PCONC(:,8)+0.65*TPK%K006(:)*PCONC(:,9)+TPK%K008(:)*P&
&CONC(:,5)+TPK%K017(:)*PCONC(:,26)+TPK%K029(:)*TPK%O3P(:)*PCONC(:,3)+TPK%K034(:&
&)*PCONC(:,14)*PCONC(:,5)+TPK%K035(:)*PCONC(:,15)*PCONC(:,3)+TPK%K037(:)*PCONC(&
&:,9)+0.7*TPK%K038(:)*PCONC(:,15)*PCONC(:,5)+TPK%K039(:)*PCONC(:,14)*PCONC(:,7)&
&+TPK%K041(:)*PCONC(:,14)*PCONC(:,9)+TPK%K042(:)*PCONC(:,1)*PCONC(:,3)+TPK%K044&
&(:)*PCONC(:,3)*PCONC(:,3)*TPK%O2(:)+TPK%K044(:)*PCONC(:,3)*PCONC(:,3)*TPK%O2(:&
&)+TPK%K045(:)*PCONC(:,5)*PCONC(:,3)+TPK%K045(:)*PCONC(:,5)*PCONC(:,3)+TPK%K046&
&(:)*PCONC(:,5)*PCONC(:,4)+TPK%K048(:)*PCONC(:,6)+TPK%K049(:)*PCONC(:,5)*PCONC(&
&:,5)+TPK%K049(:)*PCONC(:,5)*PCONC(:,5)+TPK%K071(:)*PCONC(:,26)*PCONC(:,14)+0.1&
&0530*TPK%K074(:)*PCONC(:,25)*PCONC(:,5)+0.40*TPK%K078(:)*PCONC(:,27)*PCONC(:,5&
&)+0.70*TPK%K082(:)*PCONC(:,27)*PCONC(:,1)+TPK%K089(:)*PCONC(:,27)+TPK%K090(:)*&
&PCONC(:,32)*PCONC(:,3)+0.91541*TPK%K091(:)*PCONC(:,33)*PCONC(:,3)+TPK%K092(:)*&
&PCONC(:,34)*PCONC(:,3)+0.84700*TPK%K093(:)*PCONC(:,35)*PCONC(:,3)+0.95115*TPK%&
&K094(:)*PCONC(:,38)*PCONC(:,3)+TPK%K095(:)*PCONC(:,39)*PCONC(:,3)+1.81599*TPK%&
&K096(:)*PCONC(:,40)*PCONC(:,3)+0.32440*TPK%K110(:)*PCONC(:,40)*PCONC(:,32)+0.0&
&0000*TPK%K116(:)*PCONC(:,40)*PCONC(:,39)+0.00000*TPK%K118(:)*PCONC(:,40)*PCONC&
&(:,40)+TPK%K119(:)*PCONC(:,32)*PCONC(:,5)+TPK%K120(:)*PCONC(:,33)*PCONC(:,5)+T&
&PK%K121(:)*PCONC(:,34)*PCONC(:,5)+TPK%K122(:)*PCONC(:,35)*PCONC(:,5)+TPK%K123(&
&:)*PCONC(:,38)*PCONC(:,5)+TPK%K124(:)*PCONC(:,39)*PCONC(:,5)+1.74072*TPK%K125(&
&:)*PCONC(:,40)*PCONC(:,5)+TPK%K130(:)*PCONC(:,41)*PCONC(:,3)+TPK%K131(:)*PCONC&
&(:,41)*PCONC(:,5)+TPK%KTC24(:)*PCONC(:,45)+TPK%KTR24(:)*PCONC(:,70)
!PLOSS(NO2) = +K001+K030*<O3P>+K031*<O3P>+K033*<OH>+K036*<HO2>+K043*<O3>+K046*<
!NO3>+K047*<NO3>+K083*<PHO>+K085*<ADD>+K088*<CARBOP>+KTC4+KTR4
 PLOSS(:,4) = +TPK%K001(:)+TPK%K030(:)*TPK%O3P(:)+TPK%K031(:)*TPK%O3P(:)+TPK%K0&
&33(:)*PCONC(:,14)+TPK%K036(:)*PCONC(:,15)+TPK%K043(:)*PCONC(:,1)+TPK%K046(:)*P&
&CONC(:,5)+TPK%K047(:)*PCONC(:,5)+TPK%K083(:)*PCONC(:,36)+TPK%K085(:)*PCONC(:,3&
&7)+TPK%K088(:)*PCONC(:,39)+TPK%KTC4(:)+TPK%KTR4(:)
!
!PPROD(NO3) = +0.35*K006*<HNO4>+K031*<O3P>*<NO2>+K040*<OH>*<HNO3>+K043*<O3>*<NO
!2>+K048*<N2O5>+0.71893*K070*<PAN>*<OH>+0.60*K078*<PAN>*<NO3>+KTC25*<WC_NO3>+KT
!R25*<WR_NO3>
 PPROD(:,5) = +0.35*TPK%K006(:)*PCONC(:,9)+TPK%K031(:)*TPK%O3P(:)*PCONC(:,4)+TP&
&K%K040(:)*PCONC(:,14)*PCONC(:,8)+TPK%K043(:)*PCONC(:,1)*PCONC(:,4)+TPK%K048(:)&
&*PCONC(:,6)+0.71893*TPK%K070(:)*PCONC(:,27)*PCONC(:,14)+0.60*TPK%K078(:)*PCONC&
&(:,27)*PCONC(:,5)+TPK%KTC25(:)*PCONC(:,46)+TPK%KTR25(:)*PCONC(:,71)
!PLOSS(NO3) = +K007+K008+K034*<OH>+K038*<HO2>+K045*<NO>+K046*<NO2>+K047*<NO2>+K
!049*<NO3>+K049*<NO3>+K072*<HCHO>+K073*<ALD>+K074*<CARBO>+K075*<ARO>+K076*<ALKE
!>+K077*<BIO>+K078*<PAN>+K119*<MO2>+K120*<ALKAP>+K121*<ALKEP>+K122*<BIOP>+K123*
!<AROP>+K124*<CARBOP>+K125*<OLN>+K131*<XO2>+KTC5+KTR5
 PLOSS(:,5) = +TPK%K007(:)+TPK%K008(:)+TPK%K034(:)*PCONC(:,14)+TPK%K038(:)*PCON&
&C(:,15)+TPK%K045(:)*PCONC(:,3)+TPK%K046(:)*PCONC(:,4)+TPK%K047(:)*PCONC(:,4)+T&
&PK%K049(:)*PCONC(:,5)+TPK%K049(:)*PCONC(:,5)+TPK%K072(:)*PCONC(:,22)+TPK%K073(&
&:)*PCONC(:,23)+TPK%K074(:)*PCONC(:,25)+TPK%K075(:)*PCONC(:,21)+TPK%K076(:)*PCO&
&NC(:,19)+TPK%K077(:)*PCONC(:,20)+TPK%K078(:)*PCONC(:,27)+TPK%K119(:)*PCONC(:,3&
&2)+TPK%K120(:)*PCONC(:,33)+TPK%K121(:)*PCONC(:,34)+TPK%K122(:)*PCONC(:,35)+TPK&
&%K123(:)*PCONC(:,38)+TPK%K124(:)*PCONC(:,39)+TPK%K125(:)*PCONC(:,40)+TPK%K131(&
&:)*PCONC(:,41)+TPK%KTC5(:)+TPK%KTR5(:)
!
!PPROD(N2O5) = +K047*<NO3>*<NO2>+KTC26*<WC_N2O5>+KTR26*<WR_N2O5>
 PPROD(:,6) = +TPK%K047(:)*PCONC(:,5)*PCONC(:,4)+TPK%KTC26(:)*PCONC(:,47)+TPK%K&
&TR26(:)*PCONC(:,72)
!PLOSS(N2O5) = +K048+KTC6+KTR6
 PLOSS(:,6) = +TPK%K048(:)+TPK%KTC6(:)+TPK%KTR6(:)
!
!PPROD(HONO) = +K032*<OH>*<NO>+K085*<ADD>*<NO2>+KTC27*<WC_HONO>+KTR27*<WR_HONO>
 PPROD(:,7) = +TPK%K032(:)*PCONC(:,14)*PCONC(:,3)+TPK%K085(:)*PCONC(:,37)*PCONC&
&(:,4)+TPK%KTC27(:)*PCONC(:,48)+TPK%KTR27(:)*PCONC(:,73)
!PLOSS(HONO) = +K004+K039*<OH>+KTC7+KTR7
 PLOSS(:,7) = +TPK%K004(:)+TPK%K039(:)*PCONC(:,14)+TPK%KTC7(:)+TPK%KTR7(:)
!
!PPROD(HNO3) = +K033*<OH>*<NO2>+0.3*K038*<HO2>*<NO3>+K072*<HCHO>*<NO3>+K073*<AL
!D>*<NO3>+0.91567*K074*<CARBO>*<NO3>+K075*<ARO>*<NO3>+KTC28*<WC_HNO3>+KTR28*<WR
!_HNO3>
 PPROD(:,8) = +TPK%K033(:)*PCONC(:,14)*PCONC(:,4)+0.3*TPK%K038(:)*PCONC(:,15)*P&
&CONC(:,5)+TPK%K072(:)*PCONC(:,22)*PCONC(:,5)+TPK%K073(:)*PCONC(:,23)*PCONC(:,5&
&)+0.91567*TPK%K074(:)*PCONC(:,25)*PCONC(:,5)+TPK%K075(:)*PCONC(:,21)*PCONC(:,5&
&)+TPK%KTC28(:)*PCONC(:,49)+TPK%KTR28(:)*PCONC(:,74)
!PLOSS(HNO3) = +K005+K040*<OH>+KTC8+KTR8
 PLOSS(:,8) = +TPK%K005(:)+TPK%K040(:)*PCONC(:,14)+TPK%KTC8(:)+TPK%KTR8(:)
!
!PPROD(HNO4) = +K036*<HO2>*<NO2>+KTC29*<WC_HNO4>+KTR29*<WR_HNO4>
 PPROD(:,9) = +TPK%K036(:)*PCONC(:,15)*PCONC(:,4)+TPK%KTC29(:)*PCONC(:,50)+TPK%&
&KTR29(:)*PCONC(:,75)
!PLOSS(HNO4) = +K006+K037+K041*<OH>+KTC9+KTR9
 PLOSS(:,9) = +TPK%K006(:)+TPK%K037(:)+TPK%K041(:)*PCONC(:,14)+TPK%KTC9(:)+TPK%&
&KTR9(:)
!
!PPROD(NH3) = +KTC30*<WC_NH3>+KTR30*<WR_NH3>
 PPROD(:,10) = +TPK%KTC30(:)*PCONC(:,51)+TPK%KTR30(:)*PCONC(:,76)
!PLOSS(NH3) = +K050*<OH>+KTC10+KTR10
 PLOSS(:,10) = +TPK%K050(:)*PCONC(:,14)+TPK%KTC10(:)+TPK%KTR10(:)
!
RETURN
END SUBROUTINE SUB0
!
SUBROUTINE SUB1
!
!Indices 11 a 20
!
!
!PPROD(SO2) = +KTC34*<WC_SO2>+KTR34*<WR_SO2>
 PPROD(:,11) = +TPK%KTC34(:)*PCONC(:,55)+TPK%KTR34(:)*PCONC(:,80)
!PLOSS(SO2) = +K052*<OH>+KTC14+KTR14
 PLOSS(:,11) = +TPK%K052(:)*PCONC(:,14)+TPK%KTC14(:)+TPK%KTR14(:)
!
!PPROD(SULF) = +K052*<OH>*<SO2>+KTC35*<WC_SULF>+KTR35*<WR_SULF>
 PPROD(:,12) = +TPK%K052(:)*PCONC(:,14)*PCONC(:,11)+TPK%KTC35(:)*PCONC(:,56)+TP&
&K%KTR35(:)*PCONC(:,81)
!PLOSS(SULF) = +K132+KTC15+KTR15
 PLOSS(:,12) = +TPK%K132(:)+TPK%KTC15(:)+TPK%KTR15(:)
!
!PPROD(CO) = +K010*<HCHO>+K011*<HCHO>+K012*<ALD>+0.91924*K016*<CARBO>+0.01*K054
!*<BIO>*<O3P>+0.00878*K058*<ALKA>*<OH>+K062*<HCHO>*<OH>+1.01732*K065*<CARBO>*<O
!H>+K072*<HCHO>*<NO3>+1.33723*K074*<CARBO>*<NO3>+0.35120*K079*<ALKE>*<O3>+0.360
!00*K080*<BIO>*<O3>+0.64728*K081*<CARBO>*<O3>+0.13*K082*<PAN>*<O3>
 PPROD(:,13) = +TPK%K010(:)*PCONC(:,22)+TPK%K011(:)*PCONC(:,22)+TPK%K012(:)*PCO&
&NC(:,23)+0.91924*TPK%K016(:)*PCONC(:,25)+0.01*TPK%K054(:)*PCONC(:,20)*TPK%O3P(&
&:)+0.00878*TPK%K058(:)*PCONC(:,18)*PCONC(:,14)+TPK%K062(:)*PCONC(:,22)*PCONC(:&
&,14)+1.01732*TPK%K065(:)*PCONC(:,25)*PCONC(:,14)+TPK%K072(:)*PCONC(:,22)*PCONC&
&(:,5)+1.33723*TPK%K074(:)*PCONC(:,25)*PCONC(:,5)+0.35120*TPK%K079(:)*PCONC(:,1&
&9)*PCONC(:,1)+0.36000*TPK%K080(:)*PCONC(:,20)*PCONC(:,1)+0.64728*TPK%K081(:)*P&
&CONC(:,25)*PCONC(:,1)+0.13*TPK%K082(:)*PCONC(:,27)*PCONC(:,1)
!PLOSS(CO) = +K053*<OH>
 PLOSS(:,13) = +TPK%K053(:)*PCONC(:,14)
!
!PPROD(OH) = +K004*<HONO>+K005*<HNO3>+0.35*K006*<HNO4>+K009*<H2O2>+K009*<H2O2>+
!K013*<OP1>+K014*<OP2>+K022*<O1D>*<H2O>+K022*<O1D>*<H2O>+K024*<O3>*<HO2>+K035*<
!HO2>*<NO>+0.7*K038*<HO2>*<NO3>+0.02*K054*<BIO>*<O3P>+0.00878*K058*<ALKA>*<OH>+
!0.35*K068*<OP1>*<OH>+0.44925*K069*<OP2>*<OH>+0.39435*K079*<ALKE>*<O3>+0.28000*
!K080*<BIO>*<O3>+0.20595*K081*<CARBO>*<O3>+0.036*K082*<PAN>*<O3>+K087*<ADD>*<O3
!>+KTC31*<WC_OH>+KTR31*<WR_OH>
 PPROD(:,14) = +TPK%K004(:)*PCONC(:,7)+TPK%K005(:)*PCONC(:,8)+0.35*TPK%K006(:)*&
&PCONC(:,9)+TPK%K009(:)*PCONC(:,2)+TPK%K009(:)*PCONC(:,2)+TPK%K013(:)*PCONC(:,2&
&8)+TPK%K014(:)*PCONC(:,29)+TPK%K022(:)*TPK%O1D(:)*TPK%H2O(:)+TPK%K022(:)*TPK%O&
&1D(:)*TPK%H2O(:)+TPK%K024(:)*PCONC(:,1)*PCONC(:,15)+TPK%K035(:)*PCONC(:,15)*PC&
&ONC(:,3)+0.7*TPK%K038(:)*PCONC(:,15)*PCONC(:,5)+0.02*TPK%K054(:)*PCONC(:,20)*T&
&PK%O3P(:)+0.00878*TPK%K058(:)*PCONC(:,18)*PCONC(:,14)+0.35*TPK%K068(:)*PCONC(:&
&,28)*PCONC(:,14)+0.44925*TPK%K069(:)*PCONC(:,29)*PCONC(:,14)+0.39435*TPK%K079(&
&:)*PCONC(:,19)*PCONC(:,1)+0.28000*TPK%K080(:)*PCONC(:,20)*PCONC(:,1)+0.20595*T&
&PK%K081(:)*PCONC(:,25)*PCONC(:,1)+0.036*TPK%K082(:)*PCONC(:,27)*PCONC(:,1)+TPK&
&%K087(:)*PCONC(:,37)*PCONC(:,1)+TPK%KTC31(:)*PCONC(:,52)+TPK%KTR31(:)*PCONC(:,&
&77)
!PLOSS(OH) = +K023*<O3>+K025*<HO2>+K026*<H2O2>+K032*<NO>+K033*<NO2>+K034*<NO3>+
!K039*<HONO>+K040*<HNO3>+K041*<HNO4>+K050*<NH3>+K051*<H2>+K052*<SO2>+K053*<CO>+
!K056*<CH4>+K057*<ETH>+K058*<ALKA>+K059*<ALKE>+K060*<BIO>+K061*<ARO>+K062*<HCHO
!>+K063*<ALD>+K064*<KET>+K065*<CARBO>+K066*<ORA1>+K067*<ORA2>+K068*<OP1>+K069*<
!OP2>+K070*<PAN>+K071*<ONIT>+KTC11+KTR11
 PLOSS(:,14) = +TPK%K023(:)*PCONC(:,1)+TPK%K025(:)*PCONC(:,15)+TPK%K026(:)*PCON&
&C(:,2)+TPK%K032(:)*PCONC(:,3)+TPK%K033(:)*PCONC(:,4)+TPK%K034(:)*PCONC(:,5)+TP&
&K%K039(:)*PCONC(:,7)+TPK%K040(:)*PCONC(:,8)+TPK%K041(:)*PCONC(:,9)+TPK%K050(:)&
&*PCONC(:,10)+TPK%K051(:)*TPK%H2(:)+TPK%K052(:)*PCONC(:,11)+TPK%K053(:)*PCONC(:&
&,13)+TPK%K056(:)*PCONC(:,16)+TPK%K057(:)*PCONC(:,17)+TPK%K058(:)*PCONC(:,18)+T&
&PK%K059(:)*PCONC(:,19)+TPK%K060(:)*PCONC(:,20)+TPK%K061(:)*PCONC(:,21)+TPK%K06&
&2(:)*PCONC(:,22)+TPK%K063(:)*PCONC(:,23)+TPK%K064(:)*PCONC(:,24)+TPK%K065(:)*P&
&CONC(:,25)+TPK%K066(:)*PCONC(:,30)+TPK%K067(:)*PCONC(:,31)+TPK%K068(:)*PCONC(:&
&,28)+TPK%K069(:)*PCONC(:,29)+TPK%K070(:)*PCONC(:,27)+TPK%K071(:)*PCONC(:,26)+T&
&PK%KTC11(:)+TPK%KTR11(:)
!
!PPROD(HO2) = +0.65*K006*<HNO4>+K011*<HCHO>+K011*<HCHO>+K012*<ALD>+K013*<OP1>+0
!.96205*K014*<OP2>+0.75830*K016*<CARBO>+K017*<ONIT>+K023*<O3>*<OH>+K026*<H2O2>*
!<OH>+K034*<OH>*<NO3>+K037*<HNO4>+K051*<OH>*<H2>+K052*<OH>*<SO2>+K053*<CO>*<OH>
!+0.28*K054*<BIO>*<O3P>+0.12793*K058*<ALKA>*<OH>+0.10318*K061*<ARO>*<OH>+K062*<
!HCHO>*<OH>+0.51208*K065*<CARBO>*<OH>+K066*<ORA1>*<OH>+0.02915*K069*<OP2>*<OH>+
!0.28107*K070*<PAN>*<OH>+K072*<HCHO>*<NO3>+0.63217*K074*<CARBO>*<NO3>+0.23451*K
!079*<ALKE>*<O3>+0.30000*K080*<BIO>*<O3>+0.28441*K081*<CARBO>*<O3>+0.08*K082*<P
!AN>*<O3>+0.02*K086*<ADD>*<O2>+K090*<MO2>*<NO>+0.74265*K091*<ALKAP>*<NO>+K092*<
!ALKEP>*<NO>+0.84700*K093*<BIOP>*<NO>+0.95115*K094*<AROP>*<NO>+0.12334*K095*<CA
!RBOP>*<NO>+0.18401*K096*<OLN>*<NO>+0.66*K104*<MO2>*<MO2>+0.98383*K105*<ALKAP>*
!<MO2>+K106*<ALKEP>*<MO2>+1.00000*K107*<BIOP>*<MO2>+1.02767*K108*<AROP>*<MO2>+0
!.82998*K109*<CARBOP>*<MO2>+0.67560*K110*<OLN>*<MO2>+0.48079*K111*<ALKAP>*<CARB
!OP>+0.50078*K112*<ALKEP>*<CARBOP>+0.50600*K113*<BIOP>*<CARBOP>+K114*<AROP>*<CA
!RBOP>+0.07566*K115*<CARBOP>*<CARBOP>+0.17599*K116*<OLN>*<CARBOP>+K117*<OLN>*<O
!LN>+0.00000*K118*<OLN>*<OLN>+K119*<MO2>*<NO3>+0.81290*K120*<ALKAP>*<NO3>+K121*
!<ALKEP>*<NO3>+K122*<BIOP>*<NO3>+K123*<AROP>*<NO3>+0.04915*K124*<CARBOP>*<NO3>+
!0.25928*K125*<OLN>*<NO3>+K127*<XO2>*<MO2>+KTC32*<WC_HO2>+KTR32*<WR_HO2>
 PPROD(:,15) = +0.65*TPK%K006(:)*PCONC(:,9)+TPK%K011(:)*PCONC(:,22)+TPK%K011(:)&
&*PCONC(:,22)+TPK%K012(:)*PCONC(:,23)+TPK%K013(:)*PCONC(:,28)+0.96205*TPK%K014(&
&:)*PCONC(:,29)+0.75830*TPK%K016(:)*PCONC(:,25)+TPK%K017(:)*PCONC(:,26)+TPK%K02&
&3(:)*PCONC(:,1)*PCONC(:,14)+TPK%K026(:)*PCONC(:,2)*PCONC(:,14)+TPK%K034(:)*PCO&
&NC(:,14)*PCONC(:,5)+TPK%K037(:)*PCONC(:,9)+TPK%K051(:)*PCONC(:,14)*TPK%H2(:)+T&
&PK%K052(:)*PCONC(:,14)*PCONC(:,11)+TPK%K053(:)*PCONC(:,13)*PCONC(:,14)+0.28*TP&
&K%K054(:)*PCONC(:,20)*TPK%O3P(:)+0.12793*TPK%K058(:)*PCONC(:,18)*PCONC(:,14)+0&
&.10318*TPK%K061(:)*PCONC(:,21)*PCONC(:,14)+TPK%K062(:)*PCONC(:,22)*PCONC(:,14)&
&+0.51208*TPK%K065(:)*PCONC(:,25)*PCONC(:,14)+TPK%K066(:)*PCONC(:,30)*PCONC(:,1&
&4)+0.02915*TPK%K069(:)*PCONC(:,29)*PCONC(:,14)+0.28107*TPK%K070(:)*PCONC(:,27)&
&*PCONC(:,14)+TPK%K072(:)*PCONC(:,22)*PCONC(:,5)+0.63217*TPK%K074(:)*PCONC(:,25&
&)*PCONC(:,5)+0.23451*TPK%K079(:)*PCONC(:,19)*PCONC(:,1)+0.30000*TPK%K080(:)*PC&
&ONC(:,20)*PCONC(:,1)+0.28441*TPK%K081(:)*PCONC(:,25)*PCONC(:,1)+0.08*TPK%K082(&
&:)*PCONC(:,27)*PCONC(:,1)+0.02*TPK%K086(:)*PCONC(:,37)*TPK%O2(:)+TPK%K090(:)*P&
&CONC(:,32)*PCONC(:,3)+0.74265*TPK%K091(:)*PCONC(:,33)*PCONC(:,3)+TPK%K092(:)*P&
&CONC(:,34)*PCONC(:,3)+0.84700*TPK%K093(:)*PCONC(:,35)*PCONC(:,3)+0.95115*TPK%K&
&094(:)*PCONC(:,38)*PCONC(:,3)+0.12334*TPK%K095(:)*PCONC(:,39)*PCONC(:,3)+0.184&
&01*TPK%K096(:)*PCONC(:,40)*PCONC(:,3)+0.66*TPK%K104(:)*PCONC(:,32)*PCONC(:,32)&
&+0.98383*TPK%K105(:)*PCONC(:,33)*PCONC(:,32)+TPK%K106(:)*PCONC(:,34)*PCONC(:,3&
&2)+1.00000*TPK%K107(:)*PCONC(:,35)*PCONC(:,32)+1.02767*TPK%K108(:)*PCONC(:,38)&
&*PCONC(:,32)+0.82998*TPK%K109(:)*PCONC(:,39)*PCONC(:,32)+0.67560*TPK%K110(:)*P&
&CONC(:,40)*PCONC(:,32)+0.48079*TPK%K111(:)*PCONC(:,33)*PCONC(:,39)+0.50078*TPK&
&%K112(:)*PCONC(:,34)*PCONC(:,39)+0.50600*TPK%K113(:)*PCONC(:,35)*PCONC(:,39)+T&
&PK%K114(:)*PCONC(:,38)*PCONC(:,39)+0.07566*TPK%K115(:)*PCONC(:,39)*PCONC(:,39)&
&+0.17599*TPK%K116(:)*PCONC(:,40)*PCONC(:,39)+TPK%K117(:)*PCONC(:,40)*PCONC(:,4&
&0)+0.00000*TPK%K118(:)*PCONC(:,40)*PCONC(:,40)+TPK%K119(:)*PCONC(:,32)*PCONC(:&
&,5)+0.81290*TPK%K120(:)*PCONC(:,33)*PCONC(:,5)+TPK%K121(:)*PCONC(:,34)*PCONC(:&
&,5)+TPK%K122(:)*PCONC(:,35)*PCONC(:,5)+TPK%K123(:)*PCONC(:,38)*PCONC(:,5)+0.04&
&915*TPK%K124(:)*PCONC(:,39)*PCONC(:,5)+0.25928*TPK%K125(:)*PCONC(:,40)*PCONC(:&
&,5)+TPK%K127(:)*PCONC(:,41)*PCONC(:,32)+TPK%KTC32(:)*PCONC(:,53)+TPK%KTR32(:)*&
&PCONC(:,78)
!PLOSS(HO2) = +K024*<O3>+K025*<OH>+K027*<HO2>+K027*<HO2>+K028*<HO2>*<H2O>+K028*
!<HO2>*<H2O>+K035*<NO>+K036*<NO2>+K038*<NO3>+K084*<PHO>+K097*<MO2>+K098*<ALKAP>
!+K099*<ALKEP>+K0100*<BIOP>+K0101*<AROP>+K0102*<CARBOP>+K103*<OLN>+K126*<XO2>+K
!TC12+KTR12
 PLOSS(:,15) = +TPK%K024(:)*PCONC(:,1)+TPK%K025(:)*PCONC(:,14)+TPK%K027(:)*PCON&
&C(:,15)+TPK%K027(:)*PCONC(:,15)+TPK%K028(:)*PCONC(:,15)*TPK%H2O(:)+TPK%K028(:)&
&*PCONC(:,15)*TPK%H2O(:)+TPK%K035(:)*PCONC(:,3)+TPK%K036(:)*PCONC(:,4)+TPK%K038&
&(:)*PCONC(:,5)+TPK%K084(:)*PCONC(:,36)+TPK%K097(:)*PCONC(:,32)+TPK%K098(:)*PCO&
&NC(:,33)+TPK%K099(:)*PCONC(:,34)+TPK%K0100(:)*PCONC(:,35)+TPK%K0101(:)*PCONC(:&
&,38)+TPK%K0102(:)*PCONC(:,39)+TPK%K103(:)*PCONC(:,40)+TPK%K126(:)*PCONC(:,41)+&
&TPK%KTC12(:)+TPK%KTR12(:)
!
!PPROD(CH4) = +0.04300*K079*<ALKE>*<O3>
 PPROD(:,16) = +0.04300*TPK%K079(:)*PCONC(:,19)*PCONC(:,1)
!PLOSS(CH4) = +K056*<OH>
 PLOSS(:,16) = +TPK%K056(:)*PCONC(:,14)
!
!PPROD(ETH) = +0.03196*K079*<ALKE>*<O3>
 PPROD(:,17) = +0.03196*TPK%K079(:)*PCONC(:,19)*PCONC(:,1)
!PLOSS(ETH) = +K057*<OH>
 PLOSS(:,17) = +TPK%K057(:)*PCONC(:,14)
!
!PPROD(ALKA) = 0.0
 PPROD(:,18) = 0.0
!PLOSS(ALKA) = +K058*<OH>
 PLOSS(:,18) = +TPK%K058(:)*PCONC(:,14)
!
!PPROD(ALKE) = +0.91868*K054*<BIO>*<O3P>+0.00000*K079*<ALKE>*<O3>+0.37388*K080*
!<BIO>*<O3>+0.37815*K093*<BIOP>*<NO>+0.48074*K107*<BIOP>*<MO2>+0.24463*K113*<BI
!OP>*<CARBOP>+0.42729*K122*<BIOP>*<NO3>
 PPROD(:,19) = +0.91868*TPK%K054(:)*PCONC(:,20)*TPK%O3P(:)+0.00000*TPK%K079(:)*&
&PCONC(:,19)*PCONC(:,1)+0.37388*TPK%K080(:)*PCONC(:,20)*PCONC(:,1)+0.37815*TPK%&
&K093(:)*PCONC(:,35)*PCONC(:,3)+0.48074*TPK%K107(:)*PCONC(:,35)*PCONC(:,32)+0.2&
&4463*TPK%K113(:)*PCONC(:,35)*PCONC(:,39)+0.42729*TPK%K122(:)*PCONC(:,35)*PCONC&
&(:,5)
!PLOSS(ALKE) = +K059*<OH>+K076*<NO3>+K079*<O3>
 PLOSS(:,19) = +TPK%K059(:)*PCONC(:,14)+TPK%K076(:)*PCONC(:,5)+TPK%K079(:)*PCON&
&C(:,1)
!
!PPROD(BIO) = 0.0
 PPROD(:,20) = 0.0
!PLOSS(BIO) = +K054*<O3P>+K060*<OH>+K077*<NO3>+K080*<O3>
 PLOSS(:,20) = +TPK%K054(:)*TPK%O3P(:)+TPK%K060(:)*PCONC(:,14)+TPK%K077(:)*PCON&
&C(:,5)+TPK%K080(:)*PCONC(:,1)
!
RETURN
END SUBROUTINE SUB1
!
SUBROUTINE SUB2
!
!Indices 21 a 30
!
!
!PPROD(ARO) = +0.10670*K083*<PHO>*<NO2>+1.06698*K084*<PHO>*<HO2>+K085*<ADD>*<NO
!2>+0.02*K086*<ADD>*<O2>+K087*<ADD>*<O3>
 PPROD(:,21) = +0.10670*TPK%K083(:)*PCONC(:,36)*PCONC(:,4)+1.06698*TPK%K084(:)*&
&PCONC(:,36)*PCONC(:,15)+TPK%K085(:)*PCONC(:,37)*PCONC(:,4)+0.02*TPK%K086(:)*PC&
&ONC(:,37)*TPK%O2(:)+TPK%K087(:)*PCONC(:,37)*PCONC(:,1)
!PLOSS(ARO) = +K061*<OH>+K075*<NO3>
 PLOSS(:,21) = +TPK%K061(:)*PCONC(:,14)+TPK%K075(:)*PCONC(:,5)
!
!PPROD(HCHO) = +K013*<OP1>+0.06517*K016*<CARBO>+0.05*K054*<BIO>*<O3P>+0.00140*K
!058*<ALKA>*<OH>+0.00000*K065*<CARBO>*<OH>+0.35*K068*<OP1>*<OH>+0.02915*K069*<O
!P2>*<OH>+0.57839*K070*<PAN>*<OH>+0.40*K078*<PAN>*<NO3>+0.48290*K079*<ALKE>*<O3
!>+0.90000*K080*<BIO>*<O3>+0.00000*K081*<CARBO>*<O3>+0.70*K082*<PAN>*<O3>+K090*
!<MO2>*<NO>+0.03002*K091*<ALKAP>*<NO>+1.39870*K092*<ALKEP>*<NO>+0.60600*K093*<B
!IOP>*<NO>+0.05848*K095*<CARBOP>*<NO>+0.23419*K096*<OLN>*<NO>+1.33*K104*<MO2>*<
!MO2>+0.80556*K105*<ALKAP>*<MO2>+1.42894*K106*<ALKEP>*<MO2>+1.09000*K107*<BIOP>
!*<MO2>+K108*<AROP>*<MO2>+0.95723*K109*<CARBOP>*<MO2>+0.88625*K110*<OLN>*<MO2>+
!0.07600*K111*<ALKAP>*<CARBOP>+0.68192*K112*<ALKEP>*<CARBOP>+0.34000*K113*<BIOP
!>*<CARBOP>+0.03432*K115*<CARBOP>*<CARBOP>+0.13414*K116*<OLN>*<CARBOP>+0.00000*
!K118*<OLN>*<OLN>+K119*<MO2>*<NO3>+0.03142*K120*<ALKAP>*<NO3>+1.40909*K121*<ALK
!EP>*<NO3>+0.68600*K122*<BIOP>*<NO3>+0.03175*K124*<CARBOP>*<NO3>+0.20740*K125*<
!OLN>*<NO3>+K127*<XO2>*<MO2>+KTC36*<WC_HCHO>+KTR36*<WR_HCHO>
 PPROD(:,22) = +TPK%K013(:)*PCONC(:,28)+0.06517*TPK%K016(:)*PCONC(:,25)+0.05*TP&
&K%K054(:)*PCONC(:,20)*TPK%O3P(:)+0.00140*TPK%K058(:)*PCONC(:,18)*PCONC(:,14)+0&
&.00000*TPK%K065(:)*PCONC(:,25)*PCONC(:,14)+0.35*TPK%K068(:)*PCONC(:,28)*PCONC(&
&:,14)+0.02915*TPK%K069(:)*PCONC(:,29)*PCONC(:,14)+0.57839*TPK%K070(:)*PCONC(:,&
&27)*PCONC(:,14)+0.40*TPK%K078(:)*PCONC(:,27)*PCONC(:,5)+0.48290*TPK%K079(:)*PC&
&ONC(:,19)*PCONC(:,1)+0.90000*TPK%K080(:)*PCONC(:,20)*PCONC(:,1)+0.00000*TPK%K0&
&81(:)*PCONC(:,25)*PCONC(:,1)+0.70*TPK%K082(:)*PCONC(:,27)*PCONC(:,1)+TPK%K090(&
&:)*PCONC(:,32)*PCONC(:,3)+0.03002*TPK%K091(:)*PCONC(:,33)*PCONC(:,3)+1.39870*T&
&PK%K092(:)*PCONC(:,34)*PCONC(:,3)+0.60600*TPK%K093(:)*PCONC(:,35)*PCONC(:,3)+0&
&.05848*TPK%K095(:)*PCONC(:,39)*PCONC(:,3)+0.23419*TPK%K096(:)*PCONC(:,40)*PCON&
&C(:,3)+1.33*TPK%K104(:)*PCONC(:,32)*PCONC(:,32)+0.80556*TPK%K105(:)*PCONC(:,33&
&)*PCONC(:,32)+1.42894*TPK%K106(:)*PCONC(:,34)*PCONC(:,32)+1.09000*TPK%K107(:)*&
&PCONC(:,35)*PCONC(:,32)+TPK%K108(:)*PCONC(:,38)*PCONC(:,32)+0.95723*TPK%K109(:&
&)*PCONC(:,39)*PCONC(:,32)+0.88625*TPK%K110(:)*PCONC(:,40)*PCONC(:,32)+0.07600*&
&TPK%K111(:)*PCONC(:,33)*PCONC(:,39)+0.68192*TPK%K112(:)*PCONC(:,34)*PCONC(:,39&
&)+0.34000*TPK%K113(:)*PCONC(:,35)*PCONC(:,39)+0.03432*TPK%K115(:)*PCONC(:,39)*&
&PCONC(:,39)+0.13414*TPK%K116(:)*PCONC(:,40)*PCONC(:,39)+0.00000*TPK%K118(:)*PC&
&ONC(:,40)*PCONC(:,40)+TPK%K119(:)*PCONC(:,32)*PCONC(:,5)+0.03142*TPK%K120(:)*P&
&CONC(:,33)*PCONC(:,5)+1.40909*TPK%K121(:)*PCONC(:,34)*PCONC(:,5)+0.68600*TPK%K&
&122(:)*PCONC(:,35)*PCONC(:,5)+0.03175*TPK%K124(:)*PCONC(:,39)*PCONC(:,5)+0.207&
&40*TPK%K125(:)*PCONC(:,40)*PCONC(:,5)+TPK%K127(:)*PCONC(:,41)*PCONC(:,32)+TPK%&
&KTC36(:)*PCONC(:,57)+TPK%KTR36(:)*PCONC(:,82)
!PLOSS(HCHO) = +K010+K011+K062*<OH>+K072*<NO3>+KTC16+KTR16
 PLOSS(:,22) = +TPK%K010(:)+TPK%K011(:)+TPK%K062(:)*PCONC(:,14)+TPK%K072(:)*PCO&
&NC(:,5)+TPK%KTC16(:)+TPK%KTR16(:)
!
!PPROD(ALD) = +0.96205*K014*<OP2>+0.20*K017*<ONIT>+K055*<CARBO>*<O3P>+0.08173*K
!058*<ALKA>*<OH>+0.06253*K065*<CARBO>*<OH>+0.07335*K069*<OP2>*<OH>+0.05265*K074
!*<CARBO>*<NO3>+0.51468*K079*<ALKE>*<O3>+0.00000*K080*<BIO>*<O3>+0.15692*K081*<
!CARBO>*<O3>+0.33144*K091*<ALKAP>*<NO>+0.42125*K092*<ALKEP>*<NO>+0.00000*K093*<
!BIOP>*<NO>+0.07368*K095*<CARBOP>*<NO>+1.01182*K096*<OLN>*<NO>+0.56070*K105*<AL
!KAP>*<MO2>+0.46413*K106*<ALKEP>*<MO2>+0.00000*K107*<BIOP>*<MO2>+0.08295*K109*<
!CARBOP>*<MO2>+0.41524*K110*<OLN>*<MO2>+0.71461*K111*<ALKAP>*<CARBOP>+0.68374*K
!112*<ALKEP>*<CARBOP>+0.00000*K113*<BIOP>*<CARBOP>+0.06969*K115*<CARBOP>*<CARBO
!P>+0.42122*K116*<OLN>*<CARBOP>+0.00000*K118*<OLN>*<OLN>+0.33743*K120*<ALKAP>*<
!NO3>+0.43039*K121*<ALKEP>*<NO3>+0.00000*K122*<BIOP>*<NO3>+0.02936*K124*<CARBOP
!>*<NO3>+0.91850*K125*<OLN>*<NO3>
 PPROD(:,23) = +0.96205*TPK%K014(:)*PCONC(:,29)+0.20*TPK%K017(:)*PCONC(:,26)+TP&
&K%K055(:)*PCONC(:,25)*TPK%O3P(:)+0.08173*TPK%K058(:)*PCONC(:,18)*PCONC(:,14)+0&
&.06253*TPK%K065(:)*PCONC(:,25)*PCONC(:,14)+0.07335*TPK%K069(:)*PCONC(:,29)*PCO&
&NC(:,14)+0.05265*TPK%K074(:)*PCONC(:,25)*PCONC(:,5)+0.51468*TPK%K079(:)*PCONC(&
&:,19)*PCONC(:,1)+0.00000*TPK%K080(:)*PCONC(:,20)*PCONC(:,1)+0.15692*TPK%K081(:&
&)*PCONC(:,25)*PCONC(:,1)+0.33144*TPK%K091(:)*PCONC(:,33)*PCONC(:,3)+0.42125*TP&
&K%K092(:)*PCONC(:,34)*PCONC(:,3)+0.00000*TPK%K093(:)*PCONC(:,35)*PCONC(:,3)+0.&
&07368*TPK%K095(:)*PCONC(:,39)*PCONC(:,3)+1.01182*TPK%K096(:)*PCONC(:,40)*PCONC&
&(:,3)+0.56070*TPK%K105(:)*PCONC(:,33)*PCONC(:,32)+0.46413*TPK%K106(:)*PCONC(:,&
&34)*PCONC(:,32)+0.00000*TPK%K107(:)*PCONC(:,35)*PCONC(:,32)+0.08295*TPK%K109(:&
&)*PCONC(:,39)*PCONC(:,32)+0.41524*TPK%K110(:)*PCONC(:,40)*PCONC(:,32)+0.71461*&
&TPK%K111(:)*PCONC(:,33)*PCONC(:,39)+0.68374*TPK%K112(:)*PCONC(:,34)*PCONC(:,39&
&)+0.00000*TPK%K113(:)*PCONC(:,35)*PCONC(:,39)+0.06969*TPK%K115(:)*PCONC(:,39)*&
&PCONC(:,39)+0.42122*TPK%K116(:)*PCONC(:,40)*PCONC(:,39)+0.00000*TPK%K118(:)*PC&
&ONC(:,40)*PCONC(:,40)+0.33743*TPK%K120(:)*PCONC(:,33)*PCONC(:,5)+0.43039*TPK%K&
&121(:)*PCONC(:,34)*PCONC(:,5)+0.00000*TPK%K122(:)*PCONC(:,35)*PCONC(:,5)+0.029&
&36*TPK%K124(:)*PCONC(:,39)*PCONC(:,5)+0.91850*TPK%K125(:)*PCONC(:,40)*PCONC(:,&
&5)
!PLOSS(ALD) = +K012+K063*<OH>+K073*<NO3>
 PLOSS(:,23) = +TPK%K012(:)+TPK%K063(:)*PCONC(:,14)+TPK%K073(:)*PCONC(:,5)
!
!PPROD(KET) = +0.80*K017*<ONIT>+0.03498*K058*<ALKA>*<OH>+0.00853*K065*<CARBO>*<
!OH>+0.37591*K069*<OP2>*<OH>+0.00632*K074*<CARBO>*<NO3>+0.07377*K079*<ALKE>*<O3
!>+0.00000*K080*<BIO>*<O3>+0.54531*K091*<ALKAP>*<NO>+0.05220*K092*<ALKEP>*<NO>+
!0.00000*K093*<BIOP>*<NO>+0.37862*K096*<OLN>*<NO>+0.09673*K105*<ALKAP>*<MO2>+0.
!03814*K106*<ALKEP>*<MO2>+0.00000*K107*<BIOP>*<MO2>+0.09667*K110*<OLN>*<MO2>+0.
!18819*K111*<ALKAP>*<CARBOP>+0.06579*K112*<ALKEP>*<CARBOP>+0.00000*K113*<BIOP>*
!<CARBOP>+0.02190*K115*<CARBOP>*<CARBOP>+0.10822*K116*<OLN>*<CARBOP>+0.00000*K1
!18*<OLN>*<OLN>+0.62978*K120*<ALKAP>*<NO3>+0.02051*K121*<ALKEP>*<NO3>+0.00000*K
!122*<BIOP>*<NO3>+0.34740*K125*<OLN>*<NO3>
 PPROD(:,24) = +0.80*TPK%K017(:)*PCONC(:,26)+0.03498*TPK%K058(:)*PCONC(:,18)*PC&
&ONC(:,14)+0.00853*TPK%K065(:)*PCONC(:,25)*PCONC(:,14)+0.37591*TPK%K069(:)*PCON&
&C(:,29)*PCONC(:,14)+0.00632*TPK%K074(:)*PCONC(:,25)*PCONC(:,5)+0.07377*TPK%K07&
&9(:)*PCONC(:,19)*PCONC(:,1)+0.00000*TPK%K080(:)*PCONC(:,20)*PCONC(:,1)+0.54531&
&*TPK%K091(:)*PCONC(:,33)*PCONC(:,3)+0.05220*TPK%K092(:)*PCONC(:,34)*PCONC(:,3)&
&+0.00000*TPK%K093(:)*PCONC(:,35)*PCONC(:,3)+0.37862*TPK%K096(:)*PCONC(:,40)*PC&
&ONC(:,3)+0.09673*TPK%K105(:)*PCONC(:,33)*PCONC(:,32)+0.03814*TPK%K106(:)*PCONC&
&(:,34)*PCONC(:,32)+0.00000*TPK%K107(:)*PCONC(:,35)*PCONC(:,32)+0.09667*TPK%K11&
&0(:)*PCONC(:,40)*PCONC(:,32)+0.18819*TPK%K111(:)*PCONC(:,33)*PCONC(:,39)+0.065&
&79*TPK%K112(:)*PCONC(:,34)*PCONC(:,39)+0.00000*TPK%K113(:)*PCONC(:,35)*PCONC(:&
&,39)+0.02190*TPK%K115(:)*PCONC(:,39)*PCONC(:,39)+0.10822*TPK%K116(:)*PCONC(:,4&
&0)*PCONC(:,39)+0.00000*TPK%K118(:)*PCONC(:,40)*PCONC(:,40)+0.62978*TPK%K120(:)&
&*PCONC(:,33)*PCONC(:,5)+0.02051*TPK%K121(:)*PCONC(:,34)*PCONC(:,5)+0.00000*TPK&
&%K122(:)*PCONC(:,35)*PCONC(:,5)+0.34740*TPK%K125(:)*PCONC(:,40)*PCONC(:,5)
!PLOSS(KET) = +K015+K064*<OH>
 PLOSS(:,24) = +TPK%K015(:)+TPK%K064(:)*PCONC(:,14)
!
!PPROD(CARBO) = +0.13255*K054*<BIO>*<O3P>+0.00835*K058*<ALKA>*<OH>+0.16919*K065
!*<CARBO>*<OH>+0.21863*K070*<PAN>*<OH>+0.10530*K074*<CARBO>*<NO3>+0.00000*K076*
!<ALKE>*<NO3>+0.91741*K077*<BIO>*<NO3>+0.00000*K079*<ALKE>*<O3>+0.39754*K080*<B
!IO>*<O3>+1.07583*K081*<CARBO>*<O3>+0.03407*K091*<ALKAP>*<NO>+0.45463*K093*<BIO
!P>*<NO>+2.06993*K094*<AROP>*<NO>+0.08670*K095*<CARBOP>*<NO>+0.07976*K105*<ALKA
!P>*<MO2>+0.56064*K107*<BIOP>*<MO2>+1.99461*K108*<AROP>*<MO2>+0.15387*K109*<CAR
!BOP>*<MO2>+0.06954*K111*<ALKAP>*<CARBOP>+0.78591*K113*<BIOP>*<CARBOP>+1.99455*
!K114*<AROP>*<CARBOP>+0.10777*K115*<CARBOP>*<CARBOP>+0.03531*K120*<ALKAP>*<NO3>
!+0.61160*K122*<BIOP>*<NO3>+2.81904*K123*<AROP>*<NO3>+0.03455*K124*<CARBOP>*<NO
!3>
 PPROD(:,25) = +0.13255*TPK%K054(:)*PCONC(:,20)*TPK%O3P(:)+0.00835*TPK%K058(:)*&
&PCONC(:,18)*PCONC(:,14)+0.16919*TPK%K065(:)*PCONC(:,25)*PCONC(:,14)+0.21863*TP&
&K%K070(:)*PCONC(:,27)*PCONC(:,14)+0.10530*TPK%K074(:)*PCONC(:,25)*PCONC(:,5)+0&
&.00000*TPK%K076(:)*PCONC(:,19)*PCONC(:,5)+0.91741*TPK%K077(:)*PCONC(:,20)*PCON&
&C(:,5)+0.00000*TPK%K079(:)*PCONC(:,19)*PCONC(:,1)+0.39754*TPK%K080(:)*PCONC(:,&
&20)*PCONC(:,1)+1.07583*TPK%K081(:)*PCONC(:,25)*PCONC(:,1)+0.03407*TPK%K091(:)*&
&PCONC(:,33)*PCONC(:,3)+0.45463*TPK%K093(:)*PCONC(:,35)*PCONC(:,3)+2.06993*TPK%&
&K094(:)*PCONC(:,38)*PCONC(:,3)+0.08670*TPK%K095(:)*PCONC(:,39)*PCONC(:,3)+0.07&
&976*TPK%K105(:)*PCONC(:,33)*PCONC(:,32)+0.56064*TPK%K107(:)*PCONC(:,35)*PCONC(&
&:,32)+1.99461*TPK%K108(:)*PCONC(:,38)*PCONC(:,32)+0.15387*TPK%K109(:)*PCONC(:,&
&39)*PCONC(:,32)+0.06954*TPK%K111(:)*PCONC(:,33)*PCONC(:,39)+0.78591*TPK%K113(:&
&)*PCONC(:,35)*PCONC(:,39)+1.99455*TPK%K114(:)*PCONC(:,38)*PCONC(:,39)+0.10777*&
&TPK%K115(:)*PCONC(:,39)*PCONC(:,39)+0.03531*TPK%K120(:)*PCONC(:,33)*PCONC(:,5)&
&+0.61160*TPK%K122(:)*PCONC(:,35)*PCONC(:,5)+2.81904*TPK%K123(:)*PCONC(:,38)*PC&
&ONC(:,5)+0.03455*TPK%K124(:)*PCONC(:,39)*PCONC(:,5)
!PLOSS(CARBO) = +K016+K055*<O3P>+K065*<OH>+K074*<NO3>+K081*<O3>
 PLOSS(:,25) = +TPK%K016(:)+TPK%K055(:)*TPK%O3P(:)+TPK%K065(:)*PCONC(:,14)+TPK%&
&K074(:)*PCONC(:,5)+TPK%K081(:)*PCONC(:,1)
!
!PPROD(ONIT) = +0.60*K078*<PAN>*<NO3>+K083*<PHO>*<NO2>+0.08459*K091*<ALKAP>*<NO
!>+0.15300*K093*<BIOP>*<NO>+0.04885*K094*<AROP>*<NO>+0.18401*K096*<OLN>*<NO>+K1
!03*<OLN>*<HO2>+0.67560*K110*<OLN>*<MO2>+0.66562*K116*<OLN>*<CARBOP>+2.00*K117*
!<OLN>*<OLN>+0.00000*K118*<OLN>*<OLN>+0.25928*K125*<OLN>*<NO3>
 PPROD(:,26) = +0.60*TPK%K078(:)*PCONC(:,27)*PCONC(:,5)+TPK%K083(:)*PCONC(:,36)&
&*PCONC(:,4)+0.08459*TPK%K091(:)*PCONC(:,33)*PCONC(:,3)+0.15300*TPK%K093(:)*PCO&
&NC(:,35)*PCONC(:,3)+0.04885*TPK%K094(:)*PCONC(:,38)*PCONC(:,3)+0.18401*TPK%K09&
&6(:)*PCONC(:,40)*PCONC(:,3)+TPK%K103(:)*PCONC(:,40)*PCONC(:,15)+0.67560*TPK%K1&
&10(:)*PCONC(:,40)*PCONC(:,32)+0.66562*TPK%K116(:)*PCONC(:,40)*PCONC(:,39)+2.00&
&*TPK%K117(:)*PCONC(:,40)*PCONC(:,40)+0.00000*TPK%K118(:)*PCONC(:,40)*PCONC(:,4&
&0)+0.25928*TPK%K125(:)*PCONC(:,40)*PCONC(:,5)
!PLOSS(ONIT) = +K017+K071*<OH>
 PLOSS(:,26) = +TPK%K017(:)+TPK%K071(:)*PCONC(:,14)
!
!PPROD(PAN) = +0.28107*K070*<PAN>*<OH>+0.40000*K078*<PAN>*<NO3>+0.30000*K082*<P
!AN>*<O3>+1.00000*K088*<CARBOP>*<NO2>
 PPROD(:,27) = +0.28107*TPK%K070(:)*PCONC(:,27)*PCONC(:,14)+0.40000*TPK%K078(:)&
&*PCONC(:,27)*PCONC(:,5)+0.30000*TPK%K082(:)*PCONC(:,27)*PCONC(:,1)+1.00000*TPK&
&%K088(:)*PCONC(:,39)*PCONC(:,4)
!PLOSS(PAN) = +K070*<OH>+K078*<NO3>+K082*<O3>+K089
 PLOSS(:,27) = +TPK%K070(:)*PCONC(:,14)+TPK%K078(:)*PCONC(:,5)+TPK%K082(:)*PCON&
&C(:,1)+TPK%K089(:)
!
!PPROD(OP1) = +K097*<MO2>*<HO2>+KTC40*<WC_OP1>+KTR40*<WR_OP1>
 PPROD(:,28) = +TPK%K097(:)*PCONC(:,32)*PCONC(:,15)+TPK%KTC40(:)*PCONC(:,61)+TP&
&K%KTR40(:)*PCONC(:,86)
!PLOSS(OP1) = +K013+K068*<OH>+KTC20+KTR20
 PLOSS(:,28) = +TPK%K013(:)+TPK%K068(:)*PCONC(:,14)+TPK%KTC20(:)+TPK%KTR20(:)
!
!PPROD(OP2) = +0.10149*K081*<CARBO>*<O3>+1.00524*K098*<ALKAP>*<HO2>+1.00524*K09
!9*<ALKEP>*<HO2>+1.00524*K0100*<BIOP>*<HO2>+1.00524*K0101*<AROP>*<HO2>+0.80904*
!K0102*<CARBOP>*<HO2>+1.00524*K126*<XO2>*<HO2>
 PPROD(:,29) = +0.10149*TPK%K081(:)*PCONC(:,25)*PCONC(:,1)+1.00524*TPK%K098(:)*&
&PCONC(:,33)*PCONC(:,15)+1.00524*TPK%K099(:)*PCONC(:,34)*PCONC(:,15)+1.00524*TP&
&K%K0100(:)*PCONC(:,35)*PCONC(:,15)+1.00524*TPK%K0101(:)*PCONC(:,38)*PCONC(:,15&
&)+0.80904*TPK%K0102(:)*PCONC(:,39)*PCONC(:,15)+1.00524*TPK%K126(:)*PCONC(:,41)&
&*PCONC(:,15)
!PLOSS(OP2) = +K014+K069*<OH>
 PLOSS(:,29) = +TPK%K014(:)+TPK%K069(:)*PCONC(:,14)
!
!PPROD(ORA1) = +0.00878*K058*<ALKA>*<OH>+0.15343*K079*<ALKE>*<O3>+0.15000*K080*
!<BIO>*<O3>+0.10788*K081*<CARBO>*<O3>+0.11*K082*<PAN>*<O3>+KTC37*<WC_ORA1>+KTR3
!7*<WR_ORA1>
 PPROD(:,30) = +0.00878*TPK%K058(:)*PCONC(:,18)*PCONC(:,14)+0.15343*TPK%K079(:)&
&*PCONC(:,19)*PCONC(:,1)+0.15000*TPK%K080(:)*PCONC(:,20)*PCONC(:,1)+0.10788*TPK&
&%K081(:)*PCONC(:,25)*PCONC(:,1)+0.11*TPK%K082(:)*PCONC(:,27)*PCONC(:,1)+TPK%KT&
&C37(:)*PCONC(:,58)+TPK%KTR37(:)*PCONC(:,83)
!PLOSS(ORA1) = +K066*<OH>+KTC17+KTR17
 PLOSS(:,30) = +TPK%K066(:)*PCONC(:,14)+TPK%KTC17(:)+TPK%KTR17(:)
!
RETURN
END SUBROUTINE SUB2
!
SUBROUTINE SUB3
!
!Indices 31 a 40
!
!
!PPROD(ORA2) = +0.08143*K079*<ALKE>*<O3>+0.00000*K080*<BIO>*<O3>+0.20595*K081*<
!CARBO>*<O3>+0.17307*K0102*<CARBOP>*<HO2>+0.13684*K109*<CARBOP>*<MO2>+0.49810*K
!111*<ALKAP>*<CARBOP>+0.49922*K112*<ALKEP>*<CARBOP>+0.49400*K113*<BIOP>*<CARBOP
!>+0.09955*K115*<CARBOP>*<CARBOP>+0.48963*K116*<OLN>*<CARBOP>+KTC38*<WC_ORA2>+K
!TR38*<WR_ORA2>
 PPROD(:,31) = +0.08143*TPK%K079(:)*PCONC(:,19)*PCONC(:,1)+0.00000*TPK%K080(:)*&
&PCONC(:,20)*PCONC(:,1)+0.20595*TPK%K081(:)*PCONC(:,25)*PCONC(:,1)+0.17307*TPK%&
&K0102(:)*PCONC(:,39)*PCONC(:,15)+0.13684*TPK%K109(:)*PCONC(:,39)*PCONC(:,32)+0&
&.49810*TPK%K111(:)*PCONC(:,33)*PCONC(:,39)+0.49922*TPK%K112(:)*PCONC(:,34)*PCO&
&NC(:,39)+0.49400*TPK%K113(:)*PCONC(:,35)*PCONC(:,39)+0.09955*TPK%K115(:)*PCONC&
&(:,39)*PCONC(:,39)+0.48963*TPK%K116(:)*PCONC(:,40)*PCONC(:,39)+TPK%KTC38(:)*PC&
&ONC(:,59)+TPK%KTR38(:)*PCONC(:,84)
!PLOSS(ORA2) = +K067*<OH>+KTC18+KTR18
 PLOSS(:,31) = +TPK%K067(:)*PCONC(:,14)+TPK%KTC18(:)+TPK%KTR18(:)
!
!PPROD(MO2) = +K012*<ALD>+0.03795*K014*<OP2>+K056*<CH4>*<OH>+0.65*K068*<OP1>*<O
!H>+0.13966*K079*<ALKE>*<O3>+0.03000*K080*<BIO>*<O3>+0.09016*K091*<ALKAP>*<NO>+
!0.78134*K095*<CARBOP>*<NO>+0.01390*K105*<ALKAP>*<MO2>+0.56031*K109*<CARBOP>*<M
!O2>+0.51480*K111*<ALKAP>*<CARBOP>+0.50078*K112*<ALKEP>*<CARBOP>+0.50600*K113*<
!BIOP>*<CARBOP>+K114*<AROP>*<CARBOP>+1.66702*K115*<CARBOP>*<CARBOP>+0.51037*K11
!6*<OLN>*<CARBOP>+0.09731*K120*<ALKAP>*<NO3>+0.91910*K124*<CARBOP>*<NO3>+K128*<
!XO2>*<CARBOP>+KTC39*<WC_MO2>+KTR39*<WR_MO2>
 PPROD(:,32) = +TPK%K012(:)*PCONC(:,23)+0.03795*TPK%K014(:)*PCONC(:,29)+TPK%K05&
&6(:)*PCONC(:,16)*PCONC(:,14)+0.65*TPK%K068(:)*PCONC(:,28)*PCONC(:,14)+0.13966*&
&TPK%K079(:)*PCONC(:,19)*PCONC(:,1)+0.03000*TPK%K080(:)*PCONC(:,20)*PCONC(:,1)+&
&0.09016*TPK%K091(:)*PCONC(:,33)*PCONC(:,3)+0.78134*TPK%K095(:)*PCONC(:,39)*PCO&
&NC(:,3)+0.01390*TPK%K105(:)*PCONC(:,33)*PCONC(:,32)+0.56031*TPK%K109(:)*PCONC(&
&:,39)*PCONC(:,32)+0.51480*TPK%K111(:)*PCONC(:,33)*PCONC(:,39)+0.50078*TPK%K112&
&(:)*PCONC(:,34)*PCONC(:,39)+0.50600*TPK%K113(:)*PCONC(:,35)*PCONC(:,39)+TPK%K1&
&14(:)*PCONC(:,38)*PCONC(:,39)+1.66702*TPK%K115(:)*PCONC(:,39)*PCONC(:,39)+0.51&
&037*TPK%K116(:)*PCONC(:,40)*PCONC(:,39)+0.09731*TPK%K120(:)*PCONC(:,33)*PCONC(&
&:,5)+0.91910*TPK%K124(:)*PCONC(:,39)*PCONC(:,5)+TPK%K128(:)*PCONC(:,41)*PCONC(&
&:,39)+TPK%KTC39(:)*PCONC(:,60)+TPK%KTR39(:)*PCONC(:,85)
!PLOSS(MO2) = +K090*<NO>+K097*<HO2>+K104*<MO2>+K104*<MO2>+K105*<ALKAP>+K106*<AL
!KEP>+K107*<BIOP>+K108*<AROP>+K109*<CARBOP>+K110*<OLN>+K119*<NO3>+K127*<XO2>+KT
!C19+KTR19
 PLOSS(:,32) = +TPK%K090(:)*PCONC(:,3)+TPK%K097(:)*PCONC(:,15)+TPK%K104(:)*PCON&
&C(:,32)+TPK%K104(:)*PCONC(:,32)+TPK%K105(:)*PCONC(:,33)+TPK%K106(:)*PCONC(:,34&
&)+TPK%K107(:)*PCONC(:,35)+TPK%K108(:)*PCONC(:,38)+TPK%K109(:)*PCONC(:,39)+TPK%&
&K110(:)*PCONC(:,40)+TPK%K119(:)*PCONC(:,5)+TPK%K127(:)*PCONC(:,41)+TPK%KTC19(:&
&)+TPK%KTR19(:)
!
!PPROD(ALKAP) = +1.00000*K015*<KET>+K057*<ETH>*<OH>+0.87811*K058*<ALKA>*<OH>+0.
!40341*K069*<OP2>*<OH>+1.00000*K071*<ONIT>*<OH>+0.09815*K079*<ALKE>*<O3>+0.0000
!0*K080*<BIO>*<O3>+0.08187*K091*<ALKAP>*<NO>+0.00385*K105*<ALKAP>*<MO2>+0.00828
!*K111*<ALKAP>*<CARBOP>+0.08994*K120*<ALKAP>*<NO3>
 PPROD(:,33) = +1.00000*TPK%K015(:)*PCONC(:,24)+TPK%K057(:)*PCONC(:,17)*PCONC(:&
&,14)+0.87811*TPK%K058(:)*PCONC(:,18)*PCONC(:,14)+0.40341*TPK%K069(:)*PCONC(:,2&
&9)*PCONC(:,14)+1.00000*TPK%K071(:)*PCONC(:,26)*PCONC(:,14)+0.09815*TPK%K079(:)&
&*PCONC(:,19)*PCONC(:,1)+0.00000*TPK%K080(:)*PCONC(:,20)*PCONC(:,1)+0.08187*TPK&
&%K091(:)*PCONC(:,33)*PCONC(:,3)+0.00385*TPK%K105(:)*PCONC(:,33)*PCONC(:,32)+0.&
&00828*TPK%K111(:)*PCONC(:,33)*PCONC(:,39)+0.08994*TPK%K120(:)*PCONC(:,33)*PCON&
&C(:,5)
!PLOSS(ALKAP) = +K091*<NO>+K098*<HO2>+K105*<MO2>+K111*<CARBOP>+K120*<NO3>
 PLOSS(:,33) = +TPK%K091(:)*PCONC(:,3)+TPK%K098(:)*PCONC(:,15)+TPK%K105(:)*PCON&
&C(:,32)+TPK%K111(:)*PCONC(:,39)+TPK%K120(:)*PCONC(:,5)
!
!PPROD(ALKEP) = +1.02529*K059*<ALKE>*<OH>
 PPROD(:,34) = +1.02529*TPK%K059(:)*PCONC(:,19)*PCONC(:,14)
!PLOSS(ALKEP) = +K092*<NO>+K099*<HO2>+K106*<MO2>+K112*<CARBOP>+K121*<NO3>
 PLOSS(:,34) = +TPK%K092(:)*PCONC(:,3)+TPK%K099(:)*PCONC(:,15)+TPK%K106(:)*PCON&
&C(:,32)+TPK%K112(:)*PCONC(:,39)+TPK%K121(:)*PCONC(:,5)
!
!PPROD(BIOP) = +0.00000*K059*<ALKE>*<OH>+1.00000*K060*<BIO>*<OH>
 PPROD(:,35) = +0.00000*TPK%K059(:)*PCONC(:,19)*PCONC(:,14)+1.00000*TPK%K060(:)&
&*PCONC(:,20)*PCONC(:,14)
!PLOSS(BIOP) = +K093*<NO>+K0100*<HO2>+K107*<MO2>+K113*<CARBOP>+K122*<NO3>
 PLOSS(:,35) = +TPK%K093(:)*PCONC(:,3)+TPK%K0100(:)*PCONC(:,15)+TPK%K107(:)*PCO&
&NC(:,32)+TPK%K113(:)*PCONC(:,39)+TPK%K122(:)*PCONC(:,5)
!
!PPROD(PHO) = +0.00276*K061*<ARO>*<OH>+K075*<ARO>*<NO3>
 PPROD(:,36) = +0.00276*TPK%K061(:)*PCONC(:,21)*PCONC(:,14)+TPK%K075(:)*PCONC(:&
&,21)*PCONC(:,5)
!PLOSS(PHO) = +K083*<NO2>+K084*<HO2>
 PLOSS(:,36) = +TPK%K083(:)*PCONC(:,4)+TPK%K084(:)*PCONC(:,15)
!
!PPROD(ADD) = +0.93968*K061*<ARO>*<OH>
 PPROD(:,37) = +0.93968*TPK%K061(:)*PCONC(:,21)*PCONC(:,14)
!PLOSS(ADD) = +K085*<NO2>+K086*<O2>+K087*<O3>
 PLOSS(:,37) = +TPK%K085(:)*PCONC(:,4)+TPK%K086(:)*TPK%O2(:)+TPK%K087(:)*PCONC(&
&:,1)
!
!PPROD(AROP) = +0.98*K086*<ADD>*<O2>
 PPROD(:,38) = +0.98*TPK%K086(:)*PCONC(:,37)*TPK%O2(:)
!PLOSS(AROP) = +K094*<NO>+K0101*<HO2>+K108*<MO2>+K114*<CARBOP>+K123*<NO3>
 PLOSS(:,38) = +TPK%K094(:)*PCONC(:,3)+TPK%K0101(:)*PCONC(:,15)+TPK%K108(:)*PCO&
&NC(:,32)+TPK%K114(:)*PCONC(:,39)+TPK%K123(:)*PCONC(:,5)
!
!PPROD(CARBOP) = +1.00000*K015*<KET>+0.69622*K016*<CARBO>+1.00000*K063*<ALD>*<O
!H>+1.00000*K064*<KET>*<OH>+0.51419*K065*<CARBO>*<OH>+0.05413*K069*<OP2>*<OH>+1
!.00000*K073*<ALD>*<NO3>+0.38881*K074*<CARBO>*<NO3>+0.05705*K079*<ALKE>*<O3>+0.
!17000*K080*<BIO>*<O3>+0.27460*K081*<CARBO>*<O3>+0.70000*K082*<PAN>*<O3>+1.0000
!0*K089*<PAN>+0.09532*K095*<CARBOP>*<NO>+0.05954*K109*<CARBOP>*<MO2>+0.05821*K1
!15*<CARBOP>*<CARBOP>+0.03175*K124*<CARBOP>*<NO3>
 PPROD(:,39) = +1.00000*TPK%K015(:)*PCONC(:,24)+0.69622*TPK%K016(:)*PCONC(:,25)&
&+1.00000*TPK%K063(:)*PCONC(:,23)*PCONC(:,14)+1.00000*TPK%K064(:)*PCONC(:,24)*P&
&CONC(:,14)+0.51419*TPK%K065(:)*PCONC(:,25)*PCONC(:,14)+0.05413*TPK%K069(:)*PCO&
&NC(:,29)*PCONC(:,14)+1.00000*TPK%K073(:)*PCONC(:,23)*PCONC(:,5)+0.38881*TPK%K0&
&74(:)*PCONC(:,25)*PCONC(:,5)+0.05705*TPK%K079(:)*PCONC(:,19)*PCONC(:,1)+0.1700&
&0*TPK%K080(:)*PCONC(:,20)*PCONC(:,1)+0.27460*TPK%K081(:)*PCONC(:,25)*PCONC(:,1&
&)+0.70000*TPK%K082(:)*PCONC(:,27)*PCONC(:,1)+1.00000*TPK%K089(:)*PCONC(:,27)+0&
&.09532*TPK%K095(:)*PCONC(:,39)*PCONC(:,3)+0.05954*TPK%K109(:)*PCONC(:,39)*PCON&
&C(:,32)+0.05821*TPK%K115(:)*PCONC(:,39)*PCONC(:,39)+0.03175*TPK%K124(:)*PCONC(&
&:,39)*PCONC(:,5)
!PLOSS(CARBOP) = +K088*<NO2>+K095*<NO>+K0102*<HO2>+K109*<MO2>+K111*<ALKAP>+K112
!*<ALKEP>+K113*<BIOP>+K114*<AROP>+K115*<CARBOP>+K115*<CARBOP>+K116*<OLN>+K124*<
!NO3>+K128*<XO2>
 PLOSS(:,39) = +TPK%K088(:)*PCONC(:,4)+TPK%K095(:)*PCONC(:,3)+TPK%K0102(:)*PCON&
&C(:,15)+TPK%K109(:)*PCONC(:,32)+TPK%K111(:)*PCONC(:,33)+TPK%K112(:)*PCONC(:,34&
&)+TPK%K113(:)*PCONC(:,35)+TPK%K114(:)*PCONC(:,38)+TPK%K115(:)*PCONC(:,39)+TPK%&
&K115(:)*PCONC(:,39)+TPK%K116(:)*PCONC(:,40)+TPK%K124(:)*PCONC(:,5)+TPK%K128(:)&
&*PCONC(:,41)
!
!PPROD(OLN) = +0.00000*K074*<CARBO>*<NO3>+0.93768*K076*<ALKE>*<NO3>+1.00000*K07
!7*<BIO>*<NO3>
 PPROD(:,40) = +0.00000*TPK%K074(:)*PCONC(:,25)*PCONC(:,5)+0.93768*TPK%K076(:)*&
&PCONC(:,19)*PCONC(:,5)+1.00000*TPK%K077(:)*PCONC(:,20)*PCONC(:,5)
!PLOSS(OLN) = +K096*<NO>+K103*<HO2>+K110*<MO2>+K116*<CARBOP>+K117*<OLN>+K117*<O
!LN>+K118*<OLN>+K118*<OLN>+K125*<NO3>
 PLOSS(:,40) = +TPK%K096(:)*PCONC(:,3)+TPK%K103(:)*PCONC(:,15)+TPK%K110(:)*PCON&
&C(:,32)+TPK%K116(:)*PCONC(:,39)+TPK%K117(:)*PCONC(:,40)+TPK%K117(:)*PCONC(:,40&
&)+TPK%K118(:)*PCONC(:,40)+TPK%K118(:)*PCONC(:,40)+TPK%K125(:)*PCONC(:,5)
!
RETURN
END SUBROUTINE SUB3
!
SUBROUTINE SUB4
!
!Indices 41 a 50
!
!
!PPROD(XO2) = +0.15*K054*<BIO>*<O3P>+0.10318*K061*<ARO>*<OH>+0.10162*K065*<CARB
!O>*<OH>+0.09333*K069*<OP2>*<OH>+K070*<PAN>*<OH>+0.10530*K074*<CARBO>*<NO3>+K07
!8*<PAN>*<NO3>+0.00000*K079*<ALKE>*<O3>+0.13000*K080*<BIO>*<O3>+0.13007*K091*<A
!LKAP>*<NO>+0.02563*K095*<CARBOP>*<NO>+0.13370*K105*<ALKAP>*<MO2>+0.02212*K109*
!<CARBOP>*<MO2>+0.11306*K111*<ALKAP>*<CARBOP>+0.01593*K115*<CARBOP>*<CARBOP>+0.
!16271*K120*<ALKAP>*<NO3>+0.01021*K124*<CARBOP>*<NO3>
 PPROD(:,41) = +0.15*TPK%K054(:)*PCONC(:,20)*TPK%O3P(:)+0.10318*TPK%K061(:)*PCO&
&NC(:,21)*PCONC(:,14)+0.10162*TPK%K065(:)*PCONC(:,25)*PCONC(:,14)+0.09333*TPK%K&
&069(:)*PCONC(:,29)*PCONC(:,14)+TPK%K070(:)*PCONC(:,27)*PCONC(:,14)+0.10530*TPK&
&%K074(:)*PCONC(:,25)*PCONC(:,5)+TPK%K078(:)*PCONC(:,27)*PCONC(:,5)+0.00000*TPK&
&%K079(:)*PCONC(:,19)*PCONC(:,1)+0.13000*TPK%K080(:)*PCONC(:,20)*PCONC(:,1)+0.1&
&3007*TPK%K091(:)*PCONC(:,33)*PCONC(:,3)+0.02563*TPK%K095(:)*PCONC(:,39)*PCONC(&
&:,3)+0.13370*TPK%K105(:)*PCONC(:,33)*PCONC(:,32)+0.02212*TPK%K109(:)*PCONC(:,3&
&9)*PCONC(:,32)+0.11306*TPK%K111(:)*PCONC(:,33)*PCONC(:,39)+0.01593*TPK%K115(:)&
&*PCONC(:,39)*PCONC(:,39)+0.16271*TPK%K120(:)*PCONC(:,33)*PCONC(:,5)+0.01021*TP&
&K%K124(:)*PCONC(:,39)*PCONC(:,5)
!PLOSS(XO2) = +K126*<HO2>+K127*<MO2>+K128*<CARBOP>+K129*<XO2>+K129*<XO2>+K130*<
!NO>+K131*<NO3>
 PLOSS(:,41) = +TPK%K126(:)*PCONC(:,15)+TPK%K127(:)*PCONC(:,32)+TPK%K128(:)*PCO&
&NC(:,39)+TPK%K129(:)*PCONC(:,41)+TPK%K129(:)*PCONC(:,41)+TPK%K130(:)*PCONC(:,3&
&)+TPK%K131(:)*PCONC(:,5)
!
!PPROD(WC_O3) = +KTC1*<O3>
 PPROD(:,42) = +TPK%KTC1(:)*PCONC(:,1)
!PLOSS(WC_O3) = +KTC21+KC6*<WC_HO2>+KC29*<WC_SO2>
 PLOSS(:,42) = +TPK%KTC21(:)+TPK%KC6(:)*PCONC(:,53)+TPK%KC29(:)*PCONC(:,55)
!
!PPROD(WC_H2O2) = +KTC2*<H2O2>+KC2*<WC_OH>*<WC_OH>+KC5*<WC_HO2>*<WC_HO2>
 PPROD(:,43) = +TPK%KTC2(:)*PCONC(:,2)+TPK%KC2(:)*PCONC(:,52)*PCONC(:,52)+TPK%K&
&C5(:)*PCONC(:,53)*PCONC(:,53)
!PLOSS(WC_H2O2) = +KTC22+KC1+KC4*<WC_OH>+KC30*<WC_SO2>
 PLOSS(:,43) = +TPK%KTC22(:)+TPK%KC1(:)+TPK%KC4(:)*PCONC(:,52)+TPK%KC30(:)*PCON&
&C(:,55)
!
!PPROD(WC_NO) = +KTC3*<NO>
 PPROD(:,44) = +TPK%KTC3(:)*PCONC(:,3)
!PLOSS(WC_NO) = +KTC23
 PLOSS(:,44) = +TPK%KTC23(:)
!
!PPROD(WC_NO2) = +KTC4*<NO2>+KC8*<WC_HONO>*<WC_OH>+KC10*<WC_HNO4>+KC13*<WC_HNO3
!>
 PPROD(:,45) = +TPK%KTC4(:)*PCONC(:,4)+TPK%KC8(:)*PCONC(:,48)*PCONC(:,52)+TPK%K&
&C10(:)*PCONC(:,50)+TPK%KC13(:)*PCONC(:,49)
!PLOSS(WC_NO2) = +KTC24+KC9*<WC_HO2>
 PLOSS(:,45) = +TPK%KTC24(:)+TPK%KC9(:)*PCONC(:,53)
!
!PPROD(WC_NO3) = +KTC5*<NO3>
 PPROD(:,46) = +TPK%KTC5(:)*PCONC(:,5)
!PLOSS(WC_NO3) = +KTC25+KC15*<WC_SULF>+KC16*<WC_SO2>
 PLOSS(:,46) = +TPK%KTC25(:)+TPK%KC15(:)*PCONC(:,56)+TPK%KC16(:)*PCONC(:,55)
!
!PPROD(WC_N2O5) = +KTC6*<N2O5>
 PPROD(:,47) = +TPK%KTC6(:)*PCONC(:,6)
!PLOSS(WC_N2O5) = +KTC26+KC14
 PLOSS(:,47) = +TPK%KTC26(:)+TPK%KC14(:)
!
!PPROD(WC_HONO) = +KTC7*<HONO>+KC11*<WC_HNO4>
 PPROD(:,48) = +TPK%KTC7(:)*PCONC(:,7)+TPK%KC11(:)*PCONC(:,50)
!PLOSS(WC_HONO) = +KTC27+KC8*<WC_OH>
 PLOSS(:,48) = +TPK%KTC27(:)+TPK%KC8(:)*PCONC(:,52)
!
!PPROD(WC_HNO3) = +KTC8*<HNO3>+KC12*<WC_HNO4>*<WC_SO2>+KC14*<WC_N2O5>+KC14*<WC_
!N2O5>+KC15*<WC_NO3>*<WC_SULF>+KC16*<WC_NO3>*<WC_SO2>
 PPROD(:,49) = +TPK%KTC8(:)*PCONC(:,8)+TPK%KC12(:)*PCONC(:,50)*PCONC(:,55)+TPK%&
&KC14(:)*PCONC(:,47)+TPK%KC14(:)*PCONC(:,47)+TPK%KC15(:)*PCONC(:,46)*PCONC(:,56&
&)+TPK%KC16(:)*PCONC(:,46)*PCONC(:,55)
!PLOSS(WC_HNO3) = +KTC28+KC13
 PLOSS(:,49) = +TPK%KTC28(:)+TPK%KC13(:)
!
!PPROD(WC_HNO4) = +KTC9*<HNO4>+KC9*<WC_NO2>*<WC_HO2>
 PPROD(:,50) = +TPK%KTC9(:)*PCONC(:,9)+TPK%KC9(:)*PCONC(:,45)*PCONC(:,53)
!PLOSS(WC_HNO4) = +KTC29+KC10+KC11+KC12*<WC_SO2>
 PLOSS(:,50) = +TPK%KTC29(:)+TPK%KC10(:)+TPK%KC11(:)+TPK%KC12(:)*PCONC(:,55)
!
RETURN
END SUBROUTINE SUB4
!
SUBROUTINE SUB5
!
!Indices 51 a 60
!
!
!PPROD(WC_NH3) = +KTC10*<NH3>
 PPROD(:,51) = +TPK%KTC10(:)*PCONC(:,10)
!PLOSS(WC_NH3) = +KTC30
 PLOSS(:,51) = +TPK%KTC30(:)
!
!PPROD(WC_OH) = +KTC11*<OH>+KC1*<WC_H2O2>+KC1*<WC_H2O2>+KC6*<WC_O3>*<WC_HO2>+KC
!13*<WC_HNO3>+KC28*<WC_ASO4>
 PPROD(:,52) = +TPK%KTC11(:)*PCONC(:,14)+TPK%KC1(:)*PCONC(:,43)+TPK%KC1(:)*PCON&
&C(:,43)+TPK%KC6(:)*PCONC(:,42)*PCONC(:,53)+TPK%KC13(:)*PCONC(:,49)+TPK%KC28(:)&
&*PCONC(:,63)
!PLOSS(WC_OH) = +KTC31+KC2*<WC_OH>+KC2*<WC_OH>+KC3*<WC_HO2>+KC4*<WC_H2O2>+KC7*<
!WC_SO2>+KC8*<WC_HONO>+KC19*<WC_HCHO>+KC20*<WC_ORA1>+KC23*<WC_AHMS>
 PLOSS(:,52) = +TPK%KTC31(:)+TPK%KC2(:)*PCONC(:,52)+TPK%KC2(:)*PCONC(:,52)+TPK%&
&KC3(:)*PCONC(:,53)+TPK%KC4(:)*PCONC(:,43)+TPK%KC7(:)*PCONC(:,55)+TPK%KC8(:)*PC&
&ONC(:,48)+TPK%KC19(:)*PCONC(:,57)+TPK%KC20(:)*PCONC(:,58)+TPK%KC23(:)*PCONC(:,&
&66)
!
!PPROD(WC_HO2) = +KTC12*<HO2>+KC4*<WC_H2O2>*<WC_OH>+KC10*<WC_HNO4>+2.00*KC17*<W
!C_MO2>*<WC_MO2>+KC19*<WC_HCHO>*<WC_OH>+KC20*<WC_ORA1>*<WC_OH>+KC23*<WC_AHMS>*<
!WC_OH>
 PPROD(:,53) = +TPK%KTC12(:)*PCONC(:,15)+TPK%KC4(:)*PCONC(:,43)*PCONC(:,52)+TPK&
&%KC10(:)*PCONC(:,50)+2.00*TPK%KC17(:)*PCONC(:,60)*PCONC(:,60)+TPK%KC19(:)*PCON&
&C(:,57)*PCONC(:,52)+TPK%KC20(:)*PCONC(:,58)*PCONC(:,52)+TPK%KC23(:)*PCONC(:,66&
&)*PCONC(:,52)
!PLOSS(WC_HO2) = +KTC32+KC3*<WC_OH>+KC5*<WC_HO2>+KC5*<WC_HO2>+KC6*<WC_O3>+KC9*<
!WC_NO2>+KC25*<WC_ASO5>
 PLOSS(:,53) = +TPK%KTC32(:)+TPK%KC3(:)*PCONC(:,52)+TPK%KC5(:)*PCONC(:,53)+TPK%&
&KC5(:)*PCONC(:,53)+TPK%KC6(:)*PCONC(:,42)+TPK%KC9(:)*PCONC(:,45)+TPK%KC25(:)*P&
&CONC(:,64)
!
!PPROD(WC_CO2) = +KTC13*<CO2>+KC20*<WC_ORA1>*<WC_OH>
 PPROD(:,54) = +TPK%KTC13(:)*TPK%CO2(:)+TPK%KC20(:)*PCONC(:,58)*PCONC(:,52)
!PLOSS(WC_CO2) = +KTC33
 PLOSS(:,54) = +TPK%KTC33(:)
!
!PPROD(WC_SO2) = +KTC14*<SO2>+KC22*<WC_AHMS>+KC23*<WC_AHMS>*<WC_OH>
 PPROD(:,55) = +TPK%KTC14(:)*PCONC(:,11)+TPK%KC22(:)*PCONC(:,66)+TPK%KC23(:)*PC&
&ONC(:,66)*PCONC(:,52)
!PLOSS(WC_SO2) = +KTC34+KC7*<WC_OH>+KC12*<WC_HNO4>+KC16*<WC_NO3>+KC18*<WC_MO2>+
!KC21*<WC_HCHO>+KC27*<WC_AHSO5>+KC29*<WC_O3>+KC30*<WC_H2O2>
 PLOSS(:,55) = +TPK%KTC34(:)+TPK%KC7(:)*PCONC(:,52)+TPK%KC12(:)*PCONC(:,50)+TPK&
&%KC16(:)*PCONC(:,46)+TPK%KC18(:)*PCONC(:,60)+TPK%KC21(:)*PCONC(:,57)+TPK%KC27(&
&:)*PCONC(:,65)+TPK%KC29(:)*PCONC(:,42)+TPK%KC30(:)*PCONC(:,43)
!
!PPROD(WC_SULF) = +KTC15*<SULF>+KC12*<WC_HNO4>*<WC_SO2>+2.00*KC27*<WC_AHSO5>*<W
!C_SO2>+KC28*<WC_ASO4>+KC29*<WC_SO2>*<WC_O3>+KC30*<WC_SO2>*<WC_H2O2>
 PPROD(:,56) = +TPK%KTC15(:)*PCONC(:,12)+TPK%KC12(:)*PCONC(:,50)*PCONC(:,55)+2.&
&00*TPK%KC27(:)*PCONC(:,65)*PCONC(:,55)+TPK%KC28(:)*PCONC(:,63)+TPK%KC29(:)*PCO&
&NC(:,55)*PCONC(:,42)+TPK%KC30(:)*PCONC(:,55)*PCONC(:,43)
!PLOSS(WC_SULF) = +KTC35+KC15*<WC_NO3>
 PLOSS(:,56) = +TPK%KTC35(:)+TPK%KC15(:)*PCONC(:,46)
!
!PPROD(WC_HCHO) = +KTC16*<HCHO>+2.00*KC17*<WC_MO2>*<WC_MO2>+KC22*<WC_AHMS>
 PPROD(:,57) = +TPK%KTC16(:)*PCONC(:,22)+2.00*TPK%KC17(:)*PCONC(:,60)*PCONC(:,6&
&0)+TPK%KC22(:)*PCONC(:,66)
!PLOSS(WC_HCHO) = +KTC36+KC19*<WC_OH>+KC21*<WC_SO2>
 PLOSS(:,57) = +TPK%KTC36(:)+TPK%KC19(:)*PCONC(:,52)+TPK%KC21(:)*PCONC(:,55)
!
!PPROD(WC_ORA1) = +KTC17*<ORA1>+KC19*<WC_HCHO>*<WC_OH>+KC23*<WC_AHMS>*<WC_OH>
 PPROD(:,58) = +TPK%KTC17(:)*PCONC(:,30)+TPK%KC19(:)*PCONC(:,57)*PCONC(:,52)+TP&
&K%KC23(:)*PCONC(:,66)*PCONC(:,52)
!PLOSS(WC_ORA1) = +KTC37+KC20*<WC_OH>
 PLOSS(:,58) = +TPK%KTC37(:)+TPK%KC20(:)*PCONC(:,52)
!
!PPROD(WC_ORA2) = +KTC18*<ORA2>
 PPROD(:,59) = +TPK%KTC18(:)*PCONC(:,31)
!PLOSS(WC_ORA2) = +KTC38
 PLOSS(:,59) = +TPK%KTC38(:)
!
!PPROD(WC_MO2) = +KTC19*<MO2>
 PPROD(:,60) = +TPK%KTC19(:)*PCONC(:,32)
!PLOSS(WC_MO2) = +KTC39+KC17*<WC_MO2>+KC17*<WC_MO2>+KC18*<WC_SO2>
 PLOSS(:,60) = +TPK%KTC39(:)+TPK%KC17(:)*PCONC(:,60)+TPK%KC17(:)*PCONC(:,60)+TP&
&K%KC18(:)*PCONC(:,55)
!
RETURN
END SUBROUTINE SUB5
!
SUBROUTINE SUB6
!
!Indices 61 a 70
!
!
!PPROD(WC_OP1) = +KTC20*<OP1>+KC18*<WC_MO2>*<WC_SO2>
 PPROD(:,61) = +TPK%KTC20(:)*PCONC(:,28)+TPK%KC18(:)*PCONC(:,60)*PCONC(:,55)
!PLOSS(WC_OP1) = +KTC40
 PLOSS(:,61) = +TPK%KTC40(:)
!
!PPROD(WC_ASO3) = +KC7*<WC_OH>*<WC_SO2>+KC16*<WC_NO3>*<WC_SO2>+KC18*<WC_MO2>*<W
!C_SO2>
 PPROD(:,62) = +TPK%KC7(:)*PCONC(:,52)*PCONC(:,55)+TPK%KC16(:)*PCONC(:,46)*PCON&
&C(:,55)+TPK%KC18(:)*PCONC(:,60)*PCONC(:,55)
!PLOSS(WC_ASO3) = +KC24*<W_O2>
 PLOSS(:,62) = +TPK%KC24(:)*TPK%W_O2(:)
!
!PPROD(WC_ASO4) = +KC15*<WC_NO3>*<WC_SULF>+KC26*<WC_ASO5>*<WC_ASO5>+KC26*<WC_AS
!O5>*<WC_ASO5>
 PPROD(:,63) = +TPK%KC15(:)*PCONC(:,46)*PCONC(:,56)+TPK%KC26(:)*PCONC(:,64)*PCO&
&NC(:,64)+TPK%KC26(:)*PCONC(:,64)*PCONC(:,64)
!PLOSS(WC_ASO4) = +KC28
 PLOSS(:,63) = +TPK%KC28(:)
!
!PPROD(WC_ASO5) = +KC24*<WC_ASO3>*<W_O2>
 PPROD(:,64) = +TPK%KC24(:)*PCONC(:,62)*TPK%W_O2(:)
!PLOSS(WC_ASO5) = +KC25*<WC_HO2>+KC26*<WC_ASO5>+KC26*<WC_ASO5>
 PLOSS(:,64) = +TPK%KC25(:)*PCONC(:,53)+TPK%KC26(:)*PCONC(:,64)+TPK%KC26(:)*PCO&
&NC(:,64)
!
!PPROD(WC_AHSO5) = +KC25*<WC_ASO5>*<WC_HO2>
 PPROD(:,65) = +TPK%KC25(:)*PCONC(:,64)*PCONC(:,53)
!PLOSS(WC_AHSO5) = +KC27*<WC_SO2>
 PLOSS(:,65) = +TPK%KC27(:)*PCONC(:,55)
!
!PPROD(WC_AHMS) = +KC21*<WC_SO2>*<WC_HCHO>
 PPROD(:,66) = +TPK%KC21(:)*PCONC(:,55)*PCONC(:,57)
!PLOSS(WC_AHMS) = +KC22+KC23*<WC_OH>
 PLOSS(:,66) = +TPK%KC22(:)+TPK%KC23(:)*PCONC(:,52)
!
!PPROD(WR_O3) = +KTR1*<O3>
 PPROD(:,67) = +TPK%KTR1(:)*PCONC(:,1)
!PLOSS(WR_O3) = +KTR21+KR6*<WR_HO2>+KR29*<WR_SO2>
 PLOSS(:,67) = +TPK%KTR21(:)+TPK%KR6(:)*PCONC(:,78)+TPK%KR29(:)*PCONC(:,80)
!
!PPROD(WR_H2O2) = +KTR2*<H2O2>+KR2*<WR_OH>*<WR_OH>+KR5*<WR_HO2>*<WR_HO2>
 PPROD(:,68) = +TPK%KTR2(:)*PCONC(:,2)+TPK%KR2(:)*PCONC(:,77)*PCONC(:,77)+TPK%K&
&R5(:)*PCONC(:,78)*PCONC(:,78)
!PLOSS(WR_H2O2) = +KTR22+KR1+KR4*<WR_OH>+KR30*<WR_SO2>
 PLOSS(:,68) = +TPK%KTR22(:)+TPK%KR1(:)+TPK%KR4(:)*PCONC(:,77)+TPK%KR30(:)*PCON&
&C(:,80)
!
!PPROD(WR_NO) = +KTR3*<NO>
 PPROD(:,69) = +TPK%KTR3(:)*PCONC(:,3)
!PLOSS(WR_NO) = +KTR23
 PLOSS(:,69) = +TPK%KTR23(:)
!
!PPROD(WR_NO2) = +KTR4*<NO2>+KR8*<WR_HONO>*<WR_OH>+KR10*<WR_HNO4>+KR13*<WR_HNO3
!>
 PPROD(:,70) = +TPK%KTR4(:)*PCONC(:,4)+TPK%KR8(:)*PCONC(:,73)*PCONC(:,77)+TPK%K&
&R10(:)*PCONC(:,75)+TPK%KR13(:)*PCONC(:,74)
!PLOSS(WR_NO2) = +KTR24+KR9*<WR_HO2>
 PLOSS(:,70) = +TPK%KTR24(:)+TPK%KR9(:)*PCONC(:,78)
!
RETURN
END SUBROUTINE SUB6
!
SUBROUTINE SUB7
!
!Indices 71 a 80
!
!
!PPROD(WR_NO3) = +KTR5*<NO3>
 PPROD(:,71) = +TPK%KTR5(:)*PCONC(:,5)
!PLOSS(WR_NO3) = +KTR25+KR15*<WR_SULF>+KR16*<WR_SO2>
 PLOSS(:,71) = +TPK%KTR25(:)+TPK%KR15(:)*PCONC(:,81)+TPK%KR16(:)*PCONC(:,80)
!
!PPROD(WR_N2O5) = +KTR6*<N2O5>
 PPROD(:,72) = +TPK%KTR6(:)*PCONC(:,6)
!PLOSS(WR_N2O5) = +KTR26+KR14
 PLOSS(:,72) = +TPK%KTR26(:)+TPK%KR14(:)
!
!PPROD(WR_HONO) = +KTR7*<HONO>+KR11*<WR_HNO4>
 PPROD(:,73) = +TPK%KTR7(:)*PCONC(:,7)+TPK%KR11(:)*PCONC(:,75)
!PLOSS(WR_HONO) = +KTR27+KR8*<WR_OH>
 PLOSS(:,73) = +TPK%KTR27(:)+TPK%KR8(:)*PCONC(:,77)
!
!PPROD(WR_HNO3) = +KTR8*<HNO3>+KR12*<WR_HNO4>*<WR_SO2>+KR14*<WR_N2O5>+KR14*<WR_
!N2O5>+KR15*<WR_NO3>*<WR_SULF>+KR16*<WR_NO3>*<WR_SO2>
 PPROD(:,74) = +TPK%KTR8(:)*PCONC(:,8)+TPK%KR12(:)*PCONC(:,75)*PCONC(:,80)+TPK%&
&KR14(:)*PCONC(:,72)+TPK%KR14(:)*PCONC(:,72)+TPK%KR15(:)*PCONC(:,71)*PCONC(:,81&
&)+TPK%KR16(:)*PCONC(:,71)*PCONC(:,80)
!PLOSS(WR_HNO3) = +KTR28+KR13
 PLOSS(:,74) = +TPK%KTR28(:)+TPK%KR13(:)
!
!PPROD(WR_HNO4) = +KTR9*<HNO4>+KR9*<WR_NO2>*<WR_HO2>
 PPROD(:,75) = +TPK%KTR9(:)*PCONC(:,9)+TPK%KR9(:)*PCONC(:,70)*PCONC(:,78)
!PLOSS(WR_HNO4) = +KTR29+KR10+KR11+KR12*<WR_SO2>
 PLOSS(:,75) = +TPK%KTR29(:)+TPK%KR10(:)+TPK%KR11(:)+TPK%KR12(:)*PCONC(:,80)
!
!PPROD(WR_NH3) = +KTR10*<NH3>
 PPROD(:,76) = +TPK%KTR10(:)*PCONC(:,10)
!PLOSS(WR_NH3) = +KTR30
 PLOSS(:,76) = +TPK%KTR30(:)
!
!PPROD(WR_OH) = +KTR11*<OH>+KR1*<WR_H2O2>+KR1*<WR_H2O2>+KR6*<WR_O3>*<WR_HO2>+KR
!13*<WR_HNO3>+KR28*<WR_ASO4>
 PPROD(:,77) = +TPK%KTR11(:)*PCONC(:,14)+TPK%KR1(:)*PCONC(:,68)+TPK%KR1(:)*PCON&
&C(:,68)+TPK%KR6(:)*PCONC(:,67)*PCONC(:,78)+TPK%KR13(:)*PCONC(:,74)+TPK%KR28(:)&
&*PCONC(:,88)
!PLOSS(WR_OH) = +KTR31+KR2*<WR_OH>+KR2*<WR_OH>+KR3*<WR_HO2>+KR4*<WR_H2O2>+KR7*<
!WR_SO2>+KR8*<WR_HONO>+KR19*<WR_HCHO>+KR20*<WR_ORA1>+KR23*<WR_AHMS>
 PLOSS(:,77) = +TPK%KTR31(:)+TPK%KR2(:)*PCONC(:,77)+TPK%KR2(:)*PCONC(:,77)+TPK%&
&KR3(:)*PCONC(:,78)+TPK%KR4(:)*PCONC(:,68)+TPK%KR7(:)*PCONC(:,80)+TPK%KR8(:)*PC&
&ONC(:,73)+TPK%KR19(:)*PCONC(:,82)+TPK%KR20(:)*PCONC(:,83)+TPK%KR23(:)*PCONC(:,&
&91)
!
!PPROD(WR_HO2) = +KTR12*<HO2>+KR4*<WR_H2O2>*<WR_OH>+KR10*<WR_HNO4>+2.00*KR17*<W
!R_MO2>*<WR_MO2>+KR19*<WR_HCHO>*<WR_OH>+KR20*<WR_ORA1>*<WR_OH>+KR23*<WR_AHMS>*<
!WR_OH>
 PPROD(:,78) = +TPK%KTR12(:)*PCONC(:,15)+TPK%KR4(:)*PCONC(:,68)*PCONC(:,77)+TPK&
&%KR10(:)*PCONC(:,75)+2.00*TPK%KR17(:)*PCONC(:,85)*PCONC(:,85)+TPK%KR19(:)*PCON&
&C(:,82)*PCONC(:,77)+TPK%KR20(:)*PCONC(:,83)*PCONC(:,77)+TPK%KR23(:)*PCONC(:,91&
&)*PCONC(:,77)
!PLOSS(WR_HO2) = +KTR32+KR3*<WR_OH>+KR5*<WR_HO2>+KR5*<WR_HO2>+KR6*<WR_O3>+KR9*<
!WR_NO2>+KR25*<WR_ASO5>
 PLOSS(:,78) = +TPK%KTR32(:)+TPK%KR3(:)*PCONC(:,77)+TPK%KR5(:)*PCONC(:,78)+TPK%&
&KR5(:)*PCONC(:,78)+TPK%KR6(:)*PCONC(:,67)+TPK%KR9(:)*PCONC(:,70)+TPK%KR25(:)*P&
&CONC(:,89)
!
!PPROD(WR_CO2) = +KTR13*<CO2>+KR20*<WR_ORA1>*<WR_OH>
 PPROD(:,79) = +TPK%KTR13(:)*TPK%CO2(:)+TPK%KR20(:)*PCONC(:,83)*PCONC(:,77)
!PLOSS(WR_CO2) = +KTR33
 PLOSS(:,79) = +TPK%KTR33(:)
!
!PPROD(WR_SO2) = +KTR14*<SO2>+KR22*<WR_AHMS>+KR23*<WR_AHMS>*<WR_OH>
 PPROD(:,80) = +TPK%KTR14(:)*PCONC(:,11)+TPK%KR22(:)*PCONC(:,91)+TPK%KR23(:)*PC&
&ONC(:,91)*PCONC(:,77)
!PLOSS(WR_SO2) = +KTR34+KR7*<WR_OH>+KR12*<WR_HNO4>+KR16*<WR_NO3>+KR18*<WR_MO2>+
!KR21*<WR_HCHO>+KR27*<WR_AHSO5>+KR29*<WR_O3>+KR30*<WR_H2O2>
 PLOSS(:,80) = +TPK%KTR34(:)+TPK%KR7(:)*PCONC(:,77)+TPK%KR12(:)*PCONC(:,75)+TPK&
&%KR16(:)*PCONC(:,71)+TPK%KR18(:)*PCONC(:,85)+TPK%KR21(:)*PCONC(:,82)+TPK%KR27(&
&:)*PCONC(:,90)+TPK%KR29(:)*PCONC(:,67)+TPK%KR30(:)*PCONC(:,68)
!
RETURN
END SUBROUTINE SUB7
!
SUBROUTINE SUB8
!
!Indices 81 a 90
!
!
!PPROD(WR_SULF) = +KTR15*<SULF>+KR12*<WR_HNO4>*<WR_SO2>+2.00*KR27*<WR_AHSO5>*<W
!R_SO2>+KR28*<WR_ASO4>+KR29*<WR_SO2>*<WR_O3>+KR30*<WR_SO2>*<WR_H2O2>
 PPROD(:,81) = +TPK%KTR15(:)*PCONC(:,12)+TPK%KR12(:)*PCONC(:,75)*PCONC(:,80)+2.&
&00*TPK%KR27(:)*PCONC(:,90)*PCONC(:,80)+TPK%KR28(:)*PCONC(:,88)+TPK%KR29(:)*PCO&
&NC(:,80)*PCONC(:,67)+TPK%KR30(:)*PCONC(:,80)*PCONC(:,68)
!PLOSS(WR_SULF) = +KTR35+KR15*<WR_NO3>
 PLOSS(:,81) = +TPK%KTR35(:)+TPK%KR15(:)*PCONC(:,71)
!
!PPROD(WR_HCHO) = +KTR16*<HCHO>+2.00*KR17*<WR_MO2>*<WR_MO2>+KR22*<WR_AHMS>
 PPROD(:,82) = +TPK%KTR16(:)*PCONC(:,22)+2.00*TPK%KR17(:)*PCONC(:,85)*PCONC(:,8&
&5)+TPK%KR22(:)*PCONC(:,91)
!PLOSS(WR_HCHO) = +KTR36+KR19*<WR_OH>+KR21*<WR_SO2>
 PLOSS(:,82) = +TPK%KTR36(:)+TPK%KR19(:)*PCONC(:,77)+TPK%KR21(:)*PCONC(:,80)
!
!PPROD(WR_ORA1) = +KTR17*<ORA1>+KR19*<WR_HCHO>*<WR_OH>+KR23*<WR_AHMS>*<WR_OH>
 PPROD(:,83) = +TPK%KTR17(:)*PCONC(:,30)+TPK%KR19(:)*PCONC(:,82)*PCONC(:,77)+TP&
&K%KR23(:)*PCONC(:,91)*PCONC(:,77)
!PLOSS(WR_ORA1) = +KTR37+KR20*<WR_OH>
 PLOSS(:,83) = +TPK%KTR37(:)+TPK%KR20(:)*PCONC(:,77)
!
!PPROD(WR_ORA2) = +KTR18*<ORA2>
 PPROD(:,84) = +TPK%KTR18(:)*PCONC(:,31)
!PLOSS(WR_ORA2) = +KTR38
 PLOSS(:,84) = +TPK%KTR38(:)
!
!PPROD(WR_MO2) = +KTR19*<MO2>
 PPROD(:,85) = +TPK%KTR19(:)*PCONC(:,32)
!PLOSS(WR_MO2) = +KTR39+KR17*<WR_MO2>+KR17*<WR_MO2>+KR18*<WR_SO2>
 PLOSS(:,85) = +TPK%KTR39(:)+TPK%KR17(:)*PCONC(:,85)+TPK%KR17(:)*PCONC(:,85)+TP&
&K%KR18(:)*PCONC(:,80)
!
!PPROD(WR_OP1) = +KTR20*<OP1>+KR18*<WR_MO2>*<WR_SO2>
 PPROD(:,86) = +TPK%KTR20(:)*PCONC(:,28)+TPK%KR18(:)*PCONC(:,85)*PCONC(:,80)
!PLOSS(WR_OP1) = +KTR40
 PLOSS(:,86) = +TPK%KTR40(:)
!
!PPROD(WR_ASO3) = +KR7*<WR_OH>*<WR_SO2>+KR16*<WR_NO3>*<WR_SO2>+KR18*<WR_MO2>*<W
!R_SO2>
 PPROD(:,87) = +TPK%KR7(:)*PCONC(:,77)*PCONC(:,80)+TPK%KR16(:)*PCONC(:,71)*PCON&
&C(:,80)+TPK%KR18(:)*PCONC(:,85)*PCONC(:,80)
!PLOSS(WR_ASO3) = +KR24*<W_O2>
 PLOSS(:,87) = +TPK%KR24(:)*TPK%W_O2(:)
!
!PPROD(WR_ASO4) = +KR15*<WR_NO3>*<WR_SULF>+KR26*<WR_ASO5>*<WR_ASO5>+KR26*<WR_AS
!O5>*<WR_ASO5>
 PPROD(:,88) = +TPK%KR15(:)*PCONC(:,71)*PCONC(:,81)+TPK%KR26(:)*PCONC(:,89)*PCO&