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

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 25341|回复: 52

对多时次大区域(多省份)maskout的小小改进

[复制链接]

新浪微博达人勋

发表于 2013-4-16 22:59:38 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 平流层的萝卜 于 2013-4-16 23:12 编辑

图形美化与边界控制,多时次MASKOUT命令的小小改进(2013.4.16)之前在http://bbs.06climate.com/forum.php?mod=viewthread&tid=13422 里,我就只想出华北的区域的能见度,这几天通过学习传说中的谁http://bbs.06climate.com/forum.php?mod=viewthread&tid=893&extra=&page=1http://bbs.06climate.com/forum.php?mod=viewthread&tid=8338还有clarehttp://bbs.06climate.com/forum.php?mod=viewthread&tid=7276关于maskout命令的说明,将图进一步做了美观,做到了只显示北京、河北、天津、山东、河南和山西6个地区的能见度图。基本做到了多时次mask出华北地区的任务因为严格意义上的华北地区仅仅包括山东和河南的黄河以北部分提出了多时次MASKOUT的小小改进。在自己原帖的基础上,用整合版grads的调色盘重做了能见度的色标。并对原始数据进行了筛选,只选择加密观测的08时、14时、20时数据进行出图。
话不多说,先上几张效果图:
1.png 2.png 3.png 4.png

关于多时次MASKOUT的小小改进:
关于maskout命令区域边界控制,上述两位大神都已经阐述很清楚,他们都提到了多时次maskout命令需要将mask的ctl文件的时间设置改成多时次,但我怎么运行都没弄出来,可能太笨了,给传说请教还暂未得回复~~。
最后通过尝试,我找出了一种出图的方法,不需要修改时间设置,tdef一个时次即可,也不用设置很多个ctl文件,只要在gs画图语句里,mask变量后边加(t=1),就大功告成了。如下:
'open E:\gdd\plot\nengjiandu_new\vis_02_08_14_20.ctl'
'open E:\gdd\plot\huabei.ctl'
'enable print E:\gdd\plot\nengjiandu_new\new_vis.gmf'
i=1
while(i<=124)
if(math_fmod(i,4)!=1)
'set grads off'
'set lon 110 120'
'set lat 33 43'
'set lev 1000'
'set t 'i''
'set mpdset cnriver cnworld '
'set map 1 1 9'
'set xlopts 1 6 0.2'
'set ylopts 1 6 0.2'
'set gxout shaded'
'run E:\gdd\huimai_color.gs'
'set clevs 0.1 0.3 1 10 15 20 25'
'set ccols 58 57 56 55 54 53 52 51 '
'd maskout(oacres(mask.2(t=1),vi),mask.2(t=1)-1.0)'
'cbarn.gs'
'q time'
x=subwrd(result,3)
'draw title 'x''
'print'
'c'
endif
i=i+1
endwhile
'disable print'
'reinit'
;

这是我的huabei.ctl:
DSET E:\gdd\plot\huabei.dat
TITLE Mask data
UNDEF -9999
XDEF 131 LINEAR 110 0.1
YDEF 121 LINEAR 31 0.1
ZDEF 1 LINEAR 1 1
TDEF 1 LINEAR 02z01jan2013 3hr
VARS 1
mask 0 99  background mask data
ENDVARS

关于是maskout底图如何制作,请参见传说中的谁大神的帖子http://bbs.06climate.com/forum.php?mod=viewthread&tid=893&extra=&page=1






评分

参与人数 2威望 +5 金钱 +40 贡献 +15 体力 +200 收起 理由
mofangbao + 5 + 20 + 10 + 200 平流层很冷的吧
topmad + 20 + 5 保卫萝卜

查看全部评分

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-4-16 23:21:43 | 显示全部楼层
十分给力。有一次推进了一步。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-4-16 23:22:04 | 显示全部楼层
十分给力。有一次推进了一步。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-4-17 06:47:56 | 显示全部楼层
很不错,学习了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-4-17 08:39:13 | 显示全部楼层
不错,不错!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-4-17 09:02:16 | 显示全部楼层
帖子总结的很全面 支持一下萝卜兄
密码修改失败请联系微信:mofangbao

新浪微博达人勋

0
早起挑战累计收入
发表于 2013-4-17 09:32:14 | 显示全部楼层
谢谢萝卜,笨笨这两天总是想保卫你
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-4-17 10:02:49 | 显示全部楼层
谢谢大家啊
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-4-17 18:35:11 | 显示全部楼层
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-4-17 22:45:19 | 显示全部楼层
mofangbao 发表于 2013-4-17 09:32
谢谢萝卜,笨笨这两天总是想保卫你

关于传说提出的修改mask描述ctl中的时间设置,我还是不懂。因为mask文件的大小对应的就是一个时次的,如果多个时次的话,mask文件本身的长度是否也应该相应地增大呢?(好比站点的格点数据,若是多时次,fortran生成的时候要加一层时间循环)。可是底图是从meteoinfo里截出来的,好像没法设置多个时次的,也就是说正常可用的mask及其maskctl应该只能是一个时次的吧?请教风神
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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