读取文件中的特定值
我要读取的数据文件是这个样子的No Well Type Description Oil Rate Water Rate Gas Rate FBHP
(m3 /day) (m3 /day) (m3 /day) (kPa)
---------------------------------------------------------------------------------------------
No Well Type Description Oil Rate Water Rate Gas Rate FBHP
(m3 /day) (m3 /day) (m3 /day) (kPa)
---------------------------------------------------------------------------------------------
No Well Type Description Oil Rate Water Rate Gas Rate FBHP
(m3 /day) (m3 /day) (m3 /day) (kPa)
---------------------------------------------------------------------------------------------
No Well Type Description Oil Rate Water Rate Gas Rate FBHP
(m3 /day) (m3 /day) (m3 /day) (kPa)
---------------------------------------------------------------------------------------------
Time = 32.00000 days
1 P1 6 Liquid Production RateSrf 14.81 15.19 0.000 0.000 1.026 50.63 1.1939E+04
Time = 32.00000 days
2 P2 6 Liquid Production RateSrf 14.90 15.10 0.000 0.000 1.013 50.32 1.1984E+04
Time = 32.00000 days
3 P3 6 Liquid Production RateSrf 14.86 15.14 0.000 0.000 1.019 50.47 1.1902E+04
Time = 32.00000 days
4 P4 6 Liquid Production RateSrf 14.96 15.04 0.000 0.000 1.005 50.13 1.1943E+04
Time = 32.00000 days
5 P5 6 Liquid Production RateSrf 14.65 15.35 0.000 0.000 1.048 51.18 1.2031E+04
No Well Type Description Oil Rate Water Rate Gas Rate FBHP
(m3 /day) (m3 /day) (m3 /day) (kPa)
---------------------------------------------------------------------------------------------后面都是重复这样的格式,不同的是下面这段
No Well Type Description Oil Rate Water Rate Gas Rate FBHP
(m3 /day) (m3 /day) (m3 /day) (kPa)
---------------------------------------------------------------------------------------------出现的次数不确定.
需要读取的数据是下面这段
Time = 32.00000 days
1 P1 6 Liquid Production Rate Srf 14.81 15.19 0.000 0.000 1.026 50.63 1.1939E+04中的Time后的取值,以及第二行的P1字符串和Srf后面的实数.
我的思路是首先检查某行的开头字符是否是Time,如果不是跳过这行检查下一行.如果是,则进行本行和下一行的读取操作.
我的问题是读取文件时的读取位置怎么控制,这个检查操作改怎么写?
谢谢大神!
Open 打开文件
循环 Read ，每次 read 读一行，如果不满足条件就继续下一次循环。
Backspace 回退一行。
Rewind 直接回退到起始位置。
顺序有格式文件，大概就这些控制方法。
你的思路非常正确。按照这个思路去做就可以了。
good luck!
fcode 发表于 2015-4-14 17:16
Open 打开文件
循环 Read ，每次 read 读一行，如果不满足条件就继续下一次循环。
Backspace 回退一行。
我还有一个问题想问，在读的时候，Srf后面的值可能是小数表示，也可能是科学计数E表示，写读入格式的时候怎么样才能两种数据表现都能读入。 fcode 发表于 2015-4-14 17:16
Open 打开文件
循环 Read ，每次 read 读一行，如果不满足条件就继续下一次循环。
Backspace 回退一行。
是用G编辑符吗？ 找到 Srf 的位置，然后对后面的数据 read( str ,*) 我写了循环来读取数据，但是总是读取出错，error值大于0，进入死循环，求帮忙看下，辛苦大神。
PROGRAM READPLOT
IMPLICIT NONE
INTEGER :: ERROR
CHARACTER(4) :: CTIME ! 检查是否为Time
DOUBLE PRECISION :: DTIME ! 临时存储Time后的值
INTEGER, PARAMETER:: MAXTN = 1000, MAXWN = 20 ! 最大存储1000个时间步，20口井的数据
DOUBLE PRECISION :: D(MAXTN) ! 存储Time后的所有值
CHARACTER(10) :: WELLNAME(MAXWN),WELL ! 存储井名，最大十个字符
INTEGER :: IT,IWELL ! 时间步数，井数
DOUBLE PRECISION :: OILRATE(MAXWN,MAXTN),WATERRATE(MAXWN,MAXTN) ! 油水产量
DOUBLE PRECISION :: GASRATE(MAXWN,MAXTN),GOR(MAXWN,MAXTN) ! 气产量，气油比
DOUBLE PRECISION :: WOR(MAXWN,MAXTN),WCUT(MAXWN,MAXTN),FBHP(MAXWN,MAXTN) ! 油水比，含水率，井底压力
DOUBLE PRECISION :: OR,WR,GR,G,WO,WC,FB ! 临时存储以上数据
INTEGER :: MXWN,MXTN ! 存入的最大数据个数
INTEGER :: I,J! 循环变量
! 占位变量，不进行存储
CHARACTER(4) :: AA
INTEGER :: BB
CHARACTER(27) :: CC
! 初始化
IT = 1
D = 0.0
ERROR = 0
WELLNAME = ' '
! 读取数据
OPEN(UNIT = 415, FILE = 'plot.out')
DO WHILE(ERROR >= 0) ! 检查读取文件是否结束
READ(415,IOSTAT = ERROR) CTIME,AA,DTIME,AA,BB,AA,BB
IF(ERROR == 0) THEN ! 如果文件读取正常进行以下操作
IF(CTIME == 'Time') THEN
IF(IT == 1) THEN ! 如果是第一个时间步数据
D(IT) = DTIME
READ(415,IOSTAT = ERROR) IWELL,WELL,BB,CC,OR,WR,GR,G,WO,WC,FB
WELLNAME(IWELL) = WELL
OILRATE(IWELL,IT) = OR
WATERRATE(IWELL,IT) = WR
GASRATE(IWELL,IT) = GR
GOR(IWELL,IT) = G
WOR(IWELL,IT) = WO
WCUT(IWELL,IT) = WC
FBHP(IWELL,IT) = FB
ELSE
IF(D(IT) == DTIME) THEN ! 是否相同的时间步不同的井数据
READ(415,IOSTAT = ERROR) IWELL,WELL,BB,CC,OR,WR,GR,G,WO,WC,FB
WELLNAME(IWELL) = WELL
OILRATE(IWELL,IT) = OR
WATERRATE(IWELL,IT) = WR
GASRATE(IWELL,IT) = GR
GOR(IWELL,IT) = G
WOR(IWELL,IT) = WO
WCUT(IWELL,IT) = WC
FBHP(IWELL,IT) = FB
ELSE ! 下一个时间步第一口井数据
IT = IT + 1
D(IT) = DTIME
READ(415,IOSTAT = ERROR) IWELL,WELL,BB,CC,OR,WR,GR,G,WO,WC,FB
WELLNAME(IWELL) = WELL
OILRATE(IWELL,IT) = OR
WATERRATE(IWELL,IT) = WR
GASRATE(IWELL,IT) = GR
GOR(IWELL,IT) = G
WOR(IWELL,IT) = WO
WCUT(IWELL,IT) = WC
FBHP(IWELL,IT) = FB
END IF
END IF
END IF
END IF
END DO
CLOSE(415)
! 输出数据
OPEN(UNIT = 416, FILE = 'plot.dat')
MXWN = 1
MXTN = 1
DO I = 1,MAXTN
IF (D(MXTN) == 0) THEN
EXIT
ELSE
MXTN = MXTN + 1
END IF
END DO
DO I = 1, MAXWN
IF (WELLNAME(I) == ' ') THEN
EXIT
ELSE
MXWN = MXWN + 1
END IF
END DO
DO I = 1,MXWN
WRITE(416,1003) WELLNAME(I)
WRITE(416,1004)
DO J = 1, MXTN
WRITE(416,1005) D(J),OILRATE(I,J),WATERRATE(I,J),GASRATE(I,J), &
&GOR(I,J),WOR(I,J),WCUT(I,J),FBHP(I,J)
END DO
END DO
1003 FORMAT(A10)
1004 FORMAT('TIME',9X,'OILRATE',7X,'WATERRATE',5X,'GASRATE',7X,'GOR',11X, &
& 'WOR',11X,'WCUT',10X,'FBHP')
1005 FORMAT(F10.5,3X,6(G11.4,3X),G11.4)
CLOSE(416)
END PROGRAM
fcode 发表于 2015-4-14 17:16
Open 打开文件
循环 Read ，每次 read 读一行，如果不满足条件就继续下一次循环。
Backspace 回退一行。
求大神帮忙看下我写的程序，多谢多谢 我只修改一段代码，剩下的没动。
加了一个变量：
Character(Len=512) :: cStr !// 读取用
! 读取数据
OPEN(UNIT = 415, FILE = 'plot.out')
Do !// 一直读
READ(415 ,IOSTAT = ERROR) cStr !// 先读入字符串
if ( Error /= 0 ) exit !// 如果读错，则退出循环（文件结束）
read( cStr , * ) cTime
if ( cTime == "Time" ) THEN
IF(IT == 1) THEN ! 如果是第一个时间步数据
D(IT) = DTIME
READ(415,IOSTAT = ERROR) IWELL,WELL,BB,CC,OR,WR,GR,G,WO,WC,FB
WELLNAME(IWELL) = WELL
OILRATE(IWELL,IT) = OR
WATERRATE(IWELL,IT) = WR
GASRATE(IWELL,IT) = GR
GOR(IWELL,IT) = G
WOR(IWELL,IT) = WO
WCUT(IWELL,IT) = WC
FBHP(IWELL,IT) = FB
ELSE
我算出的结果也是：
TIME OILRATE WATERRATE GASRATE GOR WOR WCUT FBHP
0.00000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
楚香饭 发表于 2015-4-15 19:26
我只修改一段代码，剩下的没动。
加了一个变量：
我也是这个结果，说明没有读入数据，哪错了呢？:-dizzy: 自己 debug 跟踪一下，看是哪些判断符合被跳过了
页:
[1]
2