python画假相当位温叠加温度平流剖面图（附代码）

西红柿毛毛雨

雨一直下 发表于 2024-6-15 10:47
请问下载的era5数据是全要素全尺寸下载吗，还是说设定了范围

这个看自己需求~~一般选择需要的区域范围就行~~这个图的数据设定过范围的，level只用了17层，200-850hpa

西红柿毛毛雨

aa... 发表于 2024-6-3 19:36
哦其实是我运行出来的温度平流变成Array格式了，不应该是DataArray吗

是不是你对数据进行过处理，比如滤波或者去趋势后可能数据格式就变了~

雨一直下

西红柿毛毛雨 发表于 2024-6-18 17:27
这个看自己需求~~一般选择需要的区域范围就行~~这个图的数据设定过范围的，level只用了17层，200-850hpa

好的感谢

某某晨

感谢楼主分享

一物栽

您好，请问这行代码里取点是23.25N,99.5E，请问是如何转换成267,397的呢？t_thetase =t[:,:,267,397]，以及计算温度平流里tadv=np.zeros((36,17,361,401))，这四个数字分别是什么意思，可以解答一下吗，不胜感激

edwardli

西红柿毛毛雨 发表于 2024-1-8 15:37
主要是假相当位温涉及公式自己计算，代码里它的变量和公式定义的很清楚了，我也只是在论坛里看的帖子，也 ...

MetPy.calc.equivalent_potential_temperature()计算的就是假相当位温，无需自己重新实现公式计算。【雷小Py-033：物理量拾遗-假相当位温_CAPE_2024.11.21-哔哩哔哩】 https://b23.tv/dLJTymo

西红柿毛毛雨

edwardli 发表于 2024-12-20 06:23
MetPy.calc.equivalent_potential_temperature()计算的就是假相当位温，无需自己重新实现公式计算。【雷 ...

谢谢李老师分享

西红柿毛毛雨

一物栽 发表于 2024-12-19 16:49
您好，请问这行代码里取点是23.25N,99.5E，请问是如何转换成267,397的呢？t_thetase =t[:,:,267,397]，以及 ...

原始数据变量t_thetase有四个维度，分别是time，level，lat，lon，这里在切片的时候时间维和高度层这两个默认全取，纬度“267”切片出现就是23.25N，同理经度“397”切片出来就是99.5E，这和你下载的格点数据分辨率有关，大概算一下就是知道你要取的经纬度在切片的哪个位置。另外还可以直接用loc方法，直接选出来例如data = f['t_thetase'].loc[:,23.25,99.5]。tadv=np.zeros((36,17,361,401))36是时间维，17是高度层数，361是纬度格点数，401是经度格点数，这个和要画图的格点范围有关。

一物栽

西红柿毛毛雨 发表于 2024-12-20 08:28
原始数据变量t_thetase有四个维度，分别是time，level，lat，lon，这里在切片的时候时间维和高度层这两个 ...

非常感谢！！！！

edwardli

西红柿毛毛雨 发表于 2024-12-20 08:28
原始数据变量t_thetase有四个维度，分别是time，level，lat，lon，这里在切片的时候时间维和高度层这两个 ...

返回的数据结构还是xarray.dataset 或者dataarray，极其推荐使用ds.sel(longitude=112.5, latitude=24.25)，再也不用数下标，这也是为什么要在“光秃秃”的numpy.ndarray纯数组上打“标签”变成更高级的数据结构的初衷、优势。