Fortran Coder

查看: 19940|回复: 10
打印 上一主题 下一主题

[文件读写] 关于数值读入读出的问题

[复制链接]

26

帖子

8

主题

0

精华

熟手

F 币
123 元
贡献
78 点
跳转到指定楼层
楼主
发表于 2015-9-23 10:59:40 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
[Fortran] 纯文本查看 复制代码
program main
implicit real*8(a-z)
integer::N,i
real*8::x1(10000),px1(10000),z1(10000),pz1(10000)
real*8::x(10000),px(10000),z(10000),pz(10000)

open(unit=11,file='x.txt')
open(unit=12,file='px.txt')
open(unit=13,file='z.txt')
open(unit=14,file='pz.txt')
open(unit=15,file='cyyb.txt')

N=10000

DO i=1,N
read(11,25)x(i)
x1(i)=x(i)
read(12,25)px(i)
px1(i)=px(i)
read(13,25)z(i)
z1(i)=z(i)
read(11,25)pz(i)
pz1(i)=z(i)
write(15,26)x1(i),px1(i),z1(i),pz1(i)

25 format(1X,F20.11)
26 format(1X,F15.6,6F15.5)
end do

end program

代码如下,但是N=10000时编译上没有问题,但是运行不了,最大的运行只能N=5000可以正常读出数据,其中x,px,z,pz都是10000行的数据,按理说不应该存在越界的问题。那么请问,问题到底出在哪里呢??
分享到:  微信微信
收藏收藏 点赞点赞 点踩点踩

490

帖子

4

主题

0

精华

大宗师

F 币
3298 元
贡献
1948 点

水王勋章元老勋章热心勋章

沙发
发表于 2015-9-23 14:01:05 | 只看该作者
1、format为啥要定义在循环里面?数据的列数和格式会改变?
2、为啥有open无close?

724

帖子

4

主题

0

精华

大师

农村外出务工人员

F 币
653 元
贡献
335 点

新人勋章爱心勋章水王勋章元老勋章热心勋章

板凳
发表于 2015-9-23 17:50:48 | 只看该作者
求助时,最重要的三个信息:错误提示 > 代码 > 操作系统及编译器版本

26

帖子

8

主题

0

精华

熟手

F 币
123 元
贡献
78 点
地板
 楼主| 发表于 2015-9-24 10:45:33 | 只看该作者
楚香饭 发表于 2015-9-23 17:50
求助时,最重要的三个信息:错误提示 > 代码 > 操作系统及编译器版本

无编译错误,但是运行不了,操作系统WIN7,编译器IVF+VS2010

26

帖子

8

主题

0

精华

熟手

F 币
123 元
贡献
78 点
5#
 楼主| 发表于 2015-9-24 10:46:29 | 只看该作者
pasuka 发表于 2015-9-23 14:01
1、format为啥要定义在循环里面?数据的列数和格式会改变?
2、为啥有open无close?
...

FORMAT 可以尝试放在循环外面,但依然不行,数据列数格式不变

有OPEN 无close 是忘写了,不过应该不影响运行

490

帖子

4

主题

0

精华

大宗师

F 币
3298 元
贡献
1948 点

水王勋章元老勋章热心勋章

6#
发表于 2015-9-24 12:45:27 | 只看该作者
yifanxing1992 发表于 2015-9-24 10:46
FORMAT 可以尝试放在循环外面,但依然不行,数据列数格式不变

有OPEN 无close 是忘写了,不过应该不影响 ...

求助时,最重要的三个信息:错误提示 > 代码 > 操作系统及编译器版本
~~~~~~重要的事情说三遍~三遍~三遍
执行读写文件操作有open无close是非常糟糕的陋习,加一个close语句的成本远小于运行出错后再去找bug

2022

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
1598 元
贡献
689 点

美女勋章热心勋章星光勋章新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

7#
发表于 2015-9-24 13:09:06 | 只看该作者
运行不了总有点现象可以供你描述的吧?

26

帖子

8

主题

0

精华

熟手

F 币
123 元
贡献
78 点
8#
 楼主| 发表于 2015-9-24 21:13:28 | 只看该作者
pasuka 发表于 2015-9-24 12:45
求助时,最重要的三个信息:错误提示 > 代码 > 操作系统及编译器版本
~~~~~~重要的事情说三遍~三遍~三遍
...

恩恩,问题解决了~~谢谢好心大神!

26

帖子

8

主题

0

精华

熟手

F 币
123 元
贡献
78 点
9#
 楼主| 发表于 2015-9-24 21:14:01 | 只看该作者
fcode 发表于 2015-9-24 13:09
运行不了总有点现象可以供你描述的吧?

问题自己解决了耶,谢谢管理大大

724

帖子

4

主题

0

精华

大师

农村外出务工人员

F 币
653 元
贡献
335 点

新人勋章爱心勋章水王勋章元老勋章热心勋章

10#
发表于 2015-9-24 21:31:28 | 只看该作者
为什么不说出来让大家借鉴呢?
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-11-23 03:07

Powered by Tencent X3.4

© 2013-2024 Tencent

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