[ferret_users] Re: [ferret_users] Re: [ferret_users] [ferret

[ferret_users] Re: [ferret_users] Re: [ferret_users] [ferret_users] 20° E seems to be magic

To: baris onol <onolba@xxxxxxxxxx>

Subject: [ferret_users] Re: [ferret_users] Re: [ferret_users] [ferret_users] 20° E seems to be magic

From: Nicolas Freychet <nicolas.freychet@xxxxxxxxx>

Date: Fri, 27 Jan 2017 12:24:36 +0000

Cc: ferret_users <ferret_users@xxxxxxxx>

Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=noaa-gov.20150623.gappssmtp.com; s=20150623; h=mime-version:in-reply-to:references:from:date:message-id:subject:to :cc:x-original-sender:x-original-authentication-results:precedence :mailing-list:list-id:list-post:list-help:list-archive :list-unsubscribe; bh=Yjqu3Q1H9Si4P+gEy0r+SnAb50nfCZJJBu44jBDvY4o=; b=EozgVKpMmBdWBqhjRMaJdHbiCT6OoUxptnNUULiS+nni/YRfWGiZ0L7izH3cop3uPP ANToptDNsQXkJtIHPGVtvz/NAsfNGxg60O1RdBkgedoAYfAsFH9TGbH7ftCvFtlx22FE YKJpw4XsQ1B59EJljDIZGaVyc6TRKOMPMab6KIrCjO0t1Y1ShBV0NLnkEVEryQJA6iRD yJWPaoptYGk+LbFQlQvXVJIZq8o9NB9pLGTMWmRsmDBCD2Lgm13qxQUzk5EVGL34bavy wmk79LZ0ziOJFFfqrpowgZaa+5vVJNUXCtXQO4EqZmEGuI3pX6Dorc2zMB5hWPzo7ZxF 0BaQ==

In-reply-to: <5aa0b3c1-dcdd-22e4-324c-3077e6684f2b@itu.edu.tr>

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>

List-unsubscribe: <mailto:googlegroups-manage+809428289204+unsubscribe@googlegroups.com>, <https://groups.google.com/a/noaa.gov/group/ferret_users/subscribe>

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

References: <CAOs80Nw+8NnyzCaUrSefk5OZQ2f40VMiHbs78H+KmcwZK-=A1g@mail.gmail.com> <28cb516a-5e29-1afb-a94d-d4df96cfef9a@io-warnemuende.de> <5aa0b3c1-dcdd-22e4-324c-3077e6684f2b@itu.edu.tr>

Sender: owner-ferret_users@xxxxxxxx

Hi Martin,

Did you try to use explicitly "20E" instead of "-20"? fill/l=1/grat/x=20E:40W/y=60S:0N sst

Also, try to do a simple "go land 1" after you filled your figure, to see if the continents are placed correctly according to the lon/lat (if not, then it may be a problem of coordinates in your data)

Cheers,

Nicolas

2017-01-27 11:51 GMT+00:00 baris onol <onolba@xxxxxxxxxx>:

how about this one?

yes? shade/l=1/x=-20:40/y=-60:0/over sst
yes? fill/l=1/x=-20.5:40.5/y=-60:0/over sst
yes? shade/over/lev=(100)/pal=grey/x=-20:40/y=-60:0 missing(sst[l=1],100)

On 01/27/2017 02:09 PM, Martin Schmidt wrote:

Dear ferreters,

I had to compile my pyferret by myself, so I am not sure if the subsequent bug is a general one.

ferret
NOAA/PMEL TMAP
PyFERRET v7 (opt)
Linux 3.19.0-64-generic - 07/01/16

yes? use coads_climatology
yes? fill/l=1/grat/x=-20:40/y=-60:0 sst
yes? shade/over /lev=(100)/pal=grey missing(sst[l=1],100)

In my figure the land west of 20°E remains white.

The result for
yes? fill/l=1/grat/x=340:400/y=-60:0 sst
yes? shade/over /lev=(100)/pal=grey missing(sst[l=1],100)
is complementary.

I find also non-filled points in the function fill_xy near 20°E, where information on adjacent ocean points seems to be "blocked" at 20°E. This looks like if the same error is behind the scene, but in the above example the problem is most obvious.

Setting a region before calling the plot commands, everything is fine for this example:
yes? set region/x=-20:40/y=-60:0
yes? fill/l=1/grat sst
yes? shade/over /lev=(100)/pal=grey missing(sst[l=1],100)

But with fill_xy there is still one point not filled at 20°E.

Seems that 20°E is still a magic boundary?

So my question: For safe code, should it generally avoided to specify regions withing the plot commands? Is there another idea to get safe code.

Best,
Martin

----------------------------------------------------

Nicolas Freychet

PDRA, School of Geosciences
University of Edinburgh
----------------------------------------------------