- 积分
- 9684
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-7-14
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
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下。
万分感谢!
|
|