如何在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
答案 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解析输入行。
注意:
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
请注意,我转换了数组,以便分配到连续的内存位置。