- 积分
- 3167
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-10-5
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 oyy. 于 2021-1-8 20:57 编辑
很简单的作图,ncl出的图 在南非缺少一块,而且多了一条线。但是用ncview看数据 又很正常。
有没有大佬提供思路?
脚本如下:
f1 = addfile("fig1_all.nc","r")
com2 = f1->com2
com2@_FillValue = 9.96921e+36
;;;===============================================================================
name = "/figure"
wks= gsn_open_wks("png",name )
res = True ; plot mods desired
res@gsnDraw = False
res@gsnFrame = False
res@gsnAddCyclic = True
;res@gsnMaximize = True
res@mpGeophysicalLineColor= "grey60"
;res@mpMinLonF = -180
;res@mpMaxLonF = 180
res@mpMinLatF = -40
res@mpMaxLatF = 60
res@mpCenterLonF = 0.
res@tmXBLabelDeltaF = -0.6 ;使其更靠近坐标轴
res@tmYLLabelDeltaF = -0.6
res@tmXBMode = "Explicit"
res@tmXBTickSpacingF = 60
res@tmXBMinorOn = True
res@tmXBMinorValues = ispan(-360,360,20)
res@gsnLeftString = ""
res@gsnRightString = ""
res@gsnLeftStringFontHeightF = 0.02
res@gsnRightStringFontHeightF = 0.02
res@mpFillOn = False
res@cnFillOn = True
res@cnLinesOn = False
;***************************************
res@cnLineLabelsOn = False
res@cnInfoLabelOn = False
res@tmXBLabelFontHeightF = 0.017 ;字体大小
res@tmYLLabelFontHeightF = 0.017
res@lbLabelBarOn = True
res@pmLabelBarWidthF = 0.07 ; 色标宽度default is shorter
res@pmLabelBarHeightF = 0.22 ; 高度default is taller
res@pmLabelBarOrthogonalPosF = 0.05 ;调整色标垂直方向距离
;res@pmLabelBarParallelPosF = -0.03 ;0.03 vor ;调整色标水平方向距离
res@lbTitleOn = True ; turn on title
res@lbTitleString = "" ;"mm day~S~-1"
res@lbTitlePosition = "Right" ; title position
res@lbTitleDirection = "Across" ; title direction
res@lbTitleFontHeightF = 0.018 ; make title smaller
res@lbOrientation = "Vertical"
res@cnFillPalette = "CBR_wet" ;色标
res@cnLevelSelectionMode = "ExplicitLevels"
res@cnLevels = 0.01*ispan(0,12,2)
res@gsnLeftString = ""
base= gsn_csm_contour_map_ce(wks,com2,res)
draw(base)
frame(wks)
数据已上传,有兴趣的可以试一下! |
-
ncl出的图
-
ncview
-
-
try.nc
566.61 KB, 下载次数: 1, 下载积分: 金钱 -5
|