linux下运行fortran程序，不知道错在哪里

xiaohundanpipi

open(unit=4,file='hgt500hpa.mon.dat',status='old',access='direct',form='unformatted',recl=1085)
read(4,*)
do t=1,nt
do y=1,ny
  do x=1,nx
read(4,*) hgt(x,y,t)
  end do
end do
end do
close(1)


运行后显示错误At line 7 of file season_mean.f90 (unit = 4, file = 'hgt500hpa.mon.dat')
Fortran runtime error: Format present for UNFORMATTED data transfer
我不知道哪里错了，这个格式不能这样读取吗？
请各位帮帮我看看

bustop

你应该想直接方式读取无格式文件。首先确定recl的值，这要依据你的记录长度了，然后 read（4，irec）irec代表ith记录。

xiaohundanpipi

bustop 发表于 2012-9-13 20:50

                               
登录/注册后可看大图

你应该想直接方式读取无格式文件。首先确定recl的值，这要依据你的记录长度了，然后 read（4，irec）irec代 ...

嗯嗯，我试试，thank you～

xiaohundanpipi

program one
implicit none
integer,parameter::nx=35,ny=21,nt=776,iyear=63
real,dimension(nx,ny,nt)::hgt
real,dimension(nx,ny,iyear)::hgtwinter
integer::i,j,k,t,irec

open(1,file='hgt500hpa.mon.dat',status='old',form='unformatted',access='direct',recl=35*21)
irec=0
do t=1,nt
irec=irec+1
read(1,rec=irec)((hgt(i,j,t),i=1,nx),j=1,ny)
end do
close(1)

这次出现的错误是I/O past end of record on unformatted file   格点范围肯定都对，我编写了ctl用grads对hgt500hpa.mon.dat'画过图的，这个错误是说我数组越界吗？ 这个dat是用grads从nc文件转过来的，到底是无格式的还是有格式的。。。好纠结，好烦，求帮助

bustop

xiaohundanpipi 发表于 2012-9-14 11:08

                               
登录/注册后可看大图

program one
implicit none
integer,parameter::nx=35,ny=21,nt=776,iyear=63

给我看你的文件，xucp11@lzu.edu.cn

xiaohundanpipi

bustop 发表于 2012-9-17 10:33

                               
登录/注册后可看大图

给我看你的文件，

你好，程序已经给你发过去了。
我换到xp上运行了一下，程序就对了，不知道到底是不是linux有什么特殊的要求
这是我求季平均的fortran程序，大家共勉。
program one
implicit none
integer,parameter::nx=35,ny=21,nt=776,iyear=64
real,dimension(nx,ny,nt)::hgt
real,dimension(nx,ny,iyear)::hgtwinter
integer::i,j,k,t,irec



open(1,file='hgt500hpa.mon.dat',status='old',form='unformatted',access='direct',recl=35*21)
do t=1,nt
  read(1,rec=t)(( hgt(i,j,t),i=1,nx),j=1,ny)      
end do
close(1)
print*,hgt(35,21,776)


do j=1,ny
  do i=1,nx
    do t=12,770,12
          if (hgt(i,j,t)/=-9.96921e+36 .and. hgt(i,j,t+1)/=-9.96921e+36 .and.hgt(i,j,t+2)/=-9.96921e+36) then
         hgtwinter(i,j,(t-12)/12+1)=-9.96921e+36
          else
         hgtwinter(i,j,(t-12)/12+1)=(hgt(i,j,t)+hgt(i,j,t+1)+hgt(i,j,t+2))/3
          end if
    end do
  end do
end do


print*,hgtwinter(30,2,:)

open(2,file='hgtwinter.mean.dat',status='replace',form='unformatted',access='direct',recl=35*21)
irec=0
do t=1,64
irec=irec+1
write(2,rec=irec)((hgtwinter(i,j,t),i=1,nx),j=1,ny)
end do
close(2)
print*,'zz'

end program one

bustop

xiaohundanpipi 发表于 2012-9-17 14:43

                               
登录/注册后可看大图

你好，程序已经给你发过去了。
我换到xp上运行了一下，程序就对了，不知道到底是不是linux有什么特殊的要 ...

我也试了的确不能linux，但将recl改成35*21*4就可以。这是因为xp系统同linux系统编译器处理数据的位数不同导致。若还有问题，发送我邮箱。

xiaohundanpipi

bustop 发表于 2012-9-17 16:46

                               
登录/注册后可看大图

我也试了的确不能linux，但将recl改成35*21*4就可以。这是因为xp系统同linux系统编译器处理数据的位数不同 ...

嗯嗯 谢谢啦。。。你真是个好人～