% This fits both Z and X
data=load('data.csv');

d=data(:,3);            % the data
E(:,1)=data(:,1);       % the longitudes
E(:,2)=-data(:,2);       % the depth
E(:,3)=1;

CME=inv(E'*E);
mhat=CME*(E'*d);

dmean=E*mhat;
dprime=d-dmean;
save Tano.txt dprime -ASCII


lonmin=min(data(:,1));
lonmax=max(data(:,1));
[Z,X]=meshgrid(-[5:10:600], lonmin:0.25:lonmax );
DMEAN=X*mhat(1) + Z*mhat(2) + mhat(3);
pcolor(X,Z,DMEAN); shading flat; colorbar


% This fits both Z and X
clear
data=load('data.csv');

d=data(:,3);            % the data
E(:,1)=-data(:,2);       % the depth
E(:,2)=1;

CME=inv(E'*E);
mhat=CME*(E'*d);

dmean=E*mhat;
dprime=d-dmean;
save Tano_Zonly.txt dprime -ASCII


lonmin=min(data(:,1));
lonmax=max(data(:,1));
[Z,X]=meshgrid(-[5:10:600], lonmin:0.25:lonmax );
DMEAN= Z*mhat(1) + mhat(2);
pcolor(X,Z,DMEAN); shading flat; colorbar