[ncl-talk] Incorrect label bar
Adam Phillips
asphilli at ucar.edu
Fri Sep 23 13:01:18 MDT 2016
Hi Beata,
Try changing this:
res at cnLevels = (/1,2,3,4/)
to this:
res at cnLevels = (/1.5,2.5,3.5,4.5/)
If that does not work please send the current version of you script, along
with your aridity colormap and data files to our ftp site so someone can
take a look. Follow the ftp instructions here:
http://www.ncl.ucar.edu/report_bug.shtml
(Under "You can ftp larger datasets and files to:")
When the transfer is completed make sure you send to ncl-talk the name of
the files you transferred.
Adam
On Fri, Sep 23, 2016 at 1:56 AM, Beata Szabo <szabo.b at czechglobe.cz> wrote:
> Dear Mary,
>
> Thank you for your response. If I apply your suggestions:
>
> res at cnLevelSelectionMode = "ExplicitLevels"
> res at cnLevels = (/1,2,3,4/)
>
> the labelbar will be correct but the maps are incorrect, because this
> method shifts the colors. I attached its result. It does not show the
> hyper-arid areas (brown) and the semi-arid levels (green) are missed in
> Iberia-Peninsula. I need same maps as I attached in the previous email.
> Could you or someone help me to fix this matter.
>
> Thank you in advance!
> Beata
>
> > Beata,
> >
> > When NCL repeats colors in a labelbar like this, it's because there
> aren't
> > enough colors in the color map for the number of contour levels being
> > requested.
> >
> > It looks like your color map only has 5 colors, so you need to use one
> > fewer contour levels:
> >
> > res at cnLevelSelectionMode = "ExplicitLevels"
> > res at cnLevels = (/1,2,3,4/)
> >
> >
> > --Mary
> >
> >
> >
> >
> > On Thu, Sep 22, 2016 at 3:06 AM, Beata Szabo <szabo.b at czechglobe.cz>
> > wrote:
> >
> >> I added
> >> resP at lbBoxCount = 5
> >>
> >> but unfortunately it does not help. I created new netcdf files with same
> >> variables, but I converted the unep variable to integer and I also tried
> >> to draw the maps with those integer variables (1,2,3,4,5) with
> >>
> >> res at cnLevels = (/1,2,3,4,5/)
> >>
> >> The resulted maps are correct but the label is still incorrect both with
> >> string and without string labels. I attached the sript and the created
> >> maps. I would appreciated if someone helped me!
> >>
> >> Beata
> >>
> >> > I don't know for certain, but you might try adding
> >> >
> >> > resP at lbBoxCount = 5
> >> >
> >> > I'm not positive that will fix it, but the docs say that the sizes of
> >> > arrays like lbLabelStrings, etc, must match the lbBoxCount, and in
> >> your
> >> > case it sounds like its defaulting to 11 (?)
> >> >
> >> > Hope that helps...
> >> > Rick
> >> >
> >> >
> >> > On Wed, Sep 21, 2016 at 2:12 AM, Beata Szabo <szabo.b at czechglobe.cz>
> >> > wrote:
> >> >
> >> >> I created a panel plot with my own netCDF files and color table be
> >> >> created
> >> >> according to
> >> >> http://www.ncl.ucar.edu/Document/Graphics/create_color_table.shtml
> >> >>
> >> >> My own color table has ncolors = 5. The netCDF files contains five
> >> float
> >> >> values (1.0, 2.0, 3.0, 4.0, 5.0). After I run the below script I
> >> >> received
> >> >> correct maps with incorrect labelbar. The labelbar containes 11 box
> >> >> where
> >> >> the first two and the last two colors are two times and the third
> >> color
> >> >> roles three times. Moreover the label strings are the follows: Arid,
> >> >> Sub-humid, Label_6, Label_8.
> >> >>
> >> >> How can I create a correct labelbar?
> >> >>
> >> >> ;-----------------------------------------------------------
> -----------
> >> >> ; Main code.
> >> >> ;-----------------------------------------------------------
> -----------
> >> >> begin
> >> >> ;************************************************
> >> >> ; read in netCDF file
> >> >> ;************************************************
> >> >> a = addfile("UNEP_E-OBS_sm.nc","r")
> >> >> b = addfile("UNEP_1971-2000_sm.nc","r")
> >> >> c = addfile("UNEP_2021-2050_sm.nc","r")
> >> >> d = addfile("UNEP_2071-2100_sm.nc","r")
> >> >>
> >> >> un1 = a->unep(:,:)
> >> >> un2 = b->unep(:,:)
> >> >> un3 = c->unep(:,:)
> >> >> un4 = d->unep(:,:)
> >> >>
> >> >> minlat = 30 ; min lat to mask
> >> >> maxlat = 75 ; max lat to mask
> >> >> minlon = -20 ; min lon to mask
> >> >> maxlon = 40 ; max lon to mask
> >> >>
> >> >> ;************************************************
> >> >> ; create plot
> >> >> ;************************************************
> >> >> wks = gsn_open_wks("png","unep") ; send graphics to PNG
> >> >> file
> >> >>
> >> >> cmap = read_colormap_file("aridity")
> >> >>
> >> >>
> >> >>
> >> >> plot = new(4,graphic)
> >> >> res = True ; plot mods desired
> >> >> res at gsnMaximize = True ; enlarge plot
> >> >> res at gsnDraw = False ; Don't draw yet
> >> >> res at gsnFrame = False ; Don't advance
> frame
> >> >> yet
> >> >> res at mpProjection = "LambertConformal"; choose projection
> >> >>
> >> >> res at cnFillOn = True ; turn on color
> >> >> res at cnLinesOn = False ; turn off contour
> >> lines
> >> >> res at lbLabelBarOn = False
> >> >> res at cnFillPalette = cmap(::-1,:) ; set color map
> >> >>
> >> >> res at mpMinLatF = minlat
> >> >> res at mpMaxLatF = maxlat
> >> >>
> >> >> res at mpMinLonF = minlon
> >> >> res at mpMaxLonF = maxlon
> >> >>
> >> >> res at gsnMaskLambertConformal = True ; turn on lc masking
> >> >> res at mpGridAndLimbOn = True
> >> >> res at mpGridLatSpacingF = 10
> >> >> res at mpGridLonSpacingF = 5
> >> >> res at gsnAddCyclic = False
> >> >>
> >> >> res at tiMainString = ""
> >> >> res at tiMainOffsetYF = 0.01 ; Move title up a little
> >> >> res at gsnRightString = "A"
> >> >> res at gsnLeftString = ""
> >> >> ; pr&lon = pr&lon-180 ; make lon go -180
> >> to
> >> >> 180
> >> >>
> >> >> plot(0) =
> >> >> gsn_csm_contour_map(wks,un1({minlat:maxlat},{minlon:maxlon}),res);
> >> >> create plot
> >> >>
> >> >> resb = True ; plot mods desired
> >> >> resb at gsnMaximize = True ; enlarge plot
> >> >> resb at gsnDraw = False ; Don't draw yet
> >> >> resb at gsnFrame = False ; Don't advance
> >> frame
> >> >> yet
> >> >> resb at mpProjection = "LambertConformal"; choose projection
> >> >>
> >> >> resb at cnFillOn = True ; turn on color
> >> >> resb at cnLinesOn = False ; turn off contour
> >> >> lines
> >> >> resb at lbLabelBarOn = False
> >> >> resb at cnFillPalette = cmap(::-1,:) ; set color map
> >> >>
> >> >> resb at mpMinLatF = minlat
> >> >> resb at mpMaxLatF = maxlat
> >> >>
> >> >> resb at mpMinLonF = minlon
> >> >> resb at mpMaxLonF = maxlon
> >> >>
> >> >> resb at gsnMaskLambertConformal = True ; turn on lc
> masking
> >> >> resb at mpGridAndLimbOn = True
> >> >> resb at mpGridLatSpacingF = 10
> >> >> resb at mpGridLonSpacingF = 5
> >> >> resb at gsnAddCyclic = False
> >> >>
> >> >> resb at tiMainString = ""
> >> >> resb at tiMainOffsetYF = 0.01 ; Move title up a little
> >> >> resb at gsnRightString = "B"
> >> >> resb at gsnLeftString = ""
> >> >>
> >> >> plot(1) =
> >> >> gsn_csm_contour_map(wks,un2({minlat:maxlat},{minlon:maxlon}),resb);
> >> >> create plot
> >> >>
> >> >> resc = True ; plot mods desired
> >> >> resc at gsnMaximize = True ; enlarge plot
> >> >> resc at gsnDraw = False ; Don't draw yet
> >> >> resc at gsnFrame = False ; Don't advance
> >> frame
> >> >> yet
> >> >> resc at mpProjection = "LambertConformal"; choose projection
> >> >> resc at cnFillOn = True ; turn on color
> >> >> resc at cnLinesOn = False ; turn off contour
> >> >> lines
> >> >> resc at lbLabelBarOn = False
> >> >> resc at cnFillPalette = cmap(::-1,:) ; set color map
> >> >>
> >> >> resc at mpMinLatF = minlat
> >> >> resc at mpMaxLatF = maxlat
> >> >>
> >> >> resc at mpMinLonF = minlon
> >> >> resc at mpMaxLonF = maxlon
> >> >>
> >> >> resc at gsnMaskLambertConformal = True ; turn on lc
> masking
> >> >> resc at mpGridAndLimbOn = True
> >> >> resc at mpGridLatSpacingF = 10
> >> >> resc at mpGridLonSpacingF = 5
> >> >> resc at gsnAddCyclic = False
> >> >>
> >> >> resc at tiMainString = ""
> >> >> resc at tiMainOffsetYF = 0.01 ; Move title up a little
> >> >> resc at gsnRightString = "C"
> >> >> resc at gsnLeftString = ""
> >> >>
> >> >>
> >> >> plot(2) =
> >> >> gsn_csm_contour_map(wks,un3({minlat:maxlat},{minlon:maxlon}),resc);
> >> >> create plot
> >> >>
> >> >>
> >> >> resd = True ; plot mods desired
> >> >> resd at gsnMaximize = True ; enlarge plot
> >> >> resd at gsnDraw = False ; Don't draw yet
> >> >> resd at gsnFrame = False ; Don't advance
> >> frame
> >> >> yet
> >> >> resd at mpProjection = "LambertConformal"; choose projection
> >> >>
> >> >> resd at cnFillOn = True ; turn on color
> >> >> resd at cnLinesOn = False ; turn off contour
> >> >> lines
> >> >> resd at lbLabelBarOn = False
> >> >> resd at cnFillPalette = cmap(::-1,:) ; set color map
> >> >>
> >> >> resd at mpMinLatF = minlat
> >> >> resd at mpMaxLatF = maxlat
> >> >>
> >> >> resd at mpMinLonF = minlon
> >> >> resd at mpMaxLonF = maxlon
> >> >>
> >> >> resd at gsnMaskLambertConformal = True ; turn on lc
> masking
> >> >> resd at mpGridAndLimbOn = True
> >> >> resd at mpGridLatSpacingF = 10
> >> >> resd at mpGridLonSpacingF = 5
> >> >> resd at gsnAddCyclic = False
> >> >>
> >> >> resd at tiMainString = ""
> >> >> resd at tiMainOffsetYF = 0.01 ; Move title up a little
> >> >> resd at gsnRightString = "D"
> >> >> resd at gsnLeftString = ""
> >> >>
> >> >>
> >> >> plot(3) =
> >> >> gsn_csm_contour_map(wks,un4({minlat:maxlat},{minlon:maxlon}),resd);
> >> >> create plot
> >> >>
> >> >> ;---Attach latitude labels
> >> >> add_lc_labels(wks,plot(0),minlat,maxlat,minlon,maxlon)
> >> >> add_lc_labels(wks,plot(1),minlat,maxlat,minlon,maxlon)
> >> >> add_lc_labels(wks,plot(2),minlat,maxlat,minlon,maxlon)
> >> >> add_lc_labels(wks,plot(3),minlat,maxlat,minlon,maxlon)
> >> >>
> >> >> ;---Drawing the plot will also draw all the attached labels.
> >> >> draw(plot(0))
> >> >> draw(plot(1))
> >> >> draw(plot(2))
> >> >> draw(plot(3))
> >> >> frame(wks)
> >> >> ;-- create panel plot
> >> >> ;************************************************
> >> >> resP = True ; modify the panel
> >> >> plot
> >> >>
> >> >> resP at gsnPanelLabelBar = True ; add common
> >> colorbar
> >> >> resP at lbLabelFontHeightF = 0.007 ; make labels
> >> smaller
> >> >> resP at cnExplicitLabelBarLabelOn = True
> >> >>
> >> >> resP at lbLabelStrings = (/"","Hyper-arid","Arid", \
> >> >> "Semi-arid","Sub-humid","Humid"/)
> >> >>
> >> >> gsn_panel(wks,plot,(/2,2/),resP)
> >> >>
> >> >> end
> >> >>
> >> >> Thank you for your help in advance!
> >> >>
> >> >> Beata
> >> >>
> >> >>
> >> >> Dr. Beata Szabo-Takacs
> >> >> Global Change Research Institute CAS
> >> >> Bělidla 986/4a
> >> >> 60300 Brno
> >> >> Czech Republic
> >> >>
> >> >> _______________________________________________
> >> >> ncl-talk mailing list
> >> >> ncl-talk at ucar.edu
> >> >> List instructions, subscriber options, unsubscribe:
> >> >> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >> >>
> >> >
> >>
> >>
> >> Dr. Beata Szabo-Takacs
> >> Global Change Research Institute CAS
> >> Bělidla 986/4a
> >> 60300 Brno
> >> Czech Republic
> >>
> >> _______________________________________________
> >> ncl-talk mailing list
> >> ncl-talk at ucar.edu
> >> List instructions, subscriber options, unsubscribe:
> >> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >>
> >>
> >
>
>
> Dr. Beata Szabo-Takacs
> Global Change Research Institute CAS
> Bělidla 986/4a
> 60300 Brno
> Czech Republic
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
--
Adam Phillips
Associate Scientist, Climate and Global Dynamics Laboratory, NCAR
www.cgd.ucar.edu/staff/asphilli/ 303-497-1726
<http://www.cgd.ucar.edu/staff/asphilli>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160923/c3c87637/attachment.html
More information about the ncl-talk
mailing list