1. 项目概述从“单打独斗”到“协同作战”的海洋模拟革命如果你正在研究海岸带风暴潮、海浪对泥沙输运的影响或者想模拟海气相互作用对区域气候的反馈那你大概率已经听说过或者正在寻找一个能把这些过程“揉”在一起计算的工具。传统的海洋、大气、波浪、沉积物模型往往是独立运行的它们之间通过文件进行“离线”数据交换这种“各干各的”模式在处理海岸带这种多尺度、多过程强烈耦合的复杂环境时就显得力不从心了。COAWSTCoupled Ocean-Atmosphere-Wave-Sediment Transport耦合模式的出现就是为了解决这个核心痛点。简单来说COAWST不是一个从零开始编写的新模型而是一个精巧的“耦合框架”。它的核心思想是把几个经过长期验证、各自领域内的“明星”模型——比如海洋环流模型ROMS、大气模型WRF、波浪模型SWAN以及沉积物模型——通过一个名为MCTModel Coupling Toolkit的中间件以“在线”的方式紧密地捆绑在一起。所谓“在线”意味着这些模型在运行时是同步的可以每时每刻比如每隔几分钟或一个时间步长进行双向数据交换波浪模型把计算出的浪高、周期、破碎等信息实时传给海洋模型海洋模型据此调整海面粗糙度、混合过程海洋模型把海表温度实时传给大气模型大气模型据此调整热通量影响风场而风场又反过来驱动波浪和海洋。这种实时的、双向的互动使得模拟结果更加物理自洽尤其能捕捉到那些非线性反馈过程。我第一次接触COAWST是在做一个台风风暴潮增水项目时当时只用海洋模型总感觉模拟的极值水位和实测有系统性的偏差。后来引入COAWST让波浪模型在线参与计算后发现波浪辐射应力对水位的抬升作用在近岸非常显著一下子就把误差拉回来了。这个经历让我深刻体会到对于现代高精度的海岸带与区域海洋研究耦合模拟已经不是“锦上添花”而是“必不可少”的基础工具。然而这个工具的“安装”却成了拦住很多人的第一道门槛。它不像安装一个普通软件那样点几下“下一步”就行其过程涉及多个大型科学计算模型的源码获取、编译环境配置、库依赖解决以及耦合接口的编译堪称一场对研究者系统管理和编译能力的综合考验。接下来我就结合自己多次在超算集群和工作站上部署COAWST的经验为你拆解这个过程中的每一个关键步骤和避坑指南。2. 系统环境准备与核心依赖解析工欲善其事必先利其器。在动手下载COAWST那一堆源码之前我们必须先把它的“家”——也就是编译和运行环境——给搭建稳妥了。这一步没做好后续的编译过程会报出各种光怪陆离的错误让人无从下手。2.1 操作系统与编译器选择COAWST及其组件ROMS WRF SWAN主要是用Fortran和C语言编写的因此一个强大的Fortran/C编译器是核心。我的经验是Linux/Unix系统是唯一推荐的选择。无论是Ubuntu、CentOS、RedHat还是专业的超算集群如Slurm/PBS环境都完美支持。绝对不推荐在Windows原生环境下尝试那会陷入无尽的依赖库和路径问题地狱。对于Windows用户唯一可行的道路是使用WSL2Windows Subsystem for Linux。我曾在WSL2Ubuntu 20.04/22.04发行版上成功编译并运行过COAWST其体验与原生Linux几乎无异是个人电脑上最佳的学习和开发环境。编译器套装主流选择有两个Intel Compiler Suiteicc ifort英特尔编译器在优化科学计算代码特别是针对英特尔CPU架构时通常能产生效率更高的可执行文件。如果你的超算节点是英特尔至强Xeon系列这往往是首选。但它的许可证可能是个问题虽然学术机构常有免费版本。GNU Compiler Collectiongcc gfortran开源免费通用性最强。虽然在某些极端优化上可能稍逊于Intel编译器但对于绝大多数研究和应用来说其性能完全足够且避免了许可烦恼。我个人的建议是尤其是初次安装优先使用GNU编译器它能减少很多因编译器特性差异导致的诡异错误。你需要确保系统安装了完整的基础开发工具和编译器。在Ubuntu/Debian系系统上可以一键安装sudo apt-get update sudo apt-get install build-essential gfortran gcc g m4 cmake csh mpich libmpich-dev这里一口气安装了编译工具链、Fortran编译器、C/C编译器、必要的宏处理器m4、构建工具cmake、一个命令行shellcsh因为某些脚本依赖它以及MPI并行环境mpich。注意有些超算集群采用Intel MPI或OpenMPI你需要根据集群手册调整。2.2 并行计算环境MPI的选型与配置COAWST设计用于大规模并行计算无论是ROMS、WRF还是SWAN都支持基于MPIMessage Passing Interface的域分解。因此一个正确配置的MPI库至关重要。MPI实现选择常见的有OpenMPI、MPICH和Intel MPI。它们基本兼容但需要与编译器匹配。例如如果你用gfortran编译那么最好也用gfortran编译过的MPI库。幸运的是通过系统包管理器安装的mpich或openmpi通常已经做好了适配。关键检查点安装后在终端输入mpicc --version和mpif90 --version确认它们调用的确实是你的目标编译器gcc/gfortran或icc/ifort。如果显示“command not found”说明MPI没装好如果调用的编译器不对后续链接会出问题。超算环境注意在集群上通常需要通过模块加载module load特定的编译器和MPI组合例如module load intel/2022.1 mpi/intelmpi/2022.1。务必遵循集群管理员提供的指南这是最稳妥的路径。2.3 第三方依赖库的攻坚战这是整个环境准备中最繁琐但也最决定成败的一环。COAWST的各个组件需要一系列数学库和I/O库的支持。库名称主要用途是否必须安装建议与避坑指南NetCDF模型输入输出气象场、地形、初始场、结果的标准格式。是必须安装NetCDF-4版本并开启HDF5支持。这是最大的坑很多系统自带的或旧版的NetCDF是3.x版本不支持压缩、分块等高级特性WRF和ROMS会编译失败。务必从源码编译并先安装好HDF5和zlib。HDF5NetCDF-4的底层数据格式支持大文件、并行I/O。是为NetCDF-4建议安装1.12.x或1.10.x稳定版。编译时开启--enable-parallel如果你需要并行NetCDF。zlib, szip数据压缩库HDF5的依赖。是为HDF5通过包管理器安装通常即可sudo apt-get install zlib1g-dev libsz2Jasper用于处理GRIB2气象数据WRF输入常用。强烈推荐处理NCEP、ECMWF等机构的再分析资料时必备。LIBPNG图像处理库Jasper的依赖。是为Jasper同样通过包管理器安装。并行NetCDF (PnetCDF)在并行计算中提供高性能的NetCDF文件读写。可选但推荐对于大规模并行运行数百核以上能显著提升输出结果的效率。需要先安装MPI。我的实操心得源码编译依赖库的“黄金路径”我强烈建议将所有依赖库HDF5, NetCDF, PnetCDF等手动编译安装到一个独立的目录下比如/opt/libraries或$HOME/local。这样做的好处是环境完全可控与系统库隔离避免污染也方便在多版本间切换。编译时务必使用相同的编译器套装和MPI环境。一个典型的NetCDF编译命令如下# 假设编译器为gcc/gfortran MPI为mpich 安装路径为$HOME/local export CCmpicc export FCmpif90 export CXXmpicxx export F77mpif77 ./configure --prefix$HOME/local --enable-shared --enable-netcdf-4 --enable-dap --disable-doxygen make -j 8 # 使用8个核心并行编译加快速度 make check # 运行测试非常重要 make install安装后将$HOME/local/bin加入PATH将$HOME/local/lib加入LD_LIBRARY_PATH并在编译COAWST时通过环境变量或配置脚本指定这些库的路径。3. COAWST源码获取与框架结构解析环境准备好后我们就可以请出“主角”了。COAWST的代码托管在USGS的Git仓库中。3.1 使用Git克隆代码仓库最规范的方式是使用git进行克隆这便于后续更新和版本管理。git clone https://github.com/USGS-COAWST/COAWST.git cd COAWST克隆完成后你会看到一个结构清晰的目录树。理解这个结构对后续编译和修改至关重要COAWST/ ├── ROMS/ # 海洋模型ROMS的源代码 ├── WRF/ # 大气模型WRF的源代码注意需要额外下载 ├── SWAN/ # 波浪模型SWAN的源代码 ├── Projects/ # 预配置的多种应用案例如飓风、河口、西太区域 ├── SCRIP/ # 网格插值权重计算工具用于不同模型网格间的数据交换 ├── MCT/ # 模型耦合工具包Model Coupling Toolkit源代码 ├── coupling/ # 核心耦合器代码定义了模型间如何交换数据 └── build_scripts/ # 各种编译脚本是我们工作的重点这里有一个关键点WRF/目录初始是空的因为WRF模型的许可证限制其源码不能直接打包在COAWST仓库中。你需要单独去WRF官网下载相同版本的WRF源码并解压到这个WRF/目录下。务必注意版本兼容性COAWST的README或官网会指明支持的WRF版本例如WRF-4.x。3.2 选择与理解编译脚本进入build_scripts/目录你会看到一堆以.sh结尾的脚本例如build_roms.shbuild_swan.shbuild_wrf.sh和最终的build_coawst.sh。这些脚本是编译过程的“总指挥”。不要直接运行这些脚本是模板你需要复制一份并修改成你自己的配置。我通常的做法是cp build_coawst.sh my_build_coawst.sh chmod x my_build_coawst.sh vi my_build_coawst.sh # 使用你喜欢的文本编辑器进行编辑脚本核心配置解析打开脚本你需要关注和修改以下几个关键部分编译器与MPI设置找到类似export COMPILERintel或export USE_MPIon的行根据你的环境修改。如果是GNU编译器可能设为export COMPILERgfortran。依赖库路径这是最容易出错的地方。脚本中会有大量如export NETCDF_INCDIR/path/to/netcdf/include和export NETCDF_LIBDIR/path/to/netcdf/lib的变量。你必须将它们全部指向你之前亲手安装的那些库的路径。如果指向了系统自带的旧版库编译必败。应用配置Application Flagexport MY_CPP_FLAGS这一行定义了预处理器宏它决定了激活ROMS中的哪些物理过程。Projects/目录下的每个案例都有一个对应的宏比如-DINLET_TEST。你需要根据你要模拟的问题选择或参考一个最接近的。初次测试可以直接使用某个案例的配置。并行设置export USE_MPIon和export USE_MPIF90on必须开启。export which_MPImpich需要与你安装的MPI类型一致。4. 分步编译与耦合集成实战有了定制化的编译脚本我们就可以开始“组装”这个庞大的系统了。建议的编译顺序是先编译底层库和工具再编译各个模型最后编译耦合器。4.1 编译基础工具与耦合库编译SCRIP这是一个用于生成不同网格间插值权重的工具耦合前如果涉及不同分辨率的网格可能需要它。cd SCRIP ./configure make cd ..编译MCTModel Coupling Toolkit这是COAWST耦合的“神经系统”。进入MCT目录通常有一个Makefile或configure脚本。你需要编辑Makefile.conf或类似文件指定你的Fortran编译器mpif90和必要的编译选项。然后执行make。编译成功后会在lib/目录下生成libmct.a和libmpeu.a等静态库。4.2 编译各分量模型这是最耗时的部分但脚本已经帮我们自动化了大部分工作。编译ROMS./my_build_roms.sh # 这是你修改后的脚本这个脚本会进入ROMS/目录根据你的配置生成一个Linux/如果使用GNU或ifort/如果使用Intel的子目录并在其中进行编译。编译成功后会生成名为oceanG的可执行文件。密切关注终端输出任何“Error”或“Fatal”都会导致编译停止。最常见的错误是NetCDF库路径不对或版本不兼容。编译SWAN./my_build_swan.shSWAN的编译相对独立。成功后会在SWAN/目录下生成swan.exe等可执行文件。编译WRF这是最复杂的一步因为WRF本身就是一个庞大的系统。你需要确保已经将正确的WRF源码放入了WRF/目录。然后运行./my_build_wrf.sh这个脚本会调用WRF自带的configure和compile流程。在配置时它会交互式地让你选择并行类型如dmpar 分布式内存并行和编译器组合。你必须选择与你环境匹配的选项例如对于GNUMPICH可能选择35. Linux x86_64, gfortran (dmpar)。编译过程很长会先后生成real.exe理想数据前处理、wrf.exe主程序等。4.3 编译耦合器并生成最终可执行文件当所有分量模型都编译成功后最后一步就是将它们“粘合”起来。编译耦合器进入coupling/目录这里有一个Makefile。你需要修改这个Makefile确保它指向正确的MCT库路径、NetCDF路径以及各分量模型的目录和库文件。然后运行make。如果一切顺利会生成耦合器相关的模块文件.mod和库。编译COAWST主程序回到COAWST根目录运行最终的构建脚本./my_build_coawst.sh这个脚本会执行一个“超级编译”它链接之前编译好的ROMS、WRF、SWAN的代码对象以及MCT和耦合器库生成一个单一的可执行文件通常命名为coawstM。看到类似“Executable created: coawstM”的提示就标志着整个COAWST系统编译安装成功了5. 编译问题排查与实战经验录即便按照指南操作编译过程中也大概率会遇到各种错误。下面是我总结的常见“坑”及其解决方案。5.1 依赖库版本冲突与路径问题症状编译早期报错提示找不到netcdf.mod、nf-config命令不存在或链接时出现undefined reference to ‘nc_open_par’等。诊断这几乎100%是NetCDF库的问题。首先用nc-config --all或nf-config --all命令检查当前生效的NetCDF版本和特性。确认它是否是NetCDF-4版本并且支持并行I/O–has-parallel4 yes。解决彻底清理旧的、错误的库路径。在编译脚本和环境变量如LD_LIBRARY_PATH中确保只指向你手动编译安装的那个正确版本。重新编译NetCDF在configure时务必加上--enable-parallel-tests如果需并行并仔细查看输出摘要确认所有需要的特性NetCDF-4 HDF5 Parallel都是“yes”。使用make distclean彻底清理ROMS/WRF的编译中间文件再重新编译。5.2 MPI与编译器不匹配症状编译通过但链接时失败报MPI相关函数未定义或者运行时直接崩溃提示MPI_Init错误。诊断检查mpif90 --show或mpif90 -v看它底层调用的是gfortran还是ifort。确保这与你在编译脚本中设置的COMPILER变量一致。解决在超算上使用module purge然后重新加载一套匹配的编译器MPI模块。在个人环境考虑卸载重装或者通过设置环境变量MPIF90gfortran如果MPI是用gfortran编译的来显式指定。5.3 WRF编译特有的难题症状在编译WRF的real.exe或wrf.exe时在链接阶段失败错误信息涉及-lgfortran或-limf等。诊断WRF的configure脚本选择的编译器选项可能不完全适合你的环境。或者你之前用Intel编译器编译过依赖库但现在用GNU编译WRF反之亦然导致库不兼容。解决进入WRF目录执行./clean -a进行彻底清理。重新运行./configure谨慎选择与你整体环境匹配的选项。手动编辑WRF/configure.wrf文件。这是一个关键技巧找到类似FCBASEOPTS的行你可能需要移除或添加一些编译器标志。例如对于较新的gfortran有时需要添加-fallow-argument-mismatch来避免某些子程序参数类型警告被视为错误。保持整个工具链的一致性所有库和模型尽量用同一套编译器编译。5.4 内存不足与进程被杀症状编译过程中特别是编译WRF或最后链接coawstM时系统卡住然后进程消失终端提示Killed。诊断这通常是内存RAM不足导致的。链接一个包含数千万行代码的巨型可执行文件需要大量的内存。解决在个人电脑上关闭其他占用内存大的程序。在编译命令前添加ulimit -s unlimited解除shell的栈大小限制。如果是在虚拟内存很小的虚拟机上考虑增加交换空间swap。最根本的在性能更强的服务器或超算计算节点上进行编译。6. 首次运行测试与后续工作流编译出coawstM只是万里长征第一步。让这个庞然大物跑起来并产出合理的结果需要正确的输入文件和配置。6.1 使用预置案例进行测试Projects/目录下的每个子目录都是一个完整的测试案例包含了运行所需的所有输入文件模型网格grid.nc、初始条件和边界条件ini.ncbry.nc、大气强迫场wrfinput*、波浪边界、以及最重要的配置文件ocean.inROMS、namelist.inputWRF和swan_inputSWAN。复制案例选择一个相对简单的案例比如Inlet_test将其整个目录复制到你的运行空间。cp -r Projects/Inlet_test /your/run/directory cd /your/run/directory准备可执行文件将编译好的coawstM可执行文件复制到该目录。运行通过MPI命令运行。例如使用4个进程mpirun -np 4 ./coawstM如果配置正确你会看到终端开始滚动输出各个模型WRF ROMS SWAN的时间步信息耦合器交换数据的日志以及最终正常结束的信息。同时目录下会生成一系列输出文件ocean_his.ncwrfout*swan.*等。6.2 配置文件的深度定制成功运行测试案例后你就需要为自己的研究区域和科学问题定制这些配置文件了。这是COAWST应用的核心技能。ocean.in(ROMS)这是海洋模型的“大脑”。你需要设置VARNAME输入文件名称。Lm Mm N网格在XI- ETA-和垂直方向上的点数。DT海洋模型时间步长秒。这是一个极其关键的参数必须满足CFL稳定性条件。通常根据网格最小分辨率和流速估算。NTIMES总的时间步数决定了模拟时长。NRREC输入/输出文件的记录索引用于热启动或从历史文件继续。大量的物理参数选项如垂直混合方案、底摩擦类型、温盐方程参数等。namelist.input(WRF)控制大气模型。核心的domains部分设置网格层数、分辨率、区域大小physics部分选择积云参数化、微物理、边界层等方案dynamics部分设置时间步长time_step 通常为网格距的6倍以秒为单位。耦合频率设置在ocean.in中有一组以COUPLING开头的参数如COUPLING_MODECOUPLING_INTERVAL。COUPLING_INTERVAL定义了海洋和大气之间交换数据的时间间隔以海洋时间步为单位。这个间隔需要仔细权衡太短如每步都耦合计算开销巨大太长则失去了紧密耦合的意义。通常设置为几分钟到几十分钟物理时间需要与模型的时间步长相匹配。6.3 输入文件制备工作流要运行自己的案例你需要制备一套完整的、相互匹配的输入文件。这是一个系统性的工作网格生成使用ROMS/Grid/下的工具如seagridgridgen或更通用的工具如ESMF_RegridWeightGen结合SCRIP生成你的海洋模型网格文件grid.nc并同时考虑WRF所需的大气网格通常通过WPS地理处理系统生成。初始场与边界场从全球或区域再分析数据如HYCOM ROMS父模型输出 ERA5 CFSv2中利用ROMS的make_ini.pymake_bry.py等工具插值得到海洋初始和边界条件。大气部分则使用WPS和real.exe处理全球气象数据生成wrfinput_d01。配置文件调整根据新网格的尺寸和参数更新所有配置文件ocean.innamelist.input等中的相应设置。测试运行与调试先进行短时间如模拟几小时的测试运行检查模型是否稳定输出是否合理。查看日志文件log.ww3ocean.log等和标准错误输出排查可能存在的初始场不平衡、边界条件不合理、时间步长过大等问题。整个COAWST的安装和应用是一个典型的“系统工程”。它考验的不仅仅是领域科学知识更是系统管理、软件编译和问题排查的综合能力。每一次成功的编译和运行都是对这些技能的一次提升。当你看到耦合模型成功地模拟出台风过境时海气界面上复杂的能量交换过程时你会觉得之前所有的折腾都是值得的。这个过程没有捷径唯有多动手、多记录、多思考逐步积累起属于自己的经验库。