WRF-Chem笔记——MOZBC边界场制作

好久没有跑WRF-Chem了，诶，模式更新换代太快，不学习马上就要out了。

今天写这篇是师弟发现WRF-Chem边界场制作从原先的MOZART-4资料改用CAM-Chem边界场资料了，我就帮重新测试了一下新文件的边界场制作，顺便记个笔记。

旧版本边界场制作使用的是MOZART数据（目前停止运营了）：
MOZART：http://www.acom.ucar.edu/wrf-chem/mozart.shtml

新版本边界场制作使用的是CAM-Chem的输出文件：
MOZART：https://www2.acom.ucar.edu/gcm/cam-chem-output

相关网址：

1. WRF-Chem工具：https://www2.acom.ucar.edu/wrf-chem/wrf-chem-tools-community

2. CAM-Chem主页：https://wiki.ucar.edu/display/camchem/Home

3. CAM-Chem作为边界场：https://wiki.ucar.edu/display/camchem/CESM2.1%3ACAM-chem+as+Boundary+Conditions

4. inp文件修改参考：https://www2.acom.ucar.edu/sites/default/files/wrf-chem/CESM-WRFchem_aerosols_20190822.pdf

下面说具体步骤：

step1: 下载边界场数据

相关处理资料下载网址：https://www.acom.ucar.edu/cam-chem/cam-chem.shtml

下载注意事项：
东西南北边界要覆盖住整个domain1，如果是盖住整个中国的话，可以选择经度范围45-160，纬度范围5-65。

时间要覆盖住整个模拟时间，可以对应模拟时间前后加一天。

step2: 利用mozbc制作边界场

进入.. /MOZBC/目录下，新的CAM-Chem output其实就是对应改一下*inp文件的spc_map,参考https://www2.acom.ucar.edu/sites/default/files/wrf-chem/CESM-WRFchem_aerosols_20190822.pdf

截图不完整哈，点进链接去看后面还有很多物种。

注意：

1. 化学机制需要和chem_opt对应上，如：CBMZ对应chem_opt=9

2. 下面两行不能丢！

 moz_var_suffix = ''def_missing_var = .true.

以CBMZ为例的inp文件修改，4bin

“vim CBMZ_CAM_Chem_4bins.inp”,

