爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4529|回复: 7

[图形美化] 新人求助!GrADS将hgt.nc转为grd文件有问题

[复制链接]

新浪微博达人勋

发表于 2015-4-24 16:51:39 | 显示全部楼层 |阅读模式

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

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

x
   如题,我在http://www.esrl.noaa.gov下载的高度场月平均资料,然后用GrADS转成100hpa的高度场grd文件,GS内容如下
'reinit'
'set gxout fwrite'
'sdfopen E:\NCEP\pressurelevel\hgt.mon.mean.nc'
'set fwrite D:\bylw\zhishu\s\hgt.grd'
i=1
while(i<=807)
'set lev 100'
'set  t 'i''
'd hgt'
i=i+1
endwhile
'disable fwrite'
;
之后用这个文件计算南亚高压面积指数,发现指数很奇怪,就写了CTL文件在GrADS里画hgt.grd的图,发现几乎每个时次的图都是这样
ceshi.png
就像这样 下面有一条很粗的线 像是把上面的等值线都顶没了
求大神帮忙看看这到底是什么原因。
新人在这谢谢大家了!

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

新浪微博达人勋

 楼主| 发表于 2015-4-24 16:52:42 | 显示全部楼层
人工置顶~吧里大神们帮帮我
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-4-24 20:44:09 | 显示全部楼层
把set x 和set y 设置加上看看
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-4-24 20:44:28 | 显示全部楼层
把set x 和set y 设置加上看看
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-4-24 21:11:12 | 显示全部楼层
顶楼上的建议
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-4-24 21:11:43 | 显示全部楼层
(1)在不设置经纬度范围时,全球是按照默认
set lon 0 360
set lat -89.5 89.5
这一范围写出数据,0和360是同一个点数据写入了两次,这种情况下,hgt.ctl中应为 xdef 145 linear 0 2.5。
(2)建议将经纬度范围加上,
set lon 0 357.5
set lat -89.5 89.5
此种写法hgt.ctl中应为 xdef 144 linear 0 2.5
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-4-24 21:36:11 | 显示全部楼层
dln 发表于 2015-4-24 21:11
(1)在不设置经纬度范围时,全球是按照默认
set lon 0 360
set lat -89.5 89.5

万分感谢!用你的办法已经解决了问题 终于可以做后面的内容了 谢谢^ ^
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-4-24 21:37:33 | 显示全部楼层

谢谢 谢谢!
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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