Fortran读入gs提取的grd数据文件

lijian

gs文件
'reinit'
'open /data/era-interim/era-40-interim.ctl'
'set fwrite test4.grd'
'set gxout fwrite'
i=18381
while(i<=18505)
j=1
while(j<=32)
'set t 'i''
'set z 'j''
'set lon 0 358.5'
'set lat -90 90'
'd Uprs'
j=j+1
endwhile
j=1
while(j<=32)
'set t 'i''
'set z 'j''
'set lon 0 358.5'
'set lat -90 90'
'd Vprs'
j=j+1
endwhile
j=1
while(j<=32)
'set t 'i''
'set z 'j''
'set lon 0 358.5'
'set lat -90 90'
'd Wprs'
j=j+1
endwhile
j=1
while(j<=32)
'set t 'i''
'set z 'j''
'set lon 0 358.5'
'set lat -90 90'
'd Zprs'
j=j+1
endwhile
i=i+1
endwhile
'disable fwrite'
;

Fortran读入数据，输出高度场
  program file_1
      program file_1
      integer,parameter::nx=240,ny=121,nz=32,nt=125
      integer i,j,k
      real u(nx,ny,nz),v(nx,ny,nz),w(nx,ny,nz),z(nx,ny,nz)
      open(12,file='test4.grd',form='unformatted',access='direct' &
            ,convert='big_endian',recl=nx*ny*4)
      open(12,file='hgt.grd',form='unformatted',access='direct' &
          ,recl=nx*ny*4)
     ! total=32*3+1=97
       irec=97
       irec1=1
       do it=1,nt
          read(12,rec=irec)(((z(i,j,k),i=1,nx),j=1,ny),k=1,nz)
          write(13,rec=irec1)(((z(i,j,k),i=1,nx),j=1,ny),k=1,nz)
          irec=irec+128
          irec1=irec+1
          t=t+1
       enddo
       close(12)
       close(13)
       end

是不是记录读错了，或是数组应该定义成4维？

river

谁告诉你是big_endian？

numa

彭国伦版的Fortran书中没有convert参数的介绍，能给讲讲么？