爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 5535|回复: 4

TRMM资料画降水粒子季度平均的垂直分布图的gs文件 画出来一直是缺测值

[复制链接]
发表于 2014-4-16 18:20:02 | 显示全部楼层 |阅读模式
GrADS
系统平台: Grads
问题截图:
问题概况: 利用TRMM 3A12卫星资料画降水粒子季度平均的垂直分布图,一直出不来图,画出的一直是一条值为0的线,求指教。
我看过提问的智慧: 看过
自己思考时长(天): 2

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

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

x
程序如下
'reinit'
'open e:\Elvin\3a12\3a12spring.ctl'
'open e:\Elvin\3a12\3a12summer.ctl'
'open e:\Elvin\3a12\3a12autumn.ctl'
'open e:\Elvin\3a12\3a12winter.ctl'
'enable print e:\Elvin\3a12\stp.composite-ns.gmf'
'set grads off'
'set grid off'
'set lat 0'
'set lon 115'
'set z 1 28'
'set cstyle 1'
'set ccolor 1'
'set cthick 8'
'set cmark 1'
'd ave(cldIce.1,t=1,t=30)'

'set lat 0'
'set lon 115'
'set z 1 28'
'set cstyle 2'
'set ccolor 2'
'set cthick 8'
'set cmark 2'
'd ave(cldIce.2,t=1,t=30)'

'set lat 0'
'set lon 115'
'set z 1 28'
'set cstyle 3'
'set ccolor 3'
'set cthick 8'
'set cmark 3'
'd ave(cldIce.3,t=1,t=30)'

'set lat 0'
'set lon 115'
'set z 1 28'
'set cstyle 4'
'set ccolor 4'
'set cthick 8'
'set cmark 4'
'd ave(cldIce.4,t=1,t=30)'

'print'
'disable print'
'reinit'
求大神指教

密码修改失败请联系微信:mofangbao
发表于 2014-4-16 21:06:33 | 显示全部楼层
你那些ctl 可以正确描述你的资料吗?那些ctl能画出来图吗?那些ctl描述的资料是什么样的?
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2014-4-17 12:17:07 | 显示全部楼层
dset E:\Elvin\3a12\autumn\3A12.%y4%m201.7.HDF
dtype hdfsds d
options template
undef -9999.9
XDEF 720 linear -179.875 0.5
ydef 160 linear -39.875  0.5
zdef 28 levels 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 8.5 9 9.5 10 11 12 13 14 15 16 17 18
tdef 30 linear 00Z01Sep2001 1mo

vars 14
surfacePrecipitation=>sp 0 x,y precipitation rate mm/hr
surfaceRain=>sr 0 x,y rain rate mm/hr
convectPrecipitation=>cp 0 x,y convect precipitation rate mm/hr
cldWater=>cldWater 28 z,x,y cloud liquid water content g/m3
rainWater=>rainWater 28 z,x,y rain water content g/m3
cldIce=>cldIce 28 z,x,y cloud ice liquid water content g/m3
snow=>snow 28 z,x,y snow liquid water content g/m3
graupel=>graupel 28 z,x,y graupel liquid water content g/m3
latentHeat=>latentHeat 28 z,x,y laten heating release C/hr
npixTotal=>nt 0 x,y pixels with pixelStatus equal to zero
npixPrecipitation=>np 0 x,y pixels with surfacePrecipitation greater than 0
fractionQuality0=>f0 0 x,y pixels with qualityFlag equal to 0(high quality) %
fractionQuality1=>f1 0 x,y pixels with qualityFlag equal to 1(medium quality) %
fractionQuality2=>f2 0 x,y pixels with qualityFlag equal to 2(low quality) %
endvars

以autumn为例的ctl文件,画变量如sp、cp都可以画出来,cldWater等就画不出来了
密码修改失败请联系微信:mofangbao
发表于 2014-4-24 11:17:53 | 显示全部楼层
垂直间隔不对吧,10到18km的是隔4km,6到10km间隔2km,反正2A12是这样的,不知道3A12是不是
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2014-4-24 18:50:31 | 显示全部楼层
已出图,是数据本身的问题。
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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