- 积分
- 541
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-10-19
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
为啥我画出来的是折线图不是平滑的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"
;************************************************
begin
;---Read in data
in=addfile("/lustre/home/niehw/ERA/11.nc","r")
tcc=in->tcc ; get tcc data
wks = gsn_open_wks ("x11","wxfb") ; send graphics to PNG file
;---Set plotting parameters
res = True ; plot mods desired
res@trXMinF =1985
res@trXMaxF =2015
res@tiXAxisString = "year"
res@tiYAxisString = "Total cloud cover/%"
tcc2=tcc@add_offset+tcc*tcc@scale_factor
copy_VarMeta(tcc,tcc2)
tcc2_mean=dim_avg_n_Wrap(tcc2,(/1,2/))
;print(tcc2_mean)
y =tcc2_mean*100
;print(y)
y1 = month_to_annual(y, 0) ; prc_annual_total(30)
y2 =y1/12
;print(y2)
asciiwrite("avg.txt ",y2)
;print(x)
x1=ispan(1985,2015,5)
plot=gsn_csm_xy(wks,x1,y2,res)
end
|
|