grads提取数据出现： Low Level I/O Error: Seek error on data file

Garfield

grads提取数据出现： Low Level I/O Error: Seek error on data file

mofangbao

先自己找原因，一个问题的出现原因是很多的，你就贴一个报错信息谁能确定是什么原因呢，只能从字面给你翻译，那就是输入输出错误，数据定位错误

易小凯

是不是文件的路径有问题啊。
GIVE MORE INFORMATION, PLEASE!

韦笑鸥

我也是在用fortran读取grads分割出的dat数据时出现了楼主一样的问题
下面是源程序，麻烦各位大大帮忙看看，数据源是nc分割出的144*73*735的dat格式uwnd数据：
program www
integer i,j,t
real uwnd(144,73,735)

OPEN(14,access='direct',ACTION='READ',FILE='uwnd925.mean.dat',FORM='unformatted',RECL=1)

do t=1,735
   do j=1,73
     do  i=1,144
     read (14,*) uwnd(i,j,t)         
         end do
   end do
  end do
  write (14,*) uwnd(1,1,1)
  Close(14)
  stop
end

Garfield

韦笑鸥 发表于 2012-3-26 23:03

                               
登录/注册后可看大图

我也是在用fortran读取grads分割出的dat数据时出现了楼主一样的问题
下面是源程序，麻烦各位大大帮忙看看， ...

你打开文件时用了记录号，下面读写的时候没把记录号加1。
还有.dat的文件read\write不用写*

韦笑鸥

Garfield 发表于 2012-3-27 20:12

                               
登录/注册后可看大图

你打开文件时用了记录号，下面读写的时候没把记录号加1。
还有.dat的文件read\write不用写*

还是没弄对，应该还是open里的设置和输出语句设置的问题，请教几个问题
1、问用普通的从nc里用grads命令sdfopen分离出来的dat文件储存方式和其他asc2码的dat格式一样吗？
2、用了一个读取nc的fortran程序读取nc文件发现读出的数据和用grads读出的除了缺测处显示一致外，发现其他其他格点不一致，比如某格点grads读书-0.260002，fortran读出的是-20490，我测试过，格点是对应的
，还有时间维的计数方式也不同，这样由于存储（或显示？）方式不同造成的差异（？）会不会影响计算结果呢？（还是，两边都没错，只是grads显示的时候，自动进行了转换，而fortran读出的是最原始的储存方式？）

ps：记录号是指哪一个？（初学者，麻烦了，不胜感谢！）

韦笑鸥

明白了，是输出命令里的recl问题，open里设定了recl=1，读写命令里recl务必与之相匹配，在循环你加入recl=recl+1，就行了

花小菜

请问楼主后来怎么解决的啊，我碰到跟你一样的问题了

killerlrving

顶~~~~~~~~~~~~~~~~~~~~~

胡小顺

program main  
implicit none
parameter n=348
real x(n),y(n),c(n),d(n),t(n)
real a,b,sx,sxy,sum2,ave2,sum1,ave1
integer i
!open(1,file='d:/interannual/q/nino3.txt')
!open(2,file='d:/interannual/q/tloop79.txt')
!open(3,file='d:/interannual/q/3.txt')
!open(4,file='d:/interannual/q/4.grd',form='binary')
!open(5,file='d:/interannual/q/5.grd',form='binary')
open(1,file='f:/wwt/lr/nino1+2.txt')
open(2,file='f:/wwt/lr/ssta.grd',form='unformatted',status='old',access='direct',recl=1)
open(3,file='f:/wwt/lr/y1.grd',form='unformatted',status='replace',access='direct',recl=1)
open(4,file='f:/wwt/lr/y2.grd',form='unformatted',status='replace',access='direct',recl=1)
open(5,file='f:/wwt/lr/t.grd',form='unformatted',status='replace',access='direct',recl=1)
do i=1,n
read(1,*) x(i)
read(2,rec=i) y(i)
enddo
do i=1,n
write(3,rec=i) y(i)
enddo
sum1=0.0
sum2=0.0
sx=0.0
sxy=0.0
do i=1,n
sum1=sum1+x(i)
sum2=sum2+y(i)
enddo
ave1=sum1/n
ave2=sum2/n
do i=1,n
c(i)=x(i)-ave1
d(i)=y(i)-ave2
enddo
do i=1,n
sx=sx+c(i)**2
sxy=sxy+c(i)*d(i)
enddo
sx=sx/n
sxy=sxy/n
b=sxy/sx
a=ave2-ave1*b
do i=1,n
t(i)=a+b*x(i)
y(i)=y(i)-t(i)
enddo
do i=1,n
write(4,rec=i) y(i)
write(5,rec=i) t(i)
enddo
end
为什么这样算出来的数据在画图的时候老是出错呢
类似这种
Low Level I/O Error:  Read error on data file    Data file name = 20091012_RCWF.dat    Error reading 141 bytes at location 79524 Data Request Error:  Error for variable 'dbz'   Error ocurred at column 1  DISPLAY error:  Invalid expression   Expression = dbz
ctl：dset f:/wwt/lr/y1.grd
title y1
undef -9.99e+08
xdef 1 linear 80.0 1.0
ydef 1 linear -15.0  1.0
zdef 1 linear 0 1.0
TDEF 348 LINEAR jan1979 1mo
vars 3
y1 0 99 yuanshi
y2 0 99 huigui
t 0 99 huigui
endvars
gs：'reinit'
'open f:\wwt\lr\y1.ctl'
'set lon 80'
'set lat -15'
'set lev 0'
'set t 1 348'
'set gxout line'
'set ccolor 4'
'd y1'
'printim f:/wwt/lr/y.png white'
;