Re: Problem of area weighted sum using "wgt_areasum2" function

Jun,

Your formula

   dy(1:jlat-2) = abs(lat1(2:jlat-1)-lat1(1:jlat-2))*rr*0.5

looks suspicious to me. It looks like your differencing latitudes at
adjacent points, but dividing by two. I would either drop the 0.5, or
replace lat1(1:jlat-2) with lat1(0:jlat-3).

Keith

-- 
******************************************************************
Keith Lindsay                http://www.cgd.ucar.edu/oce/klindsay/
email: klindsay_at_ucar.edu   phone: 303-497-1722   fax: 303-497-1700
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk