- 积分
- 164
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2021-3-10
- 最后登录
- 1970-1-1
|
10金钱
资料的经度都是-180到180这样的,我想画地图里东经100-西经100那个范围的相关系数,但好像我数据画的的是西经100-东经100,如何调整过来?还有好像显著性打点也没加上去。
脚本如下:
h2=h1(0:36,5,{200},{0:90},{100:-100})
ccr=escorc_n(o3m1,h2,0,0)
copy_VarCoords(h2(0,:,:), ccr(:,:))
wks= gsn_open_wks("png","h2")
res = True ; make plot mods
res@cnFillOn = True ; turn on color
res@cnFillPalette = "BlWhRe" ; set color map
res@gsnDraw = False ; do not draw the plot
res@gsnFrame = False
res@cnLinesOn = False ; turn off contour lines
res@gsnAddCyclic=False
res@cnLevelSelectionMode = "ManualLevels" ; manually set cn levels
res@gsnAddCyclic=False
res@cnMinLevelValF = -0.7 ; min level
res@cnMaxLevelValF = 0.7 ; max level
res@cnLevelSpacingF = .1 ; contour level spacin
res@mpMinLatF =0
res@mpMaxLatF = 90.
res@mpMinLonF = 100.
res@mpMaxLonF = 300.
res@pmTickMarkDisplayMode="Always"
res@cnRasterSmoothingOn = True
res@mpFillDrawOrder="PostDraw"
plot = gsn_csm_contour_map(wks,ccr,res)
res@pmTickMarkDisplayMode="Always"
res@cnRasterSmoothingOn = True
res@mpFillDrawOrder="PostDraw"
res@lbLabelBarOn = False
sres = True
sres@gsnDraw = False ; do not draw the plot
sres@gsnFrame = False ; set up a second resource list
sres@cnLineLabelsOn = False ; do not use line labels
sres@cnLinesOn = False ; do not draw contour lines
sres@cnInfoLabelOn = False
sres@cnFillOn = True ; color fill
sres@lbLabelBarOn= False
sres@cnLevelSelectionMode = "ExplicitLevels"
sres@cnFillDrawOrder="postDraw"
sres@cnLevels = (/-0.325,0.325/)
sres@cnMonoFillPattern = False ; want multiple patterns
sres@cnFillPatterns = (/17,-1,17/) ; the patterns
sres@cnMonoFillScale = False ; want different densities
sres@cnFillScales = (/0.6,0.6,0.6/) ; change densities
sres@cnMonoFillColor =True
sres@cnFillDotSizeF = 0.0008
plot1 = gsn_csm_contour_map(wks,ccr,sres)
overlay(plot,plot1)
draw(plot)
frame(wks)
end
|
-
|