多个文件合并的问题，有意思的探讨，图+代码+思路。

FortranCoder

如图，先解释下着十二个文件的意思，十二个BIL格式的文件代表十二个月的降雨数据，是由降雨量，经度，维度插值而成的。目的：把这十二个文件的数据合并到一个PXV格式的文件数据当中，没合并之前，换句话说，一个PXV数据里面，每一个坐标点，（月份，经度，维度）的降雨量是不同的。
思路：如果文件大小相同，那么直接把数据按照一定的步长（RECL）打开，然后把十二个二进制文件看做十二个二维数组读入到一个三维数组里面，然后在把三维数组读入到PXV文件里面。
问题：由于十二个文件的大小有两种形式，一种是798KB，一种是395KB，前者是后者的两倍，所以导致了读取方面的一些问题搞不懂了，RECL怎么设置，存入的三维数组的大小又如何设置？

求大神支招，附上自己的代码。

[Fortran] 纯文本查看 复制代码

 program main

         parameter (MROW=411,MCOL=491,MXPP=12,NY=1)        
         integer*1 grid1(MCOL,MROW,MXPP)
         integer*2 grid2(MCOL,MORW,MXPP)
         character*4 yn,dn
         integer MXP,status

       do in=1,NY
         write(yn,'(i4)') in+2007
       MXP=12
 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
                 do im=1,MXP
             write(dn,'(i4)') im+1000
             status=0
           open(1,file='.\bil\pre\pre'//yn//'_'//dn//'.bil',
      +form='binary',access='direct',recl=2*MCOL,iostat=status)
             do lat=1,MROW
              read(1,rec=lat) (grid2(lon,lat,im),lon=1,MCOL)
              IF(STATUS>0)THEN
                  GOTO 100
              ELSE
              GOTO 1000
              ENDIF
             enddo

 100       CLOSE(1)
          open(1,file='.\bil\pre\pre'//yn//'_'//dn//'.bil',
      + form='binary',access='direct',recl=MCOL,iostat=status)
                   do lat=1,MROW
               read(1,rec=lat) (grid1(lon,lat,im),lon=1,MCOL)
               enddo
 1000     Close(1)
           enddo

         open(19, file='.\pxv\pre\pre'//yn//'.pxv', form='binary',
      +access='direct', recl=MXP*2)
         irec = 0
                 do lat=1,MROW
                         if (mod(lat,100).eq.0) write(*,*) ' Row ...', lat
                         do lon=1,MCOL
                         irec = irec + 1
                         write(19,rec=irec) ((grid2(lon,lat,im)),im=1,MXP)
                         enddo
                 enddo
         close(19)
         enddo

         end 

FortranCoder

vvt 发表于 2015-4-22 10:51
这个需要你上传一定量的数据文件。

意思是这些二进制文件要上传上来“？

FortranCoder

求助啊 ，有没有能支招的

FortranCoder

vvt 发表于 2015-4-22 11:00
稍安勿躁，是的，需要上传文件。

已经上传了

FortranCoder

li913 发表于 2015-4-22 15:07
需要这两种文件的存储格式。

原始文件都是BIL格式的啊，只是大小不一样，一种是798KB，另外一种是他的1/2大小，要存储的格式为PXV格式的，就是一个文件合并的问题啊。