Skip to content

Commit c51d3ab

Browse files
fvittfvitt
committed
remove duplicate N2D_O2 reaction from waccm_t4ma_mam5
1 parent 93e53b4 commit c51d3ab

12 files changed

Lines changed: 3859 additions & 3869 deletions

‎src/chemistry/pp_waccm_t4ma_mam5/chem_mech.doc‎

Lines changed: 647 additions & 648 deletions
Large diffs are not rendered by default.

‎src/chemistry/pp_waccm_t4ma_mam5/chem_mech.in‎

Lines changed: 0 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -544,7 +544,6 @@
544544
[HO2NO2_OH] HO2NO2 + OH -> H2O + NO2 + O2 ; 4.5e-13, 610
545545
[N2D_NO,cph=543.13] N2D + NO -> N2 + O ; 7e-11
546546
[N2D_O,cph=229.61] N2D + O -> N + O ; 7e-13
547-
[N2D_O2,cph=177.51] N2D + O2 -> NO + O1D ; 5e-12
548547
[N_NO,cph=313.75] N + NO -> N2 + O ; 2.1e-11, 100
549548
[N_NO2a] N + NO2 -> N2O + O ; 2.9e-12, 220
550549
[N_NO2b] N + NO2 -> 2*NO ; 1.45e-12, 220

‎src/chemistry/pp_waccm_t4ma_mam5/chem_mods.F90‎

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -6,8 +6,8 @@ module chem_mods
66
implicit none
77
save
88
integer, parameter :: phtcnt = 114, & ! number of photolysis reactions
9-
rxntot = 439, & ! number of total reactions
10-
gascnt = 325, & ! number of gas phase reactions
9+
rxntot = 438, & ! number of total reactions
10+
gascnt = 324, & ! number of gas phase reactions
1111
nabscol = 2, & ! number of absorbing column densities
1212
gas_pcnst = 151, & ! number of "gas phase" species
1313
nfs = 2, & ! number of "fixed" species
@@ -23,8 +23,8 @@ module chem_mods
2323
indexm = 1, & ! index of total atm density in invariant array
2424
indexh2o = 0, & ! index of water vapor density
2525
clsze = 1, & ! loop length for implicit chemistry
26-
rxt_tag_cnt = 439, &
27-
enthalpy_cnt = 57, &
26+
rxt_tag_cnt = 438, &
27+
enthalpy_cnt = 56, &
2828
nslvd = 13
2929
integer :: clscnt(5) = 0
3030
integer :: cls_rxt_cnt(4,5) = 0

‎src/chemistry/pp_waccm_t4ma_mam5/m_rxt_id.F90‎

Lines changed: 283 additions & 284 deletions
Large diffs are not rendered by default.

‎src/chemistry/pp_waccm_t4ma_mam5/mo_adjrxt.F90‎

