GrADS中站点分布图的画法【转】

mofangbao

转自人人网：http://blog.renren.com/blog/236180042/728802296

---

上图是黄河流域站点分布图，由于Grads中不能显示汉字，只能用拼音什么的来代替。

Grads有个指令叫做 draw string ，写字符串，还有个指令叫做draw mark，画标记，如上图的黑点。

以一个站点为例，Grads指令如下：

1. 'q w2xy  110.10  40.50'
   
2. x=subwrd(result,3)
   
3. y=subwrd(result,6)
   
4. 'draw string 'x' 'y' Baotou'
   
5. 'q w2xy  110.10  40.50'
   
6. x=subwrd(result,3)
   
7. y=subwrd(result,6)
   
8. 'draw mark 3 'x' 'y' siz '
   

'q w2xy  110.10  40.50'<br /> x=subwrd(result,3) <br /> y=subwrd(result,6) <br /> 'draw string 'x' 'y' Baotou'<br /> 'q w2xy  110.10  40.50'<br /> x=subwrd(result,3) <br /> y=subwrd(result,6) <br /> 'draw mark 3 'x' 'y' siz '<br />

*3号标记就是实心圆，siz是标记的大小

这只是一个站的情况，而全国那么多站点，不可能一个一个手打，所以这时候需要用电脑来帮你，编fortran小程序。

1. program main
   
2. integer stn
   
3. parameter(stn=160) !定义一个常数，用以存储站点数量
   
4. real sta(stn,2)    !用来存储经纬度
   
5. character(len=10) st_name(20)  !存储站名或站号
   
6. integer i , j
   
7. character(len=80) str(7)  !定义一个长度为80的数组
   
8. real tmp open(1,file="f:\test\zd.txt")     !站点文件  分别是 经纬度 站点名或站号
   
9. do i = 1 , stn
   
10. read(1,*) sta(i,1) , sta(i,2) , st_name(i)
    
11. enddo
    
12. open(2,file="f:\test\test\out1.txt")
    
13. str(1) = "'q w2xy "
    
14. str(2) = "'"
    
15. str(3) = "x=subwrd(result,3)"
    
16. str(4) = "y=subwrd(result,6)"
    
17. str(5) = "'draw mark 3 "
    
18. str(6) = " 'x' 'y' "
    
19. str(7) = "'"
    
20. tmp = 0.1              !表示标记的大小，这里设置为0.1
    
21. do i = 1, stn
    
22. write(2,'(a8\f7.2,f7.2\a1/)') str(1),sta(i,1),sta(i,2),str(2)
    
23. write(2,'(a19/a19)') str(3),str(4)
    
24. write(2,'(a13\a9\f3.1\a1/)') str(5),str(6),tmp,str(7)
    
25. enddo close(1)
    
26. close(2) end

program main<br /> integer stn<br /> parameter(stn=160) !定义一个常数，用以存储站点数量<br /> real sta(stn,2)    !用来存储经纬度<br /> character(len=10) st_name(20)  !存储站名或站号<br /> integer i , j <br /> character(len=80) str(7)  !定义一个长度为80的数组<br /> real tmp open(1,file="f:\test\zd.txt")     !站点文件  分别是 经纬度 站点名或站号<br /> do i = 1 , stn<br /> read(1,*) sta(i,1) , sta(i,2) , st_name(i)<br /> enddo<br /> open(2,file="f:\test\test\out1.txt")<br /> str(1) = "'q w2xy "<br /> str(2) = "'"<br /> str(3) = "x=subwrd(result,3)"<br /> str(4) = "y=subwrd(result,6)"<br /> str(5) = "'draw mark 3 "<br /> str(6) = " 'x' 'y' "<br /> str(7) = "'"<br /> tmp = 0.1              !表示标记的大小，这里设置为0.1<br /> do i = 1, stn<br /> write(2,'(a8\f7.2,f7.2\a1/)') str(1),sta(i,1),sta(i,2),str(2)<br /> write(2,'(a19/a19)') str(3),str(4)<br /> write(2,'(a13\a9\f3.1\a1/)') str(5),str(6),tmp,str(7)<br /> enddo close(1)<br /> close(2) end

程序运行完之后，在txt文件中会得到所需grads代码，之后直接复制到gs文件中，即可使用。

提示：还可以使用 stnmark达到标记的效果，当然，使用MeteoInfo实现这个功能会更轻松！

飛羽_y

其实有个简单的方法。。 把变量的颜色 调成和底图一样

nana

open(2,file="f:\test\test\out1.txt")
这步的txt文件是什么啊

nana

我按着楼主的程序编译 有错误…

夏天爱裤衩

nuister 发表于 2012-5-9 14:39

                               
登录/注册后可看大图

这个出来的gs文件只有在display一个变量之后才能用，怎样绘制只含站点分布的图，楼主知道么?

在display之前可以加上set cmax -100，就可以不显示任何变量，只有站点分布

mofangbao

麻花先生 发表于 2012-5-9 20:37

                               
登录/注册后可看大图

请问下我的数据时这样的，我想画站号：
部分数据：
126.65        51.72        50353

上面几楼已经说清楚了，需要先d一个变量才可以，不是直接运行就能出图的。
南海的问题已经在另外一个帖子中回复了你

nuister

这个出来的gs文件只有在display一个变量之后才能用，怎样绘制只含站点分布的图，楼主知道么?

mofangbao

nuister 发表于 2012-5-9 14:39

                               
登录/注册后可看大图

这个出来的gs文件只有在display一个变量之后才能用，怎样绘制只含站点分布的图，楼主知道么?

你楼下的已经说了，就是那样的~

tintin

先学习一下

麻花先生

请问下我的数据时这样的，我想画站号：
部分数据：
126.65        51.72        50353
121.92        48.77        50632
119.75        49.22        50527
121.7        50.45        50434
125.23        49.17        50557
123.92        47.38        50745
126.97        47.43        50756
131.98        47.23        50788
130.17        46.49        50873
130.95        45.28        50978
126.77        45.75        50953
129.6        44.57        54094
122.05        46.08        50838

运行出来的txt是这样的：
'q w2xy  126.65  51.72'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy  121.92  48.77'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy  119.75  49.22'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy  121.70  50.45'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'
'q w2xy  125.23  49.17'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.1'

然后将它复制到一个新建的gs文件中，画不出图呢？怎么回事啊？请楼主不吝赐教！！初学者好多不懂，还有就是怎么画南海诸岛的图，看了帖子之后没有找到说的那个grads安装目录下的lib文件夹，我用的是论坛的grads，请清风楼主帮助啊！！非常感谢！！！

麻花先生

mofangbao 发表于 2012-5-9 20:46

                               
登录/注册后可看大图

上面几楼已经说清楚了，需要先d一个变量才可以，不是直接运行就能出图的。
南海的问题已经在另外一个帖子 ...

非常感谢，站点分布图已经做出来了。感谢清风的无私帮助、、、、