剖面图使用的资料是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]如果说我们直接用range（0,8）来带进去循环的话，也不是不可以，但是需要多判断几次，很容易被绕晕，所以我选择了先画前一天再画后一天的笨办法……
for time in range(0,4):
r = obj['r'][time][:][:][:] #读取相对湿度，三维的
u = obj['u'][time][:][:][:] #读取纬向风，三维的
v = o
