手把手教你用Linux编译FLEXPART-WRF从气象数据到污染物扩散模拟的完整流程当我们需要模拟大气中污染物的传输与扩散过程时FLEXPART-WRF无疑是一个强大的工具。作为一名长期从事大气环境模拟的研究者我深知从零开始搭建这个模型的痛苦——依赖库缺失、编译报错、参数配置不当等问题层出不穷。本文将带你一步步攻克这些难关完成从WRF气象场数据准备到最终污染物扩散模拟的全流程。1. 环境准备与依赖库安装在开始编译FLEXPART-WRF前我们需要确保Linux系统具备所有必要的依赖环境。以下是经过多次实践验证的基础配置方案1.1 系统要求与基础工具推荐使用Ubuntu 20.04 LTS或CentOS 7系统至少4GB内存和50GB磁盘空间。首先安装基础开发工具# Ubuntu/Debian sudo apt-get update sudo apt-get install -y build-essential gfortran m4 csh make git # CentOS/RHEL sudo yum groupinstall -y Development Tools sudo yum install -y gcc-gfortran m4 tcsh make git1.2 关键依赖库安装FLEXPART-WRF的核心依赖包括NetCDF、HDF5和Jasper库。这里提供一个经过优化的安装脚本# 安装zlib wget https://zlib.net/zlib-1.2.13.tar.gz tar -xzf zlib-1.2.13.tar.gz cd zlib-1.2.13 ./configure --prefix/usr/local make sudo make install # 安装HDF5 wget https://support.hdfgroup.org/ftp/HDF5/releases/hdf5-1.12/hdf5-1.12.2/src/hdf5-1.12.2.tar.gz tar -xzf hdf5-1.12.2.tar.gz cd hdf5-1.12.2 ./configure --prefix/usr/local --enable-fortran --enable-hl make sudo make install # 安装NetCDF-C和NetCDF-Fortran wget https://downloads.unidata.ucar.edu/netcdf-c/4.9.0/netcdf-c-4.9.0.tar.gz tar -xzf netcdf-c-4.9.0.tar.gz cd netcdf-c-4.9.0 ./configure --prefix/usr/local make sudo make install wget https://downloads.unidata.ucar.edu/netcdf-fortran/4.6.0/netcdf-fortran-4.6.0.tar.gz tar -xzf netcdf-fortran-4.6.0.tar.gz cd netcdf-fortran-4.6.0 ./configure --prefix/usr/local make sudo make install提示编译大型科学计算软件时建议设置以下环境变量到~/.bashrc中export LD_LIBRARY_PATH/usr/local/lib:$LD_LIBRARY_PATH export CPPFLAGS-I/usr/local/include export LDFLAGS-L/usr/local/lib2. 获取与准备WRF气象数据FLEXPART-WRF需要WRF模式输出的气象场作为输入。以下是处理WRF数据的实战经验2.1 WRF数据获取与验证WRF输出通常包含多个时间步的netCDF文件。使用ncdump检查数据完整性ncdump -h wrfout_d01_2020-01-01_00:00:00关键变量必须包含U/V/W (风速)T (温度)P/PB (气压)QVAPOR (水汽混合比)PH/PHB (位势高度)2.2 时间连续性处理使用NCO工具处理时间不连续的数据# 合并多个时间步文件 ncrcat wrfout_d01_2020-01-01_* wrfout_combined.nc # 修复时间变量 ncap2 -s TimesTimes wrfout_combined.nc fixed_wrfout.nc3. FLEXPART-WRF编译实战3.1 源码获取与准备从官方仓库获取最新代码git clone https://github.com/Flexpart/flexpart-wrf.git cd flexpart-wrf3.2 编译配置调整修改makefile.mom关键参数# 设置NetCDF路径 NETCDF /usr/local # 根据系统选择编译器选项 F90 gfortran F90_OPTS -O3 -fconvertbig-endian -fPIC LD_OPTS -lnetcdff -lnetcdf -lhdf5_hl -lhdf5 -lz3.3 常见编译问题解决下表总结了典型错误及解决方案错误类型表现特征解决方案NetCDF缺失netcdf.inc: No such file检查NETCDF路径确认netcdf-fortran安装内存不足segmentation fault增加swap空间或优化网格分辨率时间格式错误Time stamps not monotonic使用ncrcat重新处理输入文件编译成功后你会得到可执行文件FLEXPART_WRF可通过以下命令验证./FLEXPART_WRF --version4. 运行污染物扩散案例4.1 案例配置详解创建一个简单的点源排放案例编辑pathnames文件/home/user/flexpart-wrf/options/RELEASES /home/user/flexpart-wrf/options/COMMAND /home/user/flexpart-wrf/options/OUTGRIDRELEASES文件配置示例 RELEASE SPECIFICATIONS: YYYYMMDD HHMMSS YYYYMMDD HHMMSS Lon Lat Height #part 20200101 000000 20200101 010000 116.391 39.907 100.0 100004.2 运行与监控启动模拟并监控资源使用./FLEXPART_WRF run.log 21 tail -f run.log # 监控内存使用 watch -n 1 free -h4.3 结果可视化使用Python进行后处理import xarray as xr import matplotlib.pyplot as plt ds xr.open_dataset(grid_conc_000.nc) ds[spec001_mr].isel(time0).plot() plt.savefig(pollution_dist.png)5. 性能优化技巧经过多次项目实践我总结出以下提升FLEXPART-WRF效率的方法并行计算使用OpenMP加速export OMP_NUM_THREADS4 ./FLEXPART_WRF内存管理对于大型模拟调整COMMAND中的参数COMMAND lsubgrid .TRUE., lconvection .TRUE., ipout 1, maxpart 1000000, /输出优化只保存必要的变量和频率减少IO压力6. 实际应用案例北京地区PM2.5扩散模拟配置一个真实的污染事件模拟源参数设置位置39.9°N, 116.4°E高度50-200m建筑排放持续时间24小时排放强度10kg/h气象场处理# 使用WRF-Chem输出增强化学过程 ncks -v U,V,W,T,P,PB,PH,PHB,QVAPOR,PM2_5 wrfchemi_d01_20200101 wrfinput_d01敏感性分析对比不同边界层参数化方案测试干湿沉积参数影响
手把手教你用Linux编译FLEXPART-WRF：从气象数据到污染物扩散模拟的完整流程
手把手教你用Linux编译FLEXPART-WRF从气象数据到污染物扩散模拟的完整流程当我们需要模拟大气中污染物的传输与扩散过程时FLEXPART-WRF无疑是一个强大的工具。作为一名长期从事大气环境模拟的研究者我深知从零开始搭建这个模型的痛苦——依赖库缺失、编译报错、参数配置不当等问题层出不穷。本文将带你一步步攻克这些难关完成从WRF气象场数据准备到最终污染物扩散模拟的全流程。1. 环境准备与依赖库安装在开始编译FLEXPART-WRF前我们需要确保Linux系统具备所有必要的依赖环境。以下是经过多次实践验证的基础配置方案1.1 系统要求与基础工具推荐使用Ubuntu 20.04 LTS或CentOS 7系统至少4GB内存和50GB磁盘空间。首先安装基础开发工具# Ubuntu/Debian sudo apt-get update sudo apt-get install -y build-essential gfortran m4 csh make git # CentOS/RHEL sudo yum groupinstall -y Development Tools sudo yum install -y gcc-gfortran m4 tcsh make git1.2 关键依赖库安装FLEXPART-WRF的核心依赖包括NetCDF、HDF5和Jasper库。这里提供一个经过优化的安装脚本# 安装zlib wget https://zlib.net/zlib-1.2.13.tar.gz tar -xzf zlib-1.2.13.tar.gz cd zlib-1.2.13 ./configure --prefix/usr/local make sudo make install # 安装HDF5 wget https://support.hdfgroup.org/ftp/HDF5/releases/hdf5-1.12/hdf5-1.12.2/src/hdf5-1.12.2.tar.gz tar -xzf hdf5-1.12.2.tar.gz cd hdf5-1.12.2 ./configure --prefix/usr/local --enable-fortran --enable-hl make sudo make install # 安装NetCDF-C和NetCDF-Fortran wget https://downloads.unidata.ucar.edu/netcdf-c/4.9.0/netcdf-c-4.9.0.tar.gz tar -xzf netcdf-c-4.9.0.tar.gz cd netcdf-c-4.9.0 ./configure --prefix/usr/local make sudo make install wget https://downloads.unidata.ucar.edu/netcdf-fortran/4.6.0/netcdf-fortran-4.6.0.tar.gz tar -xzf netcdf-fortran-4.6.0.tar.gz cd netcdf-fortran-4.6.0 ./configure --prefix/usr/local make sudo make install提示编译大型科学计算软件时建议设置以下环境变量到~/.bashrc中export LD_LIBRARY_PATH/usr/local/lib:$LD_LIBRARY_PATH export CPPFLAGS-I/usr/local/include export LDFLAGS-L/usr/local/lib2. 获取与准备WRF气象数据FLEXPART-WRF需要WRF模式输出的气象场作为输入。以下是处理WRF数据的实战经验2.1 WRF数据获取与验证WRF输出通常包含多个时间步的netCDF文件。使用ncdump检查数据完整性ncdump -h wrfout_d01_2020-01-01_00:00:00关键变量必须包含U/V/W (风速)T (温度)P/PB (气压)QVAPOR (水汽混合比)PH/PHB (位势高度)2.2 时间连续性处理使用NCO工具处理时间不连续的数据# 合并多个时间步文件 ncrcat wrfout_d01_2020-01-01_* wrfout_combined.nc # 修复时间变量 ncap2 -s TimesTimes wrfout_combined.nc fixed_wrfout.nc3. FLEXPART-WRF编译实战3.1 源码获取与准备从官方仓库获取最新代码git clone https://github.com/Flexpart/flexpart-wrf.git cd flexpart-wrf3.2 编译配置调整修改makefile.mom关键参数# 设置NetCDF路径 NETCDF /usr/local # 根据系统选择编译器选项 F90 gfortran F90_OPTS -O3 -fconvertbig-endian -fPIC LD_OPTS -lnetcdff -lnetcdf -lhdf5_hl -lhdf5 -lz3.3 常见编译问题解决下表总结了典型错误及解决方案错误类型表现特征解决方案NetCDF缺失netcdf.inc: No such file检查NETCDF路径确认netcdf-fortran安装内存不足segmentation fault增加swap空间或优化网格分辨率时间格式错误Time stamps not monotonic使用ncrcat重新处理输入文件编译成功后你会得到可执行文件FLEXPART_WRF可通过以下命令验证./FLEXPART_WRF --version4. 运行污染物扩散案例4.1 案例配置详解创建一个简单的点源排放案例编辑pathnames文件/home/user/flexpart-wrf/options/RELEASES /home/user/flexpart-wrf/options/COMMAND /home/user/flexpart-wrf/options/OUTGRIDRELEASES文件配置示例 RELEASE SPECIFICATIONS: YYYYMMDD HHMMSS YYYYMMDD HHMMSS Lon Lat Height #part 20200101 000000 20200101 010000 116.391 39.907 100.0 100004.2 运行与监控启动模拟并监控资源使用./FLEXPART_WRF run.log 21 tail -f run.log # 监控内存使用 watch -n 1 free -h4.3 结果可视化使用Python进行后处理import xarray as xr import matplotlib.pyplot as plt ds xr.open_dataset(grid_conc_000.nc) ds[spec001_mr].isel(time0).plot() plt.savefig(pollution_dist.png)5. 性能优化技巧经过多次项目实践我总结出以下提升FLEXPART-WRF效率的方法并行计算使用OpenMP加速export OMP_NUM_THREADS4 ./FLEXPART_WRF内存管理对于大型模拟调整COMMAND中的参数COMMAND lsubgrid .TRUE., lconvection .TRUE., ipout 1, maxpart 1000000, /输出优化只保存必要的变量和频率减少IO压力6. 实际应用案例北京地区PM2.5扩散模拟配置一个真实的污染事件模拟源参数设置位置39.9°N, 116.4°E高度50-200m建筑排放持续时间24小时排放强度10kg/h气象场处理# 使用WRF-Chem输出增强化学过程 ncks -v U,V,W,T,P,PB,PH,PHB,QVAPOR,PM2_5 wrfchemi_d01_20200101 wrfinput_d01敏感性分析对比不同边界层参数化方案测试干湿沉积参数影响
