可以在ncl脚本里加上system的命令,也可以在正常画图所有图片输出后再运行convert命令
do i=0,ntimes-1
res@tiMainString = times(i)
plot = gsn_csm_vector_scalar_map(wks,uu(i,0,:,:),vv(i,0,:,:),pm(i,0,:,:),res)
end do
delete(wks)
system("convert -delay 50 wrf_pm2.5_vec*.png wrf_pm2.5_vec_hourly.gif")