wrf.interplevel and wrf.vinterp results in very high values for vertical temperature gradients
Chưa có ai nhận issue này.
Đánh giá
- Độ khó
- 5/5
- Thời gian dự kiến
- Hơn một tuần
- Mức phù hợp với người mới
- 25/100
Hướng nghiên cứu
Bắt đầu bằng cách tái hiện các phép tính interplevel và vinterp được báo cáo trong issue, sau đó đọc tài liệu wrf.vinterp được liên kết và so sánh các đầu vào về độ cao và nhiệt độ. Hoàn thành có nghĩa là xác định liệu các gradient cao có phải là điều được mong đợi hay xác định một lỗi nội suy cụ thể cùng một hướng khắc phục có thể tái hiện.
Do mô hình lập chỉ mục viết ra từ nội dung của issue.
Mô tả
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.
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.
- Ngôn ngữ chính
- Python
- Star
- 498
- Fork
- 178
- Chỉ số merge pull request
- Không có pull request nào được merge trong 30 ngày
Hướng dẫn đóng góp
Bắt đầu từ đâu
- Đọc hết issue, rồi đọc hướng dẫn đóng góp của dự án.
- Bình luận trên issue rằng bạn sẽ nhận — tránh hai người làm cùng một việc.
- Fork repository và làm thay đổi trên một nhánh.
- Mở pull request có tham chiếu số hiệu của issue.
Issue khác của NCAR/wrf-python
-
NCAR/wrf-python#309 · 1 người được giao ·
-
Windows builds on conda-forge Đang mở
Độ khó 4/5 3-5 ngày Mức phù hợp với người mới 38/100
NCAR/wrf-python#307 · 5 bình luận ·
-
Độ khó 4/5 3-5 ngày Mức phù hợp với người mới 35/100
NCAR/wrf-python#286 · 3 bình luận ·
-
Windows builds and testing Đang mở
Độ khó 5/5 Hơn một tuần Mức phù hợp với người mới 25/100
NCAR/wrf-python#282 · 1 bình luận ·
-
PyPI releases Đang mở
NCAR/wrf-python#274 · 5 bình luận · 1 reaction · 1 người được giao ·
Tất cả issue của NCAR/wrf-python
Issue tương tự
-
bug
Độ khó 2/5 1-3 giờ Mức phù hợp với người mới 75/100
stephrobert/dsoxlab#238 ·
-
Độ khó 2/5 1-3 giờ Mức phù hợp với người mới 75/100
-
Độ khó 2/5 1-3 giờ Mức phù hợp với người mới 75/100
sublimehq/package_control#1780 ·
-
Độ khó 2/5 1-3 giờ Mức phù hợp với người mới 65/100
-
Độ khó 2/5 1-3 giờ Mức phù hợp với người mới 70/100
nwg-piotr/nwg-displays#145 ·