Hi David,
There's no automatic way to add lat/lon labels to a masked LC plot.
I created an example showing how you might do this using gsn_add_text. Please see example 4 at:
http://www.ncl.ucar.edu/Applications/lcmask.shtml
--Mary
On Jan 23, 2011, at 11:07 PM, David Small wrote:
> Hi everyone,
>
> I'm plotting vectors/scalar and then laying two other contours over with the overlay function. The masked lambert conformal map doesn't display labels for the latitude and longitude lines. Can anyone tell me how to turn them on?
>
> Thanks.
>
> David Small
>
> wks = gsn_open_wks("eps" ,"TROP."+year(j)+"."+month(j)+"."+day(j)+"."+hour(j)) ; open a ps file
> gsn_define_colormap(wks,"ViBlGrWhYeOrRe") ; choose colormap
>
> res2 = True
> res2@mpProjection = proj
> res2@gsnMaskLambertConformal = True
> res2@cnFillOn = True ; color fill
> res2@cnLinesOn = False ; no contour lines
> res2@gsnSpreadColors = True ; use total colormap
> res2@mpGridAndLimbOn = True
> res2@pmTickMarkDisplayMode = "Always" ; turn on tickmarks
>
> res2@gsnAddCyclic = True ; regional data
> res2@mpMinLatF = lato ; range to zoom in on
> res2@mpMaxLatF = latf
> res2@mpMinLonF = lono
> res2@mpMaxLonF = lonf
> res2@mpCenterLonF = lonc
> res2@lbOrientation = "horizontal"
> res2@gsnMaximize = True ; must include w/ Paper Orientation
>
>
> res2@gsnScalarContour = True ; vectors over scalars
>
> res2@cnLevelSelectionMode = "ManualLevels"
> res2@cnMinLevelValF = 260 ; min lev
> res2@cnMaxLevelValF = 380 ; max lev
> res2@cnLevelSpacingF = 5
> res2@gsnMaximize = True ; enlarge plot
> res2@cnLineLabelsOn = False
> res2@gsnStringFont = "times-roman"
> res2@mpFillOn = True
> res2@mpOutlineBoundarySets = "Geophysical"
> res2@mpDataResolution = "Fine"
> res2@mpFillOn = False
>
> res2@vcRefMagnitudeF = 10.0 ; define vector ref mag
> res2@vcRefLengthF = 0.020 ; define length of vec ref
> ; res2@vcGlyphStyle = "CurlyVector" ; turn on curly vectors
> res2@vcGlyphStyle = "WindBarb" ; turn on curly vectors
> res2@vcMinDistanceF = 0.02 ; thin vectors
> res2@vcRefAnnoOrthogonalPosF = .1 ; move ref vector down
> res2@mpDataBaseVersion = "MediumRes"
> res2@mpGeophysicalLineThicknessF = gthick
> res2@vcRefAnnoOn = False
> ; res2@vcMinDistanceF = 0.025
> res2@gsnSpreadColorStart = 2 ; Don't start at index 2.
>
>
> res2@gsnDraw = False ; don't draw
> res2@gsnFrame = False ; don't advance frame
> res2@cnLineDashPattern = 12 ; dashed lines
> res2@cnLineThicknesses = 1.0
> res2@cnLineLabelFontHeightF = .0045
>
> res2@cnInfoLabelFontHeightF = 0.0
>
> res2@mpFillOn = False
> res2@mpGeophysicalLineThicknessF = gthick
>
> res2@mpGridLineThicknessF = 0.5
>
>
> res2@gsnLeftString = ""
> res2@gsnRightString = nme
>
> dres@gsnLeftString = ""
> dres@gsnRightString = ""
>
> res2@lbLabelAutoStride = True
> res2@tmLabelAutoStride = True
> res2@tmXTOn = True
> res2@tmYROn = True
>
> res2@mpOutlineDrawOrder = "PostDraw"
>
> res2@tmXTLabelsOn = True
> res2@tmYRLabelsOn = True
> res2@tmXBLabelFont="times-roman"
>
> plot = gsn_csm_vector_scalar_map(wks,u(j,:,:),v(j,:,:),f(j,:,:),res2)
> plot3 = gsn_csm_contour(wks,VORT850(j,:,:),vres)
> overlay(plot,plot3)
> plot2 = gsn_csm_contour(wks,dv(j,:,:),dres)
> overlay(plot,plot2)
>
> draw(plot)
> frame(wks)
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Jan 24 12:41:36 2011
This archive was generated by hypermail 2.1.8 : Tue Jan 25 2011 - 14:22:15 MST