NCL处理气溶胶AOD数据的程序

kongfeng0824

NCL的官方网站的CALIPSO版块的数据不是CAPLIPSO官方网站的HDF数据，而是经过处理的nc数据。
我的这个程序修改后可以直接对CALIPSO的AOD数据进行读取画图。
对于研究气溶胶和云的同学应该会有帮助。


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 "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/skewt_func.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/wind_rose.ncl"
begin
fi = addfile("D:\dropbox\calipso data\MYD04_L2.A2008031.0800.051.2009109172932.hdf","r")
print(fi)
data =new(dimsizes(fi->Optical_Depth_Land_And_Ocean),float,-9999)
data = (fi->Optical_Depth_Land_And_Ocean)* 0.001000000047497451 + 0
printMinMax(data,False)
lat = fi->Latitude
lat@units ="degrees_north"
lon = fi->Longitude
lon@units ="degrees_east"
data@lat2d = lat
data@lon2d = lon
wks = gsn_open_wks("pdf" ,"Optical_Depth_Land_And_Ocean")
gsn_define_colormap(wks,"Rainbow")
res                        = True
  res@mpFillOn              = False
  res@cnFillOn              = True
  res@cnFillMode            = "RasterFill"
  res@cnLinesOn             = False
  res@cnLineLabelsOn        = False
  res@cnLevelSelectionMode  = "ManualLevels"
  res@cnMinLevelValF        = 0
  res@cnMaxLevelValF        = 1.5
  res@cnLevelSpacingF       = 0.1
  res@gsnSpreadColors       = True
  res@lbLabelBarOn          = True
  res@lbLabelAutoStride     = True
  res@lbBoxLinesOn          = False

plot  = gsn_csm_contour_map(wks,data,res)
end

另外附上数据信息、程序代码、样图。

又是那隻貓

给力，谢谢分享

liutaoz21

感谢楼主分享！

ayzqs

太给力了！！！！！

加多了个宝

这是什么东西啊？在哪个平台上运行的啊。。

kongfeng0824

加多了个宝 发表于 2013-8-3 23:44

                               
登录/注册后可看大图

这是什么东西啊？在哪个平台上运行的啊。。

NCL版块当然是NCL下的程序了。

加多了个宝

kongfeng0824 发表于 2013-8-4 04:05

                               
登录/注册后可看大图

NCL版块当然是NCL下的程序了。

好吧，从来没有接触过，，嘿嘿。。。有下载链接么

qq79242943

好东西  学习了

guoli15

谢谢分享。。。

小小阿包

太好了，赶紧去试试。