[MITgcm-support] matlab script to internal_wave example , Anyone have?

[brainf]

Hi all,
   How to plot the internal_wave figure like manual.(page 20 ,Figure 1.9).
I just started to use Mitgcm,my interest is on internal_wave.
I have got a  internal_wave example result.(..\verification\internal_wave),
But I  have problem to see this result by matlab ,
Anyone have the matlab script to internal_wave example ?    Thx very.

or I have write some, who can help me continue it....

%%%%%%%%%%%%%%%% matlab script for internal_wave example start %%%%%%%%%%%%%%%%
H=rdmds('Depth');
T=rdmds('T');
T2(60,20)=0;

for j=1 : 20
T2(:,j)=T(:,1,j);
end

S=rdmds('S');
T2(60,20)=0;

for j=1 : 20
S2(:,j)=S(:,1,j);
end

rho(60,20)=0;
for j=1 : 20
rho(:,j) = densmdjwf(S2(:,j),T2(:,j).*10,j);
end

imagesc(T2');colorbar;
%imagesc(rho');colorbar;
figure;
contour(rho')
%%%%%%%%%%%%%%%%%%%%%%%%             end                 %%%%%%

Best regards. 

brainf
brainf at 163.com
2006-05-21