Re: [ferret_users] Fill and

To: William Kessler - NOAA Federal <william.s.kessler@xxxxxxxx>, ferret_users@xxxxxxxx

Subject: Re: [ferret_users] Fill and

From: Serena Illig <serena.illig.ird@xxxxxxxxx>

Date: Tue, 21 Jun 2022 18:25:18 +0200

Arc-authentication-results: i=4; mx.google.com; dkim=pass header.i=@gmail.com header.s=20210112 header.b=CkNwd8B6; arc=pass (i=2 dkim=pass dkdomain=gmail.com); spf=pass (google.com: domain of serena.illig.ird@xxxxxxxxx designates 209.85.219.71 as permitted sender) smtp.mailfrom=serena.illig.ird@xxxxxxxxx

Arc-authentication-results: i=3; mx.google.com; dkim=pass header.i=@gmail.com header.s=20210112 header.b=CkNwd8B6; arc=pass (i=2 dkim=pass dkdomain=gmail.com); spf=pass (google.com: domain of serena.illig.ird@xxxxxxxxx designates 209.85.219.71 as permitted sender) smtp.mailfrom=serena.illig.ird@xxxxxxxxx

Arc-authentication-results: i=2; mx.google.com; dkim=pass header.i=@gmail.com header.s=20210112 header.b=CkNwd8B6; arc=pass (i=1 spf=pass spfdomain=gmail.com dkim=pass dkdomain=gmail.com dmarc=pass fromdomain=gmail.com); spf=softfail (google.com: domain of transitioning serena.illig.ird@xxxxxxxxx does not designate 65.153.203.174 as permitted sender) smtp.mailfrom=serena.illig.ird@xxxxxxxxx

Arc-authentication-results: i=1; mx.google.com; dkim=pass header.i=@gmail.com header.s=20210112 header.b=CkNwd8B6; spf=pass (google.com: domain of serena.illig.ird@xxxxxxxxx designates 209.85.220.41 as permitted sender) smtp.mailfrom=serena.illig.ird@xxxxxxxxx; dmarc=pass (p=NONE sp=QUARANTINE dis=NONE) header.from=gmail.com

Arc-message-signature: i=4; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=list-archive:list-help:list-post:list-id:mailing-list:precedence:to :subject:message-id:date:from:in-reply-to:references:mime-version :sender:dkim-signature; bh=vjMgRcmdZpFjpQxWSinZStxqLn8A+MMVx3w9uFwXh2E=; b=o8Xoo6K9aYQrjwgo61MzSyjEDlRzX4Z1i0li8xlGMoMZZuOz/1OlbVxyoBsqfeFXbn Pb9MZkhFqDgLI5dwWVJVWfSaeqkdYkxrsqrGyp4oVGcTWitRwh9o54x3o5FFgwvYMaEh 1cr94SYeLF5dcOgfMAMdk86Mbaa6zCyuKOmthxHaHVfwZZLxsOfNieilSeqKJJWl3smN u7R3yB7MToGlf0W/xokWvdOtWTePES70fNOss+X4udvo5lCbq1Ouc2CPLNq7Tf9JX/Wy uxDOcxhXSqXDKs8soY/G3vqbn2XRsxLG39RPygmML3A/GZwqvZHJhWVcahh4HDRyUop5 WZgw==

Arc-message-signature: i=3; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=to:subject:message-id:date:from:in-reply-to:references:mime-version :dkim-signature; bh=vjMgRcmdZpFjpQxWSinZStxqLn8A+MMVx3w9uFwXh2E=; b=yV+yLd+7tTrM0fVl4c7njFWAIRdhjOmNtSl1lBvjY0VTkrmraAYOnTUmYbSUBpZM4k igb1iwAeHUwSthOI41mysU7UdRewGRpdR6PjXrEnsFUuP8Rtva7DKbVYw9aWCn769Xh1 InjDxtYYeQkR/eRwnWiJb5M4JLDNdT77K37He1no/x/DQM/kLpcJI19gAYkTbJZw+ULt S+1II1pzSJNe3SvEwd06OluKHtMPzel8VY5Y4bfr4TpSuDMTtH4m0ec8sC1lOSJwMimY 6C0pZoPPROHrRmdA6HrkGuTMGNww9dBPx+mwl/0F2zx3ZImOnyzYFkuxb22QGevKFq9Q o5Wg==

