fortran打开grd文件"access not allowed"

chaseechoes

为了计算北极海冰面积，先将海冰密集度的nc文件转换成grd文件
nc_grd.gs:

'reinit'
'sdfopen f:\HadISST_ICE\HadISST_ice.nc'
'set gxout fwrite'
'set fwrite f:\HadISST_ICE\area.grd'
i=1
while(i<=1701)
'set x 1 360'
'set y 91 180'
'set t 'i
'set z 1'
'd sea_ice_area_fr'
i=i+1
endwhile
'close 1'
'disable fwrite'

用fortran计算海冰总面积：
IMPLICIT NONE
INTEGER,PARAMETER::NX=360,NY=90,NT=1701
REAL FR(NX,NY,NT)
REAL S(NT)
INTEGER I,J,K,P

OPEN(10,FILE='area.grd',FORM='UNFORMATTED',STATUS='REPLACE',ACCESS='DIRECT',RECL=360*90)
OPEN(20,FILE='tot_area.grd',FORM='UNFORMATTED',STATUS='REPLACE',ACCESS='DIRECT',RECL=360*90)
OPEN(30,FILE='area.txt',FORM='FORMATTED')

P=1
DO K=1,NT
DO J=1,NY
DO I=1,NX
  READ(10,REC=P) FR(I,J,K)
  P=P+1
ENDDO
ENDDO
ENDDO

DO K=1,NT
  S(K)=0.0
  DO J=1,NY
    DO I=1,NX
        IF (FR(I,J,K)>0.0) THEN                                                                                       
      S(K)=S(K)+12351.77853*COS((J-0.5)*3.14159/180)*FR(I,J,K)
        ENDIF
        ENDDO
  ENDDO
ENDDO

P=1
DO K=1,NT
  WRITE(30)S(K)
  WRITE(20,REC=P)S(K)
  P=P+1
ENDDO

CLOSE(20)
CLOSE(30)
STOP
END


运行后fortran报错
run-time error F6414:open(area.grd) access not allowed

以为是area.grd中数据太大（360*90*1701），就把时次1701改成了10，但还是同样的错误。。
请教是为什么啊？

mofangbao

贴这么一大段代码没有用的，先用grads画图对比，确定你转出的文件是对的再继续下一步。
然后用逐步排除的方法自己先定位出fortran的错误是在哪里产生的，贴这么一大段出来这不是让人看得都眼花么

言深深

看看是不是fortran语句read的不对，我记得需要一直都用隐循环的

做个霸气的木头

打开个grd文件整的那么复杂，只要 open(10,file='area.grd',form='binary')就行了，你改改看是不是这个问题~

chaseechoes

把read改成隐循环后还是报同样的错误诶。。。

chaseechoes

做个霸气的木头 发表于 2012-7-26 16:46

                               
登录/注册后可看大图

打开个grd文件整的那么复杂，只要 open(10,file='area.grd',form='binary')就行了，你改改看是不是这个问题 ...

还是一样诶。。。

chaseechoes

言深深 发表于 2012-7-26 16:25

                               
登录/注册后可看大图

看看是不是fortran语句read的不对，我记得需要一直都用隐循环的

把read改成隐循环后还是报同样的错误诶。。。

做个霸气的木头

你nc文件的ctl，查询下（q ctlinfo），截个图上来看看，可能是精度问题在捣鬼~

chaseechoes

做个霸气的木头 发表于 2012-7-26 17:04

                               
登录/注册后可看大图

你nc文件的ctl，查询下（q ctlinfo），截个图上来看看，可能是精度问题在捣鬼~

这个吗

做个霸气的木头

chaseechoes 发表于 2012-7-26 17:16

                               
登录/注册后可看大图

这个吗

看来不是精度的问题，你这样吧，你一步步测试程序哪里出了问题
open部分改成我的，读数据的整个过程改成这样：
DO K=1,NT
   DO J=1,NY
      READ(10) (FR(I,J,K),I=1,NX)
      WRITE(*,*) FR(1,J,K)
   ENDDO
ENDDO
ENDDO
如果能在界面上输出数，而且能输完。那么就是这里问题了，如果在这部分出错了，你要试着把gs里的x、y格点也改成循环试试~