1. 项目概述：COAWST耦合模式是什么？

如果你在海洋、大气或海岸工程领域摸爬滚打过，肯定对WRF、ROMS、SWAN这些名字不陌生。WRF负责大气模拟，ROMS专攻海洋环流，SWAN则擅长计算海浪。过去，我们常常是“单打独斗”，用WRF算完风场，再把结果手动喂给ROMS去算海流，过程繁琐不说，最关键的是忽略了系统间实时的、双向的相互作用。比如，海面温度变化会立刻影响上空的大气边界层，而风应力的改变又会反过来驱动海浪和海流，这种耦合效应在台风、风暴潮等极端事件模拟中至关重要。COAWST（The Coupled Ocean-Atmosphere-Wave-Sediment Transport）系统就是为了解决这个问题而生的。

简单来说，COAWST不是一个全新的模型，而是一个精巧的“耦合框架”。它的核心价值在于，通过一套称为MCT（Model Coupling Toolkit）的库，将WRF、ROMS、SWAN这三个各自领域的“王牌”模型，像拼乐高一样紧密地连接在一起，让它们能够在运行时实时交换数据（如风场、海表温度、波浪场、流场等），实现真正意义上的“海洋-大气-波浪”双向耦合模拟。这对于研究海气相互作用、台风强度与路径预报、海岸带侵蚀与沉积物输运等复杂过程，是一个质的飞跃。

我第一次接触COAWST是在一个风暴潮增水的项目里，当时用单独的大气和海洋模型，总是对峰值水位的模拟差那么一点。后来切换到COAWST耦合模式，模拟结果与实测数据的吻合度显著提升，尤其是台风眼壁附近极端风应力与波浪破碎共同作用的区域，耦合模拟的优势体现得淋漓尽致。这套系统虽然安装和配置的门槛不低，但一旦跑通，它所带来的物理过程完整性和模拟精度，是传统离线耦合方式难以比拟的。接下来，我就结合自己多次“踩坑”的经验，带你从零开始，理清COAWST的框架，并一步步完成它的安装与编译。

2. 核心架构与耦合原理拆解

在动手安装之前，我们必须先理解COAWST是怎么工作的。知其然更要知其所以然，这能帮你在后续遇到各种编译错误或运行时异常时，快速定位问题根源。

2.1 组件模型与角色分工

COAWST耦合系统的三大核心支柱，每一个都是久经沙场的开源社区明星：

1. WRF（Weather Research and Forecasting Model）：负责大气部分。它接收来自海洋表面的温度、粗糙度（受波浪影响）等信息，计算出风场、气压、热通量等，再将这些结果传递给海洋和波浪模型。
2. ROMS（Regional Ocean Modeling System）：负责海洋环流部分。它接收来自大气的风应力、热通量、淡水通量，以及来自波浪模型的波浪参数（如斯托克斯漂流、波浪破碎引起的混合增强），计算出三维的海流、温度、盐度分布，并将海表温度等信息反馈给大气。
3. SWAN（Simulating WAves Nearshore）：负责海浪部分。它接收来自大气的风场和来自海洋的流场、水深，计算出波浪谱、波高、波向、波长等，并将波浪引起的海面粗糙度、辐射应力反馈给大气和海洋。

除了这“三巨头”，COAWST还整合了沉积物输运模块（通常基于ROMS的沉积物模块），用于模拟泥沙的侵蚀、输运和沉积过程，这对于研究海岸地貌演变至关重要。

2.2 耦合器MCT的核心作用

那么，这三个独立的模型如何“对话”呢？这就是MCT（Model Coupling Toolkit）的功劳。你可以把MCT想象成一个高效的“数据交换中心”或“通信总线”。它的核心功能包括：

• 并行数据路由：在分布式内存（MPI）并行环境下，每个模型可能运行在不同的处理器核上，甚至物理分布在不同计算节点上。MCT负责将来自WRF某个进程的数据，精准地路由到需要该数据的ROMS或SWAN的对应进程上。
• 网格插值：WRF、ROMS、SWAN通常使用不同的计算网格（水平网格、垂直层、网格 staggering方式都可能不同）。MCT集成了网格重映射算法（如双线性插值、守恒性插值），在交换数据前自动进行空间插值，确保物理量从一个模型的网格准确地映射到另一个模型的网格上。
• 时间同步：各个模型有自己的时间积分步长（比如WRF可能60秒，ROMS10秒，SWAN20秒）。MCT管理着耦合的时间步长，控制着数据交换的频率。通常，耦合步长是各模型步长的公约数或根据物理过程最快时间尺度来确定，确保相互作用能以足够的精度被解析。

