grads时间维度问题

WOT

这是我的脚本，我想绘制1990——2020年850hpa的水汽散度图像，但是我现在只能绘制某一年一个月的图像，数据时间是1948年1月1日到2022年11月1日
求哪位大神帮我改写一下脚本呢，求求了，救救孩子吧
'reinit'
'sdfopen E:\GrADS\testwww\uwnd.mon.mean.nc'
'sdfopen E:\GrADS\testwww\vwnd.mon.mean.nc'
'sdfopen E:\GrADS\testwww\rhum.mon.mean.nc'
'sdfopen E:\GrADS\testwww\shum.mon.mean.nc'
'sdfopen E:\GrADS\testwww\air.mon.mean.nc'
'set mpdset cnworld'
'set parea 1 9.5 0.5 8'
'set  grid off'
'set  grads off'
'set  lat   30  60'
'set  lon  70  100'
'set  t    876'
'set lev 850'
'define p=lev'
'define e=exp(17.67*(air.5-273.15)/(air.5-29.65))*(rhum.3/100)*6.112'
'define q=622*e/(p-0.378*e)'
'define q1=uwnd.1*q/9.8'
'define q2=vwnd.2*q/9.8'
'define p=lev'
'define e2=(6.112*exp((17.67*(air.5-273.16))/(air.5-29.65)))'
'define qs=(0.622*e2/(p-0.378*e2))'
'define q3=qs*rhum.3/100'
'define qdiv=hdivg(q3*uwnd.1/9.8,q3*vwnd.2/9.8)'
'set xlpos -10'
'set ylpos -10'
'set gxout  shaded'
'd qdiv'
'run  cbarn.gs'
'd mag(q1,q2)'
'run C:\OpenGrADS\Classic\scripts\axis.gs -type b  -position m -label on -interval 10 -sinterval 2 -size 0.1 -voffset -0.1 -suffix `3.'
'run C:\OpenGrADS\Classic\scripts\axis.gs -type L  -position o -start 30 -end 60 -suffix `3.'
'set gxout barb'
'set digsiz 0.1'
'set cthick 6'
'set ccolor 1'
'd uwnd.1;vwnd.2'
'printim  E:\GrADS\testwww\6.png white '
;