简介：Flexpart 10.4的源码压缩包，面向大气科学、环境模拟方向的研究者，也适合需要自行编译科学计算软件的Linux用户。包内以Flexpart主程序源码为核心，同时整合libemos、HDF5、NetCDF、eccodes、grib_api、zlib、szip、jasper等关键依赖库源码，对照Ubuntu 18.04环境整理，重点解决源码编译过程中依赖繁多、版本匹配难、接口配置复杂等实际问题。资源共718个文件，压缩包约328.8MB，文件类型涵盖Fortran源程序、Shell编译脚本、HTML文档、PNG示意图、CSV及GRB气象数据、多个Makefile变体等，目录结构完整，便于按依赖解析和编译顺序逐项查阅。目前已有608人学习下载，特别适合具备一定Linux基础、希望彻底掌握Flexpart编译流程并理解依赖关系的研究者。包内含大量CONTROL示例控制文件、不同分辨率运行配置及历史运行日志，可辅助排查编译后运行错误，缩短环境搭建与调试周期。若只依赖现成二进制包，往往无法灵活调整编译选项，这份源码包正是解决该痛点的完整素材。
1. flexpart.tar.gz 是什么：一个拉格朗日粒子扩散模拟程序的源码包
拿到 flexpart.tar.gz 这个压缩包时，第一件事不是急着解压，而是先认清它是什么：一个拉格朗日粒子扩散模式的源代码包。它不会像桌面软件那样双击运行，需要自己完成编译、气象数据准备和参数配置，才能输出粒子轨迹、浓度分布这类结果。做污染溯源、事故释放评估、烟羽传播模拟的人会经常撞见这个名字。
它最有价值的点是能跑正向和反向模拟：正向问污染物接下来往哪飘，反向问观测点某时刻的气团来源在哪。这个包适合那些愿意把时间花在调试而不是黑匣子上的人。编译不复杂，但气象数据和时间参数一旦没对齐，结果就会毫无意义。下面我从解包到运行一步步写清楚，包括容易翻车的地方。
2. 从压缩包到可执行程序：解压、环境检查与编译
2.1 解压前先看包里有什么
tar.gz 这类文件最常见的翻车方式就是直接tar -xzf，解完发现目录层级不对，或者包根本不是源码而是数据，白白浪费几分钟。我一般先列目录，再决定怎么解压：
# 只列出压缩包内容，不解压 tar -tzf flexpart.tar.gz | head -n 30-t是 list，-z表示 gzip 压缩，-f指定文件名，head -n 30只看前 30 条。这样能快速确认里面有没有src/、docs/、examples/这些典型源码目录。看到确实有 Fortran 源文件和 Makefile 后，再真正解压。
# 创建目标目录，解压到固定路径并去掉顶层目录层级 mkdir -p /opt/models/flexpart tar -xzf flexpart.tar.gz -C /opt/models/flexpart --strip-components=1 ls -la /opt/models/flexpart-C指定解压目标目录，--strip-components=1是把包里的顶层目录剥掉，这样/opt/models/flexpart直接就是源文件根目录。如果你不确定有没有顶层目录，可以先跑上面第一条tar -tzf看一下；有了一层目录时加--strip-components=1，没有就不加，加多了会报路径错误。
解压后重点看两个东西：一是src/下有没有 Makefile，二是docs/或examples/里有没有现成的输入示例。很多老版本包自带的示例输入文件格式和源码是匹配的，比你自己从零写省事得多。
2.2 编译环境检查：编译器、netcdf、气象数据解码库
flexpart 的主体是 Fortran 代码，这些年版本虽多，依赖却相对固定：Fortran 编译器、netcdf 格式读写库、气象数据解码库。这三样缺一样，编译阶段就会卡住。我习惯先集中检查，而不是等 make 报错再一个个装：
# 检查编译器 which gfortran which mpif90 # 检查 netcdf 库，nc-config 对应 C 库，nf-config 对应 Fortran 库 which nc-config which nf-config # 检查气象数据解码库 which grib_api_config which eccodes-configwhich能直接看出命令在不在 PATH 里。注意 netcdf 分 C 和 Fortran 两层，fortran 层如果没装，make 时会出现NetCDF.h或netcdf.mod找不到，这是最常见的问题。
如果你用的是 Ubuntu/Debian 系，常见做法是apt install gfortran libnetcdf-dev libnetcdff-dev；CentOS/RHEL 系则对应dnf install gcc-gfortran netcdf-fortran-devel。具体包名跟发行版版本有关，但思路一致：C 库和 Fortran 库都要装。气象数据解码库在多数发行版里叫 libeccodes-dev，老版本包也可能依赖 grib_api。先装上，再走编译。
这里有个容易忽略的点：不同 flexpart 版本对编译器版本敏感。老代码里常带有 Fortran 77 风格的语法，新版 gfortran 默认编译选项更严格，可能报一堆 warning 甚至 error。我的建议是先按默认参数编，warning 能忽略就忽略，error 再处理，不要一开始就开-Wall -Wextra自找麻烦。
2.3 编译最小可执行程序
编译这一步，不同版本的 Makefile 命名很不统一：有的叫Makefile，有的叫makefile_gfortran，还有的是 CMake 工程。先进src/目录看一眼再动手：
cd /opt/models/flexpart/src ls Makefile* makefile* # 以 makefile_gfortran 为例，编译并保存日志 make -f makefile_gfortran 2>&1 | tee build.logmake -f是指定用哪个 Makefile。2>&1把标准错误也重定向到标准输出，方便日志捕获。tee build.log既把日志写到文件，也实时打印到屏幕。编译输出比较长时，日志是排查的救命稻草，不用在终端里翻来翻去。
如果包是新版本的 CMake 工程，就不要走 Makefile，直接按 CMake 的常规流程：
mkdir -p build && cd build cmake .. -DCMAKE_Fortran_COMPILER=gfortran make -j4 2>&1 | tee build.log-j4是并行编译，新机器可以先从 4 开始，编译出错时并行日志会比较乱，可以去掉-j4串行重试，错误信息更干净。
编译结束后，在src/或build/下会多出一个可执行文件，不同版本叫法不同，常见的是flexpart或flexpart_*。看到这个可执行文件，环境这一步就算过了。下一步不是急着跑大案例，而是配最小输入。
3. 跑通最小模拟：气象数据、pathnames 与 COMMAND 配置
3.1 气象数据：模型的口粮
flexpart 自己不会生成风场，它必须吃外部气象数据，常见格式是 GRIB 和 NetCDF。很多人第一次跑失败，不是代码编错，是气象数据没喂对。先不要考虑高精度数据，找一个覆盖模拟时段、包含 U/V/W 风场、温度和湍流的 NetCDF 文件就能起步。
拿到数据后，我第一件事是用工具看变量和时间维，而不是直接配置进模型：
ncdump -h met_2010_01.nc | head -n 40ncdump -h只打印 NetCDF 文件的头信息，不导出全部数据。重点看两点：时间单位是不是 UTC 标准时间；垂直层有没有覆盖你要模拟的高度。如果气象数据的顶层只有 3000 米，你却把释放源放在 5000 米，模拟不会报错，但粒子一上去就飞出数据范围，结果全是 0 或空轨迹。
分辨率上，全球再分析资料常见 0.25° 和 1° 两种。0.25° 的近场细节更好，计算量也更大；1° 数据适合跑长距离输送或反向溯源。第一次验证流程，用粗分辨率数据先把链路跑通，再回来换细分辨率，这个顺序能避免把参数错误和分辨率问题混在一起。
3.2 配置 pathnames 和 COMMAND 文件
flexpart 启动时依赖几个文本配置文件，核心是pathnames和COMMAND。pathnames告诉程序输入输出目录在哪儿，路径写错的表现是程序秒退或找不到文件。一个最小示例长这样：
/opt/models/flexpart/input /opt/models/flexpart/output第一行输入目录，第二行输出目录。注意这两行不能有空格开头，也不能带结尾多余的空格，最容易踩的就是从网页复制配置时混进不可见字符。我一般用cat -A pathnames检查隐藏符号，看到行尾只有$才算干净。
COMMAND 文件放的是模拟参数，不同版本变量名有差异，但基本逻辑不变。一个很朴素的配置是：
&COMMAND ibts = 20100115000000 icts = 20100115120000 iout = 1800 ts = 3600 ldirect = 1 maxpart = 50000 /ibts和icts是模拟起止时间，格式是 YYYYMMDDHHMISS，按 UTC 写，不是本地时间。iout是输出时间步长，单位秒，1800 就是每半小时输出一次。ts是粒子释放的时间间隔，单位秒，分析快速扩散过程时收窄，长距离慢过程可以放宽。maxpart是同时存在的最大粒子数，小测试案例 5 万就够，正式模拟再提到几十万。
气象数据和 COMMAND 之间的时间对齐，是第一次跑通最常见的问题。比如数据只有 00 时和 12 时两个时次，模拟却从 06 时开始，程序会用线性插值，但头尾缺少边界时次会直接报错。我一般在模拟开始前和结束后各多留一个时次的余量，稳妥很多。
3.3 启动运行并确认第一阶段正常
输入配置没问题，就可以直接跑。最小命令不需要开 MPI，先验证流程：
cd /opt/models/flexpart/input /opt/models/flexpart/src/flexpart > run.log 2>&1 tail -n 20 run.log ls -l /opt/models/flexpart/output> run.log 2>&1把所有输出重定向到文件，避免终端刷屏。跑完后看tail能确认已经执行到最后一步。输出目录里如果出现header.txt、partposit_end、grid_conc之类的文件，说明链路已经通了。
这里要多说一句：输出文件名和格式会因为编译选项和版本不同而变化。有的版本输出grid_conc二进制文件，有的版本输出 NetCDF，还有的会按模拟时间生成多个文件。第一次跑通不要纠结格式，先确认有文件生成并且文件大小不为 0，再进入下一步看内容。
4. 读输出：浓度网格、粒子位置与可信度检查
4.1 认识输出目录里的结果文件
flexpart 的结果不像普通地图软件那样输出一个导出图，它输出的是原始模拟数据，需要自己作图或做统计。第一次看到一堆没有扩展名的文件时容易懵，按常见类型分就清楚了：
| 文件 | 内容 | 典型用途 |
|---|---|---|
| header.txt | 网格范围、分辨率、释放源信息 | 确认模拟配置是否和你预期一致 |
| partposit_end | 模拟结束时每个粒子的位置、高度、状态 | 看最终轨迹落点 |
| grid_conc | 每个网格上的浓度场，二进制格式 | 做浓度分布图 |
| grid_time | 网格上的平均时间或停留时间 | 做迁移路径与扩散快慢分析 |
打开header.txt是每次跑完的习惯动作。这里记录着输出网格的经度范围、纬度范围、垂直层数，如果和你设想的模拟区域不一致，后面的结果全部不用看。
4.2 用 Python 把 NetCDF 输出读出来
如果你编译时开了 NetCDF 输出，最省事的后处理方式是直接用 xarray。这个库能延迟加载大文件，不会一下子把几十 GB 的浓度场压爆内存：
import xarray as xr # 打开 NetCDF 输出文件，只加载元数据 ds = xr.open_dataset("grid_conc.nc") print(ds) # 选取某个时刻的浓度场 c = ds["conc"].sel(time="2010-01-15T06:00:00", method="nearest") print(c.max())open_dataset默认是延迟加载，只有真正计算时才会从磁盘读数组，这对大文件很有用。sel(time=..., method="nearest")是找离目标时刻最近的时间层，避免因为时间步长和输出步长不一致选不到值。后面用to_netcdf或to_dataframe都可以继续导数据。
如果输出是老的 Fortran 二进制格式，我一般不手写struct.unpack去解析，因为字节序和自定义头结构的坑太多。常见做法是找包里自带的转换工具，生成 NetCDF 或文本格式后再用 Python 处理。
4.3 可信度检查：先判断结果能不能信
模拟能跑完，不等于结果能信。我每次拿到新案例，先做三个快速检查：粒子数守恒、浓度极大值位置、释放源附近形态。
粒子数守恒看partposit_end里有效粒子数是否接近释放粒子数。粒子飞出模拟区域外被标记为离开可以接受，但大量粒子莫名消失就不对。浓度极大值位置：正向模拟里，浓度最高区通常还在释放源附近或沿主风向一条带，如果极大值出现在释放源上风向几百公里外，多半是气象风场方向符号理解反了。释放源附近的浓度形态：没有风的情况下浓度向四周均匀扩散，有风时形成烟羽状，如果形态完全违背常识，先回去查气象数据。
5. flexpart 使用避坑：5 个常见的编译与运行问题
5.1 报错 netcdf.mod 或 NetCDF.h 找不到
现象：编译进行到一半，报错显示找不到 netcdf 相关头文件或模块文件。原因：系统只装了 netcdf-c，没装 netcdf-fortran。解决：安装 Fortran 接口库，发行版对应的是libnetcdff-dev或netcdf-fortran-devel。装完后用nf-config --includedir看看头文件路径，再把这个路径加到 Makefile 的FCFLAGS里。这是“库没装全”而不是“代码有问题”，别去改源码自找麻烦。
5.2 程序秒退，日志没有任何有效输出
现象：运行命令后终端过两秒就返回到提示符，run.log里只有一两行，输出目录为空。原因：pathnames文件里写的输入输出目录路径不对，或者输出目录没有写权限。解决：先ls -ld /opt/models/flexpart/output确认目录存在且可写，然后cat -A /opt/models/flexpart/input/pathnames检查两行路径有没有空格。路径末尾多一个空格都会让程序读错文件。
5.3 跑到某一时刻报 meteorological data 找不到
现象：模拟前面几小时正常，到某个时间突然报错，说找不到对应时间的气象场。原因：模拟窗口超出了气象数据覆盖范围，或者数据时间单位不是模型默认的 UTC。解决：用ncdump -v time查看气象数据时间范围，把ibts和icts改成在数据范围内，并且前后各留一个时次。做跨天模拟时尤其注意，很多再分析资料的日边界是 00 时，不是 08 时或 09 时。
5.4 浓度结果全为 0，图像上一片空白
现象：程序正常运行结束，但浓度场所有值都是 0。原因：释放高度在气象数据垂直层范围外，或释放源所在网格没有气象数据覆盖。解决：先ncdump -h看垂直层坐标最大值，把释放高度降到模式顶层以下。另一个常见原因是释放源经纬度写反了，比如把经度写成了负值，源位置落到模拟区域外。检查释放文件里的坐标，再看header.txt里的网格范围是否覆盖这个点。
5.5 并行运行结果不稳定，每次都不一样
现象：用 MPI 跑 8 进程，同一配置跑两次，粒子位置不完全一致。原因：粒子运动带有随机性，并行时每个进程的随机数序列和串行时不同，这是模型本身的行为，不是 bug。解决：单看粒子路径不要指望两条轨迹完全重合，要看浓度场和统计量。浓度分布如果整体稳定，只在小尺度上有差异，就可以接受。如果多进程和单进程结果差到离谱，优先检查气象数据读取是否按进程切分正确，再检查 MPI 部署的通信库版本。
6. 让 flexpart 跑得更贴近实际：释放层、粒子数与输出频率的调参习惯
配置跑通之后，真正决定结果可用性的是几个调参细节。我通常不会一上来就用大配置，而是先做小规模测试，看结果的稳定性和合理性，再放大量粒子。
释放高度方面，不要只做一个单点高度。真实烟囱排放有浮升高度，事故释放也有垂直范围，我的习惯是至少拆成 3 层：底部、中心高度、顶部，每层分配对应质量的粒子。这样做的好处是，风场随高度变化时，扩散形态能保留垂直结构的真实感。只放一个高度的典型结果是轨迹被拉成一条线，看着奇怪，其实不一定是模型错了，是释放源设置太粗暴。
粒子数方面，第一次跑通用 5 万，正式案例我一般用 50 万到 100 万。粒子数不是越多越好，超过一定量后浓度场变化很小，计算时间却成倍涨。一个判断方法：固定其他参数，把粒子数从 20 万提到 50 万，如果模拟区域内最大浓度变化小于 5%，说明粒子数已经够用。
输出频率iout按你的分析需要设。做小时级扩散分析，1800 秒输出可以兼顾平滑和文件体积；做逐日污染过程分析，3600 秒或更长也够。输出频率太密会让文件体积迅速增大，特别是跑 100 万粒子长时间模拟时，一个小失误就能写出几十 GB 中间结果。
我自己的固定流程是：先ncdump -h看气象数据和垂直层，再改pathnames和COMMAND，跑完先看header.txt的网格范围，最后把浓度场画出来和观测对比。有一次我把释放高度写到了 5000 米，结果粒子全被高空风切变拉成一条线，直到查气象数据才发现顶层只有 3000 米。自此每次跑之前都会花三十秒确认垂直层范围，这个习惯帮我避开了很多后续的返工。希望帮到你。
本文还有配套的精品资源，点击获取
