求助！多层次站点资料的读取

whmooncool

看了一些关于站点资料的读取的帖子，但是多个层次的读取方法基本没有涉及到，自己试着按照家园整理的《站点作图详解》写了程序，但是仍然有问题，求指点！！！
station资料的读取：
program station
      parameter (mt=427,nz=4)    ! mt 站点数量 ，nz层次
      character stid(mt)*8
      integer flag,nlev
      real lat(mt),lon(mt),tim,high(mt),number(mt),lev(nz)
     real height(mt,nz),temp(mt,nz),ttd(mt,nz),ud(mt,nz),us(mt,nz)

      character stationfile(4)*8
   data stationfile/'850','700','500','200'/
  
   do k=1,nz
   open(k,file=''//trim(stationfile(k))//'_20.000')
   read(k,*)
   read(k,*)
   do i=1,mt
   
      read(k,*)stid(i),lon(i),lat(i),high(i),number(i),height(i,k),temp(i,k),ttd(i,k),ud(i,k),us(i,k)
   end do   

end do  !结束k的循环

      tim=0.0
      flag=1
      nlev=4   
   open(40,file='20.grd',form='binary')  
   do i=1,mt
  write(40)stid(i),lat(i),lon(i),tim,nlev,flag,high(i)
  
  do k=1,nz
  lev(k)=k
  write(40)lev(k),height(i,k),temp(i,k)
  write(*,*)lev(k),height(i,k)  
  end do
end do
  
nlev=0
    write(40)stid(mt-1),lat(mt-1),lon(mt-1),tim,nlev,flag
close(40)
end
显示运行无误
ctl描述文件如下：
dset f:\20.grd
dtype station
stnmap f:\20.map
undef 9999.0
title 500hpa 20
zdef 4 levels 850 700 500 200
tdef 1  linear  29jul2012 1DY
vars 3
high   0    99   high
h      4    99      height
t      4    99      temp
endvars
随后在grads中输入！stnmap -i f:/20.ctl,运行出map文件;
最后用grads画出站点资料，只有high显示有数据，而变量h和t显示空白。
请问问题出在哪？请大家指点下，谢谢~~

vera

write 语句没有将3个变量都输出，只输出了2个变量