COAWST的耦合是双向、同步的。在一个耦合时间步内，大致流程如下：首先，各模型并行运行自己的一小段时间；然后，在约定的耦合时刻暂停，通过MCT交换最新的边界场数据；接着，各模型接收来自其他模型的新数据，并将其作为边界条件或强迫场，继续下一段的积分。如此循环往复。

2.3 安装前的关键决策点

理解架构后，安装前有几个关键选择需要你根据自身研究需求和计算环境来决定：

• 耦合范围：COAWST支持灵活的耦合配置。你可以只耦合大气-海洋（WRF+ROMS），或大气-波浪（WRF+SWAN），或海洋-波浪（ROMS+SWAN），当然也可以三者全耦合。安装时需要编译你所需的组件。
• 编译方式：COAWST主要提供两种构建方式：
  1. 使用build_roms.csh脚本：这是最传统也是目前最主流的方式。该脚本是一个功能强大的Tcsh脚本，通过交互式或配置文件的方式，让你选择编译器、库路径、激活哪些组件和功能，然后自动完成所有子模型的下载、打补丁和编译。这是新手推荐的首选路径。
  2. 使用CMake（较新，仍在完善）：社区也在推动CMake的构建方式，这对于管理复杂依赖和跨平台编译更有优势。但对于COAWST这种高度定制化的系统，CMake的支持可能不如传统脚本成熟，遇到奇怪问题的概率稍高。除非你对CMake非常熟悉且有特定需求，否则建议先使用脚本方式。
• 依赖库：这是安装过程中最大的“拦路虎”。COAWST编译依赖于一系列第三方库，主要包括：
  • MPI（如OpenMPI, MPICH）：用于并行计算，必须安装。
  • NetCDF库（包括C和Fortran接口）：用于模型的输入输出，必须安装。特别注意：ROMS/COAWST对NetCDF的版本，特别是Fortran接口的兼容性非常敏感。通常推荐使用NetCDF-C 4.8.x 配合 NetCDF-Fortran 4.5.x 或 4.6.x 的稳定组合。较新的版本（如NetCDF-Fortran 4.6.0+）可能需要额外的编译选项。
  • HDF5库：NetCDF4格式的底层依赖，通常与NetCDF一起安装。
  • 其他可选库：如JASPER（用于GRIB2支持）、PGI/NVIDIA HPC SDK（如果使用PGI或NV Fortran编译器）等。

我的经验是，在开始安装COAWST主体之前，花足够的时间确保这些依赖库（尤其是NetCDF-Fortran）在你的系统上被正确安装和配置，可以节省后面大量的调试时间。一个常见的做法是使用Spack、EasyBuild等HPC包管理工具来安装这些依赖，或者在自己的家目录下用源码编译一套统一版本的库。

3. 系统环境准备与依赖库安装

这一节是实战的开始，也是最容易出错的地方。我将以在Linux系统（如Ubuntu 20.04/22.04或CentOS 7/8）上，使用GNU编译器套件（gcc, gfortran, g++）和OpenMPI为例，详细说明环境搭建。假设你已经在高性能计算集群或个人工作站上拥有了一个干净的Linux环境。

3.1 基础开发环境与编译器

首先，通过系统包管理器安装必要的编译工具和基础库。

对于Ubuntu/Debian系统：

sudo apt-get update sudo apt-get install -y build-essential gfortran gcc g++ m4 cmake csh tcsh git curl wget

对于CentOS/RHEL系统：

sudo yum groupinstall -y "Development Tools" sudo yum install -y gcc-gfortran gcc-c++ m4 cmake tcsh git curl wget

注意：COAWST的构建脚本主要使用csh或tcsh，所以必须安装。如果你习惯bash，在运行脚本时使用tcsh script.csh即可。

