用NCL读取WRF结果并输出文件

无敌泡

求助各位大神，如何把下面的程序改成可以输出时间序列的WRF结果的SWDOWN参数？我这样画出来的是曲线图，但是我还想输出数据。
谢谢大家
1 ; Example series of plotting meteograms with WRF ARW model data
  2
  3
  4 ;***********************************************
  5 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
  6 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
  7 load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
  8 ;***********************************************
  9 begin
10 ;***********************************************
11
12 a = addfile("wrfout_d01_2014-04-30_18:00:00.nc","r")
13
14   swdown   = wrf_user_getvar(a,"SWDOWN",-1)    ; get swdown for all times
15   swdown_point  = swdown(:,18,48)             ; extract a time series at a point
16   taus = (/ 1., 2., 3., 4., 5.,6.,7.,8.,9.,10.,11.,12.,13.,14.,15.,16.,17.,18.,19.,20.,21.,22.,23.,24.,25.,26.,27.,28.,29.,30.,31.,32.,33.,34.,35.,36.,37.,38.,39.,4    0.,41.,41.,43.,44.,45.,46.,47.,48. /)       ; create a time reference
17
18   wks = gsn_open_wks("pdf","meteo1")    ; open a workstation
19
20   swdown_res = True
21   swdown_plot  = gsn_csm_xy(wks,taus,swdown_point,swdown_res)
22
23 end
24
~