Fortran读取二进制dat文件出错

ningm3

Fortran代码如下
      
      parameter(iix=300,jjx=300,TX=3)         
      parameter(ix=300, jx=300,kx=37)           
      parameter(nx=300, ny=300,nz=37,cf=1.e30)  

        iuni_r1=11
        iuni_r2=22
        ire_r1=1
        ire_r2=1

        open(iuni_r1,FILE='./qva.dat',
&STATUS='UNKNOWN',FORM='unformatted',
&ACCESS='DIRECT',RECL=nx*ny*4)

             open( iuni_r2,file='./out.dat',
&status='unknown',form='unformatted',
&access='direct',recl=nx*ny*4)

                        do l=1,tx

                               print*,'readtime=',l

                              do k=1,kx      
                                    read(iuni_r1,rec=ire_r1) ((qvapor(i,j,k),i=1,ix),j=1,jx)
                                    ire_r1=ire_r1+1
                                    enddo


                       print*,'qvapor=',qvapor(20,20,4)

                          do k=1,nz
                          write(iuni_r2,rec=ire_r2) ((qvapor(i,j,k),i=1,nx),j=1,ny)
                          ire_r2=ire_r2+1
                            enddo

        enddo

        END

其中的 qva.dat 文件是用ARWpost处理的wrfout中的QVAPOR，使用这个dat文件画出的水汽图是没有问题的。
但是读入再读出后得到的out.dat中的数据非常奇怪，量级能达到10^30，（在print输出的时候就发现数据很奇怪）。
查阅了很久的资料，也搜了很多帖子，但是都解决不了。
PS：form='binary' 和 recl=nx*ny 都尝试过，但是全部报错
请求各位大神解答，感激不尽

小小书童

open(iuni_r1,FILE='./qva.dat',FORM='binary')
open( iuni_r2,file='./out.dat',form='binary')
写入是这样写：
do k=1,nz
         write(iuni_r2) ((qvapor(i,j,k),i=1,nx),j=1,ny)
   enddo
试试吧

小小书童

open(iuni_r1,FILE='./qva.dat',FORM='binary')
open( iuni_r2,file='./out.dat',form='binary')
写入是这样写：
do k=1,nz
         write(iuni_r2) ((qvapor(i,j,k),i=1,nx),j=1,ny)
   enddo
试试吧

ningm3

您好~刚刚试过了，还是有问题，数据还是错的。
程序都能运行通，但是读出的数据有问题

shayuewancom

请问楼主这个问题解决了吗？我也遇到了，想请教

Snowyea

谢谢露珠分享~

ningm3

shayuewancom 发表于 2016-6-23 10:55
请问楼主这个问题解决了吗？我也遇到了，想请教

试试是不是ctl设置的缺省值问题

shayuewancom

ningm3 发表于 2016-6-24 06:58
试试是不是ctl设置的缺省值问题

非常感谢，我试试看！

吴阿好大人

请问楼主解决这个问题了没

qq349086723

{:eb513:}{:eb513:}{:eb513:}