文件中有空行的fortran处理方法

mofangbao

今天1群有朋友问了这么个问题：
    现有一个数据文件，文件的列数固定，但是不知道会在哪一行出现空行，默认情况下，fortran读取到空行（该朋友的文件中空行中是有空格的）会自动跳过直到读取到数据为止，这样的话如果使用的是固定循环数量的话，就会出现实际循环数小于指定循环数的情况，最终就是 -end of file encountered

这个是他的示例数据：
53694 1960 1 1 -24 79 -89
53694 1960 1 2 15 97 -56
53694 1960 1 3 0 79 -41
53694 1960 1 4 40 101 -29
53694 1960 1 5 -13 15 -44
53694 1960 1 6 -41 37 -107
                                                   
0        0     0       0        0       1         2

他的要求是：空出的行跳过去，该行的数组赋值为0，如果用fortran默认的方法读取，也就是下面这个他给我的程序：
      program yy12
      integer i,titol,station,year(8),month(8),day(8)
     &,htmin(8),htmax(8),a(8)
open(222,FILE='test.txt',FORM='FORMATTED')
do i=1,8
      read(222,*) station,year(i),month(i),day(i),htmin(i),htmax(i),a(i)
print*
      print*,i
      print*,station,year(i),month(i),day(i),htmin(i),htmax(i),a(i)
enddo
700   format(i5,3x,i5,x,i4,4x,i2,3x,i2,3x,i5)
      close(222)
      end
结果将会是这样：


原因就是上面说的fortran默认读取的问题，因此要解决就需要先判断该行是否为空行，我的方法也不是很完美，只是根据第一列进行了判断，程序如下：
      program yy12
      integer i,titol,year(8),month(8),day(8)
     &,htmin(8),htmax(8),a(8)
character*5 station
open(222,FILE='test.txt')
i=1
do while(1)
read(222,'(a)',end=100,advance='no')station
if(trim(station)=="")then
  read(222,*)
  year(i)=0
  month(i)=0
  day(i)=0
  htmin(i)=0
  htmax(i)=0
  a(i)=0
  i=i+1
  cycle
endif
      read(222,*)year(i),month(i),day(i),htmin(i),htmax(i),a(i)
      !print*,station,year(i),month(i),day(i),htmin(i),htmax(i),a(i)
i=i+1
enddo
      100 close(222)
      end
我把整型的站号修改为字符型，然后来判断是不是为空，这样算是解决了他的数据格式问题，但是还有很多毛病：
1、如果空行什么数据都没有，只有一个回车，那么将会认为文件到空行就结束了；
2、如果正常情况下站号之前就存在空格，那么将会不读取任何数据
有没有朋友有更好的算法呢？

传说中的谁

fortran我不懂，我以前用VB读的时候，跳过空行的一个思路是取得改行的字符长度，如果为0就可以判定是空的。

mofangbao

传说中的谁 发表于 2011-7-27 22:48

                               
登录/注册后可看大图

fortran我不懂，我以前用VB读的时候，跳过空行的一个思路是取得改行的字符长度，如果为0就可以判定是空的。

fortran对字符串处理能力不足，很纠结啊

godenflame135

请看看我的方法可以不：
设计思想：先整体读入一行（79个字符），循环判断非空格就判定为有数据，并接着将该行分散读取到年月日等变量。
下面是我修改的程序，请参考,实际情况相应修改。
！############################################
前提：空行中有空格。
注意：1.在有可能读入数据为负数时，读取方式要稍加修改。
           2.如果单个文件的行数特别多，建议将所有数据一次性读入字符型数组中，再分行判断，这样运行速度可在一定程度
              上提升。
!#############################################

如有不当，请回帖讨论交流。

!-------------------------------------------------------------------------------------------------------------------------------------------
  program yy12
      implicit none
      integer i,titol,station,year(8),month(8),day(8)
     &,htmin(8),htmax(8),a(8)
    character*79      :: c_line    !for read one row
    integer               :: i,j,k      !for loop if need
    logical                 :: L_have_number
