求助，grads将nc写成grd出现问题

huisumin

这是我的gs
*原数据是1981.1开始的到2010年12月的平均值，数据里只有12个月
*提取冬季（12月-2月）的平均值
'SDFOPEN D:\fortran\guanxi\sstjuping\sst.mon.ltm.1981-2010.nc'
'set gxout fwrite'
'set fwrite D:\fortran\guanxi\sstjuping\sstpj.grd'
'set x 130 220'
'set y 132 170'
'define h=(sst(t=1)+sst(t=2)+sst(t=12))/3'
'd h'
'disable fwrite'
'reinit'
;

运行之后出现这样的错误
Data Request Warning:Request is completely outside file limits
请会的大神帮我看看，什么地方出了错

CDH

没遇到过这样的问题，但是你这两行set x 130 220和set y 132 170感觉有问题，因该是这里超出了你原始数据的范围了。

huisumin

CDH 发表于 2015-8-24 14:27
没遇到过这样的问题，但是你这两行set x 130 220和set y 132 170感觉有问题，因该是这里超出了你原始数据的 ...

全球的范围是set x 0 360 和set y 0 180，我想提取42-80N和130E-140W的，所以就写成set x 130 220和set y 132 170，应该是没有错吧

守望之月/;-D

huisumin 发表于 2015-8-25 08:49
全球的范围是set x 0 360 和set y 0 180，我想提取42-80N和130E-140W的，所以就写成set x 130 220和set y ...

这种情况直接写纬度和经度范围更不容易出错吧。
set lat 42 80
set lon 130 220?

huisumin

守望之月/;-D 发表于 2015-8-25 08:51
这种情况直接写纬度和经度范围更不容易出错吧。
set lat 42 80
set lon 130 220?

如果是这样的话，那经度set lon 130 220就不对了吧，应该是东经130到西经140，如果是220 的话，画出的图经度就没有东西经之分了

CDH

huisumin 发表于 2015-8-25 08:49
全球的范围是set x 0 360 和set y 0 180，我想提取42-80N和130E-140W的，所以就写成set x 130 220和set y ...

你可以用sdfopen命令，把你的数据在grads窗口打开，然后把set x 130 220 和 set y 132 170输入进去，看一下grads给出的经纬度范围是不是符合你的要求，就可以看出是不是这里超出数组范围了