- 积分
- 1885
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2018-10-7
- 最后登录
- 1970-1-1
|
发表于 2019-4-3 22:29:07
|
显示全部楼层
这是我的程序
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
setfileoption("nc","MissingToFillValue",False)
div=addfile("198812/198812div.nc","r")
dd=short2flt(div->d)
d=dd(27:30,25,:,:)
delete(div)
delete(dd)
wks=gsn_open_wks("pdf","700hPa_div1228-31")
gsn_define_colormap( wks ,"BlWhRe")
cnres=True
cnres@gsnDraw = False
cnres@gsnFrame = False
cnres@gsnAddCyclic=False
cnres@mpLimitMode="LatLon"
cnres@mpMinLatF=18
cnres@mpMaxLatF=50
cnres@mpMinLonF=75
cnres@mpMaxLonF=120
cnres@mpDataBaseVersion="Ncarg4_1"
cnres@mpDataSetName="Earth..4"
;cnres@mpOutlineBoundarySets="National"
cnres@mpOutlineSpecifiers=(/"China:states"/)
cnres@pmTickMarkDisplayMode="Always"
cnres@cnLevelSelectionMode = "ManualLevels"
cnres@cnMinLevelValF = -0.00006
cnres@cnMaxLevelValF = 0.00006
cnres@cnLevelSpacingF = 0.00001
cnres@cnFillOn = True
cnres@cnLinesOn = False
cnres@cnInfoLabelOn = False
cnres@gsnSpreadColors = True
cnres@lbLabelAutoStride = True
plot1=new(4,"graphic")
plot1(0)=gsn_csm_contour_map(wks,d(0,:,:),cnres)
plot1(1)=gsn_csm_contour_map(wks,d(1,:,:),cnres)
plot1(2)=gsn_csm_contour_map(wks,d(2,:,:),cnres)
plot1(3)=gsn_csm_contour_map(wks,d(3,:,:),cnres)
lnres=True
lnres@gsLineColor = "white"
lnres@gsLineThicknessF = 6.0
map = gsn_add_shapefile_polylines(wks,plot1,"home/Qinghai-Tibet_Plateau.shp",lnres)
resPanel = True
resPanel@gsnPanelLabelBar= True
gsn_panel(wks,plot1,(/2,2/),resPanel)
|
|