fortran读micaps第四类数据

边缘人糊涂乐

fortran程序如下：
program main
integer,parameter::x=121,y=81,z=7,g=14
real r(x,y,z,g)
character(8)::day
character(4)::cen(z)
character(4)::h(g)
character(6)::typ(1)
data cen/'1000','925','850','700','500','200','100'/
data h/'.000','.006','.012','.018','.024','.030','.036','.042','.048','.054','.060','.066','.072','.078'/
data typ/'R'/  


  do k=1,z
   do l=1,g
        open(15,file='z:\ecmwfthin\'//trim(typ(1))//'\'//trim(cen(k))//'\'//trim(init_time)//h(l),status='old')
    read(15,*)
    read(15,*)
    read(15,*)
    do i=1,y
    read(15,*) (r(j,i,k,l),j=1,x)
    enddo
    close(15)
   enddo
  enddo

  open(12,file='D:\ectin\d\grd\r.grd',status='unknown',form='binary')
  do l=1,g
   do k=1,z
        do i=y,1,-1
    write(12) (r(j,i,k,l),j=1,x)
        !print*, r(1,2,1,1)
    enddo
   enddo
  enddo
  close(12)

end

这个程序在读数据时，有些可以读出来，有些读不出来，但是数据格式都是一样的，在读 z:\ecmwfthin\r\1000\18061120.000等数据时可以读取，但是在读 z:\ecmwfthin\r\850\18061120.000时就读不出来，程序报读取有问题，但是检查过数据，都是一样的，不知道问题出在哪里？
数据截图如下：

850数据

1000数据

边缘人糊涂乐

有人吗？{:eb303:}{:eb303:}{:eb303:}{:eb303:}{:eb303:}

我的明天

请问楼主解决了吗，我也遇到了一样的问题