关于据四个时刻的nc雪深资料求逐日平均的问题~

逸儿~~

有四个文件netcdf00.nc,netcdf06.nc,netcdf12.nc,netcdf18.nc,分别为1957年9月1号至2002年8月31号每日00时，06时，12时，18时的积雪深度，想根据这四个时刻的雪深求平均作为日平均积雪深度，从而得到1957年9月1号至2002年8月31号逐日平均积雪深度~

四个nc资料的维度维数是一样的，只是观测时间不同。如下：

以下是求逐日平均的gs文件，
'reinit'
'sdfopen d:\plan4\netcdf00.nc'
'sdfopen d:\plan4\netcdf06.nc'
'sdfopen d:\plan4\netcdf12.nc'
'sdfopen d:\plan4\netcdf18.nc'
'set lon 0 360'
'set lat -90 90'
'set lev 0 0'
'set t 1 16436'
'define sdsum=sd.1+sd.2+sd.3+sd.4'
'define sddy=sdsum/4'
*'以上是根据每日00、06、12、18四个时刻的数据在全球每个格点上求日平均'

'set fwrite d:\plan4\netsddy.grd'
'set gxout fwrite'
'set lon 70 105'
'set lat 25 40'
'set lev 0 0'
'set t 1 16436'
'd sddy'
*'将青藏高原的逐日平均数据存到netsddy.grd'
'disable fwrite'
'reinit'


运行错误提示如下：


而且运行完之后grads就没反应了。。。

想问下各位大神这是为什么呢，纬向、经向格点数都一样，垂直方向也是在一层上，只因为不是同一时刻的就不能求平均吗？如果能的话，该怎样做呢？先多谢各位啦~

mofangbao

逸儿~~ 发表于 2014-11-15 22:40
又写了一个程序，可是这个得到的逐日平均数据都相同。。。
'reinit'
'sdfopen d:\plan4\netcdf00.nc'

disable frwite放到循环外面去

lqouc

逐个打开然后define再关闭，最后把define的四个平均。
不知道可不可以哈

逸儿~~

又写了一个程序，可是这个得到的逐日平均数据都相同。。。
'reinit'
'sdfopen d:\plan4\netcdf00.nc'
'sdfopen d:\plan4\netcdf06.nc'
'sdfopen d:\plan4\netcdf12.nc'
'sdfopen d:\plan4\netcdf18.nc'

'set fwrite d:\plan4\sdaa.grd'
'set gxout fwrite'
'set lon 70 105'
'set lat 25 40'
'set lev 0 0'
i=1
while(i<=16436)
'set t 'i''
'define sda=sd+sd.2(time+6hr)+sd.3(time+12hr)+sd.4(time+18hr)'
'sdaa=sda/4'
'd sdaa'
'disable fwrite'
i=i+1
endwhile
;
应该是循环的问题，到底 是哪儿错了呢？急求大神帮忙~~~

逸儿~~

mofangbao 发表于 2014-11-15 23:03
disable frwite放到循环外面去

好了耶！！太好啦！！多谢~~

逸儿~~

lqouc 发表于 2014-11-15 17:26
逐个打开然后define再关闭，最后把define的四个平均。
不知道可不可以哈

嘿嘿，谢谢哦，试过这个方法只是时间维度还是不一致，不过正因此才想到啦循环的方法！！非常感谢，嘿嘿~~

测绘

楼主，能否告诉您的雪深数据从哪里下载的？多谢啦

逸儿~~

寒区旱区科学数据中心http://westdc.westgis.ac.cn/

jolincai

想问一下 define那一句 为什么不加time+06hr等括号里的语句就会错误呢？？？