接下来，安装MPI库。这里选择OpenMPI：

# Ubuntu/Debian sudo apt-get install -y openmpi-bin libopenmpi-dev # CentOS/RHEL sudo yum install -y openmpi openmpi-devel

安装后，将MPI的库路径加入环境变量。通常可以添加以下内容到你的~/.bashrc（如果你用bash）或~/.cshrc（如果你用csh）文件中：

# 对于OpenMPI，查找mpif90的路径，通常是 /usr/lib/openmpi 或 /usr/local/openmpi export MPI_ROOT=/usr/lib/openmpi # 请根据实际路径修改 export PATH=$MPI_ROOT/bin:$PATH export LD_LIBRARY_PATH=$MPI_ROOT/lib:$LD_LIBRARY_PATH

然后执行source ~/.bashrc使环境变量生效。验证安装：which mpif90和mpif90 --version应能正确显示信息。

3.2 NetCDF和HDF5库的源码编译安装

这是最关键的一步。我强烈建议你自己用源码编译NetCDF-C和NetCDF-Fortran，而不是完全依赖系统包。这能确保版本匹配，并且将所有库安装到同一个前缀路径下，方便管理。

我们创建一个专门的目录来存放源码和安装库，例如/opt/coawst_libs（需要sudo权限）或$HOME/local/coawst_libs（用户目录下，推荐）。

mkdir -p $HOME/src $HOME/local/coawst_libs cd $HOME/src

第一步：安装HDF5HDF5是NetCDF4的依赖。下载稳定版本（如1.12.x或1.14.x）：

wget https://support.hdfgroup.org/ftp/HDF5/releases/hdf5-1.14/hdf5-1.14.4/src/hdf5-1.14.4.tar.gz tar -xzvf hdf5-1.14.4.tar.gz cd hdf5-1.14.4

配置和编译。这里启用并行I/O支持（--enable-parallel），因为COAWST模型通常并行运行。--prefix指定安装路径。

./configure --prefix=$HOME/local/coawst_libs --enable-parallel --enable-shared=no --enable-hl CC=mpicc CXX=mpicxx FC=mpif90 make -j 8 # 使用8个核心并行编译，根据你的CPU调整 make check # 运行测试（可选，但推荐） make install

第二步：安装NetCDF-CNetCDF-C依赖HDF5和zlib等。确保上一步的HDF5安装路径已加入环境变量。

cd $HOME/src wget https://github.com/Unidata/netcdf-c/archive/refs/tags/v4.9.2.tar.gz -O netcdf-c-4.9.2.tar.gz tar -xzvf netcdf-c-4.9.2.tar.gz cd netcdf-c-4.9.2

配置时，通过CPPFLAGS和LDFLAGS指定HDF5的头文件和库路径。

export HDF5_DIR=$HOME/local/coawst_libs ./configure --prefix=$HOME/local/coawst_libs --disable-dap --enable-parallel-tests --enable-shared=no CPPFLAGS="-I${HDF5_DIR}/include" LDFLAGS="-L${HDF5_DIR}/lib" make -j 8 make check make install

第三步：安装NetCDF-Fortran这是最容易出错的环节。必须使用与NetCDF-C完全匹配的版本。对于NetCDF-C 4.9.2，推荐使用NetCDF-Fortran 4.6.1。

cd $HOME/src wget https://github.com/Unidata/netcdf-fortran/archive/refs/tags/v4.6.1.tar.gz -O netcdf-fortran-4.6.1.tar.gz tar -xzvf netcdf-fortran-4.6.1.tar.gz cd netcdf-fortran-4.6.1

配置时，必须指定NetCDF-C的安装路径。

export NETCDF=$HOME/local/coawst_libs export LD_LIBRARY_PATH=$NETCDF/lib:$LD_LIBRARY_PATH ./configure --prefix=$HOME/local/coawst_libs --enable-shared=no CPPFLAGS="-I${NETCDF}/include" LDFLAGS="-L${NETCDF}/lib" make -j 8 make check # Fortran的测试有时会因为MPI环境而失败，如果只是少数测试失败，可以尝试忽略，但最好确保大部分通过 make install

第四步：设置环境变量将库的路径永久添加到你的shell配置文件中。

