run /home/ruoying/Software/opnml_matlab5/startup

dir='../obs/';

Sdat = read_s3r([dir 'ICS.S.s3r']);
Tdat = read_s3r([dir 'ICS.T.s3r']);
gdat = loadgrid('../mesh/stw2');

S = reshape(Sdat(:,3),[21 6672]);
T = reshape(Tdat(:,3),[21 6672]);

%== Surface 
subplot(2,2,1)
colormesh2d(gdat,T(1,:));
colorbar; title('Surface Temperature')

subplot(2,2,2)
colormesh2d(gdat,T(21,:));
colorbar; title('Bottom Temperature')

%== Bottom
subplot(2,2,3)
colormesh2d(gdat,S(1,:));
colorbar; title('Surface Salinity')

subplot(2,2,4)
colormesh2d(gdat,S(21,:));
colorbar; title('Bottom Salinity')

