Skip to content

Commit

Permalink
Merge pull request #176 from DirkOlivie/PR_soayield
Browse files Browse the repository at this point in the history 
noresm2_3_develop : modified SOA yield factors for isoprene
  • Loading branch information
@DirkOlivie
DirkOlivie authored Nov 19, 2024
2 parents bb9f01f + 410fcb0 commit ee8a7c6
Show file tree
Hide file tree
Showing 8 changed files with 25 additions and 25 deletions.
10 changes: 5 additions & 5 deletions src/chemistry/pp_trop_mam_oslo/chem_mech.doc
View file
Original file line number Diff line number Diff line change
Expand Up @@ -109,9 +109,9 @@ Class List
( 7) monoterp + O3 -> .15*SOA_LV rate = 8.05E-16*exp( -640./t) ( 8)
( 8) monoterp + OH -> .15*SOA_SV rate = 1.20E-11*exp( 440./t) ( 9)
( 9) monoterp + NO3 -> .15*SOA_SV rate = 1.20E-12*exp( 490./t) ( 10)
( 10) isoprene + O3 -> .05*SOA_SV rate = 1.03E-14*exp( -1995./t) ( 11)
( 11) isoprene + OH -> .05*SOA_SV rate = 2.70E-11*exp( 390./t) ( 12)
( 12) isoprene + NO3 -> .05*SOA_SV rate = 3.15E-12*exp( -450./t) ( 13)
( 10) isoprene + O3 -> .005*SOA_SV rate = 1.03E-14*exp( -1995./t) ( 11)
( 11) isoprene + OH -> .005*SOA_SV rate = 2.70E-11*exp( 390./t) ( 12)
( 12) isoprene + NO3 -> .005*SOA_SV rate = 3.15E-12*exp( -450./t) ( 13)

Heterogeneous loss species

Expand Down Expand Up @@ -155,8 +155,8 @@ Extraneous prod/loss species
d(SOA_NA)/dt = 0
d(SOA_A1)/dt = 0
d(SOA_LV)/dt = .029*r6*OH*DMS + .15*r7*O3*monoterp
d(SOA_SV)/dt = .114*r6*OH*DMS + .15*r8*OH*monoterp + .15*r9*NO3*monoterp + .05*r10*O3*isoprene
+ .05*r11*OH*isoprene + .05*r12*NO3*isoprene
d(SOA_SV)/dt = .114*r6*OH*DMS + .15*r8*OH*monoterp + .15*r9*NO3*monoterp + .005*r10*O3*isoprene
+ .005*r11*OH*isoprene + .005*r12*NO3*isoprene
d(monoterp)/dt = - r7*O3*monoterp - r8*OH*monoterp - r9*NO3*monoterp
d(isoprene)/dt = - r10*O3*isoprene - r11*OH*isoprene - r12*NO3*isoprene
d(H2O)/dt = r2*OH*H2O2
6 changes: 3 additions & 3 deletions src/chemistry/pp_trop_mam_oslo/chem_mech.in
View file
Original file line number Diff line number Diff line change
Expand Up @@ -72,9 +72,9 @@ BEGSIM
monoterp + O3 -> .15*SOA_LV ; 8.05e-16, -640.
monoterp + OH -> .15*SOA_SV ; 1.2e-11, 440.
monoterp + NO3 -> .15*SOA_SV ; 1.2e-12, 490.
isoprene + O3 -> .05*SOA_SV ; 1.03e-14, -1995.
isoprene + OH -> .05*SOA_SV ; 2.7e-11, 390.
isoprene + NO3 -> .05*SOA_SV ; 3.15e-12, -450.
isoprene + O3 -> .005*SOA_SV ; 1.03e-14, -1995.
isoprene + OH -> .005*SOA_SV ; 2.7e-11, 390.
isoprene + NO3 -> .005*SOA_SV ; 3.15e-12, -450.
End Reactions

