Fortran Coder

查看: 101|回复: 3

[文件读写] forrt1 severe<24>读取文件错误怎么搞呀

[复制链接]

3

帖子

2

主题

0

精华

入门

F 币
34 元
贡献
13 点
发表于 2021-1-7 11:12:51 | 显示全部楼层 |阅读模式
本帖最后由 xiao. 于 2021-1-7 11:24 编辑

[Fortran] 纯文本查看 复制代码

program co2h20n2
implicit none
integer,parameter::M=2966296,N=2,L=1227852
integer::cob(60000000,2)
real::co2(L,N)
real::h2o(M,N)
real::k
integer::i,j,q,p
open(unit=10,file='T_1800_P_25.0_CO2Mol_0.1_N2mol_0.9.data')
 do i=1,1227852
   do j=1,2
    read(10,*)co2(i,j)
   end do
 end do
 close(10)
 open(unit=11,file='T_1800_P_25.0_H2OMol_0.1_N2mol_0.9.data')
 do i=1,2966296
   do j=1,2
    read(11,*)h2o(i,j)
   end do
 end do
 close(11)
 q=1
 do k=0,6000,0.000001
     do p=1,1227852
       if (abs(co2(p,1)-k)<=0.0001)then
        cob(q,2)=co2(p,2)
       end if
      end do
    q=q+1
 end do
  q=1
 do k=0,6000,0.000001
  do p=1,2966296
   if (abs(h2o(p,1)-k)<=0.0001)then
     cob(q,2)=cob(q,2)+h2o(p,2)
   end if
  end do
  q=q+1
 end do
 open(unit=12,file='hecheng.txt')
  do i=1,60000000
   do j=1,2
    write (12,*)cob(i,j)
   end do
  end do
  close(12)
  end  




谢谢各位大神帮忙!!!!!
1609989696(1).jpg
1609989679(1).jpg
1609989039(1).jpg
回复

使用道具 举报

577

帖子

1

主题

0

精华

大宗师

F 币
2873 元
贡献
1815 点
发表于 2021-1-7 11:33:05 | 显示全部楼层
FAQ之 常见错误分析
http://fcode.cn/guide-36-1.html

854

帖子

0

主题

0

精华

宗师

F 币
1454 元
贡献
807 点

规矩勋章元老勋章新人勋章水王勋章

QQ
发表于 2021-1-7 11:57:38 | 显示全部楼层
   do j=1,2
    read(10,*)co2(i,j)
   end do
改为
  read(10,*)co2(i,:)


   do j=1,2
    read(11,*)h2o(i,j)
   end do
改为
    read(11,*)h2o(i,:)

3

帖子

2

主题

0

精华

入门

F 币
34 元
贡献
13 点
 楼主| 发表于 2021-1-13 10:44:42 | 显示全部楼层
vvt 发表于 2021-1-7 11:57
do j=1,2
    read(10,*)co2(i,j)
   end do

好的谢谢!
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

QQ|捐赠本站|Archiver|关于我们 About Us|群聊|Fcode

GMT+8, 2021-1-25 02:53

Powered by Discuz! X3.2

© 2001-2017 Comsenz Inc.

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