Re: [ferret_users] Counting Data-points

[Ansley Manke <ansley.b.manke@xxxxxxxx>]

Hi Luke,
It's the @NGD transformation:

Here is a mask: 1 if the variable is greater than 11, otherwise missing.

   yes? let var_hi = IF 
var[x=170E:170W,y=30S:40S,t=10-Jan-2013:10-Dec-2013] GT 11 then 1
yes? list var_hi[x=@ngd,y=@ngd,t=@ngd]

Because of the way I defined that mask, the @SUM transformation would 
give the same results.

You might want to see how the number greater than some value is changing 
over time:

yes? let var_hi = IF var[x=170E:170W,y=30S:40S] GT 11 then 1
yes? list/t=10-Jan-2013:10-Dec-2013 var_hi[x=@ngd,y=@ngd]

(Actually the data you point to looks like it's always greater than 11 
in that time and region.  Try it with, say, 18 instead.)

Ansley

On 12/4/2014 6:10 PM, Luke M wrote:
> Dear all
>
> I'm trying to count the number of data points that are above a certain 
> value for a given region and time period, but I am a bit stuck.
>
> Does anyone have any suggestions on how I could find what are, for 
> each month, the total number of data points (pixels) that have a value 
> higher than 11 degrees for a given small area & time period:
>
> x=170E:170W/y=30S:40S/t=10-Jan-2013:10-Dec-2013
>
> Using the sst dataset:
> set data "http://oceanwatch.pifsc.noaa.gov/thredds/dodsC/pfgac/monthly";
>
> Thanks!
>
>
>
>
> -