Hacktoberfest 2026: le issue che i maintainer hanno segnato per ottobre, aperte e adatte ai principianti. Sfoglia le issue Hacktoberfest

wrf.interplevel and wrf.vinterp results in very high values for vertical temperature gradients

Aperta
#234 0 commenti 0 reazioni 0 assegnatari Vedi su GitHub

Nessuno ha ancora preso questa issue.

Valutazione

Difficoltà
5/5
Tempo stimato
Più di una settimana
Idoneità per principianti
25/100
Tipo di issue
Bug
Chiarezza
Da chiarire
Stato di attività
Ferma
Stack tecnologico
python
Ambito
data

Direzione di ricerca

Inizia riproducendo i calcoli di interplevel e vinterp riportati nell’issue, quindi leggi la documentazione collegata di wrf.vinterp e confronta gli input di altezza e temperatura. Il lavoro è completato quando si determina se gli alti gradienti sono attesi oppure si identifica uno specifico difetto di interpolazione e un percorso di correzione riproducibile.

Scritto dal modello di indicizzazione a partire dal testo della issue.

Descrizione

Hello, I am trying to calculate the vertical temperature gradient between the surface and 1 km AGL.
I did the following:

temp2m= wrf.getvar(wrfnc,'T2') #wrfnc is the WRF output file for a specific hour
tk= wrf.getvar(wrfnc,'tk') # Temperature in Kelvin
height_agl= wrf.getvar(wrfnc,'height_agl') # Model Height for Mass Grid (AGL) in m
temp_1km = wrf.interplevel(tk, height_agl, 1000) # tempertaure interpolated to 1000 m above AGL

Then I caculate dthe temperature gradient between surace and 1 km AGL as follows:
del_T1= temp_1km-temp2m # kelvin
del_h1= 998 #1000 m- 2 m
Temp_gradient= -(del_T1/del_h1)*1000 # -dT/dz (K/km)

I am getting very high values for temperature gradient between surface and 1 km above ground level, and I am not sure if they are meteorologically correct. I have attached a screenshot of the values.
1712183790245

I also tried using wrf.vintep as follows:
temp_1km = wrf.vinterp(wrfnc, field=tk, vert_coord="ght_agl", interp_levels=[1], extrapolate=True,field_type="tk")
The document (https://wrf-python.readthedocs.io/en/main/user_api/generated/wrf.vinterp.html) says ’ght_agl’= grid point height agl [km], so I set the interp_levels as 1 km. But I similar results with slight difference.

Lingua principale
Python
Stelle
498
Fork
178
Nessuna PR unita negli ultimi 30g

Guida per i contributori

Apri la guida per i contributori

Come iniziare

  1. Leggi tutta la issue e poi la guida ai contributi del progetto.
  2. Commenta sulla issue per dire che te ne occupi tu — evita che due persone facciano lo stesso lavoro.
  3. Fai un fork del repository e lavora su un branch.
  4. Apri una pull request che faccia riferimento al numero della issue.

Altre issue di NCAR/wrf-python

Tutte le issue di NCAR/wrf-python

Issue simili

Altre issue su Python

Ricevi le nuove issue nella tua casella

Un breve riepilogo di issue GitHub adatte ai principianti.