[ncl-talk] Getting the grid point values of certain values in a masked array

[Barry Lynn]

Hi:

Have you tried using PrintMinMAX(Variable, False) to see the range of your
data?

Barry

On Thu, Sep 13, 2018 at 1:53 PM Ty Buckingham <
ty.buckingham at manchester.ac.uk> wrote:

> I'm not entirely sure of the syntax. If I try something like:
>
>   avo_masked  = where(avo_zoom.gt.0.006,1,avo_zoom at _FillValue)
>
>   avo_gridpts = where(avo_masked.eq.1,avo_masked,0)
>
>
> I still get all of the false values and their corresponding grid points,
> when I only want avo_gridpts to contain the values equal to 1 and
> their corresponding grid point values.
>
>
> Thanks
>
> Ty
> ------------------------------
> *From:* Barry Lynn [barry.h.lynn at gmail.com]
> *Sent:* 13 September 2018 10:50
> *To:* Ty Buckingham
> *Subject:* Re: [ncl-talk] Getting the grid point values of certain values
> in a masked array
>
> Hi:
>
> There should be.  If you want to find only those values, then set the
> middle part of your where statement equal to that.
>
> However, unless  the values are exact you might have a problem finding
> them.  Then, you might multiply by 1000 and then include the integer
> function.
>
> Barry
>
> On Thu, Sep 13, 2018 at 12:46 PM Ty Buckingham <
> ty.buckingham at manchester.ac.uk> wrote:
>
>> Hi Barry,
>>
>> It is indeed. Is there anyway to extract the grid points of all the
>> values equal to 0.006 and ignore the rest?
>>
>> Thanks
>> Ty
>> ________________________________________
>> From: Barry Lynn [barry.h.lynn at gmail.com]
>> Sent: 13 September 2018 06:00
>> To: Ty Buckingham
>> Cc: ncl-talk
>> Subject: Re: [ncl-talk] Getting the grid point values of certain values
>> in a masked array
>>
>> Hi:
>>
>> The where statement:
>>
>>
>> avo_masked  = where(avo_zoom.gt.0.006,0.006,avo_zoom at _FillValue)
>>
>>
>> will set every grid with a value of 0.006 or greater equal to 0.006.
>>
>>
>> Is that what you wanted to do?
>>
>>
>> Barry
>>
>> On Thu, Sep 13, 2018 at 12:18 AM Ty Buckingham <
>> ty.buckingham at manchester.ac.uk<mailto:ty.buckingham at manchester.ac.uk>>
>> wrote:
>> Hi all,
>>
>> I have a masked array and wish to get the grid points (or lat/lons) of
>> all of the values that are above a certain threshold. I've been messing
>> around with 'where' and 'if' statements but can't quite get it.
>>
>> I've attached my script.
>>
>> Any help?
>>
>> Thanks in advance!
>> Ty
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu<mailto:ncl-talk at ucar.edu>
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>> --
>> Barry H. Lynn, Ph.D
>> Senior Associate Scientist, Lecturer,
>> The Institute of the Earth Science,
>> The Hebrew University of Jerusalem,
>> Givat Ram, Jerusalem 91904, Israel
>> Tel: 972 547 231 170
>> Fax: (972)-25662581
>>
>> C.E.O, Weather It Is, LTD
>> Weather and Climate Focus
>> http://weather-it-is.com
>> Jerusalem, Israel
>> Local: 02 930 9525
>> Cell: 054 7 231 170
>> Int-IS: x972 2 930 9525
>>
>>
>
> --
> Barry H. Lynn, Ph.D
> Senior Associate Scientist, Lecturer,
> The Institute of the Earth Science,
> The Hebrew University of Jerusalem,
> Givat Ram, Jerusalem 91904, Israel
> Tel: 972 547 231 170
> Fax: (972)-25662581
>
> C.E.O, Weather It Is, LTD
> Weather and Climate Focus
> http://weather-it-is.com
> Jerusalem, Israel
> Local: 02 930 9525
> Cell: 054 7 231 170
> Int-IS: x972 2 930 9525
>
>

-- 
Barry H. Lynn, Ph.D
Senior Associate Scientist, Lecturer,
The Institute of the Earth Science,
The Hebrew University of Jerusalem,
Givat Ram, Jerusalem 91904, Israel
Tel: 972 547 231 170
Fax: (972)-25662581

C.E.O, Weather It Is, LTD
Weather and Climate Focus
http://weather-it-is.com
Jerusalem, Israel
Local: 02 930 9525
Cell: 054 7 231 170
Int-IS: x972 2 930 9525
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180913/15b30032/attachment.html>