Lines changed: 33 additions & 34 deletions
Original file line numberDiff line numberDiff line change
@@ -20,31 +20,31 @@ subroutine adjrxt( rate, inv, m, ncol, nlev )
2020
rate(:,:, 134) = rate(:,:, 134) * inv(:,:, 1)
2121
rate(:,:, 140) = rate(:,:, 140) * inv(:,:, 1)
2222
rate(:,:, 150) = rate(:,:, 150) * inv(:,:, 1)
23-
rate(:,:, 163) = rate(:,:, 163) * inv(:,:, 1)
24-
rate(:,:, 171) = rate(:,:, 171) * inv(:,:, 1)
23+
rate(:,:, 162) = rate(:,:, 162) * inv(:,:, 1)
24+
rate(:,:, 170) = rate(:,:, 170) * inv(:,:, 1)
25+
rate(:,:, 173) = rate(:,:, 173) * inv(:,:, 1)
2526
rate(:,:, 174) = rate(:,:, 174) * inv(:,:, 1)
2627
rate(:,:, 175) = rate(:,:, 175) * inv(:,:, 1)
27-
rate(:,:, 176) = rate(:,:, 176) * inv(:,:, 1)
28-
rate(:,:, 178) = rate(:,:, 178) * inv(:,:, 1)
29-
rate(:,:, 181) = rate(:,:, 181) * inv(:,:, 1)
30-
rate(:,:, 197) = rate(:,:, 197) * inv(:,:, 1)
28+
rate(:,:, 177) = rate(:,:, 177) * inv(:,:, 1)
29+
rate(:,:, 180) = rate(:,:, 180) * inv(:,:, 1)
30+
rate(:,:, 196) = rate(:,:, 196) * inv(:,:, 1)
31+
rate(:,:, 216) = rate(:,:, 216) * inv(:,:, 1)
3132
rate(:,:, 217) = rate(:,:, 217) * inv(:,:, 1)
32-
rate(:,:, 218) = rate(:,:, 218) * inv(:,:, 1)
33-
rate(:,:, 228) = rate(:,:, 228) * inv(:,:, 1)
34-
rate(:,:, 277) = rate(:,:, 277) * inv(:,:, 1)
33+
rate(:,:, 227) = rate(:,:, 227) * inv(:,:, 1)
34+
rate(:,:, 276) = rate(:,:, 276) * inv(:,:, 1)
35+
rate(:,:, 301) = rate(:,:, 301) * inv(:,:, 1)
3536
rate(:,:, 302) = rate(:,:, 302) * inv(:,:, 1)
3637
rate(:,:, 303) = rate(:,:, 303) * inv(:,:, 1)
37-
rate(:,:, 304) = rate(:,:, 304) * inv(:,:, 1)
38-
rate(:,:, 323) = rate(:,:, 323) * inv(:,:, 1)
39-
rate(:,:, 340) = rate(:,:, 340) * inv(:,:, 1)
38+
rate(:,:, 322) = rate(:,:, 322) * inv(:,:, 1)
39+
rate(:,:, 339) = rate(:,:, 339) * inv(:,:, 1)
40+
rate(:,:, 342) = rate(:,:, 342) * inv(:,:, 1)
4041
rate(:,:, 343) = rate(:,:, 343) * inv(:,:, 1)
41-
rate(:,:, 344) = rate(:,:, 344) * inv(:,:, 1)
42-
rate(:,:, 375) = rate(:,:, 375) * inv(:,:, 1)
43-
rate(:,:, 425) = rate(:,:, 425) * inv(:,:, 2)
44-
rate(:,:, 428) = rate(:,:, 428) * inv(:,:, 2)
45-
rate(:,:, 432) = rate(:,:, 432) * inv(:,:, 2)
42+
rate(:,:, 374) = rate(:,:, 374) * inv(:,:, 1)
43+
rate(:,:, 424) = rate(:,:, 424) * inv(:,:, 2)
44+
rate(:,:, 427) = rate(:,:, 427) * inv(:,:, 2)
45+
rate(:,:, 431) = rate(:,:, 431) * inv(:,:, 2)
46+
rate(:,:, 436) = rate(:,:, 436) * inv(:,:, 2)
4647
rate(:,:, 437) = rate(:,:, 437) * inv(:,:, 2)
47-
rate(:,:, 438) = rate(:,:, 438) * inv(:,:, 2)
4848
rate(:,:, 117) = rate(:,:, 117) * m(:,:)
4949
rate(:,:, 118) = rate(:,:, 118) * m(:,:)
5050
rate(:,:, 120) = rate(:,:, 120) * m(:,:)
@@ -101,9 +101,9 @@ subroutine adjrxt( rate, inv, m, ncol, nlev )
101101
rate(:,:, 174) = rate(:,:, 174) * m(:,:)
102102
rate(:,:, 175) = rate(:,:, 175) * m(:,:)
103103
rate(:,:, 176) = rate(:,:, 176) * m(:,:)
104-
rate(:,:, 177) = rate(:,:, 177) * m(:,:)
104+
rate(:,:, 178) = rate(:,:, 178) * m(:,:)
105105
rate(:,:, 179) = rate(:,:, 179) * m(:,:)
106-
rate(:,:, 180) = rate(:,:, 180) * m(:,:)
106+
rate(:,:, 181) = rate(:,:, 181) * m(:,:)
107107
rate(:,:, 182) = rate(:,:, 182) * m(:,:)
108108
rate(:,:, 183) = rate(:,:, 183) * m(:,:)
109109
rate(:,:, 184) = rate(:,:, 184) * m(:,:)
@@ -139,7 +139,7 @@ subroutine adjrxt( rate, inv, m, ncol, nlev )
139139
rate(:,:, 214) = rate(:,:, 214) * m(:,:)
140140
rate(:,:, 215) = rate(:,:, 215) * m(:,:)
141141
rate(:,:, 216) = rate(:,:, 216) * m(:,:)
142-
rate(:,:, 217) = rate(:,:, 217) * m(:,:)
142+
rate(:,:, 218) = rate(:,:, 218) * m(:,:)
143143
rate(:,:, 219) = rate(:,:, 219) * m(:,:)
144144
rate(:,:, 220) = rate(:,:, 220) * m(:,:)
145145
rate(:,:, 221) = rate(:,:, 221) * m(:,:)
@@ -217,13 +217,13 @@ subroutine adjrxt( rate, inv, m, ncol, nlev )
217217
rate(:,:, 293) = rate(:,:, 293) * m(:,:)
218218
rate(:,:, 294) = rate(:,:, 294) * m(:,:)
219219
rate(:,:, 295) = rate(:,:, 295) * m(:,:)
220-
rate(:,:, 296) = rate(:,:, 296) * m(:,:)
220+
rate(:,:, 297) = rate(:,:, 297) * m(:,:)
221221
rate(:,:, 298) = rate(:,:, 298) * m(:,:)
222222
rate(:,:, 299) = rate(:,:, 299) * m(:,:)
223223
rate(:,:, 300) = rate(:,:, 300) * m(:,:)
224224
rate(:,:, 301) = rate(:,:, 301) * m(:,:)
225225
rate(:,:, 302) = rate(:,:, 302) * m(:,:)
226-
rate(:,:, 303) = rate(:,:, 303) * m(:,:)
226+
rate(:,:, 304) = rate(:,:, 304) * m(:,:)
227227
rate(:,:, 305) = rate(:,:, 305) * m(:,:)
228228
rate(:,:, 306) = rate(:,:, 306) * m(:,:)
229229
rate(:,:, 307) = rate(:,:, 307) * m(:,:)
@@ -262,7 +262,7 @@ subroutine adjrxt( rate, inv, m, ncol, nlev )
262262
rate(:,:, 340) = rate(:,:, 340) * m(:,:)
263263
rate(:,:, 341) = rate(:,:, 341) * m(:,:)
264264
rate(:,:, 342) = rate(:,:, 342) * m(:,:)
265-
rate(:,:, 343) = rate(:,:, 343) * m(:,:)
265+
rate(:,:, 344) = rate(:,:, 344) * m(:,:)
266266
rate(:,:, 345) = rate(:,:, 345) * m(:,:)
267267
rate(:,:, 346) = rate(:,:, 346) * m(:,:)
268268
rate(:,:, 347) = rate(:,:, 347) * m(:,:)
@@ -305,15 +305,15 @@ subroutine adjrxt( rate, inv, m, ncol, nlev )
305305
rate(:,:, 384) = rate(:,:, 384) * m(:,:)
306306
rate(:,:, 385) = rate(:,:, 385) * m(:,:)
307307
rate(:,:, 386) = rate(:,:, 386) * m(:,:)
308-
rate(:,:, 387) = rate(:,:, 387) * m(:,:)
309-
rate(:,:, 396) = rate(:,:, 396) * m(:,:)
308+
rate(:,:, 395) = rate(:,:, 395) * m(:,:)
309+
rate(:,:, 400) = rate(:,:, 400) * m(:,:)
310310
rate(:,:, 401) = rate(:,:, 401) * m(:,:)
311311
rate(:,:, 402) = rate(:,:, 402) * m(:,:)
312-
rate(:,:, 403) = rate(:,:, 403) * m(:,:)
312+
rate(:,:, 405) = rate(:,:, 405) * m(:,:)
313313
rate(:,:, 406) = rate(:,:, 406) * m(:,:)
314314
rate(:,:, 407) = rate(:,:, 407) * m(:,:)
315-
rate(:,:, 408) = rate(:,:, 408) * m(:,:)
316-
rate(:,:, 411) = rate(:,:, 411) * m(:,:)
315+
rate(:,:, 410) = rate(:,:, 410) * m(:,:)
316+
rate(:,:, 414) = rate(:,:, 414) * m(:,:)
317317
rate(:,:, 415) = rate(:,:, 415) * m(:,:)
318318
rate(:,:, 416) = rate(:,:, 416) * m(:,:)
319319
rate(:,:, 417) = rate(:,:, 417) * m(:,:)
@@ -323,16 +323,15 @@ subroutine adjrxt( rate, inv, m, ncol, nlev )
323323
rate(:,:, 421) = rate(:,:, 421) * m(:,:)
324324
rate(:,:, 422) = rate(:,:, 422) * m(:,:)
325325
rate(:,:, 423) = rate(:,:, 423) * m(:,:)
326-
rate(:,:, 424) = rate(:,:, 424) * m(:,:)
326+
rate(:,:, 425) = rate(:,:, 425) * m(:,:)
327327
rate(:,:, 426) = rate(:,:, 426) * m(:,:)
328-
rate(:,:, 427) = rate(:,:, 427) * m(:,:)
328+
rate(:,:, 428) = rate(:,:, 428) * m(:,:)
329329
rate(:,:, 429) = rate(:,:, 429) * m(:,:)
330330
rate(:,:, 430) = rate(:,:, 430) * m(:,:)
331-
rate(:,:, 431) = rate(:,:, 431) * m(:,:)
331+
rate(:,:, 432) = rate(:,:, 432) * m(:,:)
332332
rate(:,:, 433) = rate(:,:, 433) * m(:,:)
333333
rate(:,:, 434) = rate(:,:, 434) * m(:,:)
334334
rate(:,:, 435) = rate(:,:, 435) * m(:,:)
335-
rate(:,:, 436) = rate(:,:, 436) * m(:,:)
336-
rate(:,:, 439) = rate(:,:, 439) * m(:,:)
335+
rate(:,:, 438) = rate(:,:, 438) * m(:,:)
337336
end subroutine adjrxt
338337
end module mo_adjrxt

‎src/chemistry/pp_waccm_t4ma_mam5/mo_indprd.F90‎

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -20,9 +20,9 @@ subroutine indprd( class, prod, nprod, y, extfrc, rxt, chnkpnts )
2020
! ... "independent" production for Explicit species
2121
!--------------------------------------------------------------------
2222
if( class == 1 ) then
23-
prod(:,1) =rxt(:,387)*y(:,111)*y(:,95) +rxt(:,390)*y(:,96)
24-
prod(:,2) = (.500_r8*rxt(:,340)*y(:,87) +.600_r8*rxt(:,360)*y(:,112)) &
25-
*y(:,111) +rxt(:,359)*y(:,112)*y(:,99)
23+
prod(:,1) =rxt(:,386)*y(:,111)*y(:,95) +rxt(:,389)*y(:,96)
24+
prod(:,2) = (.500_r8*rxt(:,339)*y(:,87) +.600_r8*rxt(:,359)*y(:,112)) &
25+
*y(:,111) +rxt(:,358)*y(:,112)*y(:,99)
2626
!--------------------------------------------------------------------
2727
! ... "independent" production for Implicit species
2828
!--------------------------------------------------------------------

0 commit comments

Comments
 (0)