站点资料stnmap出问题，求指点

youxia1873

想画站点降水图，可是Grads 的stnmap命令貌似错误。很费劲
如图[img=0,1]file:///C:\Users\Administrator\AppData\Roaming\Tencent\Users\1123506777\QQ\WinTemp\RichOle\QBVYABC[RL}9U$XSC6X)%XS.jpg[/img]所示


另外，rain.map出来以后，只有一行 “ GrADS_stnmapV002   
   
      ? ”
gs文件为

1. 'reinit''open F:\station\grid.ctl''open F:\station\140824r.ctl''set grads off''enable print f:\station\rain.gmf''set lon 73 135.5''set lat 17 55''set mpdset cn cnriver''define a=oacres(g,p.2,1.5)''define a1 = maskout(a,g-0.5)''define aa = smth9(a1)''set xlopts 1 10 0. 18''set ylopts 1 10 0. 18''set gxout shaded''set ccolor rainbow''d aa''cbarn.gs''set cthick 8''set clopts 1 6 0. 1''set gxout contour''set cint 50''d aa'pull dummy'print''disable print';

'reinit''open F:\station\grid.ctl''open F:\station\140824r.ctl''set grads off''enable print f:\station\rain.gmf''set lon 73 135.5''set lat 17 55''set mpdset cn cnriver''define a=oacres(g,p.2,1.5)''define a1 = maskout(a,g-0.5)''define aa = smth9(a1)''set xlopts 1 10 0. 18''set ylopts 1 10 0. 18''set gxout shaded''set ccolor rainbow''d aa''cbarn.gs''set cthick 8''set clopts 1 6 0. 1''set gxout contour''set cint 50''d aa'pull dummy'print''disable print';

两个ctl文件为

1. dset  F:\station\140824r.grd
   
2. dtype station
   
3. stnmap F:\station\rain.map
   
4. undef -999.0
   
5. title rain
   
6. tdef 1 linear Aug2014 1mo
   
7. vars 1
   
8. p 0 99 rainf

复制代码

1. dset  F:\station\grid.grd
   
2. undef -999.0
   
3. title GRID Data
   
4. xdef 71 linear 70 1
   
5. ydef 41 linear 15 1
   
6. zdef 1 linear 1000 1
   
7. tdef 1 linear Aug2014 1mo
   
8. vars 1
   
9. g 0 99 rainfall data
   
10. ENDVARS
    

复制代码

最终图像无个点数据，应该是rain.map出问题，不知各位大神都是用的什么版本GRADS，可否传我一个

river

应楼主注意区别error和warning。这个不用翻译了吧，根本就不是一个意思。你最后出来图没有，看看图是不是对的

传说中的谁

river 发表于 2014-9-24 07:49
应楼主注意区别error和waring。这个不用翻译了吧，根本就不是一个意思。你最后出来图没有，看看图是不是对 ...

就是，你看电影的时候“FBI WARNING”也不影响观看不是么

mofangbao

楼主从现象到解决方案都给出了，还要提问干嘛呢？是不是直接说谁给我个GrADS的新版安装包就可以了？
如果您不是想要这样，那就直接给出你最后出来的图形，和你预期想见到的图形之间有何差别，不要自己乱猜测，你用记事本打开一个电影文件看看你能看懂多少，map文件你能看懂那么多字不等于只有那么多内容，stnmap提示一个warning也不等于它就错了。
最后，请自己找一下我之前发过的 GrADS站点作图详解，你想问的里面都有答案。

river

传说中的谁 发表于 2014-9-24 08:45
就是，你看电影的时候“FBI WARNING”也不影响观看不是么

哈哈，大家都起来的挺早啊

youxia1873

mofangbao 发表于 2014-9-24 08:56
楼主从现象到解决方案都给出了，还要提问干嘛呢？是不是直接说谁给我个GrADS的新版安装包就可以了？
如果 ...

大神，我执行GS文件，出错误：Warning from OACRES: less than two stations  Grid will be all missing values

youxia1873

mofangbao 发表于 2014-9-24 08:56
楼主从现象到解决方案都给出了，还要提问干嘛呢？是不是直接说谁给我个GrADS的新版安装包就可以了？
如果 ...

我的站点处理程序如下
implicit none
integer n
parameter(n=1532)
integer i,j,k,irec,nlev,nflag
real tim
real lon(n),lat(n),rain(n)
integer high(n)
character*8 stid(n)
open(1,file='F:/station/14082408.000')

do i=1,n
   read(1,*) stid(i),lat(i),lon(i),high(i),rain(i)
enddo
close(1)
!pause

open(2,file='F:/station/140824r.grd',form='binary')
tim = 0.0
nlev = 1
nflag = 1
do i = 1,n
   write(2) stid(i),lat(i),lon(i),tim,nlev,nflag,rain(i)
!   write(2) stid(i),lat(i),lon(i),tim,nlev,nflag,high(i)
enddo
nlev = 0
write(2) stid(n-1),lat(n-1),lon(n-1),tim,nlev,nflag
close(2)

open(12,file='F:/station/140824r.ctl')
     write(12,'(a)') 'dset  F:\station\140824r.grd'
     write(12,'(a)') 'dtype station'
     write(12,'(a)') 'stnmap F:\station\rain.map'
     write(12,'(a)') 'undef -999.0'
     write(12,'(a)') 'title rain'
     write(12,'(a)') 'tdef 1 linear Aug2014 1mo'
     write(12,'(a)') 'vars 1'
     write(12,'(a)') 'p 0 99 rainfall data'
     write(12,'(a)') 'ENDVARS'
close(12)

end

youxia1873

传说中的谁 发表于 2014-9-24 08:45
就是，你看电影的时候“FBI WARNING”也不影响观看不是么

资料片段如下：

58128  117.31   33.30   24    8.3
55472   88.63   30.95 4671    2.4
53522  107.85   40.05 1185    4.9
57313  106.76   31.85  360    0.4
57646  109.76   28.56  302   15.2
57566  111.65   29.45   51   12.4
53965  111.56   35.51  762    1.4
57060  111.23   35.16  404   10.5
56580  102.80   27.71 2461   24.9
56578  102.75   27.06 2461   20.4
53573  112.70   40.00 1337    1.9
53869  111.71   36.58  568    8.5