Gfortran读无格式文件，写入有格式顺序文件的题目

1000 · 发表于 2020-6-11 23:27:12

这个题是要这样的：

读无格式文件’1948.dat’, ‘1949.dat’, …, ‘2013.dat’，把每个月的数据单独写成有格式顺序文件，以月份的前三个字母命名，分别为’Jan.dat’, ’Feb.dat’, …, ’Dec.dat’, 每个文件2列66行，第一列为年份，第二列为温度数据。（要求：写12个文件时，要使用循环结构，不要复制12段重复性代码）。

然后代码我也写出来了，但是不知道为什么读取总会结束在第七个文件中，这些文件的内容都是我用循环写入的，不应该有某个文件的读写单独有问题，（讲内容写入文件的代码如第一部分，运行正常，不用管file exist，如图一），该题的代码如第二部分，该题的运行结果如图二，实在不知道哪里错了。求大神指点！！！！
第一部分：

[Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode

program main
implicit none
integer::i,s=1948
character(len=5)::dirname
character(len=8)::filename
real,dimension(12)::a
open(10,file='TI.dat',form='unformatted',access='direct',recl=48)
do i=1,66
  read(10,rec=i)a
  print*,a
  write(dirname,'(i4,"s")')s
  write(filename,'(i4,".dat")')s
  call system("mkdir "//dirname)
  open(s,file=trim(dirname)//'/'//filename,form='unformatted',access='direct',recl=48)
  write(s,rec=1)a
  close(s)
  s=s+1
enddo
endprogram

第二部分：

[Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode

character(len=7),dimension(12)::monname
do i=1,66
  write(dirname,'(i4,"s")')dir
  write(filename,'(i4,".dat")')dir
  open(i,file=trim(dirname)//'/'//filename,form='unformatted',access='direct',recl=4)
  do j=1,12
    read(i,rec=j)a(i,j)
  enddo
  close(i)
  dir=dir+1
enddo

do j=1,12
  monname=['Jan.dat','Feb.dat','Mar,dat','Apr.dat','May.dat','Jun.dat','Jul.dat','Aug.dat','Sep.dat','Oct.dat','Nov.dat','Dec.dat']
  open(k,file='/'//monname(j),status='new')
  do i=1,66
    write(k,'(i4,f4.2)')s,a(i,j)
    s=s+1
  enddo
  close(k)
  k=k+1
  s=1948
enddo
endprogram