Arc-message-signature: i=2; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=to:subject:message-id:date:from:in-reply-to:references:mime-version :dkim-signature; bh=vjMgRcmdZpFjpQxWSinZStxqLn8A+MMVx3w9uFwXh2E=; b=JQ780BgzGSl9d0mY1puPkDFbbHl2TytvRJaMof+t1NFMqK8tITmgRWm+8WyQ/wkTCK 0ybLAqQ+dqnIZnCO2YNrsxjOU77aA2F1XsJXsDAzR0IDRdpceSziHMzjOYwJ2aU9Jk4a SRcacTsLobhZP2AdPBw+/P4djmPD4yAsESPDHIlMvRQCr5pRcsi/5mHUQ/+WWHbZCZh1 MgrvNisZMwZ5/v+SZ4JStxP8oWReXlfGwkLesg3XOyNEXushazJpcgNXll42aQaf7vMg EblVGqzxgy3gv0PgxEwNAsiDkEni23s4AMSCYIcmCLpbGacqUR18qmWZ1Q10kll9LiyV HMlg==

Arc-message-signature: i=1; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=to:subject:message-id:date:from:in-reply-to:references:mime-version :dkim-signature; bh=vjMgRcmdZpFjpQxWSinZStxqLn8A+MMVx3w9uFwXh2E=; b=j3/R/WF7R1l/d6ER9hzJBY7Iz4t9WKLlckbEPgXN4s9QzN+OqMjyCtMETPUiXSBSY5 2ENz8Sj0UEsRlrb95Na+shH+1izQ3o2kQEvuFBNswQIfSc+9IQ8lphU9Xdm5cAJ9jC+q lPFYDPbZwpjvIkTiXUiL1qwYtUPzZU3le27t1L7Tfob6116k76p4EP/0ubF8qs1McW7r /HMpcafO6p4UKrUj4YRzijSuJwr5mJ7g+zeGR0h5vaPNW5A114GRlSqCdWlTE/vOj7CZ EVHZxmRnRlNox6pq7RMoDuqfFY/Y0apb1bY+z0t8NKXWgzKXU2C+epi6oJ6sUD0i+BkJ rjCg==

Arc-seal: i=4; a=rsa-sha256; t=1655828737; cv=pass; d=google.com; s=arc-20160816; b=n99k15D/aNUUvKp5WuIjZQQwOuHVvvc+APBHaAX0+j5HCPnwYOTyXocZsCOdx8msN0 sofWQYIzbq97DBfJU1SI1aQLh74Bev1HxnZWoSAkOj/hqi5Fd0xiB14mfRg/Uj5nkUCt teHCTSXTIZT5P0Mytpd6UWVBkXucT5nIvaWXxlO8xi0/fVMD1nqmUELjT2ArzVjpMknC onV5K+xUogFfcZfIfGXYgG7FEG7p2Nwm8uxjZbMTmB0aEpmRU0oCsnLhGXIb98Qy8Btr ImIrtaIujWyg+av4dOFCeRg2i2YYE/zSpMs3SH7pjMs1ENglFO6bpYNWBydEqXXDuIn7 2ZKA==

Arc-seal: i=3; a=rsa-sha256; t=1655828735; cv=pass; d=google.com; s=arc-20160816; b=a5cEsa226Z3lSaKasskf4fkb/7nP00pBK3bILbGSXm0qbIBvbGUQTBrJXVXC5V56w2 yIa3uTAaYt9srwqG2UpFDp37cD0yKZHpiHW+d1KKAh2A1fNAt4h96oh1r4Rdp2WnhYGJ GAOHX/7+dsXcKxpJgd9JjPa4sPxgtTcUuvXHzASecD+O78TGJNsw7igktaG8grHoo+5U JjY9nnP+Uu0KNXKI0+K5o/5iNGcoHKSWUhY5BoPcWBXErsPIAh+5Od3r3LlUZ3UPNi+p 2JfFNOVnklrygEMLUE3kFWxsBGKR72faNHsrcw5STdQKwpLBkF8VI0YCDb1xaAnaMRH2 Ijyg==

Arc-seal: i=2; a=rsa-sha256; t=1655828732; cv=pass; d=google.com; s=arc-20160816; b=Ul5bJ8BKGXTWWAGtoICY1cHnZpwKgaMCM2cVh5+XZbnS/D8xP1Vct3iwr4pZ3BJUO7 f2+/EV/0Lveillrb0/p8XI0y8pm60Xp41hbWQvrEKaQPoxWZlB1YWArN0pKoeJ1Add1i qla3ThqDFhMOLjACPAK4liPO5BLVbBZP94sg4LWujrsGC+2YN7weYbtEVzE37QATLeSk JbBCq1muY6kM/95vy5r7K9A7ks/U0UyyZTSh7CbKGpracW5umRdNVhJO5UxXq3eEQ3MX 5lzcy5AQrcWu6Xfvlu+Yp0sV6hAMTfQcq+Y7/OAys59dJVmq9bsx36cgRQUN6xvahKiI OD2g==

