ctl无法打开

绯宜物语

读取了nc数据后，编写了如下CTL，但是有时可以打开有时打不开，打开时画图全是缺省值。这是怎么回事啊？大家帮我看看吧。谢谢。
'sdfopen f:\pr.nc'
'set gxout fwrite'
'set fwrite f:\pr.grd'
n = 1261
while(n<=1356)
'set x 1 720'
'set y 1 360'
'set t 'n
    'd pre'
n=n+1   
endwhile
'disable fwrite'
'reinit'
;

CTL:dset f:\pr.grd
title CRU
undef 9.96921e+36
xdef 720 linear -179.75 0.5
ydef 360 linear -89.75 0.5
zdef 1 linear 0 1
tdef 96 linear 00z16JAN2006 1mo
vars 1
pre=>pre 0 t,y,x precipitation
endvars

绯宜物语

我把输出的grd文件用Fortran程序看了一下，发现缺省值应该是-9.990000E+08。为什么会和nc文件里的缺省值不一样呢？

river

绯宜物语 发表于 2014-10-28 15:46
我把输出的grd文件用Fortran程序看了一下，发现缺省值应该是-9.990000E+08。为什么会和nc文件里的缺省值不 ...

这就是grads软件本身默认值

绯宜物语

river 发表于 2014-10-28 16:55
这就是grads软件本身默认值

意思是用grads提取数据时，它会改变缺省值？

river

绯宜物语 发表于 2014-10-28 18:52
意思是用grads提取数据时，它会改变缺省值？

是的，就是这个意思

绯宜物语

river 发表于 2014-10-28 19:26
是的，就是这个意思

哦哦，这样啊，学习了。谢谢你。

金色牛神

缺省值相当于一种标记flag,表示没有数据，为null,在不同的编译环境下取的数值不一样，但并不影响

绯宜物语

金色牛神 发表于 2014-10-30 23:50
缺省值相当于一种标记flag,表示没有数据，为null,在不同的编译环境下取的数值不一样，但并不影响

谢谢你。学习啦……