NCL画OLR随时间的变化趋势图

清晨迪声

请教各位大神，我用ncl画OLR随时间的变化趋势图，怎么在图上添加趋势线呢？
下面是我的脚本和画出的图
各位大神帮忙看一下呀

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/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/skewt_func.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/wind_rose.ncl"
begin
f=addfile("/yhhome/home/wangchendi/olr/olr.day.mean.nc","r")
print(f)
lat   = f->lat                                 
  lon   = f->lon                                   
  time  = f->time  
tt = time/24/365  
olr  = f->olr
olr1 = olr*0.01+327.65                  
copy_VarCoords(olr,olr1)               
printMinMax(olr1,False)                 
wtype = "png"
wtype@wkWidth = 3000
wtype@wkHeight = 3000
wks = gsn_open_wks(wtype,"beijing")
gsn_define_colormap(wks,"rainbow")
res = True
res@trYMinF= -34000
res@trYMaxF= -2000
res@trXMinF= min(tt)
res@trXMaxF= max(tt)
res@vpWidthF= 0.7
res@vpHeightF= 0.5
res@tmXBMode="Explicit"
res@tmXBValues= (/2004.0,2006.0,2008.0,2010.0,2012.0,2014.0/)
res@tmXBLabels= ispan(2004,2014,2)
res@tiMainString= "Variation of the OLR in Beijing"
res@tiXAxisString= "year"
res@tiYAxisString= "OLR (w/m^2)"
plot = gsn_csm_xy(wks,tt,olr(:,{39.9},{116.3}),res)

end

deeli

清晨迪声 发表于 2014-4-24 21:50
olr  = f->olr
olr1 = olr*0.01+327.65                  
copy_VarCoords(olr,olr1)                 ...

这个不是改数据型，是把olr的数据属性赋给olr1，读进来是数据都是有个类型的，转数据型是这种函数：short2flt为例的同类转换函数，你自己官网找找，再思考看看~~

krysunder

坐等高手解答

longlivehj

用regline做线性回归。先得到回归系数，然后做polyline趋势线。regline的用法看官网帮助，上面还有示例。

lleoiu

longlivehj 发表于 2014-4-24 12:08
用regline做线性回归。先得到回归系数，然后做polyline趋势线。regline的用法看官网帮助，上面还有示例。

我有一个疑问，在NCL里面，regline可以做趋势线，dtrend也可以做，那这两个做出来的有什么区别么？

deeli

我想说，亲，你的这个olr值貌似有点异常吧？读文件的时候要改下数据型吧？

沙颖凯

如果是日值的OLR在100 W/m^2已经很大了，而且为什么LZ的值全是负的，有点奇怪

清晨迪声

deeli 发表于 2014-4-24 20:45
我想说，亲，你的这个olr值貌似有点异常吧？读文件的时候要改下数据型吧？

好像已经改了吧

清晨迪声

沙颖凯 发表于 2014-4-24 21:35
如果是日值的OLR在100 W/m^2已经很大了，而且为什么LZ的值全是负的，有点奇怪

我也觉得挺奇怪的

清晨迪声

deeli 发表于 2014-4-24 20:45
我想说，亲，你的这个olr值貌似有点异常吧？读文件的时候要改下数据型吧？

olr  = f->olr
olr1 = olr*0.01+327.65                  
copy_VarCoords(olr,olr1)               
printMinMax(olr1,False)   
这个不是在改数据类型么？
刚刚接触这个，有点不懂唉