- 积分
- 544
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-4-17
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
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"
pfile=addfile("/data4/TRMM/2011/3B42.20110604.18.7.HDF","r")
prec1=pfile->precipitation
prec=prec1(:,::-1) ;
lon=fspan(-179.875,179.875,1440)
lat=fspan(49.875,-49.875,400)
lon@units="degrees-east"
lat@units="degrees-north"
lon@long_name="lon"
lat@long_name="lat"
lon@axis= "X"
lat@axis= "Y"
prec!0="lon"
prec!1="lat"
prec&lon=lon
prec&lat=lat
prec@_FillValue = -999
prec@missing_value= -999
res =True
res@cnLinesOn = False
res@cnLineLabelsOn = False
res@cnFillOn = True
res@gsnAddCyclic = False
res@mpMinLatF = -10.0
res@mpMaxLatF = 45.0
res@mpMinLonF = 60.
res@mpMaxLonF = 170.0
res@lbLabelBarOn = True
res@gsnSpreadColors = True ;use all range color
res@cnLevelSelectionMode = "ExplicitLevels" ; set explicit contour levels
res@cnLevels = (/0.,1.,2.,3.,4.,5./) ; set levels
res@cnFillColors = (/0,2,4,6,9,13,15/) ; set the colors to be used
prec2=prec(lat|:,lon|:)
wks=gsn_open_wks("X11","pre")
gsn_define_colormap(wks,"BlWhRe")
plot=gsn_csm_contour_map(wks,prec2,res)
|
-
|