Hi-
I am having some serious trouble drawing a polyline using gsn_add_polyline. The map appears but not the line. There is no error, just a map with no line.
; show map with the cross section line
resa = True
WRF_map_c(a, resa, 0) ; reads info from file
resa_at_tfDoNDCOverlay = True
resa_at_gsnMaximize = True ; maximize size
resa_at_gsnSpreadColors = True ; use full range of colormap
resa_at_gsnScalarContour = True ; contours desired
resa_at_cnFillMode = "RasterFill" ; raster
resa_at_cnLevelSelectionMode = "ExplicitLevels" ; explicit [unequal] cn levels
resa_at_vcGlyphStyle = "WindBarb"
resa_at_vcRefLengthF = 0.025 ; ref vec length
resa_at_vcMinDistanceF = 0.02 ; larger means sparser
resa_at_vcWindBarbTickLengthF = 0.3 ; default 0.3
resa_at_vcRefAnnoOn = False
qlevs = (/0,1,2,3/)
resa_at_mpGeophysicalLineThicknessF = 2.0
resa_at_mpUSStateLineThicknessF = 2.0
resa_at_gsnLeftString = "Wind Vectors (m/s), Temp, Hght"
resa_at_gsnRightString = "MYJ"
resa_at_gsnDraw = False
resa_at_gsnFrame = False
print("plotting")
plota = gsn_csm_map(wks,resa)
resp = True
resp_at_gsLineColor = "Black"
resp_at_gsLineThicknessF = 2.0
resp_at_gsLineLabelString = "A"
xp = (/plane(0),plane(2)/)
yp = (/plane(1),plane(3)/)
print("dum")
dum = gsn_add_polyline(wks,plota,xp,yp,resp)
print("frame")
draw(plota)
frame(wks)
James Correia Jr
Post Doc
PH: 372-6463
james.correia_at_pnl.gov
"Wisdom. Strength. Courage. Generosity. Each of us are born with one of these. We must find the other three inside of us."
-from "Into the West"
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed May 13 2009 - 02:45:32 MDT
This archive was generated by hypermail 2.2.0 : Wed May 13 2009 - 08:40:16 MDT