Heterogeneous
Expand Down
2 changes: 1 addition & 1 deletion src/chemistry/pp_trop_mam_oslo/mo_lin_matrix.F90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -49,7 +49,7 @@ subroutine linmat01( mat, y, rxt, het_rates )
mat(32) = -( het_rates(27) )
mat(4) = .114_r8*rxt(7)
mat(34) = .150_r8*rxt(9) + .150_r8*rxt(10)
mat(36) = .050_r8*rxt(11) + .050_r8*rxt(12) + .050_r8*rxt(13)
mat(36) = .005_r8*rxt(11) + .005_r8*rxt(12) + .005_r8*rxt(13)
mat(35) = -( rxt(8) + rxt(9) + rxt(10) + het_rates(28) )
mat(37) = -( rxt(11) + rxt(12) + rxt(13) + het_rates(29) )
mat(38) = -( het_rates(30) )
Expand Down
2 changes: 1 addition & 1 deletion src/chemistry/pp_trop_mam_oslo/mo_prod_loss.F90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -85,7 +85,7 @@ subroutine imp_prod_loss( prod, loss, y, rxt, het_rates )
loss(26) = ( + het_rates(26))* y(26)
prod(26) =.029_r8*rxt(7)*y(3) +.150_r8*rxt(8)*y(28)
loss(27) = ( + het_rates(27))* y(27)
prod(27) = (.050_r8*rxt(11) +.050_r8*rxt(12) +.050_r8*rxt(13))*y(29) &
prod(27) = (.005_r8*rxt(11) +.005_r8*rxt(12) +.005_r8*rxt(13))*y(29) &
+ (.150_r8*rxt(9) +.150_r8*rxt(10))*y(28) +.114_r8*rxt(7)*y(3)
loss(28) = ( + rxt(8) + rxt(9) + rxt(10) + het_rates(28))* y(28)
prod(28) = 0._r8
Expand Down
10 changes: 5 additions & 5 deletions src/chemistry/pp_tropstrat_mam_oslo/chem_mech.doc
View file
Original file line number Diff line number Diff line change
Expand Up @@ -1059,9 +1059,9 @@ Class List
(391) MTERP + O3 -> MTERP + O3 + .15*SOA_LV rate = 6.30E-16*exp( -580./t) (504)
(392) MTERP + OH -> MTERP + OH + .15*SOA_SV rate = 1.20E-11*exp( 440./t) (505)
(393) MTERP + NO3 -> MTERP + NO3 + .15*SOA_SV rate = 1.20E-12*exp( 490./t) (506)
(394) ISOP + O3 -> ISOP + O3 + .05*SOA_SV rate = 1.05E-14*exp( -2000./t) (507)
(395) ISOP + OH -> ISOP + OH + .05*SOA_SV rate = 2.54E-11*exp( 410./t) (508)
(396) ISOP + NO3 -> ISOP + NO3 + .05*SOA_SV rate = 3.03E-12*exp( -446./t) (509)
(394) ISOP + O3 -> ISOP + O3 + .005*SOA_SV rate = 1.05E-14*exp( -2000./t) (507)
(395) ISOP + OH -> ISOP + OH + .005*SOA_SV rate = 2.54E-11*exp( 410./t) (508)
(396) ISOP + NO3 -> ISOP + NO3 + .005*SOA_SV rate = 3.03E-12*exp( -446./t) (509)

Extraneous prod/loss species
( 1) BC_NI (dataset)
Expand Down Expand Up @@ -1102,8 +1102,8 @@ Extraneous prod/loss species
d(SOA_NA)/dt = 0
d(SOA_A1)/dt = 0
d(SOA_LV)/dt = .029*r356*DMS*OH + .15*r391*MTERP*O3
d(SOA_SV)/dt = .114*r356*DMS*OH + .15*r392*MTERP*OH + .15*r393*MTERP*NO3 + .05*r394*ISOP*O3 + .05*r395*ISOP*OH
+ .05*r396*ISOP*NO3
d(SOA_SV)/dt = .114*r356*DMS*OH + .15*r392*MTERP*OH + .15*r393*MTERP*NO3 + .005*r394*ISOP*O3
+ .005*r395*ISOP*OH + .005*r396*ISOP*NO3
d(ALKNIT)/dt = r236*ALKO2*NO
- j19*ALKNIT - r233*OH*ALKNIT
d(ALKOOH)/dt = r234*ALKO2*HO2
Expand Down
6 changes: 3 additions & 3 deletions src/chemistry/pp_tropstrat_mam_oslo/chem_mech.in
View file
Original file line number Diff line number Diff line change
Expand Up @@ -1173,9 +1173,9 @@ BEGSIM
MTERP + O3 -> MTERP + O3 + .15*SOA_LV ; 6.3e-16, -580
MTERP + OH -> MTERP + OH + .15*SOA_SV ; 1.2e-11, 440
MTERP + NO3 -> MTERP + NO3 + .15*SOA_SV ; 1.2e-12, 490
ISOP + O3 -> ISOP + O3 + .05*SOA_SV ; 1.05e-14, -2000
ISOP + OH -> ISOP + OH + .05*SOA_SV ; 2.54e-11, 410
ISOP + NO3 -> ISOP + NO3 + .05*SOA_SV ; 3.03e-12, -446
ISOP + O3 -> ISOP + O3 + .005*SOA_SV ; 1.05e-14, -2000
ISOP + OH -> ISOP + OH + .005*SOA_SV ; 2.54e-11, 410
ISOP + NO3 -> ISOP + NO3 + .005*SOA_SV ; 3.03e-12, -446
End Reactions

