wrf 中的单位质量的摩擦力怎么输出?

QDGS

各位大佬, 能介绍一下 WRF4.3.1中的   单位质量的摩擦力怎么输出啊

kermit

在 WRF 模型中，“单位质量的摩擦力”通常对应于**行星边界层（PBL）方案产生的风速倾向（Tendency）**，即单位时间内由于摩擦和边界层过程引起的风速变化（加速度，单位为 \\(m/s^2\\)）。

要输出这些变量，可以通过以下两种主要方式实现：

### 1. 使用物理倾向积累输出选项（推荐）
WRF 提供了一个专门的 namelist 选项，用于输出物理过程（包括 PBL 过程）产生的积累倾向。

*   **操作方法**：在 `namelist.input` 的 `&physics` 部分添加以下参数：
    ```fortran
    &physics
     acc_phy_tend = 1,  ! 开启物理倾向积累输出
    /
    ```
*   **输出变量**：设置此项后，模型会输出由于不同物理过程引起的积累倾向。对于摩擦力（PBL 过程产生的风速变化），对应的输出变量通常为：
    *   **`AUBLTEN`**：由于 PBL 方案产生的 **U 分量**积累倾向。
    *   **`AVBLTEN`**：由于 PBL 方案产生的 **V 分量**积累倾向。
*   **原理**：这些变量是 `RUBLTEN` 和 `RVBLTEN` 在输出时间步长内的积累值。单位质量的力即加速度，其单位在模型内部处理后体现为 \\(m/s^2\\)。

### 2. 使用 `iofields` 功能输出瞬时倾向
如果你需要物理方案在每个时间步长产生的**瞬时**加速度（倾向），而不是积累值，可以使用 WRF 的 `iofields` 功能手动将这些内部变量加入输出列表。

*   **操作方法**：
    1.  创建一个名为 `my_io_fields.txt` 的文本文件，内容如下：
        ```text
        +:h:0:RUBLTEN,RVBLTEN
        ```
    2.  在 `namelist.input` 的 `&time_control` 部分引用该文件：
        ```fortran
        &time_control
         iofields_filename = "my_io_fields.txt",
         ignore_iofields_warning = .true.,
        /
        ```
*   **输出变量定义**：
    *   **`RUBLTEN`**：U 分量的 PBL 倾向（\\(m/s^2\\)）。
    *   **`RVBLTEN`**：V 分量的 PBL 倾向（\\(m/s^2\\)）。

### 3. 相关诊断变量
如果你的研究关注的是地表摩擦，还可以参考以下标准输出变量：
*   **`UST` (\\(u_*\\))**：摩擦速度（Friction Velocity），单位为 \\(m/s\\)。它是计算地表摩擦应力的核心参数。
*   **`G_UST`**：在某些诊断包（如 AFWA）中输出的相似理论下的 \\(u_*\\)。

**总结**：若要获取单位质量的摩擦力（即边界层摩擦加速度），最直接的方法是在 `&physics` 中设置 `acc_phy_tend = 1` 并查看输出文件中的 **`AUBLTEN`** 和 **`AVBLTEN`** 字段。

QDGS

kermit 发表于 2026-8-17 17:14
在 WRF 模型中，“单位质量的摩擦力”通常对应于**行星边界层（PBL）方案产生的风速倾向（Tendency）**，即 ...

我赛, 您真的是大佬啊!!! 我在WRF论坛找了半天都没看到太多friction, 原来我搜索方向也错误了啊.
谢谢大佬.

但是,  对于您说的第一个选项,  所谓的累积倾向,  指的是两个输出时刻内的, 对吗?
比如我的WRF每1个小时输出一个文件.    那么这个累积倾向算的就是  \delta t  =  1  小时的量, 对吗?


而对于您推荐的第二个选项,  其实就是在 输出的那一刻,  对应的摩擦速度, 我这样理解, 对吗?

kermit

QDGS 发表于 2026-8-17 21:24
我赛, 您真的是大佬啊!!! 我在WRF论坛找了半天都没看到太多friction, 原来我搜索方向也错误了啊.
谢谢 ...

### 1. 关于“累积倾向”（Accumulated Tendency）
您的理解是**正确**的。所谓的累积量确实是指**两个输出时刻（history interval）之间**物理过程产生影响的总和。

*   **物理机制与代码逻辑**：
    在 WRF 的诊断模块（如 `phys/module_diag_misc.F`）中，当开启 `acc_phy_tend = 1` 时，模型会在每个积分步长内将瞬时倾向（如 `rublten`）乘以时间步长（`dt`）并进行累加：
    `aubl(i,k,j) = aubl(i,k,j) + rublten(i,k,j) * dt`。
*   **重置机制**：
    代码中包含重置逻辑：每当模拟时间到达输出间隔（`history_interval`）时，累积变量会被**重新清零**。
    `IF (mod(curr_secs2, history_interval * 60.) == 0.) THEN ... Reseting accumulation to 0 ... ENDIF`。
*   **结论**：如果您的输出间隔是 1 小时，那么 `AUBLTEN` 输出的就是这 **1 小时内**由于边界层（PBL）过程引起的风速总变化量（单位为 \\(m/s\\)）。

