任意两点间斜剖面做法总结【利用GRADS自带函数】

ml11178

OK了，我有畫出剖面，但是又遇到一點問題，我先貼出我的圖跟GS檔

GS檔:

'set grads off'
'set zlog on'
'set csmooth on'
'set x 1'
'set y 1'
'set lev 1000 700'
lon1 = 121.465886
lon2 = 121.539731
lat1 = 25.1108441
lat2 = 24.9693381
'define alfa=atan2('lat2-lat1','lon2-lon1')'
lon = lon1
'collect 1 free'
'collect 2 free'
'collect 3 free'
'collect 4 free'
while (lon <= lon2)
  lat = lat1 + (lat2-lat1)*(lon-lon1) / (lon2-lon1)
  'collect 1 gr2stn(U*cos(alfa)+V*sin(alfa),'lon','lat')'
  'collect 2 gr2stn(W*3,'lon','lat')'
  'collect 3 gr2stn(hgt,'lon','lat')'
  'collect 4 gr2stn(rh,'lon','lat')'
  lon = lon + 0.006
endwhile



'set map 15 1 9'
'set parea 1 10.2 1.5 8'
'set xlopts 1 6 0.1'
'set ylopts 1 6 0.1'
'set font 2'

'set x 83 85'
'set xaxis 'lon1' 'lon2
*'set clab on'
*'set gxout shaded'
*'set bargap 0'
*'set clab off'
*'set clevs 20'
*'set ccols 0 5 0'
*'d coll2gr(3,-u)'
*
'set rgb 16 255 255 255'
'set rgb 17 255 234 255'
'set rgb 18 255 213 255'
'set rgb 19 255 191 255'
'set rgb 20 255 170 255'
'set rgb 21 255 149 255'
'set rgb 22 212 107 234'
'set rgb 23 170 85 213'
'set rgb 24 128 64 191'
'set rgb 25 85 43 170'
'set rgb 26 42 21 149'
'set rgb 27 0 0 128'
'set ccols 16 17 18 19 20 21 22 23 24 25 26 27'
'set clevs 0.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 80.0 90.0 100.0'
'set gxout shaded'
'd maskout(coll2gr(4,-u),-coll2gr(3,-u))'
'cbarn.gs'
*

'set ccolor 0'
'set gxout vector'
'set arrscl 1 5'
'set clab on'
'd maskout(coll2gr(1,-u),-coll2gr(3,-u));maskout(coll2gr(2,-u),-coll2gr(3,-u));mag(coll2gr(1,-u),coll2gr(2,-u))'


畫出來的圖:




問題來了，我在GS檔裡面有設定vector的顏色為黑色，可是不管怎麼改都還是沒辦法弄成黑色，甚至是加粗也不行

請問大大是哪裡有問題嗎? 感謝感謝ORZ

river

ml11178 发表于 2014-3-17 15:28

                               
登录/注册后可看大图

OK了，我有畫出剖面，但是又遇到一點問題，我先貼出我的圖跟GS檔

GS檔:

'set ccolor 0'放在'set gxout vector'后面试一下

ml11178

river 发表于 2014-3-17 15:48

                               
登录/注册后可看大图

'set ccolor 0'放在'set gxout vector'后面试一下

一樣ORZ.... 我在思考是不是最後一行的display會綁住顏色，可是想不到替換的display方法...

river

ml11178 发表于 2014-3-17 15:57

                               
登录/注册后可看大图

一樣ORZ.... 我在思考是不是最後一行的display會綁住顏色，可是想不到替換的display方法...

哦，我忘记了，你把后面那个mag 的去掉就可以了，后面那个就是控制让矢量变成彩色的

藍鶄

貌似是好东西 赞一个

liusx2007

river 发表于 2014-3-15 11:53

                               
登录/注册后可看大图

没有，自己把2.5*2.5的插值一下不就可以了吗

谢river 楼主！
    请求仔细讲解。应该怎样插值。是否有现成的程序？

river

liusx2007 发表于 2014-3-20 10:14

                               
登录/注册后可看大图

谢river 楼主！
    请求仔细讲解。应该怎样插值。是否有现成的程序？

这个可以利用grads自带的插值功能re()来实现，论他有人发过这个功能的用法，你搜索一下吧

要想飞

求哪位大侠发给我一个地形数据文件！画剖面图用的，俺级别不够下不了！QQ：150341733，先道声谢谢！

river

要想飞 发表于 2014-3-20 23:12

                               
登录/注册后可看大图

求哪位大侠发给我一个地形数据文件！画剖面图用的，俺级别不够下不了！QQ：150341733，先道声谢谢！

这个只是限制了毛毛雨而已，你去看看别的帖子，回复一下，很快就能升到小雨

小虹

非常好的例子，谢谢分享