Area weighting giving too low answers

I am trying the code below in order to area weight my data. The data looks good except it is 5 orders of magnitude lower than I want or expected. I have tried modifying a few things but can't seem to get it to work. Can anyone see where I have gone wrong?

[x,y]=meshgrid(0:1:360,-66:1:66);

x=x(1:end-1,1:end-1); y=y(1:end-1,1:end-1);

    weight = cos(y(:,1));
    wrepa = repmat(weight,[1 120 663]);
    wrepa = permute(wrepa,[3 1 2]);

    aa = areaquad(y(:,1)-(1/2),y(:,1)*0,y(:,1)+(1/2),(y(:,1)*0)+1);
    repa = repmat(aa,[1 120 663]);
    repa = permute(repa,[3 1 2]);
    b = sum(aa)*360;

    AMW = AM.*wrepa;
    AMWW = AMW.*(repa/b);
    AMWW = nanmean3d(AMWW,1);

    AMWms = (smooth(AMWW,12));

AMWms = 663x132x120;

with values between 1 to -2 spatially.

thanks, Michael

[Michael]

so basically, I get answers 10^-8 rather than 10-03....

[Michael]

10^-03

[Michael]

so my aim is to weight my gridded data with respect to area. My grid resolution is 132 x 120.

[Supreeth Subbaraya]

What is the variable "AM" in your code?

[Michael]

AM is the pre-weighted data with same dimensions as AMWms

[Michael]

Can anyone help me?

Answer 1

[Michael]

Can anyone help me?