%read in Levitus94 temperature profile data: 
%Equatorial Pacific (170 W, 0 S) temperature
%from http://iridl.ldeo.columbia.edu/SOURCES/.LEVITUS94/.ANNUAL/.temp/

%plot temperature profile
%ian, 4/04

m=dlmread('levitus94.data','\t',4,0);
plot(m(:,1),m(:,2))
z=m(:,1);
t_data=m(:,2)
t_fit=1.2+16.3*exp(-z/600);
plot(t_data,z,t_fit,z)
set(gca,'YLim',[500 5500],'YDir','reverse')
xlabel('temperature (^oC)'),ylabel('depth (m)')
title('Equatorial Pacific temperature profile')
legend('data','fit from PS 1 soln',4)