Ext Forcing
Expand Down
10 changes: 5 additions & 5 deletions src/chemistry/pp_tropstrat_mam_oslo/mo_nln_matrix.F90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -27,12 +27,12 @@ subroutine nlnmat01( avec_len, mat, y, rxt )
mat(k,1634) = .150_r8*rxt(k,504)*y(k,127)
mat(k,1432) = .029_r8*rxt(k,469)*y(k,88)
mat(k,155) = .114_r8*rxt(k,469)*y(k,200)
mat(k,671) = .050_r8*rxt(k,509)*y(k,140) + .050_r8*rxt(k,507)*y(k,144) &
+ .050_r8*rxt(k,508)*y(k,200)
mat(k,671) = .005_r8*rxt(k,509)*y(k,140) + .005_r8*rxt(k,507)*y(k,144) &
+ .005_r8*rxt(k,508)*y(k,200)
mat(k,716) = .150_r8*rxt(k,506)*y(k,140) + .150_r8*rxt(k,505)*y(k,200)
mat(k,1732) = .050_r8*rxt(k,509)*y(k,116) + .150_r8*rxt(k,506)*y(k,127)
mat(k,1635) = .050_r8*rxt(k,507)*y(k,116)
mat(k,1433) = .114_r8*rxt(k,469)*y(k,88) + .050_r8*rxt(k,508)*y(k,116) &
mat(k,1732) = .005_r8*rxt(k,509)*y(k,116) + .150_r8*rxt(k,506)*y(k,127)
mat(k,1635) = .005_r8*rxt(k,507)*y(k,116)
mat(k,1433) = .114_r8*rxt(k,469)*y(k,88) + .005_r8*rxt(k,508)*y(k,116) &
+ .150_r8*rxt(k,505)*y(k,127)
mat(k,477) = -(rxt(k,346)*y(k,200))
mat(k,1502) = -rxt(k,346)*y(k,24)
Expand Down
4 changes: 2 additions & 2 deletions src/chemistry/pp_tropstrat_mam_oslo/mo_prod_loss.F90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -167,8 +167,8 @@ subroutine imp_prod_loss( avec_len, prod, loss, y, &
prod(k,22) =.029_r8*rxt(k,469)*y(k,200)*y(k,88) +.150_r8*rxt(k,504)*y(k,144) &
*y(k,127)
loss(k,23) = ( + het_rates(k,23))* y(k,23)
prod(k,23) = (.050_r8*rxt(k,507)*y(k,144) +.050_r8*rxt(k,508)*y(k,200) + &
.050_r8*rxt(k,509)*y(k,140))*y(k,116) &
prod(k,23) = (.005_r8*rxt(k,507)*y(k,144) +.005_r8*rxt(k,508)*y(k,200) + &
.005_r8*rxt(k,509)*y(k,140))*y(k,116) &
+ (.150_r8*rxt(k,505)*y(k,200) +.150_r8*rxt(k,506)*y(k,140)) &
*y(k,127) +.114_r8*rxt(k,469)*y(k,200)*y(k,88)
loss(k,108) = (rxt(k,346)* y(k,200) + rxt(k,19) + het_rates(k,24))* y(k,24)
Expand Down

0 comments on commit ee8a7c6

Please sign in to comment.