[ferret_users] Bug in regressxy.jnl

[Jian Ma <jianma@xxxxxxxxxx>]

Hi all,

The current version of regressxy.jnl did the x-average first and then y.
This has no problem if the region is lon-lat square, but in irregular
shaped regions like a continent or ocean, the weight would be incorrect
for each grid. I suggest to change the commands to do x and y averaging
at the same time, e.g.

let pave = pmasked[x=@ave,y=@ave]
let qave = qmasked[x=@ave,y=@ave]
let pqave = pq[x=@ave,y=@ave]
let ppave = pp[x=@ave,y=@ave]
let qqave = qq[x=@ave,y=@ave]

Well if these are based on my misunderstanding, pls point it out.

Thank you very much,
Tony