逐日资料挑出指定年份风场

自由行

自己编译想挑选指定年份的夏季平均资料，但是编译一直出错显示如图所示的错误，想问下大神错在哪里，比较急
'reinit'
*从逐日nc原始数据的u、v风场资料里提取夏季风强年的数据，写入新grd文件（单层,7月）
*1000 925 850 700 600 500 400 300 250 200 150 100 70
*t：  271   319   331   391   439   475   523 (783,1mo,1948,nc)
'sdfopen D:\studing\dipin\1979-2015\step6-reason\600hpawind\strongwind\1\uwnd.1999.nc'
'sdfopen D:\studing\dipin\1979-2015\step6-reason\600hpawind\strongwind\1\uwnd.2001.nc'
*'sdfopen D:\studing\dipin\1979-2015\step6-reason\600hpawind\strongwind\1\uwnd.2002.nc'
'sdfopen D:\studing\dipin\1979-2015\step6-reason\600hpawind\strongwind\1\uwnd.2005.nc'
'sdfopen D:\studing\dipin\1979-2015\step6-reason\600hpawind\strongwind\1\uwnd.2006.nc'
'sdfopen D:\studing\dipin\1979-2015\step6-reason\600hpawind\strongwind\1\uwnd.2009.nc'
'sdfopen D:\studing\dipin\1979-2015\step6-reason\600hpawind\strongwind\1\uwnd.2012.nc'
'set gxout fwrite'
'set lat 0 50'
'set lon 60 160'
'set t 1 366'
'set lev 600'
'set fwrite D:\studing\dipin\1979-2015\step6-reason\600hpawind\strongwind\1\stronguwnd.grd'
'define have1=ave(uwnd.1,t=152,t=243)'
'd uwnd'
'define have2=ave(uwnd.2,t=152,t=243)'
'd uwnd'
'define have3=ave(uwnd.3,t=152,t=243)'
'd uwnd'
'define have4=ave(uwnd.4,t=152,t=243)'
'd uwnd'
'define have5=ave(uwnd.5,t=152,t=243)'
'd uwnd'
'define have6=ave(uwnd.6,t=152,t=243)'
'd uwnd'
'define have7=ave(uwnd.7,t=152,t=244)'
'd uwnd'
'disable fwrite'
;

lht198358

{:lxm_20:}{:lxm_20:}

咕噜

你define的是have1\have2,但是却d uwind,要改吧

咕噜

还有貌似define和d之间也要set t，

自由行

咕噜 发表于 2017-8-22 21:29
你define的是have1\have2,但是却d uwind,要改吧

好的，我试试。太感谢您的解答，也想冒昧的问问：gs里面可以同时打开多个nc数据吗

自由行

咕噜 发表于 2017-8-22 21:31
还有貌似define和d之间也要set t，

我在ave里面已经做了时间之间的计算还需要再set t了吗

自由行

咕噜 发表于 2017-8-22 21:29
你define的是have1\have2,但是却d uwind,要改吧

试了一下还是不行

自由行

咕噜 发表于 2017-8-22 21:31
还有貌似define和d之间也要set t，

'reinit'
*从逐日nc原始数据的u、v风场资料里提取夏季风强年的数据，写入新grd文件（单层,7月）
*1000 925 850 700 600 500 400 300 250 200 150 100 70
*t：  271   319   331   391   439   475   523 (783,1mo,1948,nc)
'sdfopen D:\studing\dipin\1979-2015\step6-reason\600hpawind\strongwind\1\uwnd.1999.nc'
'sdfopen D:\studing\dipin\1979-2015\step6-reason\600hpawind\strongwind\1\uwnd.2001.nc'
'sdfopen D:\studing\dipin\1979-2015\step6-reason\600hpawind\strongwind\1\uwnd.2002.nc'
'sdfopen D:\studing\dipin\1979-2015\step6-reason\600hpawind\strongwind\1\uwnd.2005.nc'
'sdfopen D:\studing\dipin\1979-2015\step6-reason\600hpawind\strongwind\1\uwnd.2006.nc'
'sdfopen D:\studing\dipin\1979-2015\step6-reason\600hpawind\strongwind\1\uwnd.2009.nc'
'sdfopen D:\studing\dipin\1979-2015\step6-reason\600hpawind\strongwind\1\uwnd.2012.nc'
'set gxout fwrite'
'set lat 0 50'
'set lon 60 160'
'set t 1 366'
'set lev 600'
'set fwrite D:\studing\dipin\1979-2015\step6-reason\600hpawind\strongwind\1\stronguwnd.grd'
'define have1=ave(uwnd.1,t=152,t=243)'
'set t 1'
'd have1'
'define have2=ave(uwnd.2,t=152,t=243)'
'set t 1'
'd have2'
'define have3=ave(uwnd.3,t=152,t=243)'
'set t 1'
'd have3'
'define have4=ave(uwnd.4,t=152,t=243)'
'set t 1'
'd have4'
'define have5=ave(uwnd.5,t=152,t=243)'
'set t 1'
'd have5'
'define have6=ave(uwnd.6,t=152,t=243)'
'set t 1'
'd have6'
'define have7=ave(uwnd.7,t=152,t=244)'
'set t 1'
'd have7'
还是说数组越界

river

自由行 发表于 2017-8-23 09:05
'reinit'
*从逐日nc原始数据的u、v风场资料里提取夏季风强年的数据，写入新grd文件（单层,7月）
*1000  ...

所有的set t 1去掉，把你的set  t  1  366 改成set t 1

咕噜

自由行 发表于 2017-8-23 08:56
好的，我试试。太感谢您的解答，也想冒昧的问问：gs里面可以同时打开多个nc数据吗

可以打开多个的，不过我也最多就打开过2个，看样子不知道是不是7个超了，你先把程序简写一下，只打开一个nc试试，运行没问题之后再逐个添加文件，比较好找问题