- 积分
- 256
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-10-4
- 最后登录
- 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"
; ***********************************************
begin
data=fbindirread("E:/MATSA\FNL/initialdatauv/850.grd",0,(/248,2,181,360/),"float")
; wave1=onedtond( ndtooned(wave), (/2,11,11/))
u=data(:,0,:,:)
v=data(:,1,:,:)
u!0="time"
u&time=ispan(0,247,1)
u!1="lat"
u&lat=fspan(-90,90,181)
u!2="lon"
u&lon=fspan(0,360,360)
v!0="time"
v&time=ispan(0,247,1)
v!1="lat"
v&lat=fspan(-90,90,181)
v!2="lon"
v&lon=fspan(0,360,360)
ihp = 0 ; low_pass
sigma = 1.0 ; Lanczos sigma
nWgt = 81 ; loose 40 detat each end
fca = 1./10*4. ; 10days
wgtq = filwgts_lanczos (nWgt, ihp, fca, -999., sigma )
ntim = dimsizes( u )
yPlot = new ( (/248,2,181,360/) , typeof(u), getFillValue(u) )
yPlot(:,0,:,:) = wgt_runave ( u, wgtq, 0 )
yPlot(:,1,:,:) = wgt_runave ( v, wgtq, 0 )
fbinrecwrite("E:\MATSA\FNL\filterdatazjx\tcwithin\10days\850day.grd",0,yPlot)
end
这是我的程序,但是提示有警告是什么意思啊?有么有大神做过这个工作啊 求程序
warning:fbinrecwrite: end of file reached before record number, writing record a s last record in file
|
|