# 添加到 ~/.bashrc export COAWST_LIBS=$HOME/local/coawst_libs export PATH=$COAWST_LIBS/bin:$PATH export LD_LIBRARY_PATH=$COAWST_LIBS/lib:$LD_LIBRARY_PATH export NETCDF=$COAWST_LIBS export JASPERLIB=/usr/lib/x86_64-linux-gnu # 如果安装了JASPER export JASPERINC=/usr/include/jasper # 如果安装了JASPER

执行source ~/.bashrc。然后进行验证：

nc-config --all # 查看NetCDF-C配置 nf-config --all # 查看NetCDF-Fortran配置 which ncdump which ncgen

如果这些命令都能正确显示路径和版本信息，并且nf-config显示的NC_CONFIG指向你安装的NetCDF-C路径，那么依赖库的安装就基本成功了。

实操心得：编译NetCDF-Fortran时，如果遇到“Error: Size of ‘ptrdiff_t’ (8) is larger than ‘integer(kind=4)’”这类错误，通常是因为编译器或系统架构的差异。一个有效的解决方法是，在configure命令中添加FCFLAGS="-fallow-argument-mismatch"（对于gfortran 10+版本）或尝试使用稍旧一点的NetCDF-Fortran版本（如4.5.4）。另一个常见问题是make check失败，尤其是并行测试。如果只是少数几个并行I/O测试失败，而串行测试全部通过，且你的模型不需要NetCDF的并行I/O特性，可以暂时忽略，继续安装。但为了系统完整性，最好排查MPI环境或尝试关闭--enable-parallel-tests。

4. COAWST源码获取与构建脚本配置

当依赖库就绪后，我们就可以开始处理COAWST本体了。

4.1 获取COAWST源码

COAWST的代码托管在GitHub上。推荐使用git克隆，方便后续更新。

cd $HOME git clone https://github.com/COAWST/COAWST.git cd COAWST

克隆后，目录下会有若干子目录，其中最重要的包括ROMS、WRF、SWAN、Projects（存放各种预置的测试用例）以及核心的构建脚本build_roms.csh。

4.2 理解与配置 build_roms.csh

build_roms.csh是安装的“总指挥”。我们不需要直接修改它，而是通过创建一个配置文件来定制我们的构建选项。

首先，查看有哪些可用的预置配置：

