function [ LA,LO,Z]=loadgridmayer2(lo1,lo2, la1,la2);

%  [ LA,LO,Z]=loadgridmayer(lo1,lo2, la1,la2);
%   (min, max long, min max lat) ( SW06 long ~-73)

load NJshelfbathy lalodep

% cluster area 39-01, 39-04 ... 39.0167- 39.0667.
  %  73-02 -73.05  73.0333-73.0833 .
  
minla=min(la1 , la2); maxla=max(la1 , la2);
minlo=min(lo1 , lo2); maxlo=max(lo1 , lo2);

la=lalodep(:,1);
eps=.006;
some=find(la > minla-eps & la < maxla+eps);
  
ar=lalodep(some',:);
lo=ar(:,2);
los=find(lo> minlo-eps & lo < maxlo+eps);

areax=ar(los,:);
la=areax(:,1);
z=areax(:,3);
lo=areax(:,2);
loo=unique(lo);laa=unique(la);
[LO,LA]=meshgrid(loo,laa);
nloo=length(loo);
nlaa=length(laa);
LOO=reshape(LO,[nloo*nlaa,1]);
LAA=reshape(LA,[nloo*nlaa,1]);
disp([ 'lgm2: gridding la x lo ' num2str(nlaa) ' x ' num2str(nloo) ' = '...
      num2str(length(LAA)) ' pts'])
  %whos

Z=griddatan(areax(:,1:2),areax(:,3),[LAA LOO]);
Z=reshape(Z,[nlaa,nloo]);


