fortran 批处理读取txt格式的数据 ，急切求救！！！！

云大小子

我用来读取数据批处理的fortran程序如下
!written by asu to read pom outdata
      parameter(ix=148,iy=80,iz=1,it=1456)
real t(ix,iy,iz),u(ix,iy,iz),v(ix,iy,iz),s(ix,iy,iz)
character a1*1,a2*2,a3*3,a4*4,filename*100

  open(200,file='uvts20y.grd',form='unformatted',
     &   access='direct',recl=ix*iy)
irec=0
do iday=1,it
  
  if(iday.lt.10)then
   write(a1,'(i1)')iday
    a4='000'//a1
  elseif(iday.lt.100)then
   write(a2,'(i2)')iday
    a4='00'//a2
  elseif(iday.lt.1000)then
   write(a3,'(i3)')iday
                   a4='0'//a3
                    else
                        write(a4,'(i4)')iday
              endif
  filename='pomout0'//a4                   （ !前面这些都是在读取
（文件名  ，前面这些都没错  因为我在这加上 print*，filename 能输出一系列名称  ）                                                            
open(100,file="filename"//".txt", form='formatted',
     $ status='old')                                    
do  iii=1,2
read(100,*)
enddo

          read(100,101)t
101       format(8e12.5)
、 read(100,101)u
、 read(100,101)v
      read(100,101)s
  close(100)                      主要就是这前面一段关于打开文件，读取里面的数据感觉有问题
   write(6,*)iday
  do k=1,iz
   irec=irec+1

   write(200,rec=irec)((t(i,j,k),i=1,ix),j=1,iy)
   
   enddo
   
   do k=1,iz
    irec=irec+1

      write(200,rec=irec)((u(i,j,k),i=1,ix),j=1,iy)
   
   enddo
       do k=1,iz
     irec=irec+1

      write(200,rec=irec)((v(i,j,k),i=1,ix),j=1,iy)
    enddo
   do k=1,iz
   irec=irec+1
    write(200,rec=irec)((s(i,j,k),i=1,ix),j=1,iy)
     
   enddo

     enddo
  
   end
c        以上的程序 编译没错误，执行出项了如下的错误
  需要读取的txt数据里面的数据分布是下面的一个截图：
总体说一下吧 ，每一列代表一个变量，我要把一列的数据读成一个水平的148x80水平范围的grib数据。
各位有谁知道的 或者处理过的 请给小弟我指点指点 本人不胜感激！！

mofangbao

open(100 那一行改成
filename=trim(filename)//'.txt'
open(100,file=''//trim(filename)//'', form='formatted',
     $ status='old')
上面都是单引号

mofangbao

open(100 那一行改成
filename=trim(filename)//'.txt'
open(100,file=''//trim(filename)//'', form='formatted',
     $ status='old')
上面都是单引号

wlzhongouc

错误提示已经说明了 就是清风指出的问题 filename已经是字符型在open中就不需要在加上单引号

云大小子

mofangbao 发表于 2011-12-2 10:30

                               
登录/注册后可看大图

open(100 那一行改成
filename=trim(filename)//'.txt'
open(100,file=''//trim(filename)//'', form='fo ...

改了以后反而有错误提示，这个trim函数是能消除字符串之间的空格吧 但是好像我这个文件名输出的时候中间就没有空格

云大小子

wlzhongouc 发表于 2011-12-2 10:48

                               
登录/注册后可看大图

错误提示已经说明了 就是清风指出的问题 filename已经是字符型在open中就不需要在加上单引号

恩我改成了open(100,file=filename, form='formatted',
     $ status='old'
并且在这前面加了个print*,filename  执行后只打印了一个文件名pomout00001这个就中断了 说明一定在open这个地方执行不下去

mofangbao

云大小子 发表于 2011-12-2 10:49

                               
登录/注册后可看大图

改了以后反而有错误提示，这个trim函数是能消除字符串之间的空格吧 但是好像我这个文件名输出的时候中间就 ...

是首尾的空格 不是中间的空格

云大小子

mofangbao 发表于 2011-12-2 11:07

                               
登录/注册后可看大图

是首尾的空格 不是中间的空格

我改成下面那样了
filename=trim(filename)//'.txt'
open(100,file=trim(filename), form='formatted',
     $ status='old')
但是还是出现一样的问题 就像上面的截图一样

mofangbao

那你先改成固定文件名，看看还出不出错，然后输出你的filename到屏幕看看

言深深

“”filename试试看