24小时热门版块排行榜    

查看: 1730  |  回复: 2

huazai9090

新虫 (初入文坛)

[交流] 求助Fortran如何读入矩阵值的数据文件 已有2人参与

有一个子程序需要调用一个mesh.dat文件
       0       0           4
       1       0           4
       2       0           4
       3       0           4
       4       0           4
       5       0           4
       6       0           5
       7       0          -1
       8       0          -1
       9       0          -1
      10       0           6
.....
       0       1           1
       1       1           0
       2       1           0
       3       1           0
       4       1           0
       5       1           0
       6       1           5
       7       1          -1
       8       1          -1
       9       1          -1
      10       1           6
................

第一列是矩阵的i,第二列是矩阵的j,第三列是数组R(i,j)的值
请问怎样不过多耗内存的读入这个文件的数据并都成功赋值给R(i,j)呀?数组R(i,j)总共是有R(0:128,0:128)
自己写的程序:
open(21,file='mesh.dat')
         
            do j=0,ny+1
                do i=0,nx+1
                    read(21,) ii,jj,R(i,j)
                end do
            end do

Fortran刚刚学还不是很懂,写出的read老是出错,因为重复出现i,j的值,所以有没有必要用ii,jj去保存mesh文件里的i,j值呀?这个mesh.dat文件已经复制到路径里面了,求大神教教我!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

学虫dj

金虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
read的格式不对呀

[ 发自小木虫客户端 ]
穷则独善其事,达则兼济天下
2楼2015-07-07 23:55:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weiwei43

至尊木虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
open(21,file='mesh.dat')
           
             do j=0,ny+1
                 do i=0,nx+1
                     read(21,*) ii,jj,R(ii,jj)
                 end do
             end do
close(21)
3楼2015-07-08 09:41:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 huazai9090 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见