Fortran输出grd文件出错

tttccc

用高原季风强年平均降水量画降水实况图，站点转格点的过程中，用Fortran4.0写了程序，但是输出的不是grd文件，而是一个exe文件。大家可以帮我看看吗？多谢！
program year_dnz
implicit none
!
character(len=50), parameter ::  path='E:\lunwen\data\'
integer :: num1, i, k, end_num
integer, allocatable, dimension(:) :: stati
real, allocatable, dimension(:) :: lat, lon
real, allocatable, dimension(:) :: qn
character(len=8) :: id
real::tim
integer::nlev,nflag

!!!!! read stationdata !!!!!
open(1, file=trim(path)//'qnpj.txt', status='replace')
  read(1,*)   
  num1 = 0     ! 记录数据行数
  LOOP1 : do   ! 数数据的行数
    read(1,*,iostat=end_num)
    if( end_num == 0  ) then
       num1 = num1 + 1
    else
       exit LOOP1
    end if
  end do LOOP1
  !  为定义的各个变量（可调数组）分配空间
  allocate(stati(num1), lat(num1), lon(num1), qn(num1))
  rewind(1)    ! 回到数据第一行
  read(1,*)
  do k = 1, num1         
    read(1,*) stati(k), lat(k), lon(k), qn(k)
  end do
close(1)

!!!!! write txt data !!!!!

open(11,file=trim(path)//'qnpj.grd',status='replace',form='binary')
                         ! grd后缀，form=binary !
tim = 0.0
nlev = 1
nflag = 1
do i = 1,num1
  write(id,'(i8)') stati(i)  ! 站号由整形转化为字符型，长度为8
  write(11) id,lat(i),lon(i),tim,nlev,nflag,qn(i)
end do
nlev = 0
write(11) id,lat(num1),lon(num1),tim,nlev,nflag  ! nlev=0该时次结束
close(11)


deallocate(stati,lat,lon,qn)
end program year_dn