python绘制剖面图

剖面图使用的资料是ERA5 hourly data on pressure levels from 1979 to present的再分析资料，降水量图则用的是ERA5 hourly data on single levels from 1979 to present，精度均为0.25° * 0.25°。无论用什么资料，在复现的过程中请大家务必务必注意精度到底是怎么样的，否则抄代码都抄歪了。

垂直剖面图

因为想一次性输出很多张垂直剖面图，于是用到了函数块来实现，在画剖面图时把通用的设置放在了CutawayRH()这个函数里，不同时次则在外面的循环中控制。其实函数迭代啥的都不是很难，不过时间和arange元素的对应上有时候绕不过来就会出大问题……

以我给图片标题和图片名称的迭代为例子吧：

lat = obj['latitude'][:]

lon = obj['longitude'][:]

level = obj['level'][:] #读取气压层数，一维的

hour = [0, 6, 12, 18]

以上代码时根据我所使用的再分析资料来写的，我的资料里是每隔6小时有一份数据，而且我将两天(即21日、22日)的资料都放在一起了，即一共有8个时次的信息，如果说我们直接用range(0,8)来带进去循环的话，也不是不可以，但是需要多判断几次，很容易被绕晕，所以我选择了先画前一天再画后一天的笨办法……

for time in range(0,4):

r = obj['r'][time][: