radheat Mesosphere /Lower Thermosphere shortwave heating missing cpair factor
I maintainer di solito rispondono entro 3 giorni
Nessuno ha ancora preso questa issue.
Valutazione
- Difficoltà
- 2/5
- Tempo stimato
- 1-3 ore
- Idoneità per principianti
- 48/100
- Tipo di issue
- Bug
- Chiarezza
- Specificata chiaramente
- Stato di attività
- Attiva
- Stack tecnologico
- fortran
Direzione di ricerca
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.
Scritto dal modello di indicizzazione a partire dal testo della issue.
Descrizione
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
- Lingua principale
- Nessun dato sulla lingua
- Stelle
- 93
- Fork
- 186
- Merge medio
- 5g 15h
- PR unite (30g)
- 19
Preparare l'ambiente
Questo progetto non fornisce container di sviluppo, Dockerfile né guida per i contributori, quindi l'ambiente è a tuo carico: parti dal suo README e consulta la nostra guida al primo contributo per i passaggi generali.
Come iniziare
- Leggi tutta la issue e poi la guida ai contributi del progetto.
- Commenta sulla issue per dire che te ne occupi tu — evita che due persone facciano lo stesso lavoro.
- Fai un fork del repository e lavora su un branch.
- Apri una pull request che faccia riferimento al numero della issue.
Altre issue di ESCOMP/CAM
-
build-namelist hangs (infinite loop) on #-commented line in user_nl_camForse già presa @jtruesdal l’ha presa 24 giorni fa. Apertabug misc tag Must have CESM3
Difficoltà 2/5 1-3 ore Idoneità per principianti 84/100
ESCOMP/CAM#1701 · 2 commenti · 1 assegnatario ·
I maintainer di solito rispondono entro 3 giorni
-
CoupledEval3 Must have CESM3
Difficoltà 2/5 1-3 ore Idoneità per principianti 68/100
ESCOMP/CAM#1342 · 1 commento ·
I maintainer di solito rispondono entro 3 giorni
-
Remove for_testing_allow_non_annual_changes from user_nl_clmForse già presa Una pull request collegata a questa issue è aperta o già unita. Aperta
Difficoltà 1/5 Meno di un'ora Idoneità per principianti 25/100
I maintainer di solito rispondono entro 3 giorni
-
after CESM3_0 enhancement
Difficoltà 4/5 3-5 giorni Idoneità per principianti 45/100
I maintainer di solito rispondono entro 3 giorni
-
Bring in upcoming PUMAS bug fixes (ESCOMP/PUMAS#81); answer-changing for CAM7 defaultsForse già presa @peverwhee l’ha presa 10 giorni fa. Apertabug Must have CESM3
Difficoltà 5/5 Più di una settimana Idoneità per principianti 25/100
ESCOMP/CAM#1732 · 2 assegnatari ·
I maintainer di solito rispondono entro 3 giorni