用ERA5和ERA-interim的UV风场计算涡度，遇到了几个问题，请求解答，感谢！

在找55

大佬们好！我在用ERA5和ERA-interim的UV风场数据分别计算涡度，并绘制涡度图，遇到了几个问题，自己琢磨了好几天都没想明白，于是发帖求帮助，请求解答！非常感谢！问题：
（1）ERA5的数据是0.25*0.25的，但我需要用1.0*1.0的，数据该怎么截取呢？u = short2flt(a->u(iStrt:iLast,:,{20:50:4},{60:120:4}))，像这样每隔4个取经纬度对吗？我画图范围是20~50°N，60~120°E，但画出来的图都很奇怪，计算涡度是用的NCL函数uv2vrG_Wrap
（2）ERA-interim的UV风场数据，最初不会Python的批量下载，就手动下载了79-18年600hPa的UV数据，用MeteoInfoMap合并，数据类型是Double，画出来的涡度图挺正常的。但我还需要画500hPa和700hPa的，就学着批量下载，下载的数据和手动下载的只有数据类型有差别，是float型的，
画出的图有点奇怪。
搞了很久都没有想明白哪儿没做对，有点晕，请求好心人解答。附上图：
（1）手动下载ERA-interim数据画的图，气候平均600hPa1月涡度图，看起来挺正常
（2）Python批量下载的interim，气候平均700hPa1月涡度
（3）ERA5的气候平均600hPa1月涡度图

guoguohh

0.25*0.25这是分辨率，你可以直接用cdo给数据提前处理成1.0*1.0分辨率的，论坛里很多相关的东西你可以搜一下

在找55

guoguohh 发表于 2022-3-10 21:50
0.25*0.25这是分辨率，你可以直接用cdo给数据提前处理成1.0*1.0分辨率的，论坛里很多相关的东西你可以搜一 ...

好的！我试一下，谢谢你~

jack_ln

ncl有好几个计算涡度的函数，你可以试试uv2vr_cfd

sc0rp1o1

你的数据是全球的，并且是固定格点，或许你应该用 uv2vrf 或者 uv2vrF

在找55

jack_ln 发表于 2022-3-11 08:01
ncl有好几个计算涡度的函数，你可以试试uv2vr_cfd

好的！我试试，谢谢~

在找55

sc0rp1o1 发表于 2022-3-11 08:47
你的数据是全球的，并且是固定格点，或许你应该用 uv2vrf 或者 uv2vrF

好的，谢谢，我再试一下

massage

小伙伴能提供脚本和青藏高原的shp文件吗，一直想要

                               
登录/注册后可看大图

来自苹果APP客户端

在找55

massage 发表于 2022-3-12 09:57
小伙伴能提供脚本和青藏高原的shp文件吗，一直想要

嗯嗯，可以哦~
我在ncl板块重新发了一个帖子，脚本和shp文件都在里面

massage

在找55 发表于 2022-3-12 11:08
嗯嗯，可以哦~
我在ncl板块重新发了一个帖子，脚本和shp文件都在里面

谢谢楼主~~~~~~~~~~~~~~~~