希望高手指点，我这程序到底哪里出错了，出来的图怎么是这么的？

wxxxren

我想用这个插值程序把高斯格点资料（经度是等间距1.875，纬度是不等间距的）插值成0.5*0.5的均匀网格数据。只对中国区域进行插值，并且只是一个时次的。出来的结果是图中有很多空白区域，怎么回事呢？求指点啊！大家帮忙看看，非常感谢！


program inter
      implicit none
        integer,parameter::ex=37,ey=25,et=1,gx=140,gy=90  !!  Input values  ---------  =*=  ----------
  real,parameter::undef=-9.96921e+36
  
        real::glat(gx,gy),glon(gx,gy),grid(gx,gy),elat(ex,ey),elon(ex,ey),ease(ex,ey),lat(ey),f(ex,ey),p(gx,gy)
     
real::t,u,it,m,n

     integer::i,j,k,ii,jj,index,irec

!----*-----------------------------------------------------6---------7--
     open(10,file='d:\mm\mmm.grd', form='unformatted',access='direct',recl=ex*ey)
    read(10,rec=1) ((f(i,j),i=1,ex),j=1,ey)
    close(10)

      open(20,file='d:\mm\elat.txt',form='formatted')
     read(20,*) (lat(j),j=1,ey)
  close(20)
!pause
!----*-----------------------------------------------------6---------7--
      do j=1,ey
      do i=1,ex
  elat(i,j)=lat(j)
  elon(i,j)=71.25+1.875*(i-1)
enddo
  enddo

   do n=1,gy
      do m=1,gx
        glat(m,n)=15.25+0.5*real(n-1)
  glon(m,n)=70.25+0.5*real(m-1)
       enddo
  enddo
     

    do j=1,ey
        do i=1,ex
  if(f(i,j)==9.96921e+36)then
      ease(i,j)=0.0
    else
      ease(i,j)=f(i,j)
   end if
   
  enddo
  enddo
        call bilinear(elat,elon,ex,ey,ease,glat,glon,gx,gy,grid) ! VIP=Very Important Part
  write(*,*)

do n=1,gy
  do m=1,gx
    p(m,n)=grid(m,n)
  
     enddo
    enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     open(2,file='d:\mm\ccout.txt',form='formatted')
     ! do k=1
     write(2,*) ((p(i,j),i=1,gx),j=1,gy)
  close(2)
!----*-----------------------------------------------------6---------7--
      open(30,file='d:\mm\cm.grd',form='binary')
   
     write(30) ((p(i,j),i=1,gx),j=1,gy)
  close(30)
!----*-----------------------------------------------------6---------7--
      stop
   end
!----*-----------------------------------------------------6---------7--
!     --- convert from EASE to GRID using bilinear interpolation ---
!     Input   ex, ey, gx, gy, ease, elat, elon, glat, glon
!        ex
!        ey
!        gx
!        gy
!        elat(ex,ey)
!        elon(ex,ey)
!        elat(ex,ey)
!        glat(gx,gy)
!        glon(gx,gy)
!     Output  grid
!        grid(gx,gy)
!----*-----------------------------------------------------6---------7--
     subroutine bilinear(elat,elon,ex,ey,ease,glat,glon,gx,gy,grid)
!     --- convert from EASE to grid using bilinear interpolation ---
      integer::ex,ey,gx,gy
      dimension elat(ex,ey),elon(ex,ey),ease(ex,ey)
      dimension glat(gx,gy),glon(gx,gy),grid(gx,gy)
      
      do 141 i=1,gx
      do 142 j=1,gy
      do 143 ii=1,ex-1
      do 143 jj=1,ey-1
        if(glat(i,j).gt.elat(ii,jj) .and. glat(i,j).lt.elat(ii+1,jj+1) &   ! i/=ii,j/=jj
     &   .and. glon(i,j).gt.elon(ii,jj+1) .and. glon(i,j).lt.elon(ii+1,jj))then
        i0=ii
        j0=jj
        t=(glon(i,j)-elon(ii,jj))/(elon(ii+1,jj)-elon(ii,jj))
        u=(glat(i,j)-elat(ii,jj))/(elat(ii,jj+1)-elat(ii,jj))
        index=1
          if(index.eq.1)then
          grid(i,j)=(1.-t)*(1.-u)*ease(i0,j0)+ &
     &          (1.-t)*u*ease(i0,j0+1)+   &
     &          t*(1.-u)*ease(i0+1,j0)+   &
     &          t*u*ease(i0+1,j0+1)
          endif
        elseif(glat(i,j).eq.elat(ii,jj) &                                  !  i==ii,j/=jj
     &   .and. glon(i,j).gt.elon(ii,jj) .and. glon(i,j).lt.elon(ii,jj+1))then
     i0=ii
  j0=jj
  t=0.
  u=(glat(i,j)-elat(ii,jj))/(elat(ii,jj+1)-elat(ii,jj))
        index=2
    if(index.eq.2)then
          grid(i,j)=(1.-t)*(1.-u)*ease(i0,j0)+ &
     &          (1.-t)*u*ease(i0,j0+1)
          endif
        elseif(glon(i,j).eq.elon(ii,jj) &                                   !  i/=ii,j==jj
     &   .and. glat(i,j).gt.elat(ii,jj) .and. glat(i,j).lt.elat(ii+1,jj) )then
        i0=ii
        j0=jj
        t=(glon(i,j)-elon(ii,jj))/(elon(ii+1,jj)-elon(ii,jj))
        u=0.
        index=3
          if(index.eq.3)then
          grid(i,j)=(1.-t)*(1.-u)*ease(i0,j0)+ &
     &          t*(1.-u)*ease(i0+1,j0)
          endif
        elseif(glat(i,j).eq.elat(ii,jj)  &                                     !  i==ii,j==jj
     &   .and. glon(i,j).eq.elon(ii,jj))then
        i0=ii
        j0=jj
        t=0.
        u=0.
        index=4
          if(index.eq.4)then
          grid(i,j)=(1-t)*(1-u)*ease(i0,j0)
          endif
! elseif
    grid(i,j)=0.
endif
143   continue
142  continue
141  continue
      return
   end

wxxxren

C:\Documents and Settings\Administrator\桌面

wxxxren

C:\Documents and Settings\Administrator\桌面

mofangbao

你先把图片传到相册再贴过来吧，还有这个fortran太长，你还是先定位到出错的地方再来重新贴一下，要不然我们没有数据，直接看这个很崩溃的

wxxxren

mofangbao 发表于 2011-11-29 11:30

                               
登录/注册后可看大图

你先把图片传到相册再贴过来吧，还有这个fortran太长，你还是先定位到出错的地方再来重新贴一下，要不然我们 ...

恩，好的，嘿！刚才图片都没传明白！

topmad

外星人 你的图片挂了
顺便加点注释在程序里面 也好方便阅读啊

兰溪之水

又是一个破程序。。。

wxxxren

像兰溪说的这个破程序要让人疯了。下午发现判断条件写的还不够全！还得想想怎么加条件不能整成这样让大家帮忙改啊！感谢大家的积极帮忙！

wode_q_x

我也吸取教训，知道以后怎么发帖了，哈哈