求助！grads站点资料出图有问题

whmooncool

diamond 2 11年07月29日08时500百帕高空观测
12 07 29 08  500 428
  1152   14.37   67.27   13    1  560  -15  1.6  210   17
  1241    9.60   63.70    7    1  557  -16   27  205   20
  1415    5.63   58.88    9    1  551  -19   29  200   30
  2836   26.65   67.37  179    1  567  -15  2.1  255   17
  2963   23.50   60.82  103    1  576  -12  1.9  210   16
  3005   -1.17   60.13   84    1  544  -25  4.3  170    9
  3354   -1.24   53.00  117    1  552  -24    7  240   10
  3808   -5.31   50.22   88    1  557  -22   12  275   16
  3953  -10.24   51.93   14    1  555  -23   15  280   11
  4018  -22.59   63.97   38    1  568  -15   .4    5   22
  4220  -52.74   68.70   47    1  559  -17   15  210    9
  4270  -45.41   61.18   26    1  560  -17   10  165   13
  4320  -18.66   76.77   12    1  560  -12  2.0  285   12
使用数据如上，缺测值为9999
program station
      parameter (mt=428)    ! mt means nums of station
      character stid(mt)*8
      integer nflag,nlev
      real lat(mt),lon(mt),tim,td(mt)
   integer high(mt),number(mt),height(mt),temp(mt),ud(mt),us(mt)
      open(30,file='2908.txt')     
read(30,*)
read(30,*)
do i=1,mt
      read(30,*)stid(i),lon(i),lat(i),high(i),number(i),height(i),temp(i),td(i),ud(i),us(i)
end do   
close (30)
      tim=0.0
      nflag=1
      nlev=1
open(40,file='2908.dat',form='binary')
      do i=1,mt
      write(40)stid(i),lat(i),lon(i),tim,nlev,nflag,height(i)
   write(*,*)stid(i),lon(i),lat(i),tim,nlev,nflag,height(i)
      end do
       nlev=0
      write(40)stid(mt-1),lat(mt-1),lon(mt-1),tim,nlev,nflag
close(40)
end
读取的程序如上，运行显示无错误，图片1 .
dset f:\2908.dat
dtype station
stnmap f:\height.map
undef 9999
title 500hpa height
tdef 1  linear  29jul2012 1DY
vars 1
h      0   99   height
endvars
根据上面ctl运行出来的height.map也貌似没有错误
然后编写了gs文件
'reinit'
'open f:\height.ctl'
'set lat 20 60'
'set lon 90 140'
'd h'
'printim f:\height.png white'
；
最终出图如下
感觉数据就没有读进去，请问原因出在哪？求指点~~

做个霸气的木头

有没有人提醒过你，这里的   integer high(mt),number(mt),height(mt),temp(mt),ud(mt),us(mt)数据类型必须是实型数，我记得我以前出过这样的问题，你试试看先。

爱的侍者

en 改下实数型试一下~ 虽然整型可以读

river

你的gs很神奇啊，不是说站资料画图吗？我怎么连打开地图背景文件的命令都没看见，插值命令你也没使用啊

whmooncool

river 发表于 2013-1-4 16:04

                               
登录/注册后可看大图

你的gs很神奇啊，不是说站资料画图吗？我怎么连打开地图背景文件的命令都没看见，插值命令你也没使用啊

就是想简单说明下，插值命令不是在转换成格点数据才需要嘛》？

whmooncool

做个霸气的木头 发表于 2013-1-4 12:38

                               
登录/注册后可看大图

有没有人提醒过你，这里的   integer high(mt),number(mt),height(mt),temp(mt),ud(mt),us(mt)数据类型必须 ...

谢谢，换成实行的果然就写的这行就正确了，但是nlev和nflag还是不对啊，怎么不是1能？ 1152       14.37000       67.27000      0.0000000E+00  1.4012985E-45
  1.4012985E-45   560.0000   
1241       9.600000       63.70000      0.0000000E+00  1.4012985E-45
  1.4012985E-45   557.0000   
1415       5.630000       58.88000      0.0000000E+00  1.4012985E-45
  1.4012985E-45   551.0000   
2836       26.65000       67.37000      0.0000000E+00  1.4012985E-45
  1.4012985E-45   567.0000   
2963       23.50000       60.82000      0.0000000E+00  1.4012985E-45
  1.4012985E-45   576.0000   
3005      -1.170000       60.13000      0.0000000E+00  1.4012985E-45
  1.4012985E-45   544.0000   
3354      -1.240000       53.00000      0.0000000E+00  1.4012985E-45
  1.4012985E-45   552.0000   

whmooncool

爱的侍者 发表于 2013-1-4 14:14

                               
登录/注册后可看大图

en 改下实数型试一下~ 虽然整型可以读

谢谢，换成实行的果然就写的这行就正确了，但是nlev和nflag还是不对啊，不是1，居然是1.4012985E-45
，请问这是怎么回事？ 1152       14.37000       67.27000      0.0000000E+00  1.4012985E-45
  1.4012985E-45   560.0000   
1241       9.600000       63.70000      0.0000000E+00  1.4012985E-45
  1.4012985E-45   557.0000   
1415       5.630000       58.88000      0.0000000E+00  1.4012985E-45
  1.4012985E-45   551.0000   
2836       26.65000       67.37000      0.0000000E+00  1.4012985E-45
  1.4012985E-45   567.0000   
2963       23.50000       60.82000      0.0000000E+00  1.4012985E-45
  1.4012985E-45   576.0000   
3005      -1.170000       60.13000      0.0000000E+00  1.4012985E-45
  1.4012985E-45   544.0000   
3354      -1.240000       53.00000      0.0000000E+00  1.4012985E-45
  1.4012985E-45   552.0000   

爱的侍者

whmooncool 发表于 2013-1-4 16:26

                               
登录/注册后可看大图

谢谢，换成实行的果然就写的这行就正确了，但是nlev和nflag还是不对啊，不是1，居然是1.4012985E-45
，请 ...

       nlev=0
      write(40)stid(mt-1),lat(mt-1),lon(mt-1),tim,nlev,nflag
你这个是啥意思
不是只有一层数据吗？？

做个霸气的木头

whmooncool 发表于 2013-1-4 16:24

                               
登录/注册后可看大图

谢谢，换成实行的果然就写的这行就正确了，但是nlev和nflag还是不对啊，怎么不是1能？ 1152       14.370 ...

应该是小问题，自己检查下吧

river

whmooncool 发表于 2013-1-4 16:22

                               
登录/注册后可看大图

就是想简单说明下，插值命令不是在转换成格点数据才需要嘛》？

那你最后出图不就是要插值成个点的么