对多时次大区域（多省份）maskout的小小改进

平流层的萝卜

图形美化与边界控制，多时次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=1和http://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时数据进行出图。
话不多说，先上几张效果图：


关于多时次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。

kongfeng0824

十分给力。有一次推进了一步。

kongfeng0824

十分给力。有一次推进了一步。

river

很不错，学习了

曾记否

不错，不错！

topmad

帖子总结的很全面 支持一下萝卜兄

mofangbao

谢谢萝卜,笨笨这两天总是想保卫你

平流层的萝卜

谢谢大家啊

平流层的萝卜

mofangbao 发表于 2013-4-17 09:32

                               
登录/注册后可看大图

谢谢萝卜,笨笨这两天总是想保卫你

受宠若惊啊

平流层的萝卜

mofangbao 发表于 2013-4-17 09:32

                               
登录/注册后可看大图

谢谢萝卜,笨笨这两天总是想保卫你

关于传说提出的修改mask描述ctl中的时间设置，我还是不懂。因为mask文件的大小对应的就是一个时次的，如果多个时次的话，mask文件本身的长度是否也应该相应地增大呢？（好比站点的格点数据，若是多时次，fortran生成的时候要加一层时间循环）。可是底图是从meteoinfo里截出来的，好像没法设置多个时次的，也就是说正常可用的mask及其maskctl应该只能是一个时次的吧？请教风神