用micaps画多时次K指数场，生成.map出问题了请大神看看呢~~

311wen

用micaps提取了全国站4个时次的K指数


parameter(nt=4)                                       !!4个时次
integer::nlev,nflag
real::lon(120),lat(120),k(120),var(120)
character*8 sta(120)                                 !!120个全国站点，固定值
character*14 filename(nt)                          !!
open(1,file='k_4.txt')
do i=1,nt
read(1,*) filename(i)
print*,'Filename:',filename(i)
pause
enddo
close(1)               

open(2,file='k_4.dat',form='binary')             !读入数据
do j=1,nt
      open(3,file=filename(j))
      do i=1,3
      read(3,*)
      end do  
      do i=1,120                                             !  从第四行开始读每个站点的经纬高和K指数值   
      read(3,*) sta(i),lon(i),lat(i),var(i),k(i)
      if (k(i)==9999) then
      continue
      endif
      tim=0.0;nlev=1;nflag=1
      write(2)  sta(i),lat(i),lon(i),tim,nlev,nflag,k(i)
      enddo
     close(3)
      nlev=0
      write(2) sta(120),lat(120),lon(120),tim,nlev,nflag
        !一个时次的k场输入完毕

enddo
close(2)
end


--------------------------------------------------------------------------生成k_4.dat以后----------------------------------写的站点ctl文件以及！stnmap -i时的错误



   



但是我如果将ctl写为一个时次的如下，，就可以生成一个时次的.map,但后面时次的图画不出来。
dset    c:/msdev/projects/k_4/k_4.dat
dtype   STATION
stnmap  c:/msdev/projects/k_4/k_4.map
undef   9999
title   k ONSET DATE IN QTP
tdef   1 linear 14jun2007 12hr
vars    1
k      0  99 k date
endvars

谢谢大家指导！

311wen

啊实在没辙了。想不出来·~~麻烦大神看一下

311wen

啊啊没人么

311wen

up

造字仓颉

非常感谢楼主，希望好好学习一下，提高一下技术。

丫丫君

我也遇到同样的问题