1 |
nx=20; ny=16; nz=23; |
2 |
|
3 |
fld1=readbin('KPPviscA.0000000010.data',[nx ny]); |
4 |
fld2=readbin('KPPviscAz.0000000011.data',[nx ny]); |
5 |
clf, subplot(311),mypcolor(fld1');colorbar |
6 |
subplot(312),mypcolor(fld2');colorbar |
7 |
subplot(313),mypcolor(fld2'-fld1');colorbar |
8 |
|
9 |
fld1=readbin('KPPdiffS.0000000010.data',[nx ny]); |
10 |
fld2=readbin('KPPdiffKzS.0000000011.data',[nx ny]); |
11 |
clf, subplot(311),mypcolor(fld1');colorbar |
12 |
subplot(312),mypcolor(fld2');colorbar |
13 |
subplot(313),mypcolor(fld2'-fld1');colorbar |
14 |
|
15 |
fld1=readbin('KPPdiffT.0000000010.data',[nx ny]); |
16 |
fld2=readbin('KPPdiffKzT.0000000011.data',[nx ny]); |
17 |
clf, subplot(311),mypcolor(fld1');colorbar |
18 |
subplot(312),mypcolor(fld2');colorbar |
19 |
subplot(313),mypcolor(fld2'-fld1');colorbar |
20 |
|
21 |
fld1=readbin('KPP_ghat.0000000010.data',[nx ny]); |
22 |
fld2=readbin('KPPghat.0000000011.data',[nx ny]); |
23 |
clf, subplot(311),mypcolor(fld1');colorbar |
24 |
subplot(312),mypcolor(fld2');colorbar |
25 |
subplot(313),mypcolor(fld2'-fld1');colorbar |
26 |
|
27 |
fld1=readbin('KPP_hbl.0000000010.data',[nx ny]); |
28 |
fld2=readbin('KPPhbl.0000000011.data',[nx ny]); |
29 |
fld3=readbin('KPPmld.0000000010.data',[nx ny]); |
30 |
clf, subplot(311),mypcolor(fld1');colorbar |
31 |
subplot(312),mypcolor(fld2');colorbar |
32 |
subplot(313),mypcolor(fld3');colorbar |
33 |
|
34 |
fld1=readbin('KPPfrac.0000000010.data',[nx ny]); |
35 |
mypcolor(fld1');colorbar |
36 |
|
37 |
%%%%%%%%%%%%%%% |
38 |
delZ=[0 10 10 15 20 20 25 35 50 75 100 150 200 ... |
39 |
275 350 415 450 500 500 500 500 500 500 500]; |
40 |
top=cumsum(delZ(1:23)); |
41 |
bot=cumsum(delZ(2:24)); |
42 |
dpt=-(top+bot)/2; |
43 |
MLD=readbin('KPPmld.0000000010.data',[nx ny]); |
44 |
T=readbin('T.0000000011.data',[nx ny nz]); |
45 |
S=readbin('S.0000000011.data',[nx ny nz]); |
46 |
R=rho(S,T,0*T); |
47 |
DR=rho(S(:,:,1),T(:,:,1)-.8,0)-R(:,:,1); |
48 |
mld=0*MLD; |
49 |
for i=1:nx |
50 |
for j=1:ny |
51 |
if DR(i,j)>0 |
52 |
tmp=squeeze(R(i,j,:)-R(i,j,1)); |
53 |
in=closest(DR(i,j),tmp,0); |
54 |
mld(i,j)=(dpt(in(1))*abs(tmp(in(2))-DR(i,j)) + ... |
55 |
dpt(in(2))*abs(tmp(in(1))-DR(i,j))) / ... |
56 |
abs(R(i,j,in(2))-R(i,j,in(1))); |
57 |
end |
58 |
end |
59 |
end |
60 |
clf, subplot(311),mypcolor(fld3');caxis([0 70]),colorbar |
61 |
subplot(312),mypcolor(-mld');caxis([0 70]),colorbar |
62 |
subplot(313),mypcolor(-mld'-fld3');colorbar |
63 |
|
64 |
%%%%%%%%%%%%%%%%%%%%%%% |
65 |
ix=10; iy=8; |
66 |
s=squeeze(S(ix,iy,:)); |
67 |
t=squeeze(T(ix,iy,:)); |
68 |
p=pressure(abs(dpt),26); |
69 |
r=rho(s,t,0); |
70 |
dr= |
71 |
clf, plot(s,dpt,t,dpt,r-r(1),dpt,it,dpt,dr+dpt*0,dpt,'linewidth',3) |
72 |
legend('S','T','\sigma','dr',4) |
73 |
axis([0 37 -5450 0]), grid |
74 |
title(MLD(ix,iy)) |