Arc-seal: i=1; a=rsa-sha256; t=1655828730; cv=none; d=google.com; s=arc-20160816; b=GqtCTHgarM6EAXuM117ljAu6OVQfGtEz830epZq89z0i9yz/lUSQgTdP1hkoJ7Jgxc aPbmrgtakgY1fitNNwun/yPUFrLRrcEHLqCn4v1/oG1xkkjqyOmnUNkstXAW4/U1mcol 7ysU53RZy5OiuBxfifaWGzdwSk4lvFq5FK0lqNGMr2rgyJ0gJrCdKbBTxO9JROC7MX1F riecWkcSgmcAhvYuToUn1ggJRNo57DXTo2K+2AqIqsQGTWkSNuiVur1ThNEYLwuJKc2z y3qjbEO7C0iUYCgkuhmuxppIjOkXjEX+HF6jWhl34Wgf6A6D3gwWfXXW/Ewqy/7wwNrd zUOw==

Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=noaa.gov; s=google; h=sender:mime-version:references:in-reply-to:from:date:message-id :subject:to:x-original-sender:x-original-authentication-results :precedence:mailing-list:list-id:list-post:list-help:list-archive; bh=vjMgRcmdZpFjpQxWSinZStxqLn8A+MMVx3w9uFwXh2E=; b=kcHdfDuGwpVwFv0pXylrciggy+gBbMMSfsFUo1gSrTNBlQMEEHByH78pzytJVQvj/Y JfZoLk4Bogxsf9QXidwGpp9tBIRv+Ho4X8kInnWcPpYbQA+xQlOgHHtOWTlcNZKXXpVj GEkB054YyXEnjo37UDHOtlQW3bjICE11e2GLdDF4vcuetfS3lTvP90ziKiyG/mdNS3Wn Wew7VEtz+oF0wCTI6I7NzDsSQfIso1SxOTBWxB0AD59vv6FCuRopGUJ1OrNS/7iXgZ66 CbIPWcUJWd96hJSKfeRbQC9KEZdttJ4Or5rlD4uP69wfor5LHym9Cg1/+OeBC9NeOHZM xFSw==

In-reply-to: <CAK0Qe2PSt2_zORv5HrTkEatu+TsmVS=pCtsgTMKPUQ57HMENFw@mail.gmail.com>

List-archive: <https://groups.google.com/a/noaa.gov/group/ferret_users/>

List-help: <https://support.google.com/a/noaa.gov/bin/topic.py?topic=25838>, <mailto:ferret_users+help@noaa.gov>

List-id: <ferret_users.noaa.gov>

List-post: <https://groups.google.com/a/noaa.gov/group/ferret_users/post>, <mailto:ferret_users@noaa.gov>

Mailing-list: list ferret_users@xxxxxxxx; contact ferret_users+owners@xxxxxxxx

References: <CADczVkNmfe1Qo-piJn+R5=6Q2h2=8eRxiHfFavpDUTcYbjbaeA@mail.gmail.com> <CAK0Qe2PSt2_zORv5HrTkEatu+TsmVS=pCtsgTMKPUQ57HMENFw@mail.gmail.com>

Sender: owner-ferret_users@xxxxxxxx

Dear William,

Thanks for your answer. I definitely agree with you that Ferret is doing a wonderful job for contours, and plotting in general (especially time).

Concerning the contour, I did try CONSET at one point, but I did not play with SHAKEY because I only plot contours.

PPL CONSET did not solve my problem:

use coads_climatology.cdf
contour/x=50w:20e/y=30s:30n/lev=(20,23,1)/set SST[l=@ave]

ppl CONSET 0.08 1 1 0.04 0.04 0.0
ppl contour

Vijith indicated me a nice solution based on the ndigit inputs in /LEVELS=(lo, hi, delta, ndigits)

"ndigits (applies to CONTOUR command only) is the number of decimal places to use when labeling the level on individual contour lines as: -1 for integer format or -3 to omit numerical labels"

It works perfectly.

For the FILL/SHADE commands, I discovered it when I plotted some symmetric color scale for wind stress anomalies, which were not symmetric because of the missing last level... To obtain the right plot, I just added an extra level in the /LEVELS command, such as:

shade/x=50w:20e/y=30s:30n/levels=(0.014,0.030,0.001) SST[l=@ave]/1000

plots a color scale starting at 0.014 and ending at 0.029.

I do not understand why some specifications work and some do not. I will pay extra attention to my color scale, just in case the top-level disappears.

Many thanks,

Serena

