生成地图映射文件出错

schLltt

！！！！！！！！！！！！！！！！！！！！！！！！fortran程序！！！！！！！！！！！！！！
program main
   implicit none
   integer,parameter:: n=744
   character*8 stid
   character*12 filename,filena
   real lon,lat,win,dir,tim
   integer nlev,flag,error,i
   character*12 a(n)   
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  tim=0.0
   nlev=1
   flag=1
open (100,file='b1308.txt')
do i=1,n
   read(100,*) a(i)
end do
close(100)
  

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
do i=1,n
   write(filename,"(a12)")a(i)
   write(filena,"(a4)")a(i)

   open (101,file='G:\data\grads_station\1308\'//filename,status='old',iostat=error)
   
   if(error/=0)then
   write(*,*) 'open file failed.'
   end if
   open (102,file='G:\data\grads_station\1308\'//trim(adjustl(filena))//'.dat',access='append',form='binary')
   10 read(101,*,end=100)stid,lon,lat,dir,win              
  
   write(102)stid,lat,lon,tim,nlev,flag,dir,win
   goto 10
   100 continue
   close(101)
   nlev=0
   write(102)stid,lat,lon,tim,nlev,flag
   close(102)
  
   end do
  
  
   end
   
！！！！！！！！！！！！！！！！ctl文件！！！！！！！！！！！！！！！！！！
   dset ^1308.dat
   dtype station
   stnmap ^station.map
   title  wind stream
   undef 9999
   options template
   tdef 744 linear 00Z01Aug2013 1hr
   vars 2
   dir 0 99  surface wind-direction
   win 0 99  surface wind
   ENDVARS


！！！！！！！！！！！！！！！！！！Grads错误提示如图！！！！！！！！！！！！！
请问是哪里出现问题了？多谢

请问附件2里的两个文件到底有什么区别？为什么第二个（用的时候文件名都为13080101.dat）生成地图映射文件没问题，第一个就有问题？

river

建议还是先去搜索一下站点资料作图的帖子，学习一下吧

schLltt

river 发表于 2015-4-27 18:27
建议还是先去搜索一下站点资料作图的帖子，学习一下吧

方法应该是没错，可能问题还是出在程序那儿，提示time 1 stn  count= 0，没读进去

river

schLltt 发表于 2015-4-27 22:38
方法应该是没错，可能问题还是出在程序那儿，提示time 1 stn  count= 0，没读进去

不是可能，是一定。

schLltt

river 发表于 2015-4-28 12:15
不是可能，是一定。

您好，我把append去掉了，fortran生成的dat文件第一个时次正确了， 可是从第二个时次开始，生成地图映射文件时Grads总是提示has stn count=0，纠结了一天了。一个小时一个小时的读写没问题，fortran一循环就出现问题，能不能帮忙看下

schLltt

river 发表于 2015-4-28 12:15
不是可能，是一定。

每次运行都是第一个时次的数可以用，之后的就提示has stn count=0~~~~

river

schLltt 发表于 2015-4-28 17:07
每次运行都是第一个时次的数可以用，之后的就提示has stn count=0~~~~

我不是都说然你先搜索一下论坛看看站点作图的帖子吗，难道这也要别人帮你搜索好了么。
http://bbs.06climate.com/forum.php?mod=viewthread&tid=12626
http://bbs.06climate.com/forum.php?mod=viewthread&tid=10006

schLltt

river 发表于 2015-4-28 20:16
我不是都说然你先搜索一下论坛看看站点作图的帖子吗，难道这也要别人帮你搜索好了么。
http://bbs.06cli ...

解决了。。因为从第二次循环开始，nlev就被改为0了。。。。。打扰了。。。fortran学的不好。。。