1 |
prefix='stDiag'; |
prefix='dynStD'; |
2 |
namA='r527'; |
namA='r527'; |
3 |
Nexp=1; Nc=size(namA,2); |
Nexp=1; Nc=size(namA,2); |
4 |
%-- |
%-- |
5 |
|
|
64 |
end |
end |
65 |
if krd ~= 0, |
if krd ~= 0, |
66 |
ttA=squeeze(tiA(:,2,:)); |
ttA=squeeze(tiA(:,2,:)); |
67 |
ttA=ttA/86400; titT='days'; %ttA=ttA/30 ; titT='month'; ttA=ttA/12 ; titT='year'; |
ttA=ttA/3600; titT='hrs'; ttA=ttA/24; titT='days'; |
68 |
|
ttA=ttA/30 ; titT='month'; ttA=ttA/12 ; titT='year'; |
69 |
end |
end |
70 |
%========================================================= |
%========================================================= |
71 |
|
|
166 |
subplot(410+nv); |
subplot(410+nv); |
167 |
var=squeeze(vv1(2:nk,:,1,nv,n))'; |
var=squeeze(vv1(2:nk,:,1,nv,n))'; |
168 |
mnV=min(var(:)); MxV=max(var(:)); ccv=c_levs(mnV,MxV,-12); |
mnV=min(var(:)); MxV=max(var(:)); ccv=c_levs(mnV,MxV,-12); |
169 |
[cs,h]=contour(ttA(isA(n):ieA(n),n),[1:nk-1],var(isA(n):ieA(n),:)',ccv); |
%- pressure coordinate, from bottom (k=1) to top (k=nr=nk-1) |
170 |
|
%[cs,h]=contour(ttA(isA(n):ieA(n),n),[1:nk-1],var(isA(n):ieA(n),:)',ccv); |
171 |
|
%- z coordinate, from bottom (k=nr=nk-1) to top (k=1) |
172 |
|
[cs,h]=contour(ttA(isA(n):ieA(n),n),-[1:nk-1],var(isA(n):ieA(n),:)',ccv); |
173 |
%clabel(cs);isoline0(h); |
%clabel(cs);isoline0(h); |
174 |
colorbar; |
colorbar; |
175 |
if nv == 1, title(['Avr ',titv]); end |
if nv == 1, title(['Avr ',titv]); end |