24小时热门版块排行榜    

查看: 3064  |  回复: 14
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

david830wu

新虫 (初入文坛)

[求助] Fortran的格式化输入输出问题

下面是一个矩阵的内容,在文件里,想用Fortran读入,应该怎样做呢?
难点在于,空格数不固定,数字的长度也不固定,Fortran有没有能像C++那样的自动跳过空格,自动控制长度的智能读入呢?

Matrix.dat:
1 2  3   45 2
341 3 221 2 45
21     45     3 23 3
1           32 2     34 33
45  36   39    81 4

谢谢啦!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

★ ★
jjdg(金币+2): 感谢参与 2011-07-22 01:38:55
给一个我写的 Fortran 90 版本的程序,没有加前面的准备字符串部分,呵呵,希望不难懂理解……
CODE:
program read_hex

implicit none
character(len=80) :: line="3 5  F6F  1d 0 0"
character(len=80) :: buff, hex_str, fm
integer, dimension(3) :: hex
integer :: str_len, i, first, last

buff = adjustl(line)
str_len = len_trim(buff)

!! remove the extra space in the string
do i=1,str_len-1
   if (buff(i:i) == ' ' .and. buff(i+1:i+1) == ' ') buff(i:) = buff(i+1:)
end do

!! extract the hexadecimal constant
first = 1
do i=1,3
   first = index(buff(first:), ' ') + first - 1
   last = index(buff(first+1:), ' ') + first
   hex_str = buff(first+1:last-1)
   write(fm, '(a,i0,a)') '(Z', len_trim(hex_str), ')'
   read(hex_str, fm) hex(i)
   first = last
end do

!! hexadecimal to decimal constrant, and output...
write(*,*) hex

end program read_hex

11楼2011-07-21 11:45:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

★ ★
david830wu(金币+5): 嗯,以前陷到格式化输出的沟里了,其实用“*”就好了。非常感谢! 2011-07-17 17:24:03
xzhdty(金币+2): 欢迎常来 2011-07-17 21:51:22
你每行都是一样多的数字吗?或者说,你的 matrix 的维数是已知的吗?

直接
CODE:
read(12, *) ((a(i,j), i=1,5),j=1,5)

不行吗?

当然,这里都是以浮点数读入的,你可以转一下,呵呵……
2楼2011-07-17 17:13:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

david830wu

新虫 (初入文坛)

★ ★
dubo(金币+1): 欢迎常来程序语言版讨论 2011-07-17 22:51:32
dubo(金币+1): 欢迎常来程序语言版讨论 2011-07-17 22:52:42
再进一步问一下,也是我真实遇到的,就是GAMBIT的网格划分文件,数据文件是十六进制应该怎么办呢?
MATRIX.dat
1   3f  a6    21   4
a cc        f3      2  2
3    d0    b2    46   4
……
这样的文件用“*”好像就不行了,如果后面加option:Zn[.m]的话就得指定宽度,这应该怎么处理?
3楼2011-07-17 18:45:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)


dubo(金币+1): 欢迎常来程序语言版讨论 2011-07-17 22:51:45
引用回帖:
Originally posted by david830wu at 2011-07-17 18:45:35:
再进一步问一下,也是我真实遇到的,就是GAMBIT的网格划分文件,数据文件是十六进制应该怎么办呢?
MATRIX.dat
1   3f  a6    21   4
a cc        f3      2  2
3    d0    b2    46   4
……
这样的文件用 ...

这个没有啥好办法,将多个空格替换为一个空格,然后按字符串读入,以空格划分,再以内部文件的方式格式化读入……
4楼2011-07-17 19:44:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见