爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4789|回复: 5

[脚本编辑] grads算海温距平

[复制链接]

新浪微博达人勋

发表于 2017-12-28 09:40:15 | 显示全部楼层 |阅读模式

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

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

x
用grads求6月海温距平图所得结果如下图:
但是我明明设置缺省值了啊,为什么还会这样!求大神帮解答。
其中算距平的gs文件:
'reinit'
'sdfopen E:\statistic_work\fortran\eof\sst.mon.mean.nc'
'set fwrite E:\statistic_work\fortran\eof\sstjp.grd'
'set gxout fwrite'
'define asst=ave(sst,t=726,t=1506,12)'
tt=726
while(tt<=1506)
'set x 120 290'
'set y 63 150'
'set t 'tt
'set z 1'
'nsst=sst-asst'
'd nsst'
tt=tt+12
endwhile
'disable fwrite'
'reinit'
;

相应的ctl文件:
dset E:\statistic_work\fortran\eof\sstjp.grd
title sstjp
undef 1e+20
xdef 171 linear 119.5 1.0
ydef 88 linear -27.5 1.0
zdef 1 levels 0
tdef 66 linear 00z01jun1951 1yr
vars 1
nsst 0 99
endvars

我随便画了一年的,无论第几年的都像是没有考虑缺省值的样子。

所画出的海温距平图

所画出的海温距平图
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-12-28 10:07:37 | 显示全部楼层
undef  -9.99e8
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-12-28 10:56:14 | 显示全部楼层

超级感谢!!!终于可以了!!可是,我下的nc的数据显示它的缺测值是1e+20啊,不是-9.99e8,为什要改成这样?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-12-28 11:13:26 | 显示全部楼层
LIfighting 发表于 2017-12-28 10:56
超级感谢!!!终于可以了!!可是,我下的nc的数据显示它的缺测值是1e+20啊,不是-9.99e8,为什要改成这 ...

使用grads计算后输出的数据都是默认用这个缺省值,记住即可。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-12-28 11:40:25 | 显示全部楼层
啊啊啊啊!这样啊!谢谢啦{:5_213:}
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-12-28 11:40:55 | 显示全部楼层
lqouc 发表于 2017-12-28 11:13
使用grads计算后输出的数据都是默认用这个缺省值,记住即可。

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

本版积分规则

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

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

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