关于站点资料画格点图的一点摸索

youxia1873

经过一周的摸索，终于发现之前错误的原因。
将站点数据写为二进制数据时要特别注意一点：
例：

1. 58128  117.31   33.30   24    8.3
   
2. 55472   88.63   30.95 4671    2.4
   
3. 53522  107.85   40.05 1185    4.9
   
4. 57313  106.76   31.85  360    0.4
   
5. 57646  109.76   28.56  302   15.2
   
6. 57566  111.65   29.45   51   12.4
   
7. 53965  111.56   35.51  762    1.4
   
8. 57060  111.23   35.16  404   10.5
   
9. 56580  102.80   27.71 2461   24.9
   
10. 56578  102.75   27.06 2461   20.4
    
11. 53573  112.70   40.00 1337    1.9
    
12. 53869  111.71   36.58  568    8.5
    
13. 53865  111.55   36.65 1013   0.01

复制代码

其排列格式为：站点  纬度  经度  高度  降雨量

fortran编译读写程序时要注意，写入二进制文件时需要将经度放在纬度前面
例：

1. open(2,file='F:/station/140824r.grd',form='binary')
   
2. tim = 0.0
   
3. nlev = 1
   
4. nflag = 1
   
5. do i = 1,n
   
6.    write(2) stid(i),lon(i),lat(i),tim,nlev,nflag,rain(i)
   
7. enddo
   
8. nlev = 0
   
9. write(2) stid(n-1),<font color="#ff0000">lat(n-1),lon(n-1)</font>,tim,nlev,nflag
   
10. close(2)
    

复制代码

最终出图：
station.zip (151.74 KB, 下载次数: 16)

2014-10-9 15:38 上传

点击文件名下载附件
程序包

程序已给出，其中用到cnbasemap.gs 和rainbow.gs文件可在大神们的帖子中找到，就不再赘述了
{:soso_e183:}{:soso_e183:}

river

楼主是不是你自己记错了呢，GrADS要求的站点数据的结构必须是纬度在经度前面。这里有证据http://grads.iges.org/grads/gadoc/aboutstationdata.html。而且我自己做过很多了，肯定是纬度在前

youxia1873

scripts.zip 里只有cnbasemap.gs 和rainbow.gs ，已经有的不需要下载，自行修改其中参数即可

排山倒海

是纬度在前面的，前面站点信息的数据给错了吧

youxia1873

river 发表于 2014-10-9 17:41
楼主是不是你自己记错了呢，GrADS要求的站点数据的结构必须是纬度在经度前面。这里有证据http://grads.iges ...

记错是没有的,当我在写入时将纬度放前面时,程序运行结果是错误的...
我用的grads是咱们论坛2.01a版本的,原因目前还不知道.
不过以后再有类似问题,就知道怎么改了,多谢指正

novelbean

学习学习，谢谢

shanniu728

我也想学习学习。

差不多先生

{:eb502:}

fujianbeige

楼主辛苦了，谢谢分享。