如何在Fortran中读取我的csv文件?

时间:2017-09-27 09:36:15

标签: csv fortran

如何在Fortran中阅读以下数据。我无法读取此文件中的数据。我收到了以下错误

severe (64): input conversion error, unit 18

Image              PC                Routine            Line        Source             
a.out              00000000004734AA  Unknown               Unknown  Unknown

这是我的代码:

        program dataread
        implicit none
        character*15 ::head_1,head_2,head_3,head_4,head_5,head_6
        character*15:: B,C, head_7,head_8,head_9,head_10
        real,dimension(1:71385,1:10)::A
        integer::i,j

        open(unit=18, file='Tws15thHourlyData.csv' , status='old',
     &  access ='sequential',form='formatted')!,recl=71781*10)

        read(18,*) head_1,head_2,head_3,head_4,head_5,head_6,
     &          head_7,head_8, head_9,head_10


       do i=2,71385
       read(18,300)(A(i,j),j=1,10)
300    format(I5,A17,2F9.6,A8,5F4.1)
       end do

        print*, head_1,head_2,head_3,head_4,head_5,head_6,
     &          head_7,head_8,   head_9,head_10
       do i=2,71385
       read(18,300)(A(i,j),j=1,10)
300    format(I5,A17,2F9.6,A8,5F4.1)
       end do

        print*, head_1,head_2,head_3,head_4,head_5,head_6,
     &          head_7,head_8,   head_9,head_10


        do i=2,71385
              print*, (A(i,j),j=1,10)
        end do
        close(18)
        open(unit=28,file='14data.txt')
        write(28, 100),((A(i,j),j=1,10),i=2,71385)
100     format(10(71385(I4,A15,2F9.6,A8,5F4.1,2x),/))
        end program

此数据我试图阅读

TWSCODE,DISTRICT,LATITUDE_DD,LONGITUDE_DD,RECORDED_DATE,HOUR,TEMPERATURE,HUMIDITY,WIND_SPEED,WIND_DIRECTION

109,KALABURAGI,17.463587,77.42,14-08-17,0,26.2,79.4,0,168
109,KALABURAGI,17.463587,77.42,14-08-17,0,26.2,80,0,25

109,KALABURAGI,17.463587,77.42,14-08-17,0,26.1,80.4,0,25
109,KALABURAGI,17.463587,77.42,14-08-17,0,25.9,81,0,25

109,KALABURAGI,17.463587,77.42,14-08-17,1,25.8,81.7,0,25
109,KALABURAGI,17.463587,77.42,14-08-17,1,25.9,82,0,287
109,KALABURAGI,17.463587,77.42,14-08-17,1,25.9,82.5,0,299
109,KALABURAGI,17.463587,77.42,14-08-17,1,25.8,82.8,0,286
109,KALABURAGI,17.463587,77.42,14-08-17,2,25.6,83.5,0,254
109,KALABURAGI,17.463587,77.42,14-08-17,2,25.6,83.9,0,292
109,KALABURAGI,17.463587,77.42,14-08-17,2,25.6,84,0,299
109,KALABURAGI,17.463587,77.42,14-08-17,2,25.6,84.2,0,309
109,KALABURAGI,17.463587,77.42,14-08-17,3,25.5,84.4,0,327
109,KALABURAGI,17.463587,77.42,14-08-17,3,25.4,84.8,0,315
109,KALABURAGI,17.463587,77.42,14-08-17,3,25.3,84.8,0,305
109,KALABURAGI,17.463587,77.42,14-08-17,3,25.2,84.8,0,306
109,KALABURAGI,17.463587,77.42,14-08-17,4,25.2,84.9,0,305
109,KALABURAGI,17.463587,77.42,14-08-17,4,25.2,85.1,0,306
109,KALABURAGI,17.463587,77.42,14-08-17,4,25.1,85.3,0,305
109,KALABURAGI,17.463587,77.42,14-08-17,4,25.1,85.4,0,305
109,KALABURAGI,17.463587,77.42,14-08-17,5,25.1,85.7,0,308

2 个答案:

答案 0 :(得分:4)

我不喜欢@ agentp的回答,并且不同意这样的评论,即应该将该行读作字符串并将其解析为混乱。

所以我写了这个......

正如我在上面的评论中所指出的,OP的代码中的一个错误(可能是唯一的错误)是使用此语句,并格式化

       read(18,300)(A(i,j),j=1,10)
300    format(I5,A17,2F9.6,A8,5F4.1)

尝试将数字和字符串的混合读入reals数组。这永远不会好起来。退后一步,让事情有所考虑。输入文件包含一些结构良好的数据,那么为什么不定义一个很好的存储结构呢?首先是类型定义,如:

  TYPE :: met_record
     INTEGER :: TWSCODE
     CHARACTER(len=32) :: DISTRICT
     REAL :: LATITUDE_DD
     REAL :: LONGITUDE_DD
     CHARACTER(len=8) :: RECORDED_DATE
     INTEGER :: HOUR
     REAL :: TEMPERATURE
     REAL :: HUMIDITY
     REAL :: WIND_SPEED
     REAL :: WIND_DIRECTION
  END TYPE met_record

然后是那些

的数组
TYPE(met_record), DIMENSION(71385) :: weather_reports

现在读取数据非常容易......

do i = 1, 71385
    read(18,*) weather_reports(i)
end do

看!看不到明确的格式。正确声明变量并让Fortran解析输入行。

注意:

  • 我写的是相对现代的Fortran,并且没有固定格式源文件的卡车。
  • 我没有做任何阅读文件标题的安排,已有代码。
  • 撰写met_record可以像write(28,*) weather_reports(i)
  • 一样简单

答案 1 :(得分:0)

这是一种方法:

 integer,parameter :: n=71385
 character*15 :: strings(10,n)
 real :: a(7,n)
 integer :: tw(n)

..

 do i=1,n
   read(18,*)strings(:,i)
   read(strings(1,i),*)tw(i)
   read(strings(3:4,i),*)a(1:2,i)
   read(strings(6:10,i),*)a(3:7,i)
 end do

请注意,我转换了数组,以便分配到连续的内存位置。