Duplicate contour line labels with overlay shading

Good day,

I am getting duplicate contour line labels after overlaying a
stippling pattern for contours above a certain value of equivalent
potential temperature (eth).

As you can see from the excerpted script below, I reset the resource
(res_eth) contour line and contour line label to "False" before
calling "gsn_contour_shade". However, that does not stop the
duplication.

Your help with solving this problem will be greatly appreciated.

Thanks.
Arlene Laing
NCAR/MMM

;***********************************************
; set plot resources
   res_at_vpWidthF = .9 ; overwrite basic plot size
   res_at_vpHeightF = 1.0
   pltres = True

         res_xy = res
         res_xy_at_tiYAxisString = "Pressure (hPa)"
         res_xy_at_AspectRatio = 0.75
         res_xy_at_cnMissingValPerimOn = True
         res_xy_at_cnMissingValFillColor = 0
         res_xy_at_cnMissingValFillPattern = 11
         res_xy_at_sfYCStartV = max(pressure_levels)
         res_xy_at_sfYCEndV = min(pressure_levels)
         res_xy_at_trYReverse = True ; reverse Y axis
         res_xy_at_trYLog = True ; Log scale
         res_xy_at_tmYLMode = "Explicit"
         res_xy_at_tmYLValues = (/1000., 850., 700., 600., 500., 400., \
                               300., 250., 200./)
         res_xy_at_tmYLLabels = (/"1000","850","700","600","500","400", \
                               "300","250","200"/)
; res_xy_at_sfYArray = pressure_levels
         res_xy_at_tiXAxisFontHeightF = 0.02
         res_xy_at_tiYAxisFontHeightF = 0.02
         res_xy_at_tmXBMajorLengthF = 0.02
         res_xy_at_tmYLMajorLengthF = 0.02
         res_xy_at_tmYLLabelFontHeightF = 0.015
         res_xytxFontHeightF = 0.015
         res_xy_at_tiXAxisString = "Latitude"
         res_xy_at_sfXArray = lat2d(:,0)

         res_xy_at_gsnDraw = False ; Don't draw plots
         res_xy_at_gsnFrame = False ; Don't advance frame

        ; Plotting options for U
         res_u = res_xy
; res_u_at_cnInfoLabelOrthogonalPosF = 0.00
         res_u_at_cnInfoLabelOn = False ; turn off cn info label
         res_u_at_cnLineLabelPerimOn = False ; turn off label perimeter
         res_u_at_cnLabelMasking = True ; do not draw labels over contour
         res_u_at_cnLineLabelFontHeightF = 0.016
         res_u_at_cnLineLabelBackgroundColor = "white"
         res_u_at_cnLineThicknessF = 1.5
         res_u_at_ContourParameters = (/ -30., 20, 2. /)
         res_u_at_gsnContourNegLineDashPattern = 1

        ; Plotting options for RH
         res_rh = res_xy
         res_rh_at_pmLabelBarOrthogonalPosF = -0.07
         res_rh_at_ContourParameters = (/ 10., 90., 10. /)
         res_at_cnFillMode = "RasterFill" ; Raster Mode
         res_at_cnLinesOn = False ; Turn off contour lines
         res_at_cnLineLabelsOn = False ; Turn off contour lines
         res_rh_at_cnFillOn = True
         res_rh_at_cnFillColors = (/"White", "White","White",\
                                             "White","DarkOliveGreen1","Chartreuse","Green",
\
                                             "Green3", \
                                             "ForestGreen","PaleGreen4"/)

        ; Plotting options for ETH, Equivalent Pot Temperature
         res_eth = res_xy
         res_eth_at_cnInfoLabelOrthogonalPosF = 0.00
         res_eth_at_cnLineColor = "Blue"
         res_eth_at_cnLineThicknessF = 1.0
         res_eth_at_cnLineLabelFontHeightF = 0.020
         res_eth_at_cnLabelMasking = True ; do not draw labels over
contour
         res_eth_at_ContourParameters = (/ 320., 376., 2. /)

        ; Get the contour info (from WRF NCL example)
          contour_u = wrf_contour(a,wks,zave_u,res_u)
          contour_rh = wrf_contour(a,wks,zave_rh,res_rh)
          contour_eth = wrf_contour(a,wks,zave_eth,res_eth)

        ; MAKE PLOTS
         pltres_at_FramePlot = False ; do not frame plot - will do this
later

         plot =
wrf_overlays(a,wks,(/contour_eth,contour_u,contour_rh/),pltres)

         res_at_gsnRightString = ""
         res_at_gsnLeftString = ""
         res_at_gsnCenterString = ""

; testing shaded contour for eth > 343K
         res_eth_at_cnLinesOn = False ; Turn off contour lines
         res_eth_at_cnLineLabelsOn = False ; Turn off contour labels
         res_eth_at_cnMonoFillColor = True ; default color is fground
[black]
         res_eth_at_cnMonoFillPattern = False ; want multiple patterns
         res_eth_at_cnMonoFillScale = False ; want patterns w/ diff
densities
         res_eth_at_gsnShadeFillType = "Pattern"
         res_eth_at_gsnShadeHigh = 17 ; Use fill pattern #17 (stipple)

; plot=gsn_contour(wks,zave_eth,res_eth) ; plot zonal ave of eth
         plot = gsn_contour_shade(plot,-999.,343., res_eth) ; Pattern
shade contours above 343

  draw(plot)
  frame(wks)

;***********************************************
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Sep 02 2009 - 10:47:42 MDT