.¸. , . .·´¯`·
><(((º>    `·.¸.·´¯`·...¸><(((º>¸.¸. , . .·´¯`· ><(((º> .
`·.¸.·´¯`·...¸><(((º>        .
.
Serena ILLIG-THEVENIN LEGOS/IRD       .

Editorial Board Member for Scientific Reports .

14 Av. E. Belin .

31400 Toulouse .
FRANCE .

E-mail : serena.illig.ird@xxxxxxxxx .

Web-site : http://sillig.free.fr .

Phone : +33 5 61 33 29 26 .

.
<º)))><..·

On Tue, Jun 21, 2022 at 5:57 PM William Kessler - NOAA Federal <william.s.kessler@xxxxxxxx> wrote:

Hi Serena (also responding to your earlier message about contour levels).

First, contouring, while conceptually straightforward, is a hard thing to program. Figuring out where to put contour labels without obscuring other contours or overwriting other labels is again something that a human would do easily, but hard to specify in an algorithm. How many labels is "enough"? In my experience, Ferret does a better job of this than some others.

The two commands you need to control this are PPL CONSET and PPL SHAKEY; that is, you need to give your CONTOUR/SHADE/FILL command with /SET, then specify the two PPL commands, then draw the plot with PPL CONTOUR/SHADE/FILL.

PPL CONSET,hgt,nsig,narc,dashln,spacln,cay,nrng,dslab ! controls the contouring and labeling
PPL SHAKEY [0=nokey,1=key] [0=hor,1=vert] lab-size lab-inc lab-dig lab-len xlo xhi ylo yhi ! controls the color key

These are described in the "Enhanced commands" section, Appendix C 2.2 and C 2.4
https://ferret.pmel.noaa.gov/Ferret/documentation/users-guide/appendix-c-ferret-specific-pplus-enhancements/ENHANCED-COMMANDS-DESCRIPTION#_VPID_767
(Some of CONSET can also be controlled by qualifiers to CONTOUR, or by its /LEVELS options; see CONTOUR)

With respect to occasional missing end values in the colorbar, I think this may be the result of rounding somewhere in the calculation (???). You can often work around this by tweaking the upper or lower bound by a tiny amount. E.g.,

If the original command fails:
shade/x=50w:20e/y=30s:30n/levels=(0.014,0.029,0.001) SST[l=@ave]/1000
Try:
shade/x=50w:20e/y=30s:30n/levels=(0.013999,0.029,0.001) SST[l=@ave]/1000
or:
shade/x=50w:20e/y=30s:30n/levels=(0.014,0.029001,0.001) SST[l=@ave]/1000

I don't know why these sometimes succeed when the apparently correct specification fails!

Billy

On Tue, Jun 21, 2022 at 4:52 AM Serena Illig <serena.illig.ird@xxxxxxxxx> wrote:
Dear Ferret users,

I hope this email finds you well.

I am writing regarding some contour labels. I am using PyFerret v7.63 on Linux Ubuntu (via miniconda) but I have the same issues with PyFerret v7.62 and ferret v6.71.

I noticed a problem in the number of colors associated with the FILL or SHADE commands. When I use SHADE and I increase the number of levels, there are some cases for which the last level range is not drawn and is missing in the color scale. It happens when the quantity drawn is small (for instance wind stress in N/m2).

Here are 3 examples:
use coads_climatology.cdf
shade/x=50w:20e/y=30s:30n/levels=(14,29,1) SST[l=@ave]
---> 15 color levels, with levels 0(1), the colorscale is OK
shade/x=50w:20e/y=30s:30n/levels=(0.014,0.029,0.001) SST[l=@ave]/1000
---> 15 color levels, with levels 0(1/1000), the colorscale is missing the last level
shade/x=50w:20e/y=30s:30n/levels=(0.014,0.028,0.001) SST[l=@ave]/1000
---> 14 colors levels with levels 0(1/1000), the colorscale is OK
fill/x=50w:20e/y=30s:30n/levels=(0.0014,0.0029,0.0001) SST[l=@ave]/10000
---> 15 colors levels with levels 0(1/10000), the colorscale is missing the last level

I could not find any rules for which it happens. I thought at first that it is linked to the number of levels, but I found counter-examples (see my last example above).

Maybe this issue has already been reported?

Have a nice day,

Serena Illig

.¸. , . .·´¯`·
><(((º>    `·.¸.·´¯`·...¸><(((º>¸.¸. , . .·´¯`· ><(((º> .
`·.¸.·´¯`·...¸><(((º>        .
.
Serena ILLIG-THEVENIN LEGOS/IRD       .
                                                                              .
Editorial Board Member for Scientific Reports    .
.
                         14 Av. E. Belin          .
31400 Toulouse     .
FRANCE       .
                                           .
   E-mail : serena.illig.ird@xxxxxxxxx .
       Web-site : http://sillig.free.fr         .
   Phone : +33 5 61 33 29 26         .
                                                               .
                                           <º)))><..·

Follow-Ups:

Re: [ferret_users] Fill and
- From: William Kessler - NOAA Federal

References:

[ferret_users] Fill and
- From: Serena Illig
Re: [ferret_users] Fill and
- From: William Kessler - NOAA Federal