爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 6287|回复: 2

用ncl滤波,资料为GRD,输出也为grd格式的问题

[复制链接]
发表于 2015-4-3 08:18:19 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册

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

密码修改失败请联系微信:mofangbao
发表于 2016-11-28 23:35:08 来自手机 | 显示全部楼层
你好,我想请教一下nwgt的值是怎么确定的?
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2016-12-3 20:56:46 | 显示全部楼层
看你原始资料的长度,滤波有资料损失,你原始资料短的话滤波系数的个数就要少一点。 我取80点滤波就前后损失40个时次的数据。  在原始资料长度允许的情况下,理论上滤波系数个数越多误差越小。 所取滤波系数的最小值要看你的滤波的频率,就像我滤20天以上的 每天4个时次 也就是必须包含前后20个点 所以滤波系数个数不能小于41个点  但是一般41点误差还挺大的  所以我取81点。  这个值没有唯一的标准  自己根据需要来定。
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表