请教CMIP6垂直坐标转气压坐标

dnh

大家好！最近一个星期一直搜索怎样使用python把CNRM-ESM2-1输出数据中高度混合坐标转换为气压坐标。本人已下载的数据是concdust_Emon_CNRM-ESM2-1_amip_r1i1p1f2_gr_197901-201412，使用xr.open_dataset读取有变量:lev (lev): atmosphere_hybrid_sigma_pressure_coordinate
             long_name :atmospheric model level             units :1
             positive :down
             formula :p = ap + b*ps
             formula_term :ap: ap b: b ps: ps
             bounds :lev_bounds
ps (time, lat,lon): Surface Air Pressure
ap (lev) vertical coordinate formula term:ap (k)

b (lev): vertical coordinate formula term: b(k)
b_bnds (bnds, lev)：vertical coordinate formula term: b(k+1/2)
ap_bnds (bnds, lev): vertical coordinate formula term: ap(k+1/2)

由于本人非大气科班出生且第一次使用CMIP6数据，对这些变量有疑惑，不知道是否可以使用geocat.comp.interpolation.interp_hybrid_to_pressure函数进行转换，若可以，geocat.comp.interpolation.interp_hybrid_to_pressure(data['concdust'], ps, hyam, hybm, p0=100000.0, new_levels)，ps,hyam, hybm是数据中的ps,ap和b吗？在github和谷歌上搜索过，案例所使用数据都是明确了hyam, hybm。希望有知道的大佬能进行指点，非常感谢！

Lancelot

看ncl的对应函数，geocat和ncl函数形式一致

千与千叶

ncl有函数处理这个

dnh

千与千叶 发表于 2024-11-20 15:24
ncl有函数处理这个

我查了好像是vinth2p这个函数，参数基本相同，我看他们示例中都可以直接在数据中读取hyam和hybm，而我所用的这个数据中，只有ap和b，想着想着是不是ap和b对应着hyam和hybm参数

dnh

看数据中变量lev 中p=ap+b*ps，查阅pres_hybird_ccm函数中有描述ap(k) =a(k)*p0，所以hyam = ap(k)/p0，geocat.comp.interpolation.interp_hybrid_to_pressure(data['concdust'], ps, hyam, ap/p0, p0=100000.0, new_levels)不知是否正确

dnh

Lancelot 发表于 2024-11-20 10:27
看ncl的对应函数，geocat和ncl函数形式一致

感谢回复，看pres_hybird_ccm函数中有描述ap(k) =a(k)*p0，而我下载的数据中lev中p=ap+b*ps，所以我想着hyam = ap(k)/p0，hybm=b？

Lancelot

对的，实际上函数内部计算方法就像手册中描述的一样，需要一个hyam和hybm，然后简单计算就可以转换。不同模式输出的结果不同。从变量含义来看你的想法是对的