- 积分
- 541
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-10-19
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 weety9394 于 2016-11-29 12:05 编辑
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
time=in->time
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/%"
yyyy = cd_calendar(in->time, 4) ;处理时间变量
iStrt = ind(yyyy.eq.(1986*100+ 1))
iLast = ind(yyyy.eq.(2015*100+12))
tcc2=tcc@add_offset+tcc*tcc@scale_factor ;处理tcc变量
copy_VarMeta(tcc,tcc2)
tcc2_mean=dim_avg_n_Wrap(tcc2,(/1,2/)) ;区域平均
x =yyyy
y =tcc2_mean*100
y1 = month_to_annual(y, 0) ; prc_annual_total(30) ;年平均
y2 =y1/12
plot=gsn_csm_xy(wks,x,y2,res)
end
|
-
|