Grads1.9和Grads2.0出来的结果不同，显然Grads2.0有问题。

天道酬勤

DSET E:\MJO\OLRA-normal\pent.olra.y%y4.grd
options yrev    365_day_calendar template  
OPTIONS template                                                                                                           
undef -9999.                                                                    
TITLE Pentad outgoing longwave radiation anomalies                              
OPTIONS template                                                                          
XDEF 144 LINEAR 0. 2.5                                                                                                                                       
YDEF 73 LINEAR -90. 2.5                                                                                                                                       
ZDEF 01 LEVELS 1000                                                                                                                                          
TDEF 1971  LINEAR 1jan1979 5dy                                                                                                                              
VARS 1                                                                          
olra  1 -999    outgoing longwave radiation (W/m*m)                              
ENDVARS                                                                        
这是读取27个数据的ctl，一看就知道，我有27个数据，并且每个数据都是候资料，每年的时间是73,那么27年就是1971了。
下面就把问题说一下
楼主想把27个数据，放到一个grd里面去。下面是我的gs。问题出来了。
'reinit'
'open E:\MJO\OLRA-normal\pent.olra.ctl'
'set fwrite E:\MJO\OLRA-normal\pent-1979-2005.grd'
'set gxout fwrite'
tt=1
while (tt<=1971)
'set x 1 144'
'set y 1 73'
'set t 'tt''
'd olra'
tt=tt+1
endwhile
'disable fwrite'
'reinit'
;
当我用grads1.9运行的时候，是没有问题的。合成的数据是正确的，因为我用第一个ctl验证过。但如果用grads2.0的话， 就会出现问题，前1247时间点是正确的，但到1248开始，数据就不正常了。

我想问，有人知道是怎么个情况么？

sunflower

我的grads2.0整合版连while语句都不认识

mofangbao

从你的截图看是缺测值的问题

river

应该是缺测值的问题吧，改成grads默认的-9.99e8 试一下

天道酬勤

mofangbao 发表于 2013-11-25 21:11

                               
登录/注册后可看大图

从你的截图看是缺测值的问题

那为什么用两个版本的，结果会不一样？

天道酬勤

river 发表于 2013-11-25 21:56

                               
登录/注册后可看大图

应该是缺测值的问题吧，改成grads默认的-9.99e8 试一下

不行。1.9能出来结果，2.0出不来，是缺测的问题么?

天道酬勤

river 发表于 2013-11-25 21:56

                               
登录/注册后可看大图

应该是缺测值的问题吧，改成grads默认的-9.99e8 试一下

grads1.9提取出的，grads2.0打开是正确的。是不是两个版本提取的时候，对缺测值的处理不一样啊。grads2.0会改缺测值么？

river

天道酬勤 发表于 2013-11-25 22:26

                               
登录/注册后可看大图

grads1.9提取出的，grads2.0打开是正确的。是不是两个版本提取的时候，对缺测值的处理不一样啊。grads2.0 ...

这个我觉得两个版本应该是一致的，因为2.0出来之前我也用1.9的，没有出现你的问题。而且那个时候1.9的提取出来的资料的缺测值也都是-9.99e8。你可以把两个版本提取的同样的资料用fortran输出到屏幕上看看缺测值是哪个，是不是一样的。

mofangbao

不要太纠结为什么这个行那个不行的  从你觉得错误的那个按正常方法查错就行了

F平方

我和你遇到了一样的问题，不过我的是2.1版本的，同样的gs脚本，2.1出图不对，1.9正常。从图上看也觉得是这两个版本对缺测值的处理不一样。但究竟是怎么不一样一直没解决