ls ROMS/Compilers/*.mk

你会看到很多以.mk结尾的文件，如Linux-gfortran.mk、Linux-pgi.mk、Linux-ifort.mk等，它们对应不同的操作系统和编译器组合。

我们以GNU编译器在Linux上为例。复制一个最接近的模板配置文件到你的工作目录：

cp ROMS/Compilers/Linux-gfortran.mk my_linux.mk

现在，用文本编辑器（如vim或nano）打开my_linux.mk。我们需要修改几个关键部分：

# 指定编译器 export MPI_ROOT = /usr/lib/openmpi # 你的MPI安装路径 export FC = mpif90 export CC = mpicc export CXX = mpic++ export FFLAGS = -frepack-arrays -O3 -ffree-line-length-none -fallow-argument-mismatch export CFLAGS = -O3 export CXXFLAGS = -O3 export CPP = /lib/cpp export CPPFLAGS = -P -traditional export CLEAN = clean # 指定库路径 - 这是最关键的部分！ export USE_NETCDF4 = on # 使用NetCDF4格式 export NETCDF_INCDIR = $(HOME)/local/coawst_libs/include export NETCDF_LIBDIR = $(HOME)/local/coawst_libs/lib # 如果你编译需要HDF5（NetCDF4依赖），确保路径正确 export HDF5_INCDIR = $(HOME)/local/coawst_libs/include export HDF5_LIBDIR = $(HOME)/local/coawst_libs/lib # 如果你需要GRIB2支持（用于WRF的特定输入），并安装了JASPER库 # export USE_JASPER = on # export JASPER_INCDIR = /usr/include/jasper # export JASPER_LIBDIR = /usr/lib/x86_64-linux-gnu # 并行设置 export USE_MPI = on export USE_MPIF90 = on export which_MPI = openmpi # 根据你的MPI类型设置

注意：-fallow-argument-mismatch这个选项对于较新版本的gfortran（gcc 10+）几乎是必须的，因为新编译器对子程序参数类型检查更严格，而许多地球科学模型（包括WRF、ROMS）的历史代码中存在一些不严格的调用。加上这个标志可以允许编译通过。

接下来，我们需要创建一个build_roms.sh（或直接使用csh脚本）来调用构建系统。但更常用的方式是直接运行build_roms.csh并传递参数。不过，我们先创建一个简单的配置文件来指定我们的项目和应用设置。

在COAWST根目录下，可以创建一个my_build.csh文件：

#!/bin/tcsh # 设置环境变量 setenv MY_CPP_FLAGS "-D" setenv MY_PROJECT_DIR `pwd`/Projects/Upwelling # 我们使用一个测试案例，如上涌流案例 setenv MY_ROMS_SRC `pwd`/ROMS setenv COMPILERS_DIR `pwd`/ROMS/Compilers setenv MY_HEADER_DIR ${MY_PROJECT_DIR} setenv MY_ANALYTICAL_DIR ${MY_PROJECT_DIR} setenv BINDIR `pwd`/Bin setenv which_MPI openmpi setenv FORT gfortran setenv USE_NETCDF4 on setenv USE_MPI on setenv USE_MPIF90 on # 调用主构建脚本 ./build_roms.csh -j 8 -noclean -compiler my_linux.mk -app UPWELLING -d ${MY_PROJECT_DIR}

这个脚本做了几件事：

1. 设置了一些路径和环境变量。
2. 指定使用我们修改过的my_linux.mk编译器配置文件。
3. 指定构建UPWELLING这个应用（一个经典的海洋上涌流测试案例）。
4. -j 8表示用8个进程并行编译。
5. -noclean表示在编译前不清理旧的构建文件（第一次编译不需要，后续重新编译时可以加快速度）。
6. -d指定项目目录。

在运行之前，给脚本执行权限：chmod +x my_build.csh。

4.3 首次编译与问题排查

激动人心的时刻到了，运行编译脚本：

tcsh my_build.csh

编译过程会持续一段时间（可能10-30分钟，取决于机器性能），因为它会依次下载（如果尚未存在）WRF、SWAN的特定版本源码，打上COAWST的耦合补丁，然后进行编译。

在这个过程中，你极有可能遇到错误。不要慌张，这是学习COAWST安装的“必修课”。编译错误主要来自以下几个方面：

1. 依赖库路径错误：这是最常见的。症状是编译初期就报错，提示找不到netcdf.mod、nf-config命令不存在或mpif.h找不到。

   • 排查：仔细检查my_linux.mk文件中的NETCDF_INCDIR、NETCDF_LIBDIR、MPI_ROOT等路径是否正确。可以用ls命令确认目录和文件是否存在。
   • 解决：确保环境变量NETCDF已设置并指向包含lib和include的父目录。运行echo $NETCDF和nf-config --all来验证。

2. 编译器标志不兼容：症状是在编译WRF或ROMS的特定文件时，出现“参数不匹配”、“隐式接口”等Fortran编译错误。

   • 排查：查看完整的错误信息，定位到出错的源文件和行号。
   • 解决：在my_linux.mk的FFLAGS中添加更多的兼容性标志，如-fallow-argument-mismatch、-fallow-invalid-boz、-w（抑制所有警告，不推荐）等。对于非常顽固的错误，可能需要去对应的源码文件（通常在WRF或ROMS子目录下）进行微调，但这需要谨慎。

3. 内存不足或进程被杀死：在编译WRF这类大型模型时，如果使用-j选项指定了过多的并行任务，可能会耗尽内存。

   • 解决：减少并行编译任务数，将-j 8改为-j 4或-j 2。

4. 特定组件的补丁失败：脚本会自动从官网下载WRF、SWAN等组件的特定版本源码并打补丁。如果网络问题导致下载失败，或者补丁文件与源码版本不匹配，会报错。

   • 解决：可以尝试手动下载所需版本的源码（版本号在build_roms.csh或相关makefile中有指定），放入COAWST目录下的LIBRARIES文件夹内，然后重新运行编译脚本。脚本会检测到本地已有源码，跳过下载。

当编译最终成功时，你会在COAWST/Bin目录下看到生成的可执行文件，通常命名为oceanG或coawstM（取决于耦合的组件）。同时，在Projects/Upwelling目录下，会有该案例的输入文件（如ocean_upwelling.in）。

实操心得：第一次编译，建议从一个最简单的、非耦合的ROMS案例（如UPWELLING）开始。成功后再尝试加入WRF或SWAN进行耦合编译。build_roms.csh脚本功能强大但也很复杂，它的-noclean选项在调试时非常有用。如果你只修改了my_linux.mk编译器配置，可以带上-noclean重新编译，这样会跳过耗时的代码提取和补丁步骤，只重新编译链接，速度很快。另外，编译产生的中间文件在Build_roms目录下，如果编译彻底混乱了，可以删除这个目录和Bin目录下的可执行文件，然后不带-noclean选项重新开始。

5. 运行测试案例与结果验证

编译成功只是第一步，让模型跑起来并得到合理的结果才是最终目标。我们以UPWELLING案例为例，进行一个简短的运行测试。

5.1 准备运行环境与输入文件

UPWELLING案例已经自带了所有必要的输入文件。进入项目目录：

cd $HOME/COAWST/Projects/Upwelling ls

你会看到一些.in文件（如ocean_upwelling.in，这是主输入文件）、*.nc文件（初始条件、边界条件、强迫场等NetCDF文件）以及可能的数据文件。

在运行前，需要确保你的运行时环境能找到所有动态链接库：

export LD_LIBRARY_PATH=$HOME/local/coawst_libs/lib:$LD_LIBRARY_PATH

5.2 执行模型运行

由于我们编译的是并行版本，需要使用mpirun或mpiexec来启动。假设使用4个MPI进程运行：

cd $HOME/COAWST/Projects/Upwelling mpirun -np 4 $HOME/COAWST/Bin/oceanG ocean_upwelling.in

• -np 4：指定使用4个MPI进程。
• $HOME/COAWST/Bin/oceanG：你编译出的可执行文件的绝对路径。
• ocean_upwelling.in：主输入文件名。

运行开始后，终端会输出大量的日志信息，包括读取的输入参数、网格信息、时间步进状态等。对于UPWELLING这种小案例，可能几秒到一分钟就能跑完。

5.3 结果分析与可视化

运行结束后，在当前目录下会生成一系列输出文件，通常是NetCDF格式，命名类似ocean_avg_0001.nc（时间平均输出）、ocean_his_0001.nc（历史快照输出）、ocean_dia_0001.nc（诊断输出）等。

你可以使用多种工具来查看和验证结果：

1. 使用ncdump快速查看文件结构和变量：

   ncdump -h ocean_avg_0001.nc

   这会输出文件的维度、变量和属性信息，帮助你了解输出了哪些物理量（如温度temp、盐度salt、流速u/v等）。

2. 使用Python（xarray, matplotlib, cartopy）进行可视化： 这是最灵活的方式。一个简单的Python脚本示例：

   import xarray as xr import matplotlib.pyplot as plt import numpy as np # 打开输出文件 ds = xr.open_dataset('ocean_avg_0001.nc') # 查看变量 print(ds.data_vars) # 绘制最后一个时间步的海表温度 sst = ds.temp.isel(s_rho=-1, ocean_time=-1) # 选择最表层的最后一个时间点 plt.figure(figsize=(10,6)) sst.plot.contourf(levels=20, cmap='RdBu_r') plt.title('Sea Surface Temperature at Final Time Step') plt.xlabel('Longitude') plt.ylabel('Latitude') plt.colorbar(label='Temperature (C)') plt.savefig('sst_upwelling.png', dpi=150) plt.show()

3. 使用NCL或MATLAB：如果你熟悉这些传统的地球科学数据分析工具，它们也有强大的NetCDF读取和绘图功能。

对于UPWELLING案例，你应该能看到一个经典的上涌流温度场：近岸由于风生埃克曼输运，表层海水离岸，导致下层低温海水上涌，在海岸线附近形成一个低温区。

注意事项：第一次运行时，如果模型立即崩溃或没有任何输出文件生成，请首先检查：

1. 输入文件路径：在ocean_upwelling.in中，所有指向外部NetCDF文件（如网格文件、初始文件、边界文件）的路径是否正确。默认是相对路径，如果从其他目录运行，可能需要修改为绝对路径或调整相对路径。
2. 权限问题：确保你有在当前目录的写入权限。
3. MPI环境：mpirun命令是否正确。可以尝试先用mpirun --version测试，或者用-np 1单进程运行，排除并行问题。
4. 库链接：运行ldd $HOME/COAWST/Bin/oceanG检查可执行文件依赖的动态库是否都能找到。如果有not found，需要将对应库的路径加入LD_LIBRARY_PATH。

6. 进阶配置：启用耦合与自定义案例

成功运行了单一的海洋模型后，我们就可以向真正的耦合模拟迈进了。

6.1 编译耦合版本（WRF+ROMS+SWAN）

要编译包含WRF和SWAN的完全耦合版本，我们需要修改构建配置。关键在于build_roms.csh的-app参数和项目目录的选择。

COAWST在Projects目录下提供了一些耦合测试案例，例如Coupled目录下可能有一个Flux案例（用于测试海气通量交换）。我们可以基于这些案例进行编译。

首先，创建一个新的构建脚本my_build_coupled.csh：

#!/bin/tcsh setenv MY_CPP_FLAGS "-D" setenv MY_PROJECT_DIR `pwd`/Projects/Coupled/Flux # 指向一个耦合案例 setenv MY_ROMS_SRC `pwd`/ROMS setenv COMPILERS_DIR `pwd`/ROMS/Compilers setenv MY_HEADER_DIR ${MY_PROJECT_DIR} setenv MY_ANALYTICAL_DIR ${MY_PROJECT_DIR} setenv BINDIR `pwd`/Bin setenv which_MPI openmpi setenv FORT gfortran setenv USE_NETCDF4 on setenv USE_MPI on setenv USE_MPIF90 on # 激活耦合选项！这是关键。 setenv USE_COUPLING on setenv USE_WRF on setenv USE_SWAN on # 如果需要沉积物模块 # setenv USE_SEDIMENT on ./build_roms.csh -j 8 -compiler my_linux.mk -app COAWST -d ${MY_PROJECT_DIR}

注意这里的变化：

• MY_PROJECT_DIR指向了一个耦合案例目录。
• 设置了USE_COUPLING、USE_WRF和USE_SWAN环境变量为on。
• -app参数从UPWELLING变成了COAWST，这告诉构建系统我们要编译耦合的可执行文件。

运行这个脚本：

tcsh my_build_coupled.csh

编译过程会更长，因为需要额外编译WRF和SWAN。成功后，在Bin目录下生成的可执行文件可能名字不同，比如coawstM。

6.2 配置耦合案例的输入文件

耦合案例的输入文件比单一模型复杂得多。通常包括：

• ocean.in：ROMS的输入文件。
• wrfinput_d01和wrfbdry_d01：WRF的初始和边界条件文件（NetCDF格式）。
• swan.in或swan_input：SWAN的输入文件。
• coupling.in：一个关键的文本文件，用于配置耦合参数。它定义了：
  • 各组件模型的名称和数量。
  • 耦合器的类型（通常是MCT）。
  • 耦合网格的映射关系（有时通过单独的weights文件指定）。
  • 耦合时间步长：这是最重要的参数之一。它需要是各模型时间步长的整数倍，并且通常基于物理过程（如海洋表层混合层响应时间、大气边界层调整时间）和计算稳定性来设定。设置过大会丢失高频相互作用，设置过小会急剧增加计算成本。

一个简化的coupling.in片段示例：

&coupling_fields coupled_components = 3, component_name = 'WRF', 'ROMS', 'SWAN', ... / &coupling_schedule coupling_interval = 600.0, # 耦合间隔，单位：秒 start_coupling = 0.0, end_coupling = 86400.0, /

这个例子表示每600秒（10分钟），三个模型通过耦合器交换一次数据。

6.3 运行耦合模拟

运行耦合模拟的命令与单一模型类似，但需要确保所有必要的输入文件都在正确的位置，并且路径在各自的.in文件中配置正确。

cd $HOME/COAWST/Projects/Coupled/Flux mpirun -np 12 $HOME/COAWST/Bin/coawstM coupling.in

这里使用了12个MPI进程。耦合模拟的资源分配策略更复杂：你需要决定多少进程分配给WRF，多少给ROMS，多少给SWAN。这通常在coupling.in或通过MPI的进程映射功能来指定。对于初学者，可以先使用对称分配（例如各4个进程）。

耦合模拟的运行日志会显示来自不同模型和耦合器的信息，交织在一起。你需要仔细查看输出，确认数据交换正在按预期进行，没有出现“找不到变量”或“插值错误”等耦合器报错。

7. 常见问题、性能调优与心得

即使成功安装并运行，在实际科研应用中你还会遇到各种挑战。这里分享一些积累下来的经验和常见问题的解决方法。

7.1 编译与运行常见错误速查表

7.2 性能调优建议

1. MPI进程布局：对于耦合模拟，WRF（结构化网格）和ROMS（曲线网格）对并行域分解的偏好不同。WRF喜欢规则的矩形分解，而ROMS在负载均衡上可能更复杂。需要根据你的网格大小和计算节点核心数进行试验，找到最优的进程分配比例。有时，让WRF和ROMS使用相同数量的进程并保持网格分解一致，可以简化耦合插值并提升性能。
2. I/O优化：NetCDF输出是主要的性能瓶颈之一。在输入文件中，减少历史输出(HIS)和平均输出(AVG)的频率，或者只输出关键变量和层次。考虑使用时间切片输出（将长时间积分分成多个文件输出），避免单个文件过大。
3. 编译器优化：在my_linux.mk的FFLAGS中尝试更激进的优化选项，如-O3 -march=native -ffast-math。但要注意，-ffast-math可能牺牲一些IEEE浮点标准的严格性，在极少数情况下影响结果的逐位可重复性，但对于物理模拟通常是安全的，并能带来显著性能提升。
4. 耦合频率：在coupling.in中设置的coupling_interval是性能与精度权衡的关键。对于海气相互作用研究，通常600秒（10分钟）是一个合理的起点。对于台风等高动态过程，可能需要缩短到60-120秒。通过敏感性试验来确定不影响物理结论的最低耦合频率。

7.3 个人心得与建议

从我自己的使用经历来看，COAWST的安装确实是一道坎，但跨过去之后，其强大的模拟能力会让你觉得一切努力都是值得的。这里给几点最后的建议：

• 从小案例开始：绝对不要一开始就用自己的高分辨率区域和复杂配置去编译和运行。务必从Upwelling、Flux这样的标准测试案例开始。它们不仅能验证安装，其输入文件也是你学习如何配置自己案例的最佳模板。
• 版本控制：对my_linux.mk、coupling.in、ocean.in等所有配置文件进行版本控制（如使用git）。每次修改前都做一个备份或提交。当模拟出现奇怪行为时，能快速回退到之前可工作的状态。
• 善用社区：COAWST有活跃的Google Groups邮件列表和GitHub Issues页面。遇到错误时，首先仔细阅读错误信息，尝试搜索邮件列表的存档。提问时，务必提供完整的错误信息、你的编译器版本、库版本和相关的配置文件片段。
• 理解物理，而不仅是代码：COAWST是一个强大的工具，但工具的输出质量取决于使用者的物理认知。花时间理解你所激活的每一个参数（如湍流闭合方案、底摩擦公式、波浪破碎参数化）的物理意义，比盲目调参重要得多。耦合模拟的结果分析也更具挑战性，你需要能区分哪些现象是大气驱动的，哪些是海洋反馈的，这需要扎实的多学科基础。

安装和配置COAWST的过程，本身就是对地球系统数值模拟技术栈的一次深度遍历。这个过程会迫使你去理解从编译器、并行库、数据格式到各学科物理模型的完整链条。虽然初期会遇到不少挫折，但每一次问题的解决都会让你对这套系统的掌控力更深一分。当看到第一个耦合模拟的动画结果成功展示出海气相互作用的复杂图景时，那种成就感无疑是巨大的。希望这份基于实战经验的指南，能帮你更平稳地度过入门阶段，早日将COAWST这个利器运用到你的研究中去。