应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 怎样用read读取文件赋值二维数组
71/1返回列表
查看: 2961  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

kardrine

新虫 (初入文坛)

[求助] 怎样用read读取文件赋值二维数组 已有1人参与

用以下程序得到随机数，想保存起来，以便以后每次计算都能得到相同的结果，可是再次读取数据的时候 总是出现断点，求指教，下面把俩程序都贴出来。
      program two simension random seed
      implicit none
      real::a(300,300),b(300,300)
      real::rnd
      integer i,j,n
      real::x
      n=300
      do i=1,n
          do j=1,n
              a(i,j)=0.0
              b(i,j)=0.0
          end do
      end do
      CALL RANDOM_SEED()
      open(1,file='E:\试验程序\one dimension random seed\data1.dat')
      open(2,file='E:\试验程序\one dimension random seed\data2.dat')
      do i=1,n
          do j=1,n
              call random_number(rnd)
              x=int(rnd*10)
              x=x/10
              a(i,j)=x
          end do
      end do
      do i=1,n
          do j=1,n
              call random_number(rnd)
              x=int(rnd*10)
              x=x/10
              b(i,j)=x
          end do
      end do
      write(1,20)((a(i,j),j=1,n),i=1,n)
      write(2,20)((b(i,j),j=1,n),i=1,n)
20    FORMAT (1X,300D15.6)
      end

      program mass of disorder
      
      !1.寻找正确读取格式
      dimension fai(300,300)
      OPEN(100,FILE='E:\试验程序\one dimension random seed\data1.dat',
     1STATUS='OLD')
      open(101,file='E:\试验程序\one dimension random seed\fai.dat')
        DO i=1,300
          do j=1,300
10                READ(100,'(D15.6)')fai(i,j)
          end do
      END DO
      write(101,20)((fai(i,j),j=1,300),i=1,300)
20    FORMAT (1X,300D15.6)
      end
总是在第二个程序标着10的read行出现断点，求大神解答，初来咋到，并没有钱...
回复此楼

» 猜你喜欢
1楼 2016-11-09 17:20:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

7152075

铁杆木虫 (著名写手)

解决问题的方法比遇到的问题多

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
kardrine: 金币+5, ★★★★★最佳答案, 很感谢你的帮助，刚看到，马上好好看看，我好像有94金币，但是现在最多只能给你5个了，实在抱歉，周围遇不到会FORTRAN的人，以后可能还会麻烦您，多谢多谢 2016-11-09 19:43:21
你这个代码是有问题的，很多东西你不是太了解，相关说明我在注释中给你了，你看一下，注意一些基本的问题。加油练，慢慢写你会越来越好的，祝成功。
CODE:
  
program two_simension_random_seed !//主程序名中间不要有空格
    implicit none
    real::a( 300,300 ), b( 300,300 )
    real::rnd
    real::fai( 300,300 )
    integer i,j,n
    real::x
    n=300
    do i=1,n
        do j=1,n
            a( i,j ) =0.0
            b( i,j ) =0.0
        end do
    end do
    CALL RANDOM_SEED()
    open( 1,file='E:\试验程序\one dimension random seed\data1.dat' )
    open( 2,file='E:\试验程序\one dimension random seed\data2.dat' )!//意义不明，不知道你要干什么
    open( 101,file='E:\试验程序\one dimension random seed\fai.dat' )
    do i=1,n
        do j=1,n
            call random_number( rnd )
            x=int( rnd*10 )
            x=x/10
            a( i,j )=x
        end do
    end do
    do i=1,n
        do j=1,n
            call random_number( rnd )
            x=int( rnd*10 )
            x=x/10
            b( i,j )=x
        end do
    end do
    write( 1,20 )( ( a( i,j ),j=1,n ),i=1,n )
    write( 2,20 )( ( b( i,j ),j=1,n),i=1,n )
20  FORMAT ( 1X,300D15.6 )
    rewind( 1 ) !//将文件转回开头，从头读就行了，不必再重新open一次
    rewind( 2 ) !//另外，通道号为10以内的最好不要用，因为它们是程序预留的通道号

    DO i=1,300
        READ(1,'( 300D15.6 )' ) fai( :,i ) !//这个地方就实现了你的要求，冒号代表将从1~100，它是根据你声明部分的多少而定的
        write( 101,20 ) fai( :,i )
    END DO

    close( 1 ) !//一定要有程序结束后将文件关闭的习惯，否则以后编大程序的时候，会出各种各样的问题
    close( 2 )
    close( 101 )
    end program two_simension_random_seed
赞一下(2人) 回复此楼
Chemistry
2楼2016-11-09 19:36:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

7152075

铁杆木虫 (著名写手)

