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

maliyasunny · 发表于 2013-4-24 08:35:47

谢谢呢，我懂了，后来我也试了一下这样子的，也可以画出来呢，两种方法都可以
'open d:/k/zddata/7cl/74/znr/yue/huanan.ctl'
'define m=mask.1'
'set lon 102 125'
'set lat 18 30'
'set t 1'
'open d:/k/zddata/7cl/74/znr/yue/'mo'.ctl'
i=1
while(i<=51)
'set lon 102 125'
'set lat 18 30'
'set t 'i''
'set grads off'
'set mpdset huanan'
'define a=oacres(m,r.2,50,20,10,5,1)'
'define a1=maskout(a,m-1.0)'
'define aa=smth9(a1)'
'set gxout shaded'
'd aa'
'set gxout contour'
'd aa'
'cbarn'
'printim d:/k/zddata/7cl/74/znr/yue/tu/'mo'.'i+1959'.png white'
i=i+1
'c'
endwhile

平流层的萝卜 · 发表于 2013-4-24 12:41:08

苏绮发表于 2013-4-23 22:03

登录/注册后可看大图

为啥我总是在运行自定义色标这里出问题呢？？run green-65.gs 但是仍是彩虹图。着急死了。。。可能是哪里的 ...

gs怎么写的呢？

平流层的萝卜 · 发表于 2013-4-24 12:41:38

maliyasunny 发表于 2013-4-24 08:35

登录/注册后可看大图

谢谢呢，我懂了，后来我也试了一下这样子的，也可以画出来呢，两种方法都可以
'open d:/k/zddata/7cl/74/z ...

好的，不客气！

丹徵 · 发表于 2013-5-8 21:41:06

这是正常的图
请大神们帮看看，我是maskout高原的，有了那个高原的边界ctl，和dat ，后来编了gs，结果出了问题

这是我的gs

求指点迷津

丹徵 · 发表于 2013-5-8 21:43:00

图不见了
正常图

错误

我的gs

平流层的萝卜 · 发表于 2013-5-9 10:07:19

丹徵发表于 2013-5-8 21:43

登录/注册后可看大图

图不见了
正常图
错误

你的mm.ctl是怎么写的，是maskout的ctl吗？怎么看起来是格点的？贴出来看看吧~

丹徵 · 发表于 2013-5-9 13:24:41

平流层的萝卜发表于 2013-5-9 10:07

登录/注册后可看大图

你的mm.ctl是怎么写的，是maskout的ctl吗？怎么看起来是格点的？贴出来看看吧~

我的mm.ctl和mm.dat是高原轮廓的

mm.ctl
DSET F:\snoweof\mm.dat
TITLE Mask data
UNDEF -9999
XDEF 51 LINEAR 60 1
YDEF 31 LINEAR 20 1
ZDEF 1 LINEAR 1 1
tdef 1 linear Jun1951 1yr
VARS 1
g 0 99 background mask data
ENDVARS

复制代码

时间什么的都没有问题

丹徵 · 发表于 2013-5-9 13:26:25

这是mm.dat

mm.dat (6.18 KB, 下载次数: 0)

平流层的萝卜 · 发表于 2013-5-9 15:14:30

丹徵发表于 2013-5-9 13:26

登录/注册后可看大图

这是mm.dat

在gs里，把所有g.1改成g.1(t=1),再试试

丹徵 · 发表于 2013-5-9 18:25:08

平流层的萝卜发表于 2013-5-9 15:14

登录/注册后可看大图

在gs里，把所有g.1改成g.1(t=1),再试试

灰常感谢大神，能画出来了，就是好丑，而且没有边界，而且我想请教一下，为什么别的时次不行？我的资料是做过eof的有8个特征场的资料，正常所有图都能出，但这个maskout就只能画出第一张。。

		自动登录	找回密码
密码			立即注册