radheat Mesosphere /Lower Thermosphere shortwave heating missing cpair factor
メンテナーはふだん 3 日以内に返信
まだ誰も着手していません。
評価
- 難易度
- 2/5
- 見積もり時間
- 1〜3時間
- 初心者へのやさしさ
- 48/100
- issue の種類
- バグ
- 明瞭さ
- 明確に書かれている
- 活発さ
- 活発
- 技術スタック
- fortran
調査の方向性
The issue points to radheat.F90 and identifies the qrs_mlt conversion before merge_qrx as the change. Read the radheat_tend implementation and check how composition-dependent cpairv is used; then run CAM with an MT model_top configuration and verify the shortwave heating has consistent J/kg/s units. The reporter says they plan to address it and need help, so coordinate with them before starting.
索引モデルが issue の本文から書いたものです。
説明
What happened?
From @briandobbins
In radheat_tend, the idealized Mesosphere /Lower Thermosphere shortwave heating qrs_mlt is built from a K/day profile ÷ 86400 (so K/s), then merged via merge_qrx with the RRTMGP qrs (which is qrs_prime, in J/kg/s). The two have inconsistent units, so wherever the merge weight qrl_wt < 1 the M/LT SW heating is ~cpair (≈1004×) too small. The LW partner qrlfomichev is already in J/kg/s, so only the SW branch is affected. This only matters for MT (model_top="mt") configs, where the top ~6 levels (p ≲ 20 Pa) run cold — full Fomichev cooling with almost no SW heating. Low-top grids have qrl_wt = 1 everywhere and are unaffected.
Fix: convert qrs_mlt to J/kg/s with the composition-dependent cpairv before merging.
What are the steps to reproduce the bug?
Run CAM
What CAM tag were you using?
latest
What machine were you running CAM on?
CISL machine (e.g. derecho)
What compiler were you using?
Intel
Path to a case directory, if applicable
No response
Will you be addressing this bug yourself?
Yes, but I will need some help
Extra info
One-line code fix:
/glade/u/home/pel/src/for-cecile-radheat/radheat.F90
- 主要言語
- 言語のデータがありません
- スター
- 93
- フォーク
- 186
- 平均マージ
- 5日 15時間
- マージ済み PR（30日）
- 19
環境構築
このプロジェクトには開発コンテナ、Dockerfile、コントリビューションガイドがありません。まず README を読み、一般的な手順ははじめてのコントリビューションガイドを参照してください。
はじめの一歩
- issue を最後まで読み、次にプロジェクトのコントリビューションガイドを読みます。
- 着手することを issue にコメントします — 二人が同じ作業をするのを防げます。
- リポジトリをフォークし、ブランチを切って変更します。
- issue 番号を参照したプルリクエストを送ります。
ESCOMP/CAM のほかの issue
-
build-namelist hangs (infinite loop) on #-commented line in user_nl_cam対応中かも @jtruesdal が 23 日前に担当しました。 オープンbug misc tag Must have CESM3
難易度 2/5 1〜3時間 初心者へのやさしさ 84/100
ESCOMP/CAM#1701 · コメント 2 件 · 担当者 1 名 ·
メンテナーはふだん 3 日以内に返信
-
CoupledEval3 Must have CESM3
難易度 2/5 1〜3時間 初心者へのやさしさ 68/100
メンテナーはふだん 3 日以内に返信
-
Remove for_testing_allow_non_annual_changes from user_nl_clm対応中かも このイシューにリンクされたプルリクエストがオープン中、またはマージ済みです。 オープン
難易度 1/5 1時間未満 初心者へのやさしさ 25/100
メンテナーはふだん 3 日以内に返信
-
after CESM3_0 enhancement
難易度 4/5 3〜5日 初心者へのやさしさ 45/100
メンテナーはふだん 3 日以内に返信
-
Bring in upcoming PUMAS bug fixes (ESCOMP/PUMAS#81); answer-changing for CAM7 defaults対応中かも @peverwhee が 9 日前に担当しました。 オープンbug Must have CESM3
難易度 5/5 1週間以上 初心者へのやさしさ 25/100
メンテナーはふだん 3 日以内に返信