Sdat = read_s3r('ICS.S.s3r');
Tdat = read_s3r('ICS.T.s3r');
gdat = loadgrid('../MESH/stw2');

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

figure(1)
colormesh2d(gdat,S(21,:));
colorbar

figure(2)
colormesh2d(gdat,T(21,:));
colorbar
