一般我们使用的时间序列都是以年为单位，比如下图：
而如果要绘制月尺度的时间变化图，则需要调整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')
下面是运行结果：
❤️欢迎点赞收藏❤️