- 积分
- 47560
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-11-14
- 最后登录
- 1970-1-1
![[小淇的小呆] 粉丝数:75 微博数:55 新浪微博达人勋](source/plugin/sina_login/img/light.png)
|
发表于 2013-4-24 21:14:05
|
显示全部楼层
本帖最后由 lqouc 于 2013-4-24 21:15 编辑
Bubble_liu 发表于 2013-4-24 20:59
我检查了一下,如果把open里面的recl=1改为recl=4,将会得到unformatted和binary一致的结论,但是如果将r ...
乘不乘以4那是编译器的问题,一般fortran不需要乘以4,G77则需要(引自彭国伦的fortran教程)
至于binary,我用的时候就没加rec和recl,运行很好。
open(13,file='f:\data\new\last\visnew.grd',status='replace',form='binary')
do i=1,395
o=0.0
p=1
q=1
do j=1,m
write(stat,'(i8)')data3(1,j,i)
write(13)stat,data3(2,j,i),data3(3,j,i),o,p,q,data3(4,j,i),data3(5,j,i),data3(6,j,i),data3(7,j,i),data3(8,j,i)
end do
p=0
write(13)stat,data3(2,j,i),data3(3,j,i),o,p,q
end do
close(13)
这一段改编自清风的grads站点绘图,可以看看和你写的有什么不一样。
至于你的报错,说数据的记录长度不确定。我个人感觉是数据的问题,可能是数据里面的数有的带小数点有的不带。 |
|