对站点进行格点转换，但是出不了图，全是缺测，为什么

刘晶

parameter(m=66)
real vec(m)
     open(1,file='e:\sx\pr_ens4_evf.grd',form='binary')
     read(1)(vec(i),i=1,m)
     close(1)   
  call stntogrd(vec)
  end
   subroutine stntogrd(h)
   parameter (n=66)
         real lat(N),lon(N),h(n)
         character*8 stid(N)
   open(20,file='e:\sx\lon_lat-deal.txt')
   
!读经纬度
do  k=1,n
      read(20,100)lat(k),lon(k)
   print*,lat(k),lon(k)
enddo
100 format(f5.2,2x,f6.2)
      close(20)
!站号赋值
       do 2 i=1,n
2    stid(i)=char(i)


!cccc  Note: recl=9*4 means: stid*8: 2*4,  rlat,rlon tim,nlev,nflag,x,y:7*4
!cccc   so total recl=9*4
OPEN(9,FILE='e:\sx\66winter.grd',form='binary')
OPEN(10,FILE='e:\sx\66winter.txt')
       TIM=0.0
       NLEV=1
       NFLAG=1
       IREC=1
       DO 40 I=1,n
       WRITE(9) STID(I),LAT(I),LON(I),TIM,NLEV,NFLAG,h(i)
WRITE(10,*) STID(I),LAT(I),LON(I),TIM,NLEV,NFLAG,h(i)
40     continue
!cccc   On end of file write last time group terminator.
       NLEV = 0
       WRITE(9) STID(I-1),LAT(I-1),LON(I-1),TIM,NLEV,NFLAG
    WRITE(10,*) STID(I-1),LAT(I-1),LON(I-1),TIM,NLEV,NFLAG
  Close(9)
  close(10)
       return
end

一闭一睁

飘过，我也在学习啊！