24小时热门版块排行榜    

查看: 3066  |  回复: 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的回帖

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的回帖
查看全部 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的回帖

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的回帖

david830wu

新虫 (初入文坛)


dubo(金币+1): 欢迎常来程序语言版讨论 2011-07-31 13:21:47
嗯,我这个问题已经解决了:
Buffer = "3  2a  5  f0  0  0"  ! 就是说2a, 5, f0为三角网格的三个节点的索引,就是所要的数据,程序如下:
        !------------------------------------------------------------------------
        ! This subroutine results the indices of the three node, which
        ! forms the triangle face, givening a character string buffer.
        ! node_list is a dimension(3) integer vector which saves the indices of
        ! the three nodes.
        subroutine import_face(buffer, node_list)
                implicit none
                character(len = *), intent(in) :: buffer
                integer, dimension(3) :: node_list
                integer, dimension(4) :: hex_buffer
                integer :: i, j, m
               
                j = 3
                do i = 1, 3
                        do while(buffer(j:j) == ' ')
                                j = j + 1
                        end do
                        m = 1
                        do while(buffer(j:j) /= ' ' .and. m < 5)
                                read(buffer(j:j),"(Z1)" hex_buffer(m)
                                m = m + 1
                                j = j + 1
                        end do
                        m = m - 1
                        node_list(i) =  hex_buffer(1)*(16**(m-1))+hex_buffer(2)*(16**(m-2))+hex_buffer(3)*(16**(m-3))+hex_buffer(4)*1
                end do
        end subroutine import_face

鼓励批评改进啊!
5楼2011-07-18 13:08:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见