对CMIP6模式结果使用python插值气压坐标

dnh

在使用CMIP6数据中遇到模式垂直坐标转换气压坐标系问题，前段时间尝试将CNRM-ESM2-1数据中高度混合坐标转换为气压坐标（http://bbs.06climate.com/forum.php?mod=viewthread&tid=110716），但geocat中interp_hybrid_to_pressure函数好像并不适用将hybrid height coordinate 转换成垂直坐标系，因自身能力有限，周边认识的伙伴也都没用过CMIP6数据和垂直坐标转换问题，因此又来请教大家，不知道大家有没有什么建议，非常感谢

人可山

可以看看https://github.com/jzhe3/xcesm/blob/master/xcesm/core/xcesm.py中的hybrid_to_pressure函数。

edwardli

https://xcdat.readthedocs.io/en/ ... -to-pressure-levels

Lancelot

你给的图都告诉你答案了

dnh

edwardli 发表于 2024-11-26 09:20
https://xcdat.readthedocs.io/en/doc-485-dask-guide/examples/regridding-vertical.html#4:-Remap-cloud- ...

感谢，又学会了一种新方法，用xcdat的和geocat中函数相互验证了下，xcdat处理的更快

dnh

Lancelot 发表于 2024-11-26 12:36
你给的图都告诉你答案了

哈哈，没办法，中途出家太菜了，还得多向你们学习请假

dnh

人可山 发表于 2024-11-26 08:53
可以看看https://github.com/jzhe3/xcesm/blob/master/xcesm/core/xcesm.py中的hybrid_to_pressure函数。

感谢感谢

千与千叶

https://github.com/Yefee/xcesm

Lancelot

正确的方法是按照数据里给出的方程计算每个混合坐标对应的气压，然后对每个点用对数插值。说xcesm的怕不是都没用过

dnh

Lancelot 发表于 2024-11-26 16:09
正确的方法是按照数据里给出的方程计算每个混合坐标对应的气压，然后对每个点用对数插值。说xcesm的怕不是 ...

哦哦，这样子，以为有现成的python函数处理这个，那我沿你的这个思路先试试。感谢