如何实现不同图形使用同一比色刻度尺？

大老鹰

如附件所示，想做一个gif动画。但是图里1月和8月比色刻度尺不一致，系统默认是按照比例赋色，现在需要按照固定值赋色。
另外，矢量箭头也存在这个问题，请问应该如何处理？

1. 'reinit'
   
2. 'sdfopen D:\NOAA\ESRL.2.5\air.mon.ltm.nc'
   
3. 'sdfopen D:\NOAA\ESRL.2.5\rhum.mon.ltm.nc'
   
4. 'sdfopen D:\NOAA\ESRL.2.5\uwnd.mon.ltm.nc'
   
5. 'sdfopen D:\NOAA\ESRL.2.5\vwnd.mon.ltm.nc'
   
6. *输入文件
   
7. 'set mpdraw off'
   
8. *关闭原有底图
   
9. 'set lon 30 150'
   
10. 'set lat -10 45'
    
11. *设置经纬度
    
12. 'set lev 700'
    
13. t=1
    
14. while(t<=12)
    
15. *设置层级，时间起始
    
16. 'set t 't''
    
17. 'set grads off'
    
18. 'define p=lev'
    
19. 'define e=exp(17.67*(air-273.15)/(air-29.65))*(rhum.2/100)*6.112'
    
20. 'define q=622*e/(p-0.378*e)'
    
21. 'define q1=uwnd.3*q/9.8'
    
22. 'define q2=vwnd.4*q/9.8'
    
23. *水汽通量计算
    
24. 'set gxout shaded'
    
25. 'run D:\NOAA\200colors.gs'
    
26. 'set cint 300'
    
27. 'run cbarn.gs'
    
28. *阴影图，导入论坛GrADS色带，间隔300
    
29. 'set gxout vector'
    
30. 'd mag(q1,q2)'
    
31. 'set ccolor 4'
    
32. 'd uwnd.3*2.5;vwnd.4*2.5'
    
33. *矢量图，设置U\V分量，蓝色，U\V分量频率
    
34. 'set gxout shp'
    
35. 'set line 1'
    
36. 'draw shp D:\NOAA\country.shp'
    
37. *shp图，白色线形，画shp图 文件路径
    
38. 'q dim'
    
39. a=sublin(result,5)
    
40. a=subwrd(a,9)
    
41. 'draw title 'a
    
42. 'set rgb 50 255 255 255'
    
43. *查询维度，导入月份时间作为标题。
    
44. 'printim D:\NOAA\'%t%'.png x1300 y650 white'
    
45. 'C'
    
46. t=t+1
    
47. *输出图片，清屏，循环。
    
48. endwhile
    

复制代码

country.rar (1.12 MB, 下载次数: 0)

2018-8-5 08:53 上传

点击文件名下载附件 200colors.gs (6.1 KB, 下载次数: 0)

2018-8-5 08:52 上传

点击文件名下载附件
water_vapor_vector.gs (1.17 KB, 下载次数: 0)

2018-8-5 08:52 上传

点击文件名下载附件

river

大老鹰 发表于 2018-8-5 09:57
经过论坛里搜索
把

真逗，你得先display，才能'run cbarn.gs' 。屏幕上什么都没有，它不报错才怪呢

yuchen

set vrange 试试。

大老鹰

经过论坛里搜索
把

'set gxout shaded'
'run D:\NOAA\200colors.gs'
'set cint 300'
'run cbarn.gs'

改成

'set gxout shaded'
'run D:\NOAA\200colors.gs'
'set clevs 300 600 900 1200 1500 1800 2100 2400 2700 3000 3300'
'run cbarn.gs'

色带倒是统一了，但是色度卡无法显示。
提示
cannot plot color bar no shading information

大老鹰

把
'run cbarn.gs'
命令放到后面就好啦

大老鹰

良好脚本，色卡不一致问题解决。

'reinit'
'sdfopen D:\NOAA\ESRL.2.5\air.mon.ltm.nc'
'sdfopen D:\NOAA\ESRL.2.5\rhum.mon.ltm.nc'
'sdfopen D:\NOAA\ESRL.2.5\uwnd.mon.ltm.nc'
'sdfopen D:\NOAA\ESRL.2.5\vwnd.mon.ltm.nc'
*输入文件
'set mpdraw off'
*关闭原有底图
'set lon 30 150'
'set lat -10 45'
*设置经纬度
'set lev 700'
t=1
while(t<=12)
*设置层级，时间起始
'set t 't''
'set grads off'
'define p=lev'
'define e=exp(17.67*(air-273.15)/(air-29.65))*(rhum.2/100)*6.112'
'define q=622*e/(p-0.378*e)'
'define q1=uwnd.3*q/9.8'
'define q2=vwnd.4*q/9.8'
*水汽通量计算
'set gxout shaded'
'run D:\NOAA\200colors.gs'
'set clevs 300 600 900 1200 1500 1800 2100 2400 2700 3000 3300'

*阴影图，设置色带，间隔300
'set gxout vector'
'd mag(q1,q2)'
'set ccolor 4'
'd uwnd.3*2.5;vwnd.4*2.5'
*矢量图，设置U\V分量，蓝色，U\V分量频率
'set gxout shp'
'set line 1'
'draw shp D:\NOAA\country.shp'
*shp图，白色线形，画shp图 文件路径
'q dim'
a=sublin(result,5)
a=subwrd(a,9)
'draw title 'a
'set rgb 50 255 255 255'
'run cbarn.gs'
*查询维度，导入月份时间作为标题。
'printim D:\NOAA\'%t%'.png x1300 y650 white'
'C'
t=t+1
*输出图片，清屏，循环。
endwhile

大老鹰

最后一个问题，如何统一图中矢量箭头大小？

大老鹰

箭头长度这里有一个解释，不过看得不是很明白。我想按照风力大小给箭头长度赋值
http://gradsusr.org/pipermail/gradsusr/2014-September/037460.html

大老鹰

《实用手册》里的语法是这样的

set arrscl size <magnitude>
设置矢量箭头的长度为size 英寸（虚页英寸），选项magnitude 的值设定矢量箭头的大小，
缺省时所有矢量同长，c 或d 命令重新设置。
set arrowhead size
设置箭头大小，缺省为0.05，如取为0 不画头，如取为负与矢量值成比例。

大老鹰

回复出错，删除