grads作图数据文件转换的时候，fortran程序里的recl为什么设成nx*ny*4？

超级大菜鸟

十分不理解这个4是怎么回事？？？哪位能帮忙解答一些，谢谢了！！！！

Aires

不用设置*4吧

...该...

问问版主吧

鲸鱼便子

http://bbs.06climate.com/forum.php?mod=viewthread&tid=7773
这个教材的234页有具体说明。
希望对你有帮助~~~~~

超级大菜鸟

例子：现有ASCII码（十进制存放格式）数据资料文件u.dat、v.dat和SST.dat，其空间范围：60-150°E，0-40°N；层次：u、v为850、 200hPa；时段：1982.1-1985.12；分辨率：2.5*2.5。要求编写出将这三个文件转换成二进制无格式直接存取（Grads格式.grd）文件的Fortran程序。
程序：c   打开原始数据文件
open(1,file='u.dat')
open(2,file='v.dat')
open(3,file='sst.dat')
c   打开目标文件，recl为记录的长度，对于直接存取文件必须指定记录长度 。
open(12,file=‘mhy.grd',form='unformatted',
#         access='direct', recl=nx*ny*4)
c 将数据资料写入新文件中， irec重新进行记录累
irec=0      
do 200 it=1,nt
do 30 iz=1,nz
    irec=irec+1
   write(12,rec=irec) ((u(i,j,iz,it),i=1,nx),j=1,ny)
30 Continue
  do 31 iz=1,nz
     irec=irec+1
   write(12,rec=irec) ((v(i,j,iz,it),i=1,nx),j=1,ny)
31 continue
    irec=irec+1
   write(12,rec=irec) ((sst(i,j,it),i=1,nx),j=1,ny)
200   continue
end
有的需要，有的不需要，但是不知道为什么？？？？

超级大菜鸟

鲸鱼便子 发表于 2013-4-3 08:33

                               
登录/注册后可看大图

http://bbs.06climate.com/forum.php?mod=viewthread&tid=7773
这个教材的234页有具体说明。
希望对你有帮 ...

按教材的理论，只需设置成nx*ny 就行。。

东风急流

一般要看你用的编译器，如果是intel的，则不用乘以4。如果是在intel编译器中调通了，换到pgi编译器中则要乘以4，否则会提示超出范围。