| 
 
	积分9124贡献 精华在线时间 小时注册时间2013-7-14最后登录1970-1-1 
 | 
 
| 
研一菜鸟一枚,研究方向要处理MLS数据,刚刚接处不久,目前用NCL处理,现在想要用以下脚本读取一个MLS数据并做出卫星轨迹图,代码如下:
x
登录后查看更多精彩内容~您需要 登录 才可以下载或查看,没有帐号?立即注册 
  ;*****************************************
 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下。
 万分感谢!
 
 
 
 
 
 | 
 |