fortran读取有缺失的数据问题

daydaycode

data.txt (12.64 KB, 下载次数: 3)

2020-3-24 16:18 上传

点击文件名下载附件
样本数据

fortran小白请教下fortran处理txt文件问题。需要处理的数据为地面自动站气象数据，数据包含的数据类型主要有字符串和浮点及整数等类型格式，数据为62247*212的矩阵，但是在57418行到最后的前6列都没有值，这样在数据读取的时候就会提示读取文件末尾报错，请问想正确读入这类型的数据，该怎么编写fortran代码。自己参考别人的代码，写的fortran程序，不能处理以上的问题，代码及数据和报错信息如下，还请各位前辈指教下怎么修改

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

PROGRAM predeal_dimianzidong
IMPLICIT NONE

INTEGER                                ::i,j                !读取数据的行数
INTEGER, PARAMETER      ::row=62247,col=212
CHARACTER,DIMENSION(row,col)     ::data_read,data_1

OPEN(UNIT=3, FILE='SURF_CHN_MUL_HOR_2019080100.txt', STATUS='OLD',ACTION='READ',FORM='FORMATTED',ACCESS='sequential')        !打开文件
!READ(3,*) ((data_read(i,j),j=1,col),i=1,row)
DO i = 1,row
READ(3,*)(data_read(i,j),j=1,col)
END DO
!data_1 = transpose(data_read) 
print*,SHAPE(data_read)
!PRINT*,data_1
CLOSE(UNIT=3)        !关闭文件
STOP
END PROGRAM predeal_dimianzidong

Rookie

Rookie 发表于 2022-4-22 08:11
好的。关于字符串的缺失，我的电脑直接输出空白。所以我在派生数据的声明阶段，直接初始化了。
PS.我所编 ...

刚刚发现是txt文档数据结尾有两个空白行，才会导致iostat仍为0.

Rookie

fcode 发表于 2022-4-20 17:14
没有统一规则。看具体“保存该文件的程序”的策略。

好的。关于字符串的缺失，我的电脑直接输出空白。所以我在派生数据的声明阶段，直接初始化了。
PS.我所编译的程序，在读取外部文件的TXT文件时候，到达文档结尾的时候，IOSTAT仍为0.不知道前辈们有没有遇到过类似的情况？

fcode

Rookie 发表于 2022-4-20 14:54
请问这种缺失的数据，读取的时候会保存为什么？

没有统一规则。看具体“保存该文件的程序”的策略。

Rookie

vvt 发表于 2020-3-24 17:25
你可以看看你的数据是如何得来的，是否可以转存为其他数据格式？比如 csv，用逗号间隔开。

比如

请问这种缺失的数据，读取的时候会保存为什么？

daydaycode

li913 发表于 2020-3-25 17:50
有的，更麻烦。读取整行，以tab为分隔符，将一行分为多个记录，存入数组。 ...

好的，多谢。

li913

daydaycode 发表于 2020-3-25 13:47
十分感谢，现在问题解决了。请问，像这种情况只能把txt转成csv吗，还有别的方法吗？ ...

有的，更麻烦。读取整行，以tab为分隔符，将一行分为多个记录，存入数组。

daydaycode

vvt 发表于 2020-3-24 17:25
你可以看看你的数据是如何得来的，是否可以转存为其他数据格式？比如 csv，用逗号间隔开。

比如

好的，多谢。

daydaycode

li913 发表于 2020-3-25 10:56
看起来你的文件中数据是以tab分隔的，你可以用excel打开数据，检查一下是否有误，然后另存为csv格式。读取c ...

十分感谢，现在问题解决了。请问，像这种情况只能把txt转成csv吗，还有别的方法吗？

li913

看起来你的文件中数据是以tab分隔的，你可以用excel打开数据，检查一下是否有误，然后另存为csv格式。读取csv格式就方便了。

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

program test
  implicit none
  integer i
  INTEGER, PARAMETER      ::row=11,col=212
  character(40) dat(col-6, row)
  dat = 'undefined'
  open(10,file='1.csv')
  read(10,*)
  do i = 1, row
    read(10,*) dat(:,i)
  end do
  close(10)
end program