爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3018|回复: 2

[求助] 关于fortran的二进制输出

[复制链接]

新浪微博达人勋

发表于 2013-12-20 20:32:40 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
最近在用气候中心的环流数据来写作业。
我用fortran读取十进制的txt,再转化为二进制的dat,发现用grads读取会出现9.80909e-45这样的数。
我又输出成十进制发现没有问题。
请问,如何输出才能使grads能够正确读取?
下面是我的fortran。

Program fg
parameter(iva=74,imon=12,iyr=100,ifg=16)
integer i,j,k,irec
dimension hc068(iva,imon,iyr)
open(10,file='C:\Users\hp\Desktop\NCEP\hc068.txt',err=200)
      read(10,100) (((hc068(k,i,j),k=1,iva),i=1,imon),j=1,iyr)
100        format(37I5)
200        continue
close(10)

open(11,file='C:\Users\hp\Desktop\NCEP\fugao.txt')
write(11,"(I5)")((hc068(ifg,i,j),i=1,imon),j=1,iyr)
close(11)

open(19,file='C:\Users\hp\Desktop\NCEP\fugao.dat',access='direct',form='binary',recl=imon*iyr*4)
write(19,rec=1)((hc068(ifg,i,j),i=1,imon),j=1,iyr)
close(19)

End program fg
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-12-20 21:04:37 | 显示全部楼层
open(19,file='C:\Users\hp\Desktop\NCEP\fugao.dat',form='binary')
write(19)((hc068(ifg,i,j),i=1,imon),j=1,iyr)
close(19)
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-12-20 21:33:30 | 显示全部楼层

我已经解决问题了,是77的问题,把格式改成(,*)就可以了
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表