<div dir="ltr"><div dir="ltr"><div>Hi all,</div><div><br></div><div>I am trying to do area weighted averaging for a data on curvilinear grid.</div><div>var = f->a   ; a(time, yc,xc)<br></div><div>lat2d = f->lat    ; lat(yc,xc) <br></div><div>lon2d = f->lon  ; lon(yc,xc)</div><div><br></div><div>Is the following right way?</div><div><br></div><div>re = 6.37122e06 ; [m] average radius of earth<br>rad  = 4.0 * atan(1.0) / 180.0            <br>wgt = cos(lat2d*rad)           ; cosine of latitude</div><br></div><div dir="ltr">areaAve = wgt_areaave2(var, wgt , 1) <br><div><br></div><div>Thanks</div><div>Siva<br></div><div><br></div></div></div>