1. 项目概述从零到一掌握CMAQ如果你正在研究空气质量模拟或者你的课题、项目需要评估污染物在大气中的传输、转化和沉降过程那么CMAQCommunity Multiscale Air Quality模型几乎是一个绕不开的工具。作为一个开源的第三代空气质量模型CMAQ能够综合模拟臭氧、颗粒物PM2.5/PM10、有毒有害空气污染物等多种污染物的复杂物理化学过程。我第一次接触CMAQ是为了做一个区域性的臭氧污染来源解析项目当时面对海量的源代码、复杂的依赖库和动辄几十个的配置文件确实感到头皮发麻。网上能找到的教程要么过于简略跳过了关键步骤要么年代久远环境依赖早已更新换代照着做总是卡在某个编译错误或运行报错上。经过几个项目的反复折腾我逐渐摸清了从系统准备、源码编译、案例运行到结果解读的一整套流程。这篇内容就是把这些踩过的坑和总结的经验系统化目标很明确让你能在一台干净的Linux服务器或工作站上成功安装、配置并运行起CMAQ最后还能看懂它输出的那些数据到底在说什么。无论你是环境科学、大气科学专业的学生还是从事环境规划、评估的工程师这篇手把手的指南应该都能帮你省下大量摸索的时间。2. CMAQ模式安装与配置全解析2.1 系统环境准备与依赖库安装CMAQ的运行严重依赖一个稳定且兼容的Linux环境。我强烈推荐使用CentOS 7或Rocky Linux 8/9这类企业级稳定的发行版Ubuntu也可以但需要特别注意库版本的兼容性。我的主力环境是Rocky Linux 8它提供了长期支持软件源稳定。第一步不是急着下载CMAQ代码而是搭建好它的“地基”——编译器和基础依赖库。核心编译器套件我选择Intel oneAPI HPC Toolkit。虽然GNU编译器也能用但Intel编译器针对Intel处理器优化得更好生成的代码运行效率通常有可观的提升。你可以从Intel官网免费下载安装。安装后记得通过source /opt/intel/oneapi/setvars.sh路径可能因安装而异来初始化编译环境。接下来是一系列基础依赖库的安装它们就像是CMAQ的“砖瓦”NetCDF库这是CMAQ读写数据的基础必须安装。需要同时安装C、Fortran和C接口。我习惯从Unidata官网下载源码编译这样可以严格控制版本。通常NetCDF-C和NetCDF-Fortran需要单独编译且Fortran库依赖于C库。编译时务必指定--prefix参数安装到自定义路径如/opt/libs/netcdf并确保在环境变量中正确设置。MPI库如果你想进行并行计算对于区域稍大的模拟这是必须的需要安装MPI。Intel oneAPI自带Intel MPI直接可用。也可以选择OpenMPI或MPICH但需要确保其与编译器的兼容性。IOAPI库这是CMAQ模型输入输出专用的I/O库同样需要源码编译。它依赖于NetCDF。编译IOAPI时需要修改Makefile或Makeinclude文件正确指向你的NetCDF库路径和编译器。一个关键的实操心得是统一管理所有依赖库。我会创建一个专门的目录比如/opt/libs/将所有自行编译的库NetCDF, IOAPI等都安装在这里。然后通过~/.bashrc文件永久设置相关的环境变量export CCicc export CXXicpc export FCifort export F77ifort export CFLAGS-O3 -xHost export FFLAGS-O3 -xHost export NETCDF/opt/libs/netcdf export IOAPI/opt/libs/ioapi export PATH$NETCDF/bin:$IOAPI/bin:$PATH export LD_LIBRARY_PATH$NETCDF/lib:$IOAPI/lib:$LD_LIBRARY_PATH export LIBRARY_PATH$NETCDF/lib:$IOAPI/lib:$LIBRARY_PATH export CPATH$NETCDF/include:$IOAPI/include:$CPATH这样做的好处是环境清晰避免系统自带库的干扰并且在后续编译CMAQ时配置脚本能自动找到这些路径。完成这些后务必执行source ~/.bashrc使环境变量生效并通过nc-config --all和which mpiifort等命令验证库是否可用。2.2 CMAQ源码获取与编译构建环境搭好后就可以获取CMAQ的源代码了。最新的稳定版本通常发布在美国环保署EPA的CMAS中心网站。下载的源码包解压后你会看到一个结构清晰的目录树。核心的编译配置集中在CMAQ_REPO/config_cmaq.csh或config_cmaq.sh这个配置脚本中。编译的第一步是配置。你需要根据你的环境编辑这个配置文件。主要修改以下几项编译器选择将CMAQ_COMPILER设置为intel如果你用Intel编译器或gcc。并行模式将CMAQ_MPI设置为openmpi或mpi对于Intel MPI。如果你只是测试可以设为serial进行串行编译。库路径设置这是最关键的一步。确保NETCDF和IOAPI的路径指向你之前安装的位置。脚本里会有类似setenv NETCDF /your/path/to/netcdf的语句修改它。应用配置CMAQ包含多个可执行程序称为“应用”如气象化学耦合接口ICON/BCON、主化学传输模块CCTM等。在配置脚本中你可以选择需要编译哪些应用。配置完成后在源码根目录下执行编译脚本如./bldit_cctm.csh或./bldit_project.csh。编译过程会持续一段时间如果前面的依赖库配置正确通常能一次通过。如果报错最常见的错误信息是“找不到xxx.h文件”或“对xxx函数未定义的引用”。这几乎百分之百是库路径或环境变量设置有问题。你需要回头仔细检查NETCDF和IOAPI的路径是否包含include和lib子目录并且LD_LIBRARY_PATH是否已正确导出。编译成功后会在各个应用目录的BLD子目录下生成可执行文件例如CCTM模块的可执行文件可能是cctm_5.4_intel_mpi.exe。你可以将其拷贝到一个集中的bin目录方便后续调用。3. 运行案例与关键配置详解3.1 案例数据准备与运行目录结构CMAQ官方通常会提供测试案例Test Case这是上手最快的方式。案例包中包含了运行一次完整模拟所需的所有输入文件气象驱动数据如WRF模型的输出、初始和边界条件、排放清单、土地利用数据、光解速率文件以及各种配置表格如化学机制文件CHEMISTRY_*.nml。拿到案例数据后第一件事不是直接运行而是理解其目录结构。一个典型的CMAQ运行目录会包含以下子目录/inputs/存放所有静态和动态输入文件。例如/inputs/mcip/里是气象化学接口处理器MCIP处理后的气象文件/inputs/emis/里是各种排放源文件/inputs/surface/里是土地利用文件等。/scripts/存放运行脚本。这是控制模拟的核心。脚本会设置模拟的起止时间、网格信息、化学机制、要调用的可执行程序路径等。/outputs/或/data/运行脚本中指定的输出目录模拟结果将写入这里。你需要将编译好的可执行文件如cctm.exe也放到这个运行目录下或者修改运行脚本中的路径指向它。一个重要的注意事项是检查所有输入文件的路径引用。运行脚本通常是.csh或.sh文件里会用环境变量或绝对路径来指定输入文件的位置。你必须根据你实际存放案例数据的路径修改这些脚本中的路径设置否则模型会因为找不到输入文件而立即失败。3.2 运行脚本核心参数解析与修改运行脚本是CMAQ模拟的“大脑”。我们以最常见的CCTM运行脚本为例拆解几个你必须修改的关键参数模拟网格定义setenv GRID_NAME “案例网格名”这个名称必须与输入文件如气象文件的网格标识完全一致。setenv GRIDDESC “./GRIDDESC”指定了网格描述文件路径该文件定义了网格的投影方式、范围、分辨率等。模拟时间设置setenv START_DATE和setenv END_DATE定义了模拟的开始和结束日期格式为YYYY-MM-DD。同时setenv START_TIME和setenv END_TIME定义了每天的开始和结束小时UTC时间。模型会按照这个时间循环执行。化学机制选择setenv MECH “cb6r3_ae7_aq”指定了使用的化学机制。不同的机制包含的物种和化学反应不同这直接影响了模拟的化学过程复杂度和结果。你必须确保排放清单等输入文件中的物种与此化学机制兼容。输入文件路径一系列setenv命令定义了气象文件、初始条件文件、边界条件文件、排放文件等的具体路径。你必须逐项核对确保路径有效。输出控制setenv CTM_CONC_1 “T”控制是否输出浓度文件。setenv AVG_CONC_1 “1”表示输出1小时平均浓度。你还可以设置输出哪些物种通过setenv CONC_SPCS以控制输出文件的大小。并行设置setenv NPCOL 4和setenv NPROW 4定义了在水平方向上将网格划分为4x416个进程进行计算。这个数值需要根据你的总CPU核心数和网格大小来调整通常设置为能整除网格X和Y方向格点数的值以获得较好的负载均衡。修改完脚本后不要急于提交。先执行source命令加载脚本设置的环境变量然后尝试用echo命令打印关键变量如echo $GRID_NAME, $START_DATE确认设置无误。最后通过mpirun -np 16 ./cctm_5.4_intel_mpi.exe | tee run.log这样的命令提交并行任务并将屏幕输出重定向到日志文件run.log中便于后续排查问题。4. 结果解译与可视化分析实战4.1 理解CMAQ输出文件结构与内容CMAQ成功运行后会在输出目录生成一系列NetCDF格式的文件。文件名通常包含网格名、模拟日期、变量类型等信息例如CCTM_CONC_1_20180701.nc。理解这些文件是解译结果的第一步。使用ncdump -h filename.nc命令可以查看文件的“头信息”这里面包含了文件的维度、变量、全局属性等元数据。关键维度通常包括TSTEP时间步长维度。LAY垂直层维度。ROW和COL网格的行和列维度。VAR变量维度在某些文件格式中。关键的变量就是各种污染物的浓度单位通常是ppmV气态污染物或μg/m³颗粒物。例如变量O3代表臭氧PM25_TOT代表PM2.5总质量浓度。除了浓度文件还可能包含干湿沉降通量、能见度等诊断变量。一个核心技巧是使用Python的xarray和netCDF4库来处理这些数据。它们比传统的NCL或GrADS更灵活、更强大。你可以轻松地读取文件、选取特定时间、层次、区域或物种进行切片计算。import xarray as xr # 打开浓度文件 ds xr.open_dataset(‘CCTM_CONC_1_20180701.nc’) # 提取地表层通常LAY0的O3浓度 o3_surface ds[‘O3’].isel(LAY0) # 计算整个区域在某个时间点的平均值 o3_avg o3_surface.mean(dim(‘ROW‘ ’COL‘)) # 计算某个城市点特定行、列的时间序列 o3_city ds[‘O3’].sel(ROW100, COL150, LAY0)4.2 基础分析与可视化从数据到图表得到数据后最基本的分析就是时空分布特征。例如绘制某一天PM2.5的平均浓度空间分布图import matplotlib.pyplot as plt import cartopy.crs as ccrs # 假设我们已经提取了某一天地表PM2.5的数据pm25_data fig plt.figure(figsize(10, 8)) ax plt.axes(projectionccrs.PlateCarree()) # 绘制填色图 mesh ax.pcolormesh(lon, lat, pm25_data, transformccrs.PlateCarree(), cmap‘viridis’) ax.coastlines() # 添加海岸线 ax.gridlines() # 添加网格线 plt.colorbar(mesh, label‘PM2.5 Concentration (μg/m³)’) plt.title(‘Surface PM2.5 on YYYY-MM-DD’) plt.show()这里lon和lat需要从网格描述文件或MCIP输出的经纬度文件中获取。Cartopy库可以方便地处理地图投影。除了空间分布时间序列分析也至关重要。你可以分析某个站点对应网格点污染物浓度的日变化、过程变化或者分析整个区域平均浓度的变化趋势并与观测数据进行对比验证。这是评估模拟性能的关键步骤。如果模拟值与观测值在量级和变化趋势上差异巨大就需要回头检查排放清单、气象输入或模型配置。4.3 高级诊断与过程分析看懂浓度场只是第一步。CMAQ的强大之处在于其过程分析Process Analysis PA和敏感性分析功能这能帮你回答“为什么”浓度会这样变化。过程分析PA工具可以将污染物浓度的净变化分解为若干个物理化学过程的贡献例如水平平流输送、垂直对流、气相化学转化、气溶胶过程、干沉降、湿沉降等。输出文件中的变量如ADVx_OXX方向平流对氧化剂贡献、CHEM_OX化学转化贡献等。通过分析这些过程项你可以判断在污染事件中究竟是本地化学生成主导还是区域输送更重要。敏感性分析通过DDM-3D或HDDM-3D模块可以定量计算某个源区或某种排放类型如交通源、工业源对目标区域浓度的贡献率即敏感性系数。这对于污染来源解析、制定减排策略具有直接的指导意义。解译这些高级输出需要更专业的知识。你需要清楚每个过程项的正负号含义是增加浓度还是减少浓度并能够将多个过程项的时间序列叠加与净浓度变化进行对比验证确保能量守恒各过程贡献之和约等于净变化。这通常是深入分析空气质量问题的核心。5. 常见问题排查与性能调优5.1 编译与运行阶段典型错误即使按照指南操作遇到报错也是家常便饭。这里记录几个我反复遇到的“坑”及其解决方法编译时提示“找不到 -lnetcdff”或类似链接错误问题这明确指向库路径问题。编译器在链接阶段找不到NetCDF Fortran库。排查首先确认NETCDF环境变量指向的目录下确实有lib/libnetcdff.so文件。其次检查LD_LIBRARY_PATH是否包含了该lib目录。最彻底的方法是进入CMAQ的编译目录手动检查Makefile或编译日志看-L和-l参数是否正确指向了你的库路径。有时CMAQ配置脚本没有正确导出所有路径需要你在编译前手动执行export LD_LIBRARY_PATH$NETCDF/lib:$LD_LIBRARY_PATH。运行时立即报错“APPLICATION TERMINATED DUE TO MPI ERROR”或直接挂起问题MPI并行任务启动失败。排查首先用which mpirun和mpirun --version确认MPI环境正常。其次检查运行脚本中设置的进程数NPCOL * NPROW是否超过了实际可用的CPU核心数。最后尝试先用一个非常小的网格和单进程串行模式运行排除是否是并行域分解导致的问题。如果串行能跑并行不能很可能是网格数不能被进程数整除或者MPI库版本与编译器不兼容。运行中途报错“ERROR OPENING FILE xxx”问题模型找不到某个输入文件。排查这是最常出现的运行时错误。仔细查看错误日志中给出的完整文件路径。然后检查该路径下的文件是否存在。检查文件权限是否可读。检查文件名是否与脚本中设置的一模一样包括大小写和扩展名。如果路径是相对路径确认模型启动时的当前工作目录是否正确。模拟结果中某些区域浓度出现极不合理的极高值如NaN或10^30问题数值计算不稳定导致“爆炸”。排查这通常与物理参数化方案或化学机制在特定气象条件下的兼容性有关。可以尝试缩短输出时间步长CTM_TSTEP。调整平流方案CTM_ADVEC_METHOD例如从“YAMO”换成“HSLT”。检查气象输入数据特别是垂直风速和边界层高度是否存在异常突变。从一个“干净”的初始条件通常使用清洁背景场开始模拟并给予几天的“spin-up”时间让模型达到稳定。5.2 性能调优与资源管理CMAQ模拟通常是计算和I/O密集型任务。合理的调优能显著缩短运行时间。并行效率使用NPCOL和NPROW进行二维域分解。理想情况下每个进程负责的网格数NCOLS/NPCOL * NROWS/NPROW应大致相等且不要太小通常不少于20x20以减少进程间通信开销。可以通过工具如mpiP或Intel VTune分析负载是否均衡。I/O优化输出文件是性能瓶颈之一。尽量减少不必要的输出变量和输出频率。如果可能使用CMAQ的“拼接”功能将每小时输出文件按天或按月拼接成一个大文件这比处理成千上万个小文件效率高得多。内存使用大型模拟如全国范围、高分辨率可能消耗大量内存。在编译时可以尝试调整优化级别如-O2而非-O3以降低内存占用或使用动态内存分配如果模型支持。监控运行时的内存使用情况确保没有超出物理内存导致交换swapping否则性能会急剧下降。使用高性能存储将输入/输出目录放在高性能并行文件系统如Lustre, GPFS上能极大缓解I/O压力尤其是在多任务同时运行时。最后保持工作目录的整洁。每次运行前清理旧的输出文件和临时文件。使用脚本自动化编译、运行和后处理流程并详细记录每次运行的配置参数和关键环境变量。建立一个稳定的、可复现的模拟环境是开展长期、可靠科学研究的基础。模型运行本身只是工具如何设计科学的模拟实验如何结合观测数据深入解译模拟结果背后的物理化学机制才是真正产生价值的开始。