爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 5418|回复: 11

[脚本编辑] Grads的底图设置的求助

[复制链接]

新浪微博达人勋

发表于 2018-4-23 13:26:31 | 显示全部楼层 |阅读模式

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

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

x
刚开始接触Grads,计算研究区的水汽通量,已输出流域的水汽通量图(使用的ERA数据,附上下载地址:http://bbs.06climate.com/forum.p ... on=newthread&fid=36),
但在底图设置(参考“传说中的谁”使用meteinfo设置)时出现错误,报错未定义格网,不知道怎么回事?
请大家看看我写的“ERA数据成图的程序”对着不?
sdfopen函数是不是会自动对ERA数据(*Scale_factor+add_offect),不需要再手动进行计算?
烦请大家帮忙解答,谢谢。

ERA数据成图gs文件如下:
ERA数据成图.jpg
meteinfo输出如下:
meteinfo输出文件.jpg
报错信息如下:
报错信息.jpg


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

新浪微博达人勋

发表于 2018-4-23 18:17:25 | 显示全部楼层
1、slhall.dat的ctl里TDEF后0改成1
2、gs里d mask.2后加上(t=1)
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-4-23 20:29:17 | 显示全部楼层
balfulosa 发表于 2018-4-23 18:17
1、slhall.dat的ctl里TDEF后0改成1
2、gs里d mask.2后加上(t=1)

谢谢您的回答,我是新手,没有金币和积分奖励,不好意思。
还有一个问题就是我按照您的建议修改了,结果是重合了,但是显示的流域范围是以等高线形式显示的,我使用了“set gxout line”,但是没有效果,您看看怎么回事,谢谢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-4-23 20:29:22 | 显示全部楼层
balfulosa 发表于 2018-4-23 18:17
1、slhall.dat的ctl里TDEF后0改成1
2、gs里d mask.2后加上(t=1)

谢谢您的回答,我是新手,没有金币和积分奖励,不好意思。
还有一个问题就是我按照您的建议修改了,结果是重合了,但是显示的流域范围是以等高线形式显示的,我使用了“set gxout line”,但是没有效果,您看看怎么回事,谢谢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-4-23 20:29:29 | 显示全部楼层
balfulosa 发表于 2018-4-23 18:17
1、slhall.dat的ctl里TDEF后0改成1
2、gs里d mask.2后加上(t=1)

谢谢您的回答,我是新手,没有金币和积分奖励,不好意思。
还有一个问题就是我按照您的建议修改了,结果是重合了,但是显示的流域范围是以等高线形式显示的,我使用了“set gxout line”,但是没有效果,您看看怎么回事,谢谢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-4-23 20:29:32 | 显示全部楼层
balfulosa 发表于 2018-4-23 18:17
1、slhall.dat的ctl里TDEF后0改成1
2、gs里d mask.2后加上(t=1)

谢谢您的回答,我是新手,没有金币和积分奖励,不好意思。
还有一个问题就是我按照您的建议修改了,结果是重合了,但是显示的流域范围是以等高线形式显示的,我使用了“set gxout line”,但是没有效果,您看看怎么回事,谢谢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-4-23 20:29:45 | 显示全部楼层
balfulosa 发表于 2018-4-23 18:17
1、slhall.dat的ctl里TDEF后0改成1
2、gs里d mask.2后加上(t=1)

谢谢您的回答,我是新手,没有金币和积分奖励,不好意思。
还有一个问题就是我按照您的建议修改了,结果是重合了,但是显示的流域范围是以等高线形式显示的,我使用了“set gxout line”,但是没有效果,您看看怎么回事,谢谢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-4-24 10:08:41 | 显示全部楼层
Zhangzh 发表于 2018-4-23 20:29
谢谢您的回答,我是新手,没有金币和积分奖励,不好意思。
还有一个问题就是我按照您的建议修改了,结果 ...

因d mask是以等值线的形式把mask里的边界绘制出来了。建议绘制这种边界之前可以用meteinfo把需要的地图转换成shp格式,GrADS2.0.a8以上版本是支持打开shp的,这样绘制地图边界简便一点。具体可以参考http://bbs.06climate.com/forum.php?mod=viewthread&tid=2263
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-4-24 11:36:16 | 显示全部楼层
balfulosa 发表于 2018-4-24 10:08
因d mask是以等值线的形式把mask里的边界绘制出来了。建议绘制这种边界之前可以用meteinfo把需要的地图转 ...

明白了,原来是版本的事。谢谢您。还有一个问题是ERA数据里有一个Scale_factor和Add_offset,这个是不是得对比湿等数据进行校正,不能用Grads的sdfopen函数或者arcgis的“创建NetCDF栅格图层”直接用?
您能否加一下我qq:270990857,张智华,以后请您多多指教,谢谢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-4-24 19:50:31 | 显示全部楼层
Zhangzh 发表于 2018-4-24 11:36
明白了,原来是版本的事。谢谢您。还有一个问题是ERA数据里有一个Scale_factor和Add_offset,这个是不是 ...

您的第二个问题不是很清楚,但ERA_Interim确实有很多变量在GrADS里不能直接打开,有些是因为变量名称里包含有_所以打不开。如果是这个问题的话可以参考家园关于xdfopen用法的介绍,重新给不能直接打开的变量添加ctl,这样有些问题可以解决。
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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