<div>Brian,</div> <div> </div> <div>I would not agree with your explantion simply because you assumed implicitly that U and T (in your agrument) are 3D vectors inside the operators aave and amean. For model data, this is not applicable because model data is stored on grids, and the area-average operator is applied only on 2D data. Therefore, U and T lie on the tangent surface at each grid point and is thus orthogonal with normal vector A and cos(u) and cos(t) should be zero always (surface A as I understand in the grads user guide is defined to be the lat-lon surface). The differences between amean and aave are explained on grads website. It is due solely to the fact that the latter uses lattitudes to weight the average (a square area defined by lat-lon will shrink to zero at pole) whereas the former does not. So they can not be the same, obviously.</div> <div> </div> <div>Cheers.. </div> <div><BR><B><I>brian
vant-hull <brianvanthull@YAHOO.COM></I></B> wrote:</div> <BLOCKQUOTE class=replbq style="PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: #1010ff 2px solid">You have a mathematical misperception. I can't think of any better way to describe this then by vector mathematics, so go back to first year physics. Think of weighting as a vector dot product. Let A be the area vector, U and T be the other data vectors. Then the weighted averages are A*U/n and A*T/n while the unweighted averages are |U|/n and |T|/n. Canceling out n, you are claiming that<BR><BR>A*U/A*T = |U||T|<BR><BR>but A*U = |A||U|cos(u) and A*T = |A||T|cos(t) where u and t are the angles between the vectors A and U or T. So you are claiming that<BR><BR>|A||U|cos(u)/|A||T|cos(t) = |U|/|T|<BR><BR>This is only true if the angle between the weighting vector and both data vectors are the same. Rarely happens.<BR><BR><B><I>Heiner Körnich
<heiner@MISU.SU.SE></I></B> wrote: <BLOCKQUOTE class=replbq style="PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: rgb(16,16,255) 2px solid">Hi,<BR><BR>I am wondering about the precision of the functions amean and aave. If I<BR>take the data from the GrADS example (model.ctl and model.dat), I would<BR>assume that the following expressions should be equal.<BR><BR>ga-> d amean(u,g)/amean(T,g)<BR>Result value = 0.0550876<BR>ga-> d aave(u,g)/aave(T,g)<BR>Result value = 0.0609393<BR><BR>I thought, the difference between aave and amean is only the final weighting<BR>of the sum with the area.<BR><BR>Do I overlook something or what is wrong here?<BR><BR>Cheers,<BR>Heiner<BR></BLOCKQUOTE><BR><BR><BR>Brian Vant-Hull<BR>301-646-1149 <div> <HR SIZE=1> Want to start your own business? Learn how on <A href="http://us.rd.yahoo.com/evt=41244/*http://smallbusiness.yahoo.com/r-index">Yahoo! Small Business.</A></BLOCKQUOTE><BR><p> Send instant messages to your online
friends http://uk.messenger.yahoo.com