01zzCMAQ

CMAQ入门回顾

下图是CMAQ整体的运行流程图,这里需要注意的是SMOKE模块是没有在CMAQ程序中的
图一
整个CMAQ的运行流程如下:

  1. 由WRF或者是MM5模型输出的气象模型输出数据可以通过MCIP生成【带有坐标和网格定义信息的网格描述文件】、【单点、交差点、边界气象文件】。这是整个CMAQ所需要的最基础的一些文件。
  2. 通过SMOKE模型生成【为指定气相化学机理和 PM 模型物种化的依赖时间的 2-D或 3-D 排放文件的名称和位置】、【点源排放文件(the point emissions file for secter)】、【(stack groups files for secter)】、【物种形成廓线文件(包含用于将聚集清单污染物排放总和分为CCTM需要的模型物种排放的因子)】、【泛化生物排放用于 BEIS3 的输入】、【转换霜冻日期的文件(依赖时间的并且通常包含一整年的数据,365或366。它使用一个变量SEASON,其值为0,网格单元格应该使用当天的冬季因子,或1,网格单元格应该使用当天的夏季因子)】。以上有些文件是CCTM的可选输入文件。
  3. 通过ICON生成IC文件,这是CCTM的初始时间污染物的浓度。
  4. 通过BCON生成BC文件,这是CCTM的处在边界上污染物的浓度。

注意到还有CCTM输出到ICON和BCON的箭头,这是为了进行嵌套计算的,因为当第一次CCTM运行的时候,ICON是从一个默认的ic_profile文件中读取的初始值而生成的IC文件,BCON是从一个默认的bc_profile文件中读取的边界初始值而生成的BC文件。如果进行嵌套模拟,则需要选择CTM_CONC_1,即输出的关于污染物的名称和位置的估计,小时均值,CONC文件。

CMAQ整体核心的计算式:

∂C/∂t = Adv + Diff + Rc + Ec – Sc

其中:
∂C/∂t是每个单元格上随时间改变的浓度
Adv = 平流
Diff = 扩散
Rc = 物种 c 的化学转化
Ec = 物种 c 的排放
Sc = 物种 c 的损失过程
在 CMAQ 中,平流和排放项分别基于由气象和排放模型生成的输入文件计算; 扩散、化学转化和损失过程项在 CCTM计算。(平流和扩散的区别是什么?)

CMAQ各个模块的个人理解

以下所有的分析是参照了CMAQ中文手册和github(CMAQv5.2.1 Operational Guidance Document)上相关资料加上自己的理解,全属个人观点。

MCIP

Meteorological Mechanism Interface Program,气象化学接口程序就是使用WRF或者MM5气象模型的输出文件创建适合SMOKE(计算CMAQ的排放输入的排放处理器)和CMAQ使用的netCDF格式的输入文件。

BCON

Boundary Condition,边界条件就是考虑到在域的边界上可能有处在域外的污染物和其他化学物种向域内的平流,这里是假设不存在通过域的顶部边界进行的交换(即垂直交换)。有两种选择,一种是模型默认的污染物浓度初值(profile),另一种就是进行嵌套模拟而选择的CTM_CONC_1文件(m3conc)。

ICON

Initial Condition,初始条件就是模拟在模拟域中第一个小时的化学条件,有两种选择,一种是模型默认的污染物浓度初值(profile),另一种就是进行嵌套模拟而选择的CTM_CONC_1文件(m3conc)。

JPROC

INLINE_PHOT_PREPROC,晴空光解速率处理器就是利用制成表的吸收截面积和量子产率(CSQY)数据计算固定高度、太阳时间角以及纬度带的指定化学机理晴空光解速率。需要注意的是CCTM包含了一个使用预测的臭氧和气溶胶计算光解速率的嵌入光解选项,如果在编译CCTM的时候ModPhot 设置为 phot_inline则不需要JPROC;如果ModPhot设置为phot_table则需要编译运行JPROC生成成日光解速率查找表。

CCTM

CMAQ Chemical Transport Model,CMAQ化学传输模型就是用于对平流、扩散和对流机制的大气化学过程、传输和沉降过程进行模拟。CCTM 整合前面介绍的前处理程序(JPROC, BCON, ICON 和 MCIP)的输出结果以及为 CMAQ 准备的排放输入参数(如 SMOKE 的输出),用来模拟连续的大气化学条件。

CMAQv5.1个人说明

刚开始只是为了能够将CMAQ跑起来,并没有认真仔细的研究CMAQ中的script,即各个模块的编译选项和执行选项,然后是从官网上下载的data作为测试数据跑起来的。下面总结一下使用测试数据所遇到的keng。

  • 下面是开始从官网上下载的测试数据直接解压得到的data。
    在这里插入图片描述
  • 下面是运行CMAQ之后得到的data,多了一个cctm文件夹。
    在这里插入图片描述
  1. 参考网上的安装教程和官网的安装手册,它会说要我们去编译script文件夹中mcip,但是没说要运行mcip的可执行文件。其实根本不需要编译运行mcip,因为mcip的输入文件是wrf或者mm5的输出文件,data里面并没有;况且在给的data文件中有一个mcip文件夹,这里面其实就是运行mcip的输出文件。
  2. 这里面最大的keng就是其实在给出的data文件中有icon和bcon两个文件夹,文件夹中分别有一个文件,这其实是icon和bcon编译执行所输出的文件。那再看script中cctm的run.cctm可执行文件,其中所指定的ICFILE和BCFILE就是给定的data中icon和bcon文件夹里面的文件。因此其实也根本不需要自己再去分别编译执行icon和bcon。
    在这里插入图片描述
  • 6
    点赞
  • 19
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值