爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 7228|回复: 15

[脚本编辑] 求助 想要画冬季平均500hPa 位势高度场 如图 结果自己画不出来 成下图这样 求解释

[复制链接]
发表于 2016-3-8 11:44:25 | 显示全部楼层 |阅读模式

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

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

x
各位大神  走过路过的帮个忙  学地理的学生一枚  现在想要做冬季平均500hPa 位势高度场图  可是 不会啊 完全不会编程  然后看了一个月了 还是不会  写了一个GS文件   结果运行出来成这个样子了  求解释。我用到的是1987 1979 1999 201 2002年的nc数据  不会编写ctl文件就直接写了gs文件。如下
'sdfopen D:\hgt.mon.mean.nc'
'set lon 60 150'
'set lat 10 70'
'set lev 500'
'define h1=ave(hgt,t=372,t=374,1)'
'define h2=ave(hgt,t=612,t=614,1)'
'define h3=ave(hgt,t=636,t=639,1)'
'define h4=ave(hgt,t=648,t=650,1)'
'define h5=ave(hgt,t=708,t=710,1)'
'define ha=(h1+h2+h3+h4+h5)/5'
'set gxout contour'
'set ccolor rainbow'
'set cint 5'
'set clevs 5360'
'set cthick 6'
'set ccols 2'
'd ha'
'printim d:h500.png white'
'print'
'disable print'
真的实在是不会了   运行出来的图是这样   没人指导的文科妹子一枚   求指导   想要的图如下  大神们能不能帮我完善一下程序啊

想要的图阴影区为距平位势高度绝对值大于15 gpm 区域

想要的图阴影区为距平位势高度绝对值大于15 gpm 区域

自己做的傻图

自己做的傻图
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2016-3-8 11:45:31 | 显示全部楼层
大神们   快点帮忙指导知道啊    好多天都不能睡个安心觉了{:eb303:}{:eb303:}{:eb303:}
密码修改失败请联系微信:mofangbao
发表于 2016-3-8 14:57:42 | 显示全部楼层
'set clevs 5360'
这个语句的意思是只画出5360这条等值线,当然只能画成你那个样子
把它去掉就好了
密码修改失败请联系微信:mofangbao
发表于 2016-3-8 15:05:41 | 显示全部楼层
去掉'set clevs 5360'
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2016-3-10 21:02:07 | 显示全部楼层
river 发表于 2016-3-8 15:05
去掉'set clevs 5360'

我现在想求距平然后填色 为什么一直不行啊
'sdfopen D:\hgt.mon.mean.nc'
'set lon 60 150'
'set lat 10 70'
'set lev 500'
'set mpdset hires'
'set map color 0'
'set map thickness 3'
'set grads off'
'define h1=ave(hgt,t=372,t=374,1)'
'define h2=ave(hgt,t=612,t=614,1)'
'define h3=ave(hgt,t=636,t=639,1)'
'define h4=ave(hgt,t=648,t=650,1)'
'define h5=ave(hgt,t=708,t=710,1)'
'define h6=ave(hgt.1,t+155,t+803,12)'
'define h7=ave(hgt.1,t+156,t+804,12)'
'define h8=ave(hgt.1,t+157,t+805,12)'
'define h9=(h1+h2+h3+h4+h5)/5'
'define  hju=h9-(h6+h7+h8)/3'
'set gxout shaded'
'set rbcols 5'
'set cblack -15 15'
'run charc.gs'
'd hju'
'set gxout contour'
'set ccolor rainbow'
'set cint 30'
'set cthick 2'
'set ccols 2'
'd h9'
'printim d:/h5001.png white'
'print'
'disable print'
然后出现问题一直  实在不会了  各种设置也是在手册里看了  再看了论坛的  还是理解不透彻  望指教啊
出错点和做出的图   图不知道怎嘛上传了
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2016-3-10 21:02:46 | 显示全部楼层
C:\Users\thk229\Desktop
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2016-3-10 21:06:09 | 显示全部楼层
river 发表于 2016-3-8 15:05
去掉'set clevs 5360'

就是说 地图设置错误 还有'set cblack -15 15'
'run charc.gs'  这两句设置错误   出图全部是蓝色的
密码修改失败请联系微信:mofangbao
发表于 2016-3-10 22:07:44 | 显示全部楼层
lp0517 发表于 2016-3-10 21:02
我现在想求距平然后填色 为什么一直不行啊
'sdfopen D:\hgt.mon.mean.nc'
'set lon 60 150'

错误提示还有错误的图传上来
密码修改失败请联系微信:mofangbao
发表于 2016-3-11 08:17:51 | 显示全部楼层
粗看了一下,跳过计算部分,前面的应该没问题。从下面开始给你解释一下这些语句是什么效果,你自己看看哪里有问题,自己改。
'set gxout shaded'----------------------------------绘阴影
'set rbcols 5'----------------------------------------阴影填充色为4(蓝色)
'set cblack -15 15'--------------------------------不画-15~15之间的值(这一句没用,500hPa高度场没这种值)
'run charc.gs'------------------------------------你想干什么?是否要cbarn.gs?
'd hju'--------------------------------------画变量hju为阴影(其实只有蓝色)
'set gxout contour'----------------------设置所绘变量为等值线输出
'set ccolor rainbow'---------------------变量颜色用彩虹色表示
'set cint 30'-----------------------------变量间距30
'set cthick 2'---------------------------线条粗细值为2(很细)
'set ccols 2'----------------------------设置变量颜色为2(红色)
'd h9'------------------------------------绘制变量h9

如果你前面没问题,绘出的图至多是蓝色背景上有几条很细的红色等值线。改的思路是:
1、先set gxout shaded,画hju为阴影。
2、 run cbarn.gs,出色标。
3、set gxout contour--------------------------开始绘等值线
4、set cthick (4-6差不多),set ccolor 1(黑色等值线),set cint (***,等值线间隔)
5、d h9----------------------------画h9的等值线
6、最后出图。

Ps:如果要在留言里附图的话,360浏览器有时候会出问题,画个浏览器试试。
密码修改失败请联系微信:mofangbao
发表于 2016-3-11 08:35:34 | 显示全部楼层
唉,智能看看了
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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