一般我们使用的时间序列都是以年为单位，比如下图：

而如果要绘制月尺度的时间变化图，则需要调整X轴的标注。下面代码展示了如何绘制小时尺度的降水数据。

[sname2,lon2,lat2] = kml2xy('GZ_.kml');

nc_bound2 = [lon2,lat2];

area_ind2=inpolygon(era_tp.lon,era_tp.lat,nc_bound2(:,1),nc_bound2(:,2));

%% region area

area_scale=cal_grid_region(era_tp);

GZ_rg = zeros(size(era_tp.rg,3),1);

GZ_tt = GZ_rg;

for ii=1:size(era_tp.rg,3)

    GZ_rg(ii)=sum(sum((era_tp.rg(:,:,ii).*area_ind2.*...

        area_scale))/sum(sum(area_ind2.*area_scale)));

    GZ_tt(ii)=era_tp.tt(ii,1);

    disp(ii)

end

ttime = datenum(era_tp.tt(:,2:end));

plot(ttime,GZ_rg*1000,'-')

datetick('x','mm-dd')

xlabel('Date')

ylabel('Hourly precipitation/mm')

title('ERA5')

下面是运行结果：

❤️欢迎点赞收藏❤️