spc_map = 'o3->O3', 'n2o -> N2O', 'no -> NO','no2 -> NO2', 'nh3 -> NH3', 'hno3 -> HNO3', 'hno4 -> HO2NO2','n2o5 -> N2O5', 'h2o2 -> H2O2','ch4 -> CH4', 'co -> CO', 'ch3ooh -> CH3OOH','hcho -> CH2O', 'ch3oh -> CH3OH', 'c2h4 -> C2H4','ald -> CH3CHO', 'acet -> CH3COCH3', 'mgly -> CH3COCHO','pan -> PAN', 'mpan -> MPAN', 'macr -> MACR','mvk -> MVK', 'c2h6 -> C2H6', 'c3h6 -> C3H6', 'c3h8 -> C3H8','c2h5oh -> C2H5OH', 'c10h16 -> MTERP','isopr -> ISOP','acetol -> HYAC', 'mek -> MEK','bigene -> BIGENE', 'bigalk -> BIGALK','tol -> TOLUENE', 'benzene -> BENZENE', 'xylenes -> XYLENES','cres -> CRESOL', 'dms -> DMS', 'so2 -> SO2','oc_a03->0.1164*pom_a1+0.0000*soa1_a2+0.1164*soa1_a1+0.0000*soa2_a2+0.1164*soa2_a1+0.0000*soa3_a2+0.1164*soa3_a1+0.1164*soa4_a1+0.1164*soa5_a1;1.e9','oc_a04->0.0002*pom_a1+0.0000*soa1_a2+0.0002*soa1_a1+0.0000*soa2_a2+0.0002*soa2_a1+0.0000*soa3_a2+0.0002*soa3_a1+0.0002*soa4_a1+0.0002*soa5_a1;1.e9','bc_a01->0.0093*bc_a1+0.0093*bc_a4+0.1123*bc_a1+0.1123*bc_a4;1.e9','bc_a02->0.3835*bc_a1+0.3835*bc_a4+0.3783*bc_a1+0.3783*bc_a4;1.e9',
'bc_a03->0.1077*bc_a1+0.1077*bc_a4+0.0087*bc_a1+0.0087*bc_a4;1.e9','bc_a04->0.0002*bc_a1+0.0002*bc_a4+0.0000*bc_a1+0.0000*bc_a4;1.e9','so4_a01->0.7510*so4_a2+0.0093*so4_a1+0.0000*so4_a3+0.2376*so4_a2+0.1123*so4_a1+0.0000*so4_a3;1.e9','so4_a02->0.0113*so4_a2+0.3835*so4_a1+0.0000*so4_a3+0.0001*so4_a2+0.3783*so4_a1+0.0002*so4_a3;1.e9','so4_a03->0.0000*so4_a2+0.1077*so4_a1+0.0061*so4_a3+0.0000*so4_a2+0.0087*so4_a1+0.0934*so4_a3;1.e9','so4_a04->0.0000*so4_a2+0.0002*so4_a1+0.4020*so4_a3+0.0000*so4_a2+0.0000*so4_a1+0.4983*so4_a3;1.e9','nh4_a01->0.1410*so4_a2+0.0033*so4_a1+0.0000*so4_a3+0.0446*so4_a2+0.0017*so4_a1+0.0000*so4_a3;1.e9','nh4_a02->0.0021*so4_a2+0.0210*so4_a1+0.0000*so4_a3+0.0000*so4_a2+0.0720*so4_a1+0.0000*so4_a3;1.e9','nh4_a03->0.0000*so4_a2+0.0202*so4_a1+0.0011*so4_a3+0.0000*so4_a2+0.0001*so4_a1+0.0175*so4_a3;1.e9','nh4_a04->0.0000*so4_a2+0.0000*so4_a1+0.0755*so4_a3+0.0000*so4_a2+0.0000*so4_a1+0.0935*so4_a3;1.e9','no3_a01->0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3+0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3;1.e9','no3_a02->0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3+0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3;1.e9','no3_a03->0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3+0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3;1.e9','no3_a04->0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3+0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3;1.e9','na_a01->0.2954*ncl_a2+0.0037*ncl_a1+0.0000*ncl_a3+0.0935*ncl_a2+0.0442*ncl_a1+0.0000*ncl_a3;1.e9','na_a02->0.0045*ncl_a2+0.1509*ncl_a1+0.0000*ncl_a3+0.0000*ncl_a2+0.1488*ncl_a1+0.0000*ncl_a3;1.e9','na_a03->0.0000*ncl_a2+0.0424*ncl_a1+0.0024*ncl_a3+0.0000*ncl_a2+0.0034*ncl_a1+0.0367*ncl_a3;1.e9','na_a04->0.0000*ncl_a2+0.0000*ncl_a1+0.1582*ncl_a3+0.0000*ncl_a2+0.0000*ncl_a1+0.1960*ncl_a3;1.e9','cl_a01->0.4555*ncl_a2+0.0056*ncl_a1+0.0000*ncl_a3+0.1441*ncl_a2+0.0681*ncl_a1+0.0000*ncl_a3;1.e9',
'cl_a02->0.0068*ncl_a2+0.2326*ncl_a1+0.0000*ncl_a3+0.0000*ncl_a2+0.2295*ncl_a1+0.0000*ncl_a3;1.e9','cl_a03->0.0000*ncl_a2+0.0654*ncl_a1+0.0037*ncl_a3+0.0000*ncl_a2+0.0055*ncl_a1+0.0567*ncl_a3;1.e9','cl_a04->0.0000*ncl_a2+0.0001*ncl_a1+0.2439*ncl_a3+0.0000*ncl_a2+0.0000*ncl_a1+0.3023*ncl_a3;1.e9','oin_a01->0.7510*dst_a2+0.0093*dst_a1+0.0000*dst_a3+0.2376*dst_a2+0.1123*dst_a1+0.0000*dst_a3;1.e9','oin_a02->0.0113*dst_a2+0.3835*dst_a1+0.0000*dst_a3+0.0001*dst_a2+0.3783*dst_a1+0.0002*dst_a3;1.e9','oin_a03->0.0000*dst_a2+0.1077*dst_a1+0.0061*dst_a3+0.0000*dst_a2+0.0087*dst_a1+0.0934*dst_a3;1.e9','oin_a04->0.0000*dst_a2+0.0002*dst_a1+0.4020*dst_a3+0.0000*dst_a2+0.0000*dst_a1+0.4983*dst_a3;1.e9','num_a01->0.9502*num_a2+0.2509*num_a1+0.0000*num_a3+0.0494*num_a2+0.4626*num_a1+0.0000*num_a3;1.0','num_a02->0.0004*num_a2+0.2470*num_a1+0.0007*num_a3+0.0000*num_a2+0.0377*num_a1+0.0232*num_a3;1.0','num_a03->0.0000*num_a2+0.0016*num_a1+0.1886*num_a3+0.0000*num_a2+0.0000*num_a1+0.4372*num_a3;1.0','num_a04->0.0000*num_a2+0.0000*num_a1+0.2935*num_a3+0.0000*num_a2+0.0000*num_a1+0.0566*num_a3;1.0','oc_a01->0.1216*pom_a1+0.9886*soa1_a2+0.122*soa1_a1+0.98*soa2_a2+0.1216*soa2_a1+0.98*soa3_a2+0.1216*soa3_a1+0.98*soa4_a2+0.1216*soa4_a1+0.98*soa5_a2+0.1216*soa5_a1;1.e9','oc_a02->0.7618*pom_a1+0.0114*soa1_a2+0.762*soa1_a1+0.0114*soa2_a2+0.76*soa2_a1+0.0114*soa3_a2+0.76*soa3_a1+0.0114*soa4_a2+0.76*soa4_a1+0.0114*soa5_a2+0.76*soa5_a1;1.e9',

注意:

1. 我这里因为是4个bin，所以手动把参考文件中的对应8个bin的两两相加了。

2. oc_a01,oc_a02我放到最后，这两个略微有点偏差，因为这两行太长了超出单行字符限制，所以我把某几个四位小数约成了两位。。。没办法，我fortran太菜，不知道咋改，引号中用&换行貌似不奏效，有知道怎么处理的UU也可以指导我一下。

step3: 修改目录运行

修改对应的domain，dir_wrf，等

“./mozbc < CBMZ_CAM_Chem_4bins.inp”

domain =1,2,3…… 从外到里修改进行边界场制作

WRF-Chem笔记——MOZBC边界场制作相关推荐

区域气象-大气化学在线耦合模式（WRF/Chem）在大气环境中的应用
[查看原文]区域气象-大气化学在线耦合模式(WRF/Chem)在大气环境中的应用随着我国经济快速发展,我国面临着日益严重的大气污染问题.近年来,严重的大气污染问题已经明显影响国计民生,引起政府.学界 ...
Docker笔记三 Docker镜像制作
Docker笔记三 Docker镜像制作 1.Docker镜像制作方法: docker commit 保存当前container的状态到镜像,生成p_w_picpath. docker build 利 ...
WRF 和 WPS 3.9.1 以及 wrf chem、wrf hydro安装
WRF 的安装其实没有什么难点,主要是依赖库比较繁琐,还有两个用的比较多的 chem 和 hydro,我在安装的时候在网上找的教程都不尽如人意,所以在这里总结一下机器配置以及编译器选择 intel ...
Unity学习笔记—二次元日系游戏制作（实践篇-游戏初始化场景制作）
原教程:siki:二次元日系游戏制作工具 - live2dSDK入门教程 http://www.sikiedu.com/my/course/282 (上)Unity学习笔记-二次元日系游戏制作(理论篇 ...
Flutter学习笔记--Gridview网格组件制作相册
Flutter学习笔记–Gridview网格组件制作相册效果: GridView可以构建一个二维网格列表,其官方定义参考: https://book.flutterchina.club/chapte ...
unity3d-学习笔记10-卡牌游戏制作（实现点击卡片消失，并将后面的卡片填充的功能）
unity3d-学习笔记10-卡牌游戏制作(实现点击卡片消失,并将后面的卡片填充的功能) 实现点击卡片消失,并将后面的卡片填充的功能实现更换unity启动的图标实现点击卡片消失,并将后面的卡片填充 ...
咸鱼Maya笔记—NURBS边界
咸鱼Maya笔记-NURBS边界使用菜单栏中的"曲面>边界"命令,可以在已创建的边界曲线之间进行相应的插值填充,从而创建出三维曲面,具体操作方法如下. 打开Maya应用程序 ...
[Unity学习笔记：FPS游戏制作（3）]子弹拖尾，碰撞与枪口火焰效果
往期博客[Unity学习笔记:FPS游戏制作(2)] 发射子弹----(2021.6.20学习笔记) 文章目录一,实现思路二,粒子效果的实现 (1)子弹拖尾特效的实现 (2)枪口火焰特效的实现 ( ...
[Unity学习笔记：FPS游戏制作（2）] 发射子弹————（2021.6.20学习笔记）
往期博客 [Unity学习笔记:FPS游戏制作(1)]角色的移动,旋转与推进上升----(2021.6.13学习笔记) 文章目录一,实现思路二,实现代码三,脚本的使用方法四,最终效果一,实现 ...
Polyworks脚本开发学习笔记(九)-公差控制及制作报告
Polyworks脚本开发学习笔记(九)-公差控制及制作报告定义公差系统有默认的公差设置,可以在选项中进行系统的默认设置,但往往不是想要的.比如下图的XYZ三向都是-/+1,我想只控制Y向并且公差 ...

WRF-Chem笔记——MOZBC边界场制作

WRF-Chem笔记——MOZBC边界场制作相关推荐

最新文章

热门文章