脚本程序的错误

花楹

我有一个ncl脚本，但是循环到第四年就没了（1948、1949、1950、1951年，后面的没有了），我的脚本大概意思就是：从grd文件中读出格点数据（判断闰年与否，是闰年366天，不是闰年365天），然后插值到站点。意思很简单，但就是运行一个闰年就卡住了，研究许久也不知道哪错误，高手帮忙看看啊~多谢，有点长但是else后面和前面是一样的 ，所以不是很复杂的，大家帮帮忙啊，附件中给出ncl文件，和下面粘贴出来的一样的！ gausslat.asc (790 Bytes, 下载次数: 9)

2012-6-13 16:29 上传

点击文件名下载附件


location.asc (11.62 KB, 下载次数: 1)

2012-6-13 16:28 上传

点击文件名下载附件


interpolate_gauss.ncl (1.11 KB, 下载次数: 1)

2012-6-13 16:28 上传

点击文件名下载附件

location.asc (11.62 KB, 下载次数: 1)

2012-6-13 16:29 上传

点击文件名下载附件

load "/cygdrive/e/cygwin/usr/local/lib/ncarg/nclscripts/csm/contributed.ncl"
load "/cygdrive/e/cygwin/usr/local/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "/cygdrive/e/cygwin/usr/local/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
begin
do i=2008,2010
fili="/cygdrive/f/yin/gauss/ave/"+i+".grd"
if(mod(i,4).eq.0) then
z2=fbindirread(fili,0,(/366,94,192/),"float")
f=new((/366,753/),float)
z3=asciiread("/cygdrive/f/yin/location.asc",(/753,2/),"float")
lat1=asciiread("/cygdrive/d/essay2/gausslat.asc",-1,"float")
lon=z3(:,1)
lat=z3(:,0)
lon1=fspan(0.,360.,192)
f=linint2_points(lon1,lat1,z2,True,lon,lat,0)
;print(f)
asciiwrite("/cygdrive/f/yin/asc/ave/"+i+".txt",sprintf("%7.3f",f))
else
delete(z2)
delete(f)
delete(z3)
f=new((/365,753/),float)
z2=fbindirread(fili,0,(/365,94,192/),"float")
z3=asciiread("/cygdrive/f/yin/location.asc",(/753,2/),"float")
lat1=asciiread("/cygdrive/d/essay2/gausslat.asc",-1,"float")
lon=z3(:,1)
lat=z3(:,0)
lon1=fspan(0.,360.,192)
f=linint2_points(lon1,lat1,z2,True,lon,lat,0)
;print(f)
asciiwrite("/cygdrive/f/yin/asc/ave/"+i+".txt",sprintf("%7.3f",f))
end if
end do
end

mofangbao

这个不懂，只能帮顶了

hzwjy

首先说闰年，闰年的判断规则是：
能被4整除，但不能被100整除；
如能被100整除，需要能同时被400整除

可以使用ncl的函数 isleapyear 来判断是否闰年。
至于脚本问题，你在判断是闰年的时候是没忘记delete，但是非闰年的时候忘了

1.       delete(z2)
   
2.       delete(f)
   
3.       delete(z3)

复制代码

所以再次遇到闰年的时候就会报错。

花楹

hzwjy 发表于 2012-6-13 18:20

                               
登录/注册后可看大图

首先说闰年，闰年的判断规则是：
能被4整除，但不能被100整除；
如能被100整除，需要能同时被400整除

非常感谢，是自己粗心了~

花楹

hzwjy 发表于 2012-6-13 18:20

                               
登录/注册后可看大图

首先说闰年，闰年的判断规则是：
能被4整除，但不能被100整除；
如能被100整除，需要能同时被400整除

您好 我现在用这个脚本已经插值出来数据了，但是这个数据很明显就不对，比如说某个站1月1日和1月2日的数值相差非常大，不知道是这个脚本中哪错了呢？怎么都没检查出来 不知道能不能再帮忙检查下呢？

hzwjy

花楹 发表于 2012-6-14 14:20

                               
登录/注册后可看大图

您好 我现在用这个脚本已经插值出来数据了，但是这个数据很明显就不对，比如说某个站1月1日和1月2日的数值 ...

这个就没办法了，你可以检查下原始数据是否就差距很大。

花楹

hzwjy 发表于 2012-6-14 14:38

                               
登录/注册后可看大图

这个就没办法了，你可以检查下原始数据是否就差距很大。

这个检查不出来啊 原来的是格点数据，插值好的是站点数据