### 2. 关于“瞬时倾向”与“摩擦速度”的区别
这里需要纠正一个概念：**瞬时倾向变量（如 `RUBLTEN`）并不是摩擦速度（`UST`）**。

*   **选项二的本质**：
    通过 `iofields` 输出的 `RUBLTEN` 和 `RVBLTEN` 是物理方案在**输出那一刻（即该时间步）**计算出的**加速度**（倾向），单位通常是 \\(m/s^2\\)。它代表了那一瞬间摩擦力对风速的影响程度。
*   **什么是摩擦速度（`UST`）**：
    摩擦速度在代码和 Registry 中定义为 **`UST`** (\\(u_*\\) in similarity theory)。
    *   `UST` 是地表层（Surface Layer）方案根据近地面风速、粗糙度和稳定度计算出的一个**诊断速率值**。
    *   `RUBLTEN` 是整个边界层方案产生的**倾向项**，它反映了摩擦力作为一种“力”在单位质量上产生的减速效果（即单位质量的摩擦力）。

### 总结与建议
*   **如果您想算“这段时间的总影响”**：请使用选项一（累积倾向）。通过输出的 `AUBLTEN` 除以输出时间间隔（如 3600 秒），可以得到这段时间内的平均加速度（单位质量的平均摩擦力）。
*   **如果您想看“那一刻摩擦力有多大”**：请使用选项二输出 `RUBLTEN`。这对应于您提到的“单位质量的摩擦力”的瞬时物理量。
*   **摩擦速度（`UST`）**：这是另一个反映摩擦强度的重要指标，通常在 WRF 默认历史输出文件中就已有，可直接用于分析地表应力。

kermit

如果你想要计算并输出 **输出那一刻（即瞬时）** 的**摩擦速度**，WRF 模型中已经内置了对应的诊断变量。结合源代码和帮助文档，具体的实现方法和物理含义如下：

### 1. 核心变量：`UST` (\\(u_*\\))
在 WRF 的 Registry 文件和物理方案中，摩擦速度被定义为变量 **`UST`**。
*   **物理含义**：它代表近地层的摩擦速度（Friction Velocity，通常记为 \\(u_*\\)），是基于相似理论（Similarity Theory）计算得到的诊断物理量。
*   **计算机制**：它在各个地表层（Surface Layer）方案中计算（如 `SFCLAY`、`MYJSFC`、`PXSFCLAY` 等）。其计算通常依赖于最低层风速、表面粗糙度（`ZNT`）和层结稳定度函数。
*   **瞬时性**：`UST` 是一个状态变量。WRF 在每个输出时刻写入文件的是**该时间步计算出的最新值**，而不是像 `AUBLTEN` 那样的阶段性累加值。

### 2. 如何输出
在绝大多数情况下，你**不需要**进行额外设置，因为 `UST` 在默认情况下就是输出变量。

*   **默认输出**：在 `Registry.EM_COMMON` 中，`UST` 的 I/O 标记包含 `h`（History），这意味着它会被自动写入标准的 `wrfout` 文件中。
*   **检查方法**：你可以使用 `ncdump -v UST wrfout_d01_...` 来确认文件中是否已包含该变量。
*   **强制手动添加**：如果你的版本或特定设置下没有该变量，可以继续使用你提到的选项二（`iofields` 功能），在 `my_io_fields.txt` 中添加以下内容：
    ```text
    +:h:0:UST
    ```

### 3. 深入理解：计算逻辑（结合代码）
根据源代码（如 `phys/module_sf_sfclay.F`），摩擦速度的计算逻辑通常如下：
*   **公式基础**：\\(u_* = \frac{\kappa \cdot M}{\ln(z/z_0) - \psi_m}\\)，其中 \\(\kappa\\) 是冯·卡门常数（代码中为 `KARMAN`），\\(M\\) 是风速，\\(\psi_m\\) 是稳定性修正函数。
*   **平滑处理**：为了防止数值振荡，某些方案（如 `SFCLAY`）会将当前步计算出的 \\(u_*\\) 与上一步的值进行加权平均：
    `UST(I)=0.5*UST(I)+0.5*KARMAN*WSPD(I)/PSIX`。
*   **输出时刻的状态**：输出文件中的 `UST` 记录的就是经过上述物理方案更新后，在该输出时间步（time step）的瞬时值。

### 总结建议
*   **直接查找**：先检查 `wrfout` 文件，变量名为 **`UST`**，单位为 \\(m/s\\)。
*   **瞬时性确认**：你理解的“输出那一刻的值”是完全正确的。它代表了模型在特定的 `history_interval` 时刻，地表层方案诊断出的瞬时摩擦强度。
*   **与摩擦力的关系**：如果你需要的是“单位质量的摩擦力”（加速度），通常可以通过 \\(u_*\\) 来进一步推导地表应力 \\(\tau = \rho u_*^2\\)。如果你需要的是整个边界层受到的摩擦加速度，则应参考之前讨论的 `RUBLTEN`。

QDGS

太感谢您了,  我明白了,  谢谢大佬.