[ncl-talk] Area weighted average - cuvilinear

[Prasad, Siva]

Hi all,

I am trying to do area weighted averaging for a data on curvilinear grid.
var = f->a   ; a(time, yc,xc)
lat2d = f->lat    ; lat(yc,xc)
lon2d = f->lon  ; lon(yc,xc)

Is the following right way?

re = 6.37122e06 ; [m] average radius of earth
rad  = 4.0 * atan(1.0) / 180.0
wgt = cos(lat2d*rad)           ; cosine of latitude

areaAve = wgt_areaave2(var, wgt , 1)

Thanks
Siva
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20181113/bc9a98ac/attachment.html>