站点资料处理无法生成map文件

baby1117

处理站点资料，运行stnmap.exe，输入如下图的命令，然后回车，然后就闪退，也没有生成map文件。用grads运行也是出错，请大家帮忙看一看~

integer,parameter:: sta_num=824
character*5:: fstation(824,31)
integer:: i,j,tim
integer:: fyear(824,31)
integer:: fmon(824,31),fday(824,31)
integer:: pre1(824,31),pre2(824,31),pre_total(824,31)
integer:: c_code1(824,31),c_code2(824,31),c_code3(824,31)
real lat(824),lon(824),hgt(824)
integer:: lon1(824,31),lat1(824,31),hgt1(824,31)
character*8:: id(824)

open(unit=1,file='H:\data00\daily meteorological data_1979_2013\SURF_CLI_CHN_MUL_DAY-PRE-13011-200810.TXT') !!precipitation data

open(unit=22,file='H:\data00\daily rain\sta_decimal.txt')
open(unit=23,file='H:\data00\daily rain\200810.dat',form='binary',status='unknown')
open(unit=24,file='H:\data00\daily rain\200810.txt',form='formatted',status='unknown')

read(1,'(A5,I5,I6,I7,I5,2I3,3I7,3I2)') ((fstation(i,j),lat1(i,j),lon1(i,j),hgt1(i,j),fyear(i,j),fmon(i,j),fday(i,j),pre1(i,j),pre2(i,j),pre_total(i,j),c_code1(i,j),c_code2(i,j),c_code3(i,j),j=1,31),i=1,824)
read(22,*)(id(i),lat(i),lon(i),hgt(i),i=1,824 )
do j=1,31,1
  tim=0
   lev=1
   nflag=1
   do i=1,824,1
      write(23)id(i),lat(i),lon(i),tim,lev,nflag,pre_total(i,j)
             write(24,*)id(i),lat(i),lon(i),tim,lev,nflag,pre_total(i,j)
   enddo
enddo

end

baby1117

现在又出现这个问题