用Fwrite函数提取nc资料后，grd文件中数据排列是什么格式呢？

jinxinhao

我用的1948-2013年hgt月资料，选择500hpa，全部792个时次，分辨率2.5*2.5，我想用fortran继续处理，请问grads生成的grd文件中数组的形式是什么？

容与

'set gxout fwrite'
……………………
'set fwrite e:\***.grd'
'd hgt'
'disable fwrite'
;

维克多·蓝晶

open(10,file="e:\hgt.grd",form='binary')
  do m=1,t
     do j=1,y
        do i=1,x
            read(10) h
        end do
     end do
  end do

snowNJU

     open(12,file='hgt.grd',form='unformatted',
     +     status='old',access='direct',recl=m*n*4)
     iw=1
         do it=1,792
         read(12,rec=iw) ((hgt(i,j,it),i=1,m),j=1,n)
         iw=iw+1
         enddo

jinxinhao

维克多·蓝晶 发表于 2015-1-27 09:17
open(10,file="e:\hgt.grd",form='binary')
  do m=1,t
     do j=1,y

谢谢解答。。。

jinxinhao

snowNJU 发表于 2015-1-27 13:26
open(12,file='hgt.grd',form='unformatted',
     +     status='old',access='direct',recl=m*n*4) ...

谢谢解答，我试试。。

雨恬魅影

二进制格式的。。