- 积分
- 82
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-7-7
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
求助关于linint_Wrap 的问题,插值画后不出图,我用 linint_Wrap 这个函数想把ERA-interim 1.5*1.5 的插值成 NCEP格式2.5*2.5,结果插值好了,运行也没提示有错误或警告,但就是画不出等值线,只有底图,而色标,单位都有的,print插值后的数也是有值的,就是画不出来,请各位大虾帮忙解决!
下边是脚本:
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"
begin
latS=0.
latN=60.
lonL=60.
lonR=160.
diri="/home/clear/ncl/mon/"
h=addfile(diri+"slp.mon.mean.nc","r")
latout=h->lat
lonout=h->lon
h1=addfile(diri+"slp7-8.nc","r")
TIME=h1->time
YYYY =ut_calendar(TIME,0)
iYYYY =ind(YYYY(:,1).eq.7)
iYYYYY =ind((YYYY(:,0).ge.1979.and.YYYY(:,0).le.2011).and.(YYYY(:,1).eq.8))
h7=short2flt(h1->msl(iYYYY,:,)
h8=short2flt(h1->msl(iYYYYY,:,)
h7ave=dim_avg(h7(latitude|:,longitude|:,time|)
h8ave=dim_avg(h8(latitude|:,longitude|:,time|)
h2=(h7ave+h8ave)/200
copy_VarCoords(h7(0,:,,h2)
latin=h2&latitude
lonin=h2&longitude
var=h2
h3=linint2_Wrap(lonin,latin,var,True,lonout,latout,0)
wks=gsn_open_wks("x11","slpera")
gsn_define_colormap(wks,"gsdtol")
resCN =True
resCN@cnFillOn =False
resCN@cnLinesOn =True
resCN@gsnAddCyclic =False
resCN@mpOutlineBoundarySets = "National"
resCN@mpGeophysicalLineThicknessF =1.
resCN@mpNationalLineThicknessF =1.
resCN@gsnSpreadColors =False
resCN@mpMinLatF = latS
resCN@mpMaxLatF = latN
resCN@mpMinLonF = lonL
resCN@mpMaxLonF = lonR
resCN@mpFillOn =False
resCN@pmTickMarkDisplayMode ="Conditional"
resCN@cnLineLabelInterval = 1
resCN@cnLevelSpacingF = 2
resCN@gsnContourLineThicknessesScale=2
resCN@cnLineLabelPlacementMode = "Constant"
resCN@cnInfoLabelOn = False
plot=gsn_csm_contour_map(wks,h3(:,,resCN)
end
|
|