爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 13004|回复: 18

[其他] NCL处理MLS卫星数据读取出错,望有大神指教~

[复制链接]

新浪微博达人勋

发表于 2014-3-18 23:01:26 | 显示全部楼层 |阅读模式

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

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

x
研一菜鸟一枚,研究方向要处理MLS数据,刚刚接处不久,目前用NCL处理,现在想要用以下脚本读取一个MLS数据并做出卫星轨迹图,代码如下:
;*****************************************
diri= "./"
fili= "MLS-Aura_L2GP-O3_v03-43-c01_2013d328.he5"

pltType ="ps"
pltName ="hdfeos_3"

;*******************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "HDFEOS_LIB.ncl"
;*******************************************
;   MAIN
;*******************************************

begin
f = addfile(diri+fili,"r")
print(f)

;read some data
print(names)
lat=f->Latitude_O3
lon=f->Longitude_O3
lev=f->Pressure_O3
time=f->Time_O3
ntim=dimsizes(time)

T=f->L2gpValue_O3

printVarSummary(T)
printMinMax(T,True)
print("===============")

telapse=(time-time(0))/60
telapse@long_name="Elapsed Time(minutes)"
telapse@units="minutes since"+time(0)

;***********************************
;plot
;***********************************
wks=gsn_open_wks(pltType,pltName)
gsn_define_colormap(wks,"amwg")

dumb=NhlNewColor(wks,0.7,0.7,0.7)
res=True
res@gsnMaximize=True
res@gsnSpreadColors=True
res@gsnSpreadColorEnd=-2
;res@gsnPaperOrientation="portrait"

res@tiMainString=fili

res@cnFillOn=True
res@cnLinesOn=False
res@cnLinesLabelsOn=False
res@cnFillMode="RasterFill"
res@cnRasterSmoothingOn=True

res@lbLabelAutoStride=True
res@lbOrientation="Vertical"

res@cnLevelSelectionMode="ManualLevels"
res@cnMinLevelValF=160
res@cnMaxLevelValF=260
res@cnLevelSpacingF=10

res@trYReverse=True
T&nLevels_O3=lev
T&nTimes_O3=telapse

res@gsnYAxisIrregular2Linear=True
res@tmYLFormat="f"
res@tiYAxisString="Pressure(hPa)"
res@gsnLeftString=T@Title

plot=gsn_csm_contour(wks,O(nLevels_O3|:,nTimes_O3|:),res)

end

运行后提示错误是:An internal error has occured. The file format requested does not appea to be supported,could not open(./MLS-Aura_L2GP-O3_v03-43-c01_2013d328.he5)
并且提示错误可能发生在21行,也就是从读取数据开始的:lat=f->Latitude_O3。。。。。。

请问哪位用NCL处理过MLS数据的大神帮忙一下,为什么不支持数据格式文件打不开,难道是我的NCL安装有问题?PS:我的NCL是在Linus fedora下。
万分感谢!




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

新浪微博达人勋

发表于 2014-3-19 08:22:02 | 显示全部楼层

回帖奖励 +1 金钱

“The file format requested does not appea to be supported”
换个最新版本的NCL试试
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-3-19 09:20:57 | 显示全部楼层
MLS是关于降水的数据吗?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-3-19 09:23:07 | 显示全部楼层
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-3-19 09:38:28 | 显示全部楼层
The file format requested does not appea to be supported

打开文件出错,先检查文件格式,看看是否下载完整
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-3-19 13:44:21 | 显示全部楼层
hzwjy 发表于 2014-3-19 08:22
“The file format requested does not appea to be supported”
换个最新版本的NCL试试

嗯,实在不行我就更换一下NCL版本,但是对于安装我不太了解,不要请求别人帮忙,忧桑~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-3-19 13:46:15 | 显示全部楼层
kongfeng0824 发表于 2014-3-19 09:20
MLS是关于降水的数据吗?

MLS是NASA的Aura卫星上的临边探测探测器,没有降水的,主要是温度,臭氧,水汽等。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-3-19 13:46:43 | 显示全部楼层
freekiller 发表于 2014-3-19 09:38
The file format requested does not appea to be supported

打开文件出错,先检查文件格式,看看是否下 ...

嗯嗯,我再更换一下别的日期的数据,谢谢了~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-3-19 14:51:56 | 显示全部楼层
NCL版本至少要高于5.2.0,能用最新的建议用最新的。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-3-19 19:39:55 | 显示全部楼层
hzwjy 发表于 2014-3-19 14:51
NCL版本至少要高于5.2.0,能用最新的建议用最新的。

哦哦,原来处理MLS还有这一讲究,谢谢了!我不清楚我的是哪个版本的,是师姐帮着装的,太囧了~
请问最新的版本已经到哪个了呀?
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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