[ncl-talk] Mean with same lat and lon

[Rick Brownrigg]

I don't really know for certain, but what about something like this (i.e.,
explicitly computing the average):

data           :=  new((/nlat,mlon/),"float"
data           =  0
dataCount := new((/nlat,mlon/), "integer")
dataCount = 0

do i=.....
   data({lat(i)},{lon(i)}) = data({lat(i)},{lon(i)}) + var(i)
  dataCount({lat(i)},{lon(i)}) = data({lat(i)},{lon(i)}) + 1
end do

data = where(dataCount.ne.0, data/dataCount, data)

Hope that helps...
Rick


On Sat, Jun 27, 2020 at 6:01 PM Vanúcia Schumacher via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Hi users,
>
> I need some help to figure out how I could get the average of the
> variable when there is more than 1 value for the same pixel (lat and lon
> values that are within the same pixel).
>
> In this data example (attached), 3 values of lat and lon are very close,
> so in my script (attached) I get the sum of these values, but I would like
> to change it to get the average value when this happens (lat and lon close).
>
> Variable: var
> Type: float
> Total Size: 48 bytes
>             12 values
> Number of Dimensions: 1
> Dimensions and sizes: [12]
> Coordinates:
> (0) -64.611
> (1) -23.999
> (2) 12.76
> (3) -22.784
> (4) -16.24
> .....
>
> Variable: xvalid
> Type: float
> Total Size: 40 bytes
>             10 values
> Number of Dimensions: 1
> Dimensions and sizes: [10]
> Coordinates:
> Number Of Attributes: 1
>   _FillValue : -999
> (0) -75.85      ======>>>>> sum of  -64.611, -23.999, 12.76 values
> (1) -22.784
> .......
>
> I appreciate some help
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200629/c64a281e/attachment.html>