解决问题的方法比遇到的问题多
★
jjdg: 金币+1, 辛苦了 2016-11-11 11:45:05
怪了，这插入代码功能也不知道出什么问题了，重新贴。
    program two_simension_random_seed !//主程序名中间不要有空格
    implicit none
    real::a( 300,300 ), b( 300,300 )
    real::rnd
    real::fai( 300,300 )
    integer i,j,n
    real::x
    n=300
    do i=1,n
        do j=1,n
            a( i,j ) =0.0
            b( i,j ) =0.0
        end do
    end do
    CALL RANDOM_SEED()
    open( 1,file='E:\试验程序\one dimension random seed\data1.dat' )
    open( 2,file='E:\试验程序\one dimension random seed\data2.dat' )!//意义不明，不知道你要干什么
    open( 101,file='E:\试验程序\one dimension random seed\fai.dat' )
    do i=1,n
        do j=1,n
            call random_number( rnd )
            x=int( rnd*10 )
            x=x/10
            a( i,j )=x
        end do
    end do
    do i=1,n
        do j=1,n
            call random_number( rnd )
            x=int( rnd*10 )
            x=x/10
            b( i,j )=x
        end do
    end do
    write( 1,20 )( ( a( i,j ),j=1,n ),i=1,n )
    write( 2,20 )( ( b( i,j ),j=1,n),i=1,n )
20  FORMAT ( 1X,300D15.6 )
    rewind( 1 ) !//将文件转回开头，从头读就行了，不必再重新open一次
    rewind( 2 ) !//另外，通道号为10以内的最好不要用，因为它们是程序预留的通道号

    DO i=1,300
        READ(1,'( 300D15.6 )' ) fai( :,i ) !//这个地方就实现了你的要求，冒号代表将从1~100，它是根据你声明部分的多少而定的
        write( 101,20 ) fai( :,i )
    END DO

    close( 1 ) !//一定要有程序结束后将文件关闭的习惯，否则以后编大程序的时候，会出各种各样的问题
    close( 2 )
    close( 101 )
    end program two_simension_random_seed
赞一下(1人) 回复此楼
Chemistry
3楼2016-11-09 19:40:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

7152075

铁杆木虫 (著名写手)

解决问题的方法比遇到的问题多
内容已删除
回复此楼

» 本帖已获得的红花（最新10朵）

kardrine
Chemistry
4楼2016-11-09 20:29:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kardrine

新虫 (初入文坛)
送红花一朵
引用回帖:
4楼: Originally posted by 7152075 at 2016-11-09 20:29:00
小事，小事，还是希望你多写多练，同时希望你能加入一个名叫Fortran Coder的群，百度上直接打Fortran coder，那个网站是专门玩Fortran的，资料，教学视频很多，建议看看。同时希望你能加入那个群，群中大神如云，而 ...

前辈，你说的那网站一注册了，很有帮助。另外，还得让您帮忙看看，按照你给出的代码，我改了一下，
      program mass_of_disorder
      integer::i,j,k
      dimension::fai(300,300)
      
      open(1,file='E:\试验程序\one dimension random seed\data1.dat')
      open(101,file='E:\试验程序\one dimension random seed\fai.dat')
      DO i=1,300
      READ(1,'(300D15.6)') fai(:,i)
      end do
      do i=1,300
          do j=1,300
                  fai(i,j)=mod(fai(i,j),1.0)
          end do
      end do
      
      write(101,20) ((fai(i,j),j=1,300),i=1,300)
20    FORMAT (1X,300D15.6)
      !end do
      close(1)
      close(101)
      end program mass_of_disorder
之所以用mod，是因为1里面的都是[0,1)之间的小数，但是read（）fai（:,i)以后，有很多值都超过1，甚至10了，这一点弄不明白。不过反正是我需要的是随机值，也没太影响。以上这个程序可以，我把这个过程用到其他程序里面，比如以下
      DO i=1,300         
          READ(1,'(300D15.6)') fai(:,i)
          READ(1,'(300D15.6)') fai0(:,i)
          READ(1,'(300D15.6)') cata(:,i)
      end do

      do i=1,300
          do j=1,300
                  fai(i,j)=mod(fai(i,j),1.0)
                  fai0(i,j)=mod(fai(i,j),1.0)
                  cata(i,j)=mod(fai(i,j),1.0)
          end do
      end do
就总是在read行出现断点，跟最开始一样了，怎么破？
赞一下 回复此楼
5楼2016-11-10 17:28:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kardrine

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by 7152075 at 2016-11-09 19:36:26
你这个代码是有问题的，很多东西你不是太了解，相关说明我在注释中给你了，你看一下，注意一些基本的问题。加油练，慢慢写你会越来越好的，祝成功。
    program two_simension_random_seed !//主程序名中间不要有 ...

前辈，这一行，READ(1,'( 300D15.6 )' ) fai( :,i ) !//这个地方就实现了你的要求，冒号代表将从1~100，它是根据你声明部分的多少而定的，冒号代表的是1~300吧，我定义的     fai（i,j)是300*300的二维数组，另外有时间请看看我刚发的另一个求助帖吧，劳驾前辈。
赞一下 回复此楼
6楼2016-11-10 21:45:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

7152075

铁杆木虫 (著名写手)

