用GRADS从NC中提取GRD后数据变成原来的四倍是什么原因？

南十字星奇迹

从2000年vwnd.nc文件提取平均径向风，gs如下：'reinit''sdfopen e:\data1\v.nc'
'set fwrite e:\data1\test.grd'
'set gxout fwrite'
'set lat -90 90'
'set lon 0 360'
'set lev 1000 10'
'set t 1 366'
'define v=ave(vwnd,lon=0,lon=360)'
'define v1=ave(v,t=1,t=366)'
'set t 1'
'set lon 0'
'd v1'
'disable fwrite'
;
高度共有17层，lev有73层，应该生成的grd数据是17*73个，可实际生成的是17*73*4个数据。用grds画图没有异常，但要用fortran继续编程时因多出4倍的数据没法继续下去。不知是什么情况，和计算机有关吗？我的是win7,64位机。

Aires

@mofangbao程序看着没问题啊

南十字星奇迹

Aires 发表于 2012-10-16 12:45

                               
登录/注册后可看大图

@mofangbao程序看着没问题啊

程序没问题，不知为什么出来的数据多了4倍

mofangbao

改成用while的循环做，你这样写有点乱，而且设置的时候用xy表示范围，不要用经纬度，要不然你的新文件范围容易出现偏差

南十字星奇迹

mofangbao 发表于 2012-10-16 20:20

                               
登录/注册后可看大图

改成用while的循环做，你这样写有点乱，而且设置的时候用xy表示范围，不要用经纬度，要不然你的新文件范围容 ...

问题解决了，是我fortran编程不对，太马虎了。非常感谢您的建议！