- 积分
- 638
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2017-12-8
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
请问各位用gsn_csm_vector_scalar函数画风廓线图的时候有出现过纵坐标画不出的情况吗?我的纵坐标是自己设置的,如果使用默认设置就可以画出,但是用explicit就不行了,请问各位有遇到过这种问题吗?绘图部分的设置如下:
res = True
res@gsnDraw = False ; do not draw
res@gsnFrame = False ; do not advance the frame
res@tiMainString = "Wind Profile" ; add title
res@vcGlyphStyle = "WindBarb" ; select wind barbs
res@vcRefLengthF = 0.045 ; define length of vec ref
res@vcRefAnnoOn = False ; turn off ref wind barb
res@trYReverse = True ; reverse the Y-axis
res@vcMonoWindBarbColor = False ; color barbs by scalar
res@lbOrientation = "vertical" ; vertical labelbar
; res@vcLevelSelectionMode = "ExplicitLevels" ; will explicitly set the levels
res@vcWindBarbLineThicknessF= 1.5 ; set the wind barb thickness
; res@vcLevels = ispan(5,65,5) ; set the levels
; res@vcLevelColors = (/20,30,40,45,50,55,60,65,70,75,80,85,90,95/) ; set the colors associated with the wind speed
res@tmXTOn = False
res@tmYROn = False
res@tmYLOn = True
res@tmYLMode = "Explicit"
; res@tmYLValues = (/50.,100.,150.,200.,250.,300.,350.,400.,450.,500.,550.,600.,650.,700.,750.,800.,850.,900.,950.,1000./)
res@tmYLValues = ispan(50,1000,50)
res@tmYLLabels = ""+res@tmYLValues
; res@trYMinF = 50
; res@trYMaxF = 1000
res@tiYAxisString = "Pressure(hPa)"
res@tiXAxisString = "wind(A-B)(m/s)"
res@tmXBMode = "Explicit"
res@tmXBValues = ispan(0,3,1)
res@tmXBLabels = (/"00Z","06Z","12Z","18Z"/)
res@vcMapDirection = False ; needed to disassociate the map from the vectors
res@gsnAddCyclic = False ; data is not cyclic
res@gsnLeftString = ""
plot = gsn_csm_vector_scalar(wks,u,v,spd,res)
getvalues plot ; get the X/Y axis min/max for use in the loglin plot
"trXMinF" : trxmin
"trXMaxF" : trxmax
"trYMinF" : trymin
"trYMaxF" : trymax
end getvalues
loglin = create "logling" logLinPlotClass wks ; draw a loglin plot, with expanded X/Y axis
"trXMinF" : trxmin-1
"trXMaxF" : trxmax+1
"trYMinF" : trymin
"trYMaxF" : trymax+1
"trYReverse" : True
"vpXF" : .15 ; set the X-axis NDC starting point
"vpYF" : .8 ; set the Y-axis NDC starting point
"vpWidthF" : .7 ; set the width of the plot in NDC units
"vpHeightF" : .45 ; set the height of the plot in NDC units
end create
overlay(loglin,plot) ; overlay plot with the loglin plot
draw(loglin) ; draw the plot
frame(wks) ; advance the frame
|
|