解决问题的方法比遇到的问题多
★
jjdg: 金币+1, 谢谢 2016-11-11 11:45:31
引用回帖:
6楼: Originally posted by kardrine at 2016-11-10 21:45:50
前辈，这一行，READ(1,'( 300D15.6 )' ) fai( :,i ) !//这个地方就实现了你的要求，冒号代表将从1~100，它是根据你声明部分的多少而定的，冒号代表的是1~300吧，我定义的     fai（i,j)是300*300的二维数组，另外有 ...

不是你说的那样的，里面的随机数都是0~1之间的，你可能把其中的0.60000000D+00当成大于1的数，这是错误的，这个代表的是0.60000000x10^0(0.6乘以10的0次方)，其实是0.600000，小于1的数。
      DO i=1,300         
          READ(1,'(300D15.6)') fai(:,i)
          READ(1,'(300D15.6)') fai0(:,i)
          READ(1,'(300D15.6)') cata(:,i)
      end do

      do i=1,300
          do j=1,300
                  fai(i,j)=mod(fai(i,j),1.0)
                  fai0(i,j)=mod(fai(i,j),1.0)
                  cata(i,j)=mod(fai(i,j),1.0)
          end do
      end do
这一部分，首先第一个部分，你这么写的意思是将文件中前三行分别读完，然后赋值给了fai，fai0和cata，如果程序要是正确的话，总行数需要900行才行，但是你的文件总行数才300行，所以错了。文件读完了一行后，会直接转到下一行，并不会像你想象的那样还从头读。后面的部分不太理解你要干啥，不还是其本身嘛。前面那部分这样改：
    rewind(1) ！//别忘了让文件转到开头部分
    DO i=1,300         
        READ(1,'(300D15.6)') fai(:,i)
    end do
    rewind(1)
    DO i=1,300         
        READ(1,'(300D15.6)') fai0(:,i)
    end do
    rewind(1)
    DO i=1,300         
        READ(1,'(300D15.6)') cata(:,i)
    end do
PS：请养成一个良好的习惯，程序声明部分用implicit none。
赞一下 回复此楼
Chemistry
7楼2016-11-11 08:38:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kardrine 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 小木虫没落了，除了祈祷帖子，几乎看不到有价值的帖子 +14 botar 2026-07-14 15/750 2026-07-19 00:15 by beefly
[基金申请] 面上的水很深 +13 布布和一二 2026-07-18 18/900 2026-07-18 20:37 by 日月天地
[基金申请] 国社科通讯评审 +4 潜水冠军123 2026-07-18 5/250 2026-07-18 19:58 by chenmodeyu
[基金申请] 5个%2F，非固定段22个字母 +8 jaasxulu 2026-07-15 9/450 2026-07-18 17:40 by dxcgbm
[基金申请] 生命口C13面上会评时间 +7 luckinging 2026-07-14 11/550 2026-07-17 20:43 by 虫卜1988
[教师之家] 内心匮乏 +16 水冰月月野兔 2026-07-12 16/800 2026-07-17 19:25 by otani
[基金申请] 生命口面上会评结束了吗大家有知道结果吗？现在还没有结果正常不？ +6 碧水00 2026-07-16 9/450 2026-07-17 15:21 by jaasxulu
[论文投稿] 投英文期刊如何查重 10+5 不知道取什么名� 2026-07-15 6/300 2026-07-17 09:07 by bobvan
[基金申请] 不要再数国自然申请书的 filecode 的分隔符个数了 +14 sadpencil 2026-07-12 21/1050 2026-07-17 00:15 by wangze12014
[有机交流] 伯胺上甲基求助 100+4 002tmumu 2026-07-14 7/350 2026-07-16 19:47 by 002tmumu
[论文投稿] 跨出版社商投稿 20+5 倪好520 2026-07-13 5/250 2026-07-16 17:24 by 无视军团
[基金申请] 生生命学部会评 +4 wuke100666 2026-07-15 4/200 2026-07-16 00:12 by 啦啦小妮子
[基金申请] H口，面上，时间戳7月7日。 +3 plumage5 2026-07-15 3/150 2026-07-15 21:06 by 6543yes
[基金申请] 会评 +5 无名者登山 2026-07-15 5/250 2026-07-15 18:15 by 杠就是你对
[基金申请] E口会评 +11 布布和一二 2026-07-13 13/650 2026-07-15 16:34 by Kittylucky
[基金申请] E口青C会评结束了吗？ 40+3 小小书虫c 2026-07-14 9/450 2026-07-15 16:04 by 安静的知了
[论文投稿] 农业机械学报投稿周期 +3 G0DLIN 2026-07-14 7/350 2026-07-15 09:26 by xs74101122
[基金申请] 国自然面上D口祈祷 +4 monkeynana 2026-07-14 4/200 2026-07-15 09:14 by 白水煮
[基金申请] 没收到消息，再次陪跑 +3 天空星辰fly 2026-07-14 3/150 2026-07-14 19:14 by 6543yes
[基金申请] 明天E口面上会评 +8 布布和一二 2026-07-12 9/450 2026-07-13 22:59 by QTSorange
信息提示
关闭
请填处理意见
关闭