open(222,FILE='test.txt',FORM='FORMATTED',err=998)
j=0      !index of  station,year...

! begin loop to read------------------------------------------------
do while()     !----------exit by going to 999 when read-----
799   read(222,'(a79)',err=999)c_line
          L_have_number=    .false.
!   ---loop to judge have number or not--------------------------
    do i=1,79
          if(ichar(c_line(i:i))  /= 32      )  then     !   any one which is not space
               L_have_number  =  .true.
                go to 800
          end if
     end do
!   ----end judge have number--------------------------------------
!judge to write data to station(j),year(j)....
       j=j+1
800  if(L_have_number)then   
              write(c_line(1:5),*)      station (j)
              write(c_line(7:10),*)    year    (j)
              write(c_line(12:12),*)  month(j)
              write(c_line(13:13),*)  day     (j)
             if(ichar(15)  ==  45   )then                      !   '-' sign   ,htmin(j)<100,means have only 2 digitals.
                 write(c_line(15:17),*)htmin(j)
                 k=17
              else                                                        ! no '-'  sign,htmin(j)<100,means have only 2 digitals                 
                  write(c_line(15:16),*)htmin(j)   
                  k=16      
              end if
             if (ichar( c_line(k+4))  /=32    )     then    !htmax(j) >100,means have only 3 digitals.
                 write(c_line(k+1:k+4),*)htmax(j)  
                 k=k+4
            else                                                         !htmax(j)<100,means have only 2 digitals.   
                 write(c_line(k+1:k+3),*)htmax(j)     
                 k=K+3
             end if
          if(ichar(k+4) /= 32 )then                           !  all of a(j) have  '-' sign ,-999<a(j)<-100.
                   write(c_line(k+1:k+4),*)a(j)
                   k=k+4
          else                                                           ! all of a(j) have '-' sign ,-100<a(j)<-10.
                   write(c_line(k+1:k+3),*)a(j)
                    k=K+3
           end if
          ! -----end of write a line------------------------

           print*,station(j),year(j),month(j),day(j),htmin(j),htmax(j),a(j)

         !----end of print each number-----------------
             go to 799        !continue to read next row

    else                                                            ! have no number,no write and read next row.
          station(j)  =  0
          year(j)      =  0
          month(j)   =  0
          day(j)        =  0
          htmin(j)     =   0
          htmax(j)    =   0
          a(j)            =   0
            
             print*,station(j),year(j),month(j),day(j),htmin(j),htmax(j),a(j)

            !----end of print each number-----------------

          go to 799     !continue to read next row

        end if    !------end of judge to write  number or not-----------------
end do   !----------exit by going to 999 when read-----
go to 1000     ! ----------go to the end

!  error process--------------------------------------------------
999 print *,'read error, or EOF'
      go to 1000
998 print *,'open error  ,  or file not exist.'
1000 continue

close(222)

      end

!----------------------------------------------------------------------------------------------------------------------------------------

mofangbao

godenflame135 发表于 2011-7-31 22:48

                               
登录/注册后可看大图

请看看我的方法可以不：
设计思想：先整体读入一行（79个字符），循环判断非空格就判定为有数据，并接着将 ...

调试中发现一点问题，
1、重复定义了i变量，station不知道你是不是定义为数组的，上下有一些不一致
2、程序第38行附近出现编译错误，大概是操作的参数有问题，不知道是不是应该把：
if (ichar( c_line(k+4))  /=32    ) 改为：if (ichar( k+4)  /=32    )
3、由于没有调试成功，还没发进一步测试
呵呵 不错你的想法挺不错的 不知道你用的哪个fortran版本，出结果了没？

757084134

还好可以用

Trinity

最近也遇到个文件中有几行不需要的情况。不知上面那位的方法能用不

wangxue910402

{:eb502:}{:eb502:}{:eb502:}

黑夜的北极星

如果我的数据是这样的
12,12,12,12,12， ，12,12
11,11,11,11， ，11,11,11
这样用fortran怎么读取呢？

qswl

cycle用的不错