利用ERA5数据计算假相当位温并做垂直剖面

lc2016

最近用到了剖面图进行诊断分析，总结整理一下，方便以后拿来用。

用era5数据绘制剖面图比较简单，在绘图时选择plot=gsn_csm_contour(vars（：，{38:50}，：）,res)  ,设置合适的res基本就可以绘制一幅简单的剖面图。但是实际一次天气过程中气象要素场的分布不会是理想的水平或者垂直的，有时可能需要斜着进行剖面来分析，因此就需要插值来得到最终绘图的变量分布。

后来又对程序简单处理了一下，改成了给出任意两点的经纬度值，就可以绘制两点间的剖面。较为通用。但应该也只适用于era5这种等经纬度数据，像wrf那种数据可能需要多一步插值得到等经纬度数据，或者更改合适的插值函数。

利用NCL 处理ERA5数据计算得到多层假相当位温，作假相当位温垂直剖面。

ERA5多气压层，多气象要素，0.25°*0.25°，nc格式，数据下载：https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-pressure-levels?tab=form

需要多个气压层的比湿q和温度，计算得到假相当位温后，通过插值函数插值到设定的经纬度上，绘图即可得到剖面图。

若要叠加地形还需要下载该时刻对应的地面气压数据：https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-single-levels?tab=form。

可以对任意两点进行垂直剖面，直接修改参数起始点（lat0，lon0）——结束点（lat1，lon1）即可。

结果参考图：经度剖面：（lat0：35，lon0:110；lat1:35,lon1:120）

经度剖面

纬度剖面：（lat0：35，lon0:116；lat1：44:,lon1:116）

纬度剖面

任意两点剖面

任意两点

omak718

是不是上传错附件了，上传nc文件干啥呢，不应该是分享代码吗

小其其格

我之前都是用wrf_user_intrp2d或者wrf_user_intrp3d函数来做插值，感觉还是楼主这个理解起来更简单些

lc2016

小其其格 发表于 2024-2-25 19:01
我之前都是用wrf_user_intrp2d或者wrf_user_intrp3d函数来做插值，感觉还是楼主这个理解起来更简单些

你说的这两个函数可能适用于wrf那种不规则网格的数据，不过大致思路是差不多的

lc2016

omak718 发表于 2024-2-25 00:21
是不是上传错附件了，上传nc文件干啥呢，不应该是分享代码吗

有nc文件也有ncl文件

TingTing_Hsu

最近正在搞这个 解燃眉之急 感谢楼主分享 楼主万岁万岁万万岁

星火

我看脚本里Dxxx取的0.1，资料的分辨率应该是0.1*0.1是吧？

lc2016

星火 发表于 2024-7-2 15:03
我看脚本里Dxxx取的0.1，资料的分辨率应该是0.1*0.1是吧？

这是插值到剖面后的水平分辨率，是我们自己定义的。

maomao(陪你长大

感谢这种解释+代码的帖子，参考意义最大，谢谢!

Kminl

小其其格 发表于 2024-2-25 19:01
我之前都是用wrf_user_intrp2d或者wrf_user_intrp3d函数来做插值，感觉还是楼主这个理解起来更简单些

请问怎么计算wrfout的假相当位温呢，要把每个变量都先wrf_user_intrp3d插值，再按公式进行计算吗