求助！关于grads站点资料作图！

realzly

在论坛里找到了清风的关于将站点资料用grads作图的帖子，仔细研究了一番。现在出现以下问题：
在用grads创建.map文件时出现了错误，提示：
ga-> ! stnmap -i e:/rain/sk/rain.ctl
  Name of binary data set: e:\rain\sk\rain1.grd
  Number of times in the data set: 1
  Number of surface variables: 1
  Number of level dependent variables: 0

Starting scan of station data binary file.
Binary data file open: e:\rain\sk\rain1.grd
Processing time = 1
  Invalid station hdr found in station binary file
  Possible causes:  Invalid level count in hdr
                    Descriptor file mismatch
                    File not station data
                    Invalid relative time
    levs = 892398276  flag = 538976304  time = 2.57843e-09
后又用清风给的数据做了一下，仍然是这样的错误。
万望做出来的大侠们给予指点！

abd

清风给的数据也是这个数据啊？这个有点奇怪了。是MICAPS第一类数据吗？你是用你自己编的fortran程序转换的吗？还是又一个exe文件转换的呀？站点数据确定没问题的话，那是不是stnmap.exe有问题啊？你换一个试试吧！

mofangbao

你转换的fortran程序贴出来看一下，应该是转换的地方出错了，我给的例子有一些地方有点笔误，看的时候要加上自己的思考

wode_q_x

最好把转换的fortran程序贴一下，这样才可以找到问题所在。

yinlingr

进来学习，标记一下，以后碰到问题也有个参考。

realzly

这是fortran程序：就是清风给的程序

program skjs
      implicit none
        real lon,lat,rain,var,tim
        integer nlev,flag
        character*8 stid(160)
      tim=0.0
        nlev=1
        flag=1
        open(1,file='E:\rain\sk\raindata.txt',status='old')   
      open(2,file='E:\rain\sk\raindata.grd',status='replace'
     &,form='binary')
        
        read(1,*)
        read(1,*)
10    read(1,*,end=100)stid,lon,lat,var,rain
        write(2)stid,lat,lon,tim,nlev,flag,rain
        goto 10
100   continue
      close(1)
        nlev=0
        write(2)stid,lat,lon,tim,nlev,flag
        close(2)
        end

这是ctl:
DSET e:/rain/sk/raindata.grd
DTYPE station
stnmap e:/rain/sk/raindata.map
UNDEF -32766
TDEF 1 linear 00z27JUL2012 12hr
VARS 1
RAIN 0 99 RAINFALL DATA
ENDVARS

望大侠们给予指点！
        

姚小娟

偶试了一下多个时次的站点数据，也出现同样的错误，请教大虾了！！！

1417920146

UNDEF -32766 ,这个是不是要改数值，好多实况缺测用999或9999

TANG

求助啊，遇到同样问题

lightxx

realzly 发表于 2012-8-22 18:33
这是fortran程序：就是清风给的程序

program skjs

你的问题解决了没哦 现在遇到一样的问题了 解决的话可以的话交流下哦