CMAQ中如何修改N2O5非均相反应速率常数

congxue

  在CMAQ5.3.3中N2O5非均相反应速率常数如何修改？ 其中# 1.0~<HETERO_N2O5IJ> 等是什么定义？

931404656

在mechanism文件中，例如mech_cb6r3_ae6_aq.def

congxue

931404656 发表于 2022-9-27 10:54
在mechanism文件中，例如mech_cb6r3_ae6_aq.def

您好，非常谢谢您能够回复我的问题，我是初学很茫然，现在能够找到这个文件，但是不知道如何操作去修改这些反应的非均相反应速率常数，请问能具体说明一下吗？

931404656

<HET_NTR2> NTR2 = HNO3 # 1.4~<HETERO_NTR2>;

<HET_N2O5IJ>    N2O5 = HNO3 + H2NO3PIJ  # 1.0~<HETERO_N2O5IJ>;
<HET_N2O5K>     N2O5 = HNO3 + H2NO3PK   # 1.0~<HETERO_N2O5K>;

<HET_H2NO3PIJA> H2NO3PIJ = HNO3         # 1.0~<HETERO_H2NO3PAIJ>;
<HET_H2NO3PKA>  H2NO3PK  = HNO3         # 1.0~<HETERO_H2NO3PAK>;

<HET_H2NO3PIB>  H2NO3PIJ + ACLI = CLNO2 # 1.0~<HETERO_H2NO3PBIJ>;
<HET_H2NO3PJB>  H2NO3PIJ + ACLJ = CLNO2 # 1.0~<HETERO_H2NO3PBIJ>;
<HET_H2NO3PKB>  H2NO3PK  + ACLK = CLNO2 # 1.0~<HETERO_H2NO3PBK>;

<HET_N02> NO2 = 0.5*HONO + 0.5*HNO3 # 1.0~<HETERO_NO2>;
上述#后面的数即为反应速率常数

congxue

931404656 发表于 2022-9-27 20:15
NTR2 = HNO3 # 1.4~;

    N2O5 = HNO3 + H2NO3PIJ  # 1.0~;

您好，很感谢您再一次的回复，我还有一些不明白的地方，请问：1、修改反应速率常数是在#后面直接修改吗？还是说还有其他code也要修改 ；2、<HETERO_N2O5IJ>、<HETERO_H2NO3PAIJ> 等这些参数表示的是什么呀？为什么默认的是反应速率常数是 1.0~<HETERO_N2O5IJ>呀？

931404656

1. #后面直接修改；2.异质反应的参与物，其中I  J分别表示爱根核模和积聚模；3.按code显示反应常数是1，所以你可以改改

congxue

931404656 发表于 2022-9-28 08:35
1. #后面直接修改；2.异质反应的参与物，其中I  J分别表示爱根核模和积聚模；3.按code显示反应常数是1，所 ...

您好，感谢您的回复，就具体问题我有点困惑，我想根据公式k=(ω1&#8729;γ&#8729;S_a)/4修改其反应速率常数，请问该如何修改#后面的速率常速，其中<HETERO_N2O5IJ>是一个label还是一个特定的参数，“~”表示的是什么？谢谢！

931404656

~只是特定的筛选符，用于Fortran读取

931404656

具体参考RXNS_FUNC_MODULE.F90文件中写法

pip2013

931404656 发表于 2022-9-28 08:35
1. #后面直接修改；2.异质反应的参与物，其中I  J分别表示爱根核模和积聚模；3.按code显示反应常数是1，所 ...

您好，请问下我修改了mech*.def文件，cctm重新编译了，但是修改前后模拟出来的值完全一致，硝酸盐小时浓度保留三位小数，366个时次完全一致，可能是什么原因呢？