Documentation of nino3

Global Index (short | long) | Local contents | Local Index (short | long)

Function Synopsis

x=nn3(sst,nxa,nya,frame);

Help text

   x = nn3(sst, nxa, nya, frame);

   sst = input sst field
   frame = input frame to average over:  [minlon maxlon minlat maxlat]
           default is [180 270 -6 6]
   x = output time series of averaged sst

Cross-Reference Information

Listing of function nino3

function x=nn3(sst,nxa,nya,frame);

if nargin < 4; frame=[180 270 -6 6]; end;

xkeep=[];
for i=1:length(nxa);
  if (nxa(i) >= frame(1)) & (nxa(i) <= frame(2));
    xkeep = [xkeep i];
  end 
end
ykeep=[];
for i=1:length(nya);
  if (nya(i) >= frame(3)) & (nya(i) <= frame(4))
    ykeep = [ykeep i];
  end
end

keep1=length(nxa)*(ykeep-1);

keep=[];
for i=1:length(ykeep)
  keep=[keep (keep1(i)+xkeep)];
end

x=sst(:,keep);
[m,n] = size(x);
if n > 1
  x=mean2(x');
else
  x=x';
end