grads中的fwrite将nc格式的写成dat格式的为何出错？

姚小娟

本人不会用fortran读取nc文件，只好先使用grads写成dat文件才进行数据处理，本人将olr.day.mean.nc中2013年5月份的olr逐日资料提取出来，写入dat文件中，其gs文件如下：
'reinit'
'sdfopen E:/NCEP/olr/olr.day.mean.nc'
'set gxout fwrite'
'set fwrite E:/file/climate/data/olr/olr.daily.201305.dat'
'set x 1 144'
'set y 1 73'
it=4139
while(it<=4169)
'set t 'it
'd olr'

it=it+1
endwhile

'disable fwrite'

'reinit'


问题就来了，我写出来的dat数据里面是空的，当使用grads来dispaly的时候，结果却是-9e+08 to 0 interval 1e+08.
我检查了其他，没发现有写错漏的，请问各位大虾，这是哪出现了问题呢？？？

姚小娟

unice 发表于 2013-7-17 20:24

                               
登录/注册后可看大图

楼主问题解决了没有，我也碰到这个问题了，很是头疼啊

恩，解决了，是缺测值问题，nc文件自带的缺测值有时候不一定准确，用set gxout state可查到真正的缺测值

le_soleil

'reinit'
'sdfopen G:\by\hgt.mon.mean.nc'
'set gxout fwrite'
'set fwrite G:\by\hgt5.grd'
i=5
while(i<=780)
k=12
'set x 1 144'
'set y 1 73'
'set z 'k''
'set t 'i''
'd hgt'
i=i+12
endwhile
'disable fwrite'
;
这是一个模板，我觉得你应该是把经纬度描述放到循环之外了，此外，GRADS读取数据才是正途，很方便，我看到过一个FORTRAN读取的，特别长，看的心烦，极简主义才是好的，我也是个新手，慢慢学习吧

姚小娟

补充，用grads将olr的月平均资料的nc文件写成dat文件时没问题的。但日平均资料无论怎么修改缺测值，都没正确过。

lqouc

grads没有报错？

river

明显是你提取出来以后编的ctl有问题，缺测值肯定不对

guo2008yitao

楼主是要用全球的资料吗？如果不是的话，尽量不要这样'set x 1 144'，'set y 1 73'设置经纬度，用那块就转换那块。因为个点上1和144是一个点，读写时可能会出错。我以前也遇到过这种问题，改了经纬度设置后就好了！

蔚蓝的天

顶一个哈哦

nanqiyuan

不懂。。。不知道咋的

unice

楼主问题解决了没有，我也碰到这个问题了，很是头疼啊

nnlovess

unice 发表于 2013-7-17 20:24

                               
登录/注册后可看大图

楼主问题解决了没有，我也碰到这个问题了，很是头疼啊

应该是缺测值的问题，你看看你那个值域， -9e+08 to 0 interval 1e+08. 把缺测值换成下限这个值-9.99e8