有关用ncl绘制温度时间序列图的问题

彭雪娇

我模拟的是预测2019年9月5日-8日未来三天的气象，输出数据是每隔三小时。在用ncl处理的时候却只有5个时间点，我尝试过修改但是总是出错，请问怎样才能绘制出全部的时间序列图，脚本应该怎样修改？还希望有关前辈可以指导一下，不胜感激！另外附上我的出图和脚本。

                               

load"$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"

begin

a = addfile("wrfout_d01_2019-09-05_00:00:00","r")

t2    = wrf_user_getvar(a,"T2",-1)

t2    = t2-273.16   

t2_point  = t2(:,18,20)

taus = (/ 1., 2., 3., 4., 5. /)   

times_in_file = a->Times

dims = dimsizes(times_in_file)

times = new(dims(0),string)

do i=0,dims(0)-1

times(i) = chartostring(times_in_file(i,8:12))

end do

wks = gsn_open_wks("x11","meteo2")  

n t2_res = True

t2_res@tmXBMode       = "Explicit"

t2_res@tmXBValues     = taus  

t2_res@tmXBLabels     = times

t2_res@tmXTOn         = False  

t2_res@xyLineThicknesses = 2

t2_res@xyLineColor    =  "blue"  

t2_plot  = gsn_csm_xy(wks,taus,t2_point,t2_res)

end