剖面图使用的资料是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][: