爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 7624|回复: 12

grads写文件画不出图

[复制链接]
发表于 2014-8-30 14:27:42 | 显示全部楼层 |阅读模式
GrADS
系统平台:
问题截图: -
问题概况: 用nc资料写出500hpa1957至2008年逐年4-6月平均位势场dat文件,为什么写出来的文件画不出图,显示是都是常数0.下面是脚本gs和ctl文件
'reinit'
'sdfopen e:\0001\hgt.mon.mean.nc'
'set gxout fwrite'
'set fwrite e:\0001\anomal\hgt.all.dat'

'set lev 500'
'set lon 100 150'
'set lat 0 30'

i=112
while(i<=726)
'set t 'i''
'define hgtave=(hgt(t='i')+hgt(t='i+1')+hgt(t='i+2'))/3'
'd hgtave'
i=i+12
endwhile


'disable fwrite'
*end


ctl文件:
dset e:\0001\anomal\hgt.all.dat
title monthly longterm mean hgt from the NCEP reanalysis
undef -9.96921e+36
xdef 21 linear 100 2.5
ydef 13 linear 0 2.5
zdef 1 levels 850
tdef 52 linear 00z01Apr1957 1yr
vars 1
hgtave 1 99 monthly longterm mean of hgt
endvars
我看过提问的智慧: 看过
自己思考时长(天): 2

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

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

x

密码修改失败请联系微信:mofangbao
 楼主| 发表于 2014-8-30 15:02:19 | 显示全部楼层
求回复
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

发表于 2014-8-30 15:52:46 | 显示全部楼层
hgt.mon.mean.nc 你这个资料是哪年到哪年的?
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2014-8-30 15:57:21 | 显示全部楼层
Eegle 发表于 2014-8-30 15:52
hgt.mon.mean.nc 你这个资料是哪年到哪年的?

1948年1月到2013年11月逐月数据
密码修改失败请联系微信:mofangbao
发表于 2014-8-30 16:26:28 | 显示全部楼层
是不是缺省值有问题~~
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2014-8-30 16:39:42 | 显示全部楼层
zhouling 发表于 2014-8-30 16:26
是不是缺省值有问题~~

写出来的文件就有问题,只有56k,很奇怪,按理说很多年的数据应该很大才对
密码修改失败请联系微信:mofangbao
发表于 2014-8-30 17:17:31 | 显示全部楼层
缺测值改成-9.99e+8试试
密码修改失败请联系微信:mofangbao
发表于 2014-8-30 22:02:22 | 显示全部楼层
'set lev 850'  看看  
gs里面和ctl里面的不一致。
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2014-8-30 23:23:45 | 显示全部楼层
苏涛 发表于 2014-8-30 22:02
'set lev 850'  看看  
gs里面和ctl里面的不一致。

改了,写出来的dat文件大小还是56k,好奇怪
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2014-8-30 23:31:09 | 显示全部楼层
river 发表于 2014-8-30 17:17
缺测值改成-9.99e+8试试

改过还是一样,设置一个时间d出来显示是:constant field.value=0
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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