Replaced resolution(ipts,1)*resolution(ipts,2) with area(ipts). resolution should not be used as it is only working for lat-lon grids. See ticket #563