请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

新浪微博登陆

只需一步, 快速开始

QQ登录

只需一步,快速开始

搜索
查看: 601|回复: 11

[脚本编辑] 如何实现不同图形使用同一比色刻度尺?

[复制链接] |关注本帖

新浪微博达人勋

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
如附件所示,想做一个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)
8.png
1.png
密码修改失败请联系qq:937062711

新浪微博达人勋

发表于 2018-8-5 10:18:46 | 显示全部楼层 |取消关注该作者的回复
大老鹰 发表于 2018-8-5 09:57
经过论坛里搜索

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

评分

参与人数 1金钱 +10 贡献 +2 收起 理由
大老鹰 + 10 + 2 谢谢提醒,才学习一天,多多指教。

查看全部评分

密码修改失败请联系qq:937062711
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2018-8-5 09:45:13 | 显示全部楼层 |取消关注该作者的回复
set vrange 试试。

评分

参与人数 1金钱 +10 贡献 +2 收起 理由
大老鹰 + 10 + 2 谢谢您的热心

查看全部评分

密码修改失败请联系qq:937062711

新浪微博达人勋

 楼主| 发表于 2018-8-5 09:57:43 | 显示全部楼层 |取消关注该作者的回复
经过论坛里搜索

'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
密码修改失败请联系qq:937062711

新浪微博达人勋

 楼主| 发表于 2018-8-5 10:01:10 | 显示全部楼层 |取消关注该作者的回复
本帖最后由 大老鹰 于 2018-8-5 10:21 编辑


'run cbarn.gs'
命令放到后面就好啦
密码修改失败请联系qq:937062711

新浪微博达人勋

 楼主| 发表于 2018-8-5 10:02:19 | 显示全部楼层 |取消关注该作者的回复
本帖最后由 大老鹰 于 2018-8-5 10:23 编辑

良好脚本,色卡不一致问题解决。
'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

密码修改失败请联系qq:937062711

新浪微博达人勋

 楼主| 发表于 2018-8-5 10:02:54 | 显示全部楼层 |取消关注该作者的回复
最后一个问题,如何统一图中矢量箭头大小?
密码修改失败请联系qq:937062711

新浪微博达人勋

 楼主| 发表于 2018-8-5 10:20:57 | 显示全部楼层 |取消关注该作者的回复
箭头长度这里有一个解释,不过看得不是很明白。我想按照风力大小给箭头长度赋值
http://gradsusr.org/pipermail/gradsusr/2014-September/037460.html
密码修改失败请联系qq:937062711

新浪微博达人勋

 楼主| 发表于 2018-8-5 10:23:02 | 显示全部楼层 |取消关注该作者的回复
《实用手册》里的语法是这样的

set arrscl size <magnitude>
设置矢量箭头的长度为size 英寸(虚页英寸),选项magnitude 的值设定矢量箭头的大小,
缺省时所有矢量同长,c 或d 命令重新设置。
set arrowhead size
设置箭头大小,缺省为0.05,如取为0 不画头,如取为负与矢量值成比例。
密码修改失败请联系qq:937062711

新浪微博达人勋

 楼主| 发表于 2018-8-6 00:25:37 | 显示全部楼层 |取消关注该作者的回复
本帖最后由 大老鹰 于 2018-8-6 07:56 编辑

回复出错,删除
密码修改失败请联系qq:937062711
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表