求助，画散度剖面图出问题

其微

RT

'reinit'
'sdfopen D:\mine\nc\uwnd.2012.nc'
'sdfopen D:\mine\nc\vwnd.2012.nc'
'define div=hdivg(Uwnd.1,Vwnd.2)*100000'
'set lon 120'
'set lat 20 40'
'set lev 1000 100'
'd ave(div,t=53,t=69)'
'printim D:\mine\201202\div.jpg x1000 y800 white'

想得到的是1000hpa到100hpa的散度剖面图，结果出来的时折线图，感觉挺简单的图，可是不会画

mofangbao

'set lev 1000 100'
'd ave(div,t=53,t=69)'
不要这样写，先求出每一层的平均值之后，输出到文件，然后重新打开那个文件画图，这样你的思路会清晰一点。

兰溪之水

'reinit'
'sdfopen D:\mine\nc\uwnd.2012.nc'
'sdfopen D:\mine\nc\vwnd.2012.nc'
'set lon 110 130'
'set lat 20 40'
'set lev 1000 100'
'set t 53 69'
'define div=hdivg(Uwnd.1,Vwnd.2)*100000'
'set lon 120'
'set lat 20 40'
'set lev 1000 100'
'd ave(div,t=53,t=69)'
'printim D:\mine\201202\div.jpg x1000 y800 white'

其微

谢谢啦，已解决，只要前面定义一个范围就可以画了~

其微

╮(╯▽╰)╭后记：
按照兰溪之水给出的gs画图，可以画出来，可是回到文件夹看图片发现和画出来的不一样，图上就几个乱码似的小点点，于是按着gs文件一步步画找答案，到'd ave(div,t=53,t=69)'这一步grads上出来很多英文，然后在清楚图片之前好像什么也做不了了。

因此，懒人处理法把gs改成了这样。。。连循环都懒得设置。。。主要是没什么时间了><
'reinit'
'sdfopen D:\mine\nc\uwnd.2012.nc'
'sdfopen D:\mine\nc\vwnd.2012.nc'
'set lon 110 130'
'set lat 15 45'
'set lev 1000 100'
'set t 53'
'define div53=hdivg(Uwnd.1,Vwnd.2)*100000'
'set t 54'
'define div54=hdivg(Uwnd.1,Vwnd.2)*100000'
'set t 55'
'define div55=hdivg(Uwnd.1,Vwnd.2)*100000'
'set t 56'
'define div56=hdivg(Uwnd.1,Vwnd.2)*100000'
'set t 57'
'define div57=hdivg(Uwnd.1,Vwnd.2)*100000'
'set t 58'
'define div58=hdivg(Uwnd.1,Vwnd.2)*100000'
'set t 59'
'define div59=hdivg(Uwnd.1,Vwnd.2)*100000'
'set t 60'
'define div60=hdivg(Uwnd.1,Vwnd.2)*100000'
'set t 61'
'define div61=hdivg(Uwnd.1,Vwnd.2)*100000'
'set t 62'
'define div62=hdivg(Uwnd.1,Vwnd.2)*100000'
'set t 63'
'define div63=hdivg(Uwnd.1,Vwnd.2)*100000'
'set t 64'
'define div64=hdivg(Uwnd.1,Vwnd.2)*100000'
'set t 65'
'define div65=hdivg(Uwnd.1,Vwnd.2)*100000'
'set t 66'
'define div66=hdivg(Uwnd.1,Vwnd.2)*100000'
'set t 67'
'define div67=hdivg(Uwnd.1,Vwnd.2)*100000'
'set t 68'
'define div68=hdivg(Uwnd.1,Vwnd.2)*100000'
'set t 69'
'define div69=hdivg(Uwnd.1,Vwnd.2)*100000'
'define

div=div53+div54+div55+div56+div57+div58+div59+div60+di

v61+div62+div63+div64+div65+div66+div67+div68+div69'
'set lon 120'
'set lat 20 40'
'set lev 1000 100'
'set grads off'
'd div/17'
'draw title (a)'
'printim D:\mine\201202\div.jpg x1000 y800 white'

期待看起来美妙点的gs文件~

兰溪之水

其微 发表于 2012-7-24 08:35

                               
登录/注册后可看大图

╮(╯▽╰)╭后记：
按照兰溪之水给出的gs画图，可以画出来，可是回到文件夹看图片发现和画出来的不一样， ...

怎么不一样？图贴出来~

其微

不好意思，才刚看到


就是这样，一片都是白的，只有些点点

你可不可以自己试一下，会不会出现一样的问题？

其微

就是grads上面的图是画好了的，但是没有title，好像画了图，在c之前的操作统统失灵了。
去看保存到文件夹的图片就是上面那个白的

米米面面

。。。。好像不行

18351819511

学习了，学习了