% Thickness levels for Samoan Passage
% Have 10-m thicknesses between:
% Min depth 3792.495 (shallowest depth of 1 C isotherm)
% Max depth 5870.09814453125 (deepest bathymetry)

load ../grid/thk90

ik=find(thk90<82);
thk=thk90(ik);
m=1.05;
while sum(thk)<3792
    thk=[thk 10+(thk(end)-10)*m];
    m=m*.9905;
end
while sum(thk)<5871
    thk=[thk 10];
end
clf, plot(thk,-cumsum(thk),'.'), hold on
plot([0 100],-[3792 3792])
xlabel('thikness (m)')
ylabel('depth (m)')
title(['Numer of vertical levels: ' int2str(length(thk))])
print -djpeg SamoanDepths

dpt=(cumsum(thk)+cumsum([0 thk(1:end-1)]))/2;
save thk325 dpt thk