| 
 
	积分47561贡献 精华在线时间 小时注册时间2012-11-14最后登录1970-1-1 
 | 
 
 发表于 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站点绘图,可以看看和你写的有什么不一样。
 至于你的报错,说数据的记录长度不确定。我个人感觉是数据的问题,可能是数据里面的数有的带小数点有的不带。
 | 
 |