| 1 |
% $Header: $ |
| 2 |
% $Name: $ |
| 3 |
|
| 4 |
rDir1='./'; rDir2=rDir1; |
| 5 |
%rDir2='res_c00/'; |
| 6 |
%rDir2='./'; |
| 7 |
|
| 8 |
%kdif=0; %- do diff plot (kdif=1) or plot both exp./var/time (kdif=0) |
| 9 |
kdif=0; nf=0; k=0; ccB=[0 0]; ccD=[0 0]; |
| 10 |
|
| 11 |
%- iters to load |
| 12 |
deltaT=3600; iter=[0:5]*24; |
| 13 |
%iter=[0 5]*24; |
| 14 |
%iter=iter(2:end); %- iter=0 missing for some fields |
| 15 |
Nit=length(iter); |
| 16 |
|
| 17 |
%- select variables to plot |
| 18 |
%namv1='ice_fract'; ccB=[-.1 1.1]; |
| 19 |
namv1='ice_iceH'; ccB=[.1 1.1]/2.1; ccD=[-1 1]*0.06; %ccD=[-1 10]/2; |
| 20 |
%namv1='ice_Tsrf'; ccB=[-5 .1]; |
| 21 |
%namv1='T'; ccB=[-.1 0.8]-1.62; |
| 22 |
%namv1='AREA'; ccB=[-.1 1.1]; |
| 23 |
%namv1='HEFF'; ccB=[.1 1.1]/1.1; ccD=[-1 1]*0.06; %ccD=[-1 10]/2; |
| 24 |
%namv1='Qsw'; ccB=[-1.2 0.1]*10; |
| 25 |
%namv1='UICE'; nf=1; ccB=[-.1 0.7]; |
| 26 |
%namv1='VICE'; nf=2; ccB=[-1 1]*.21; |
| 27 |
namv2=namv1; cc2=ccB; |
| 28 |
%namv1='U'; namv2='V'; nf=1; ccB=[-1 6]*.1; cc2=[-1 1]*.2; |
| 29 |
%namv2='ice_iceH'; |
| 30 |
%-------------------------- |
| 31 |
%namv='UICE'; nf=3; |
| 32 |
%namv='U'; k=1; |
| 33 |
%namv='KPPdiffKzS'; k=2; |
| 34 |
|
| 35 |
%- select time index to plot |
| 36 |
nt1=Nit; |
| 37 |
nt2=nt1; |
| 38 |
nt1=2; nt2=Nit; |
| 39 |
|
| 40 |
%- load grid-files: |
| 41 |
G=load_grid(rDir1,0); |
| 42 |
nx=G.dims(1); ny=G.dims(2); |
| 43 |
msk1=ceil(G.hFacC(:,:,1)); |
| 44 |
|
| 45 |
%- load fields: |
| 46 |
v1=rdmds([rDir1,namv1],iter); |
| 47 |
v2=rdmds([rDir2,namv2],iter); |
| 48 |
|
| 49 |
%- kconv=1 : to convert in-situ thickness to effective thickness |
| 50 |
% kconv=-1: or the other way |
| 51 |
kconv=0; |
| 52 |
if kconv==1, |
| 53 |
namv3='ice_fract'; |
| 54 |
v3=rdmds([rDir1,namv3],iter); |
| 55 |
v4=rdmds([rDir2,namv3],iter); |
| 56 |
if strcmp(namv1,'ice_iceH'), v1=v1.*v3; namv1='Ice-Vol'; end |
| 57 |
if strcmp(namv2,'ice_iceH'), v2=v2.*v4; namv2='Ice-Vol'; end |
| 58 |
end |
| 59 |
if kconv==-1, |
| 60 |
namv3='AREA'; |
| 61 |
v3=rdmds([rDir1,namv3],iter); |
| 62 |
v4=rdmds([rDir2,namv3],iter); |
| 63 |
v3(find(v3==0))=-1; v3=1./v3; v3=max(v3,0); |
| 64 |
v4(find(v4==0))=-1; v4=1./v4; v4=max(v4,0); |
| 65 |
if strcmp(namv1,'HEFF'), v1=v1.*v3; namv1='Ice-H'; end |
| 66 |
if strcmp(namv2,'HEFF'), v2=v2.*v4; namv2='Ice-H'; end |
| 67 |
end |
| 68 |
|
| 69 |
var=abs(v2-v1); %var(1:5,:)=0; |
| 70 |
nDim=length(size(v1)); if Nit > 1, nDim=nDim-1; end |
| 71 |
if nDim == 2, ccM=max(max(var)); k=0; |
| 72 |
else ccM=max(max(max(var))); end |
| 73 |
if max(ccM) > 0, fprintf(' %e',ccM) ; fprintf('\n'); end |
| 74 |
%return |
| 75 |
|
| 76 |
titex1=rDir1(1:end-1); titex1=strrep(titex1,'_','\_'); |
| 77 |
titex2=rDir2(1:end-1); titex2=strrep(titex2,'_','\_'); |
| 78 |
titva1=namv1; titva1=strrep(titva1,'_','\_'); |
| 79 |
titva2=namv2; titva2=strrep(titva2,'_','\_'); |
| 80 |
titim1=sprintf('t= %4.1f d',iter(nt1)*deltaT/86400); |
| 81 |
titim2=sprintf('t= %4.1f d',iter(nt2)*deltaT/86400); |
| 82 |
|
| 83 |
xax=[1:nx]-.5 ; yax=[1:ny]-.5; |
| 84 |
xtxt=xax(1); ytxt=yax(1)-5*(yax(2)-yax(1)); |
| 85 |
xydp=[50 20]; xyp0=[50 20]+nf*xydp; xysp=[500 700]; xydp=[100 40]; |
| 86 |
|
| 87 |
nf=nf+1; xyp0=xyp0+xydp; |
| 88 |
figure(nf); set(nf,'position',[xyp0 xysp]);clf; |
| 89 |
subplot(211); |
| 90 |
%for nt=1:Nit, |
| 91 |
if k > 0, var=v1(:,:,k,nt1); else var=v1(:,:,nt1); end |
| 92 |
if strcmp(namv1,'U') & strcmp(namv2,'V'), |
| 93 |
var=var.*var; vtmp=var; |
| 94 |
vtmp(1:nx-1,:)=vtmp(1:nx-1,:)+var(2:nx,:); vtmp(nx,:)=vtmp(nx,:)+var(1,:); |
| 95 |
var=v2(:,:,nt2); var=var.*var; vtmp=vtmp+var; |
| 96 |
vtmp(:,1:ny-1)=vtmp(:,1:ny-1)+var(:,2:ny); |
| 97 |
var=sqrt(vtmp*0.5); titva1='|uVel|'; |
| 98 |
end |
| 99 |
var(find(msk1==0))=NaN; |
| 100 |
mnV=min(var(:)); MxV=max(var(:)); |
| 101 |
imagesc(xax,yax,var'); set(gca,'YDir','normal'); |
| 102 |
if ccB(2) > ccB(1), caxis(ccB/1); end |
| 103 |
change_colmap(-1); |
| 104 |
colorbar; |
| 105 |
grid |
| 106 |
text(xtxt,ytxt,sprintf('min,Max= %9.5g , %9.5g', mnV, MxV)) |
| 107 |
titplot=[titex1,' : ',titva1,' ; ',titim1]; |
| 108 |
title(titplot); |
| 109 |
|
| 110 |
if k > 0, var=v2(:,:,k,nt2)-kdif*v1(:,:,k,nt1); else var=v2(:,:,nt2)-kdif*v1(:,:,nt1); end |
| 111 |
var(find(msk1==0))=NaN; |
| 112 |
mnV=min(var(:)); MxV=max(var(:)); |
| 113 |
subplot(212); |
| 114 |
imagesc(xax,yax,var'); set(gca,'YDir','normal'); |
| 115 |
if cc2(2) > cc2(1) & kdif == 0, caxis(cc2); end |
| 116 |
if ccD(2) > ccD(1) & kdif == 1, caxis(ccD); end |
| 117 |
change_colmap(-1); |
| 118 |
colorbar; |
| 119 |
grid |
| 120 |
text(xtxt,ytxt,sprintf('min,Max= %9.5g , %9.5g', mnV, MxV)) |
| 121 |
if kdif == 1, |
| 122 |
%titplot=['Diff: ',titex2,' - ',titex1]; |
| 123 |
tt1=' -'; tt2=''; |
| 124 |
if ~strcmp(rDir1,rDir2), tt1=[tt1,' ',titex1]; tt2=[tt2,' ',titex2]; end |
| 125 |
if ~strcmp(namv1,namv2), tt1=[tt1,' ',titva1]; tt2=[tt2,' ',titva2]; end |
| 126 |
if nt1 ~= nt2, tt1=[tt1,' ',titim1]; tt2=[tt2,' ',titim2]; end |
| 127 |
titplot=['Diff:',tt2,tt1]; |
| 128 |
else titplot=[titex2,' : ',titva2,' ; ',titim2]; end |
| 129 |
title(titplot); |
| 130 |
put_date; |
| 131 |
%return |
| 132 |
|
| 133 |
clin='kbcmrgy'; |
| 134 |
nf=nf+1; xyp0=xyp0+[xysp(1) 0]; xysp=[400 500]; |
| 135 |
figure(nf); set(nf,'position',[xyp0 xysp]);clf; |
| 136 |
is=61; is2=21; |
| 137 |
%is=67; is=38; |
| 138 |
subplot(211); |
| 139 |
if k > 0, var=v1(is,:,k,:); else var=v1(is,:,:); end |
| 140 |
var=reshape(var,[ny Nit]); |
| 141 |
m1d=squeeze(msk1(is,:)); [J]=find(m1d==0); |
| 142 |
for nt=1:Nit, |
| 143 |
jt=1+rem(nt-1,length(clin)); |
| 144 |
if nt > length(clin); ccln=[clin(jt:jt),'--']; else ccln=[clin(jt:jt),'-'];end |
| 145 |
var(J,nt)=NaN; |
| 146 |
plot(yax,var(:,nt),ccln); |
| 147 |
if nt==1, hold on; end |
| 148 |
end |
| 149 |
hold off; |
| 150 |
AA=axis; axis([0 ny AA(3:4)]); |
| 151 |
grid |
| 152 |
titplot=[titex1,' : ',titva1,' ; is=',int2str(is)]; |
| 153 |
title(titplot); |
| 154 |
|
| 155 |
subplot(212); |
| 156 |
if kdif == 1, is2=is; end |
| 157 |
if k > 0, var=v2(is2,:,k,:)-kdif*v1(is,:,k,:); else var=v2(is2,:,:)-kdif*v1(is,:,:); end |
| 158 |
var=reshape(var,[ny Nit]); |
| 159 |
m2d=squeeze(msk1(is2,:)); [J]=find(m2d==0); |
| 160 |
for nt=1:Nit, |
| 161 |
jt=1+rem(nt-1,length(clin)); |
| 162 |
if nt > length(clin); ccln=[clin(jt:jt),'--']; else ccln=[clin(jt:jt),'-'];end |
| 163 |
var(J,nt)=NaN; |
| 164 |
plot(yax,var(:,nt),ccln); |
| 165 |
if nt==1, hold on; end |
| 166 |
end |
| 167 |
hold off; |
| 168 |
AA=axis; axis([0 ny AA(3:4)]); |
| 169 |
grid |
| 170 |
if kdif == 1, titplot=['Diff: ',titex2,' - ',titex1]; |
| 171 |
else titplot=[titex2,' : ',titva2,' ; is=',int2str(is2)]; end |
| 172 |
title(titplot); |
| 173 |
put_date; |
| 174 |
%return |
| 175 |
|
| 176 |
clin='kbcmrgy'; |
| 177 |
nf=nf+1; xyp0=xyp0+[0 xysp(2)]; |
| 178 |
figure(nf); set(nf,'position',[xyp0 xysp]);clf; |
| 179 |
js=41; js2=6; |
| 180 |
subplot(211); |
| 181 |
if k > 0, var=v1(:,js,k,:); else var=v1(:,js,:); end |
| 182 |
var=reshape(var,[nx Nit]); |
| 183 |
m1d=squeeze(msk1(:,js)); [J]=find(m1d==0); |
| 184 |
for nt=1:Nit, |
| 185 |
jt=1+rem(nt-1,length(clin)); |
| 186 |
if nt > length(clin); ccln=[clin(jt:jt),'--']; else ccln=[clin(jt:jt),'-'];end |
| 187 |
var(J,nt)=NaN; |
| 188 |
plot(xax,var(:,nt),ccln); |
| 189 |
if nt==1, hold on; end |
| 190 |
end |
| 191 |
hold off; |
| 192 |
grid |
| 193 |
titplot=[titex1,' : ',titva1,' ; js=',int2str(js)]; |
| 194 |
title(titplot); |
| 195 |
|
| 196 |
subplot(212); |
| 197 |
if kdif == 1, js2=js; end |
| 198 |
if k > 0, var=v2(:,js2,k,:)-kdif*v1(:,js,k,:); else var=v2(:,js2,:)-kdif*v1(:,js,:); end |
| 199 |
var=reshape(var,[nx Nit]); |
| 200 |
m2d=squeeze(msk1(:,js2)); [J]=find(m2d==0); |
| 201 |
for nt=1:Nit, |
| 202 |
jt=1+rem(nt-1,length(clin)); |
| 203 |
if nt > length(clin); ccln=[clin(jt:jt),'--']; else ccln=[clin(jt:jt),'-'];end |
| 204 |
var(J,nt)=NaN; |
| 205 |
plot(xax,var(:,nt),ccln); |
| 206 |
if nt==1, hold on; end |
| 207 |
end |
| 208 |
hold off; |
| 209 |
grid |
| 210 |
if kdif == 1, titplot=['Diff: ',titex2,' - ',titex1]; |
| 211 |
else titplot=[titex2,' : ',titva2,' ; js=',int2str(js2)]; end |
| 212 |
title(titplot); |
| 213 |
put_date; |
| 214 |
return |
| 215 |
|
| 216 |
for nG=1:2, |
| 217 |
if nG==1, var=v1; titplot=[titex1,' : ',titvar,' ; it=',int2str(nit)]; |
| 218 |
% else var=v2; titplot=[titex2,' : ',titvar,' ; it=',int2str(nit)]; end |
| 219 |
else var=v2-v1; titplot=['Diff: ',titex2,' - ',titex1]; end |
| 220 |
subplot(210+nG); |
| 221 |
|
| 222 |
ccM=max(abs(var)); |
| 223 |
if ccM > 0, ccE=floor(log10(ccM)); ccM=ceil(ccM*10^-ccE)*10^ccE; else ccM = 1; end |
| 224 |
%ccM=0; |
| 225 |
%if nG == 2, ccM=ccM/1000; end |
| 226 |
ccB=[-1 1]*ccM; |
| 227 |
plot(xax,var','k'); |
| 228 |
|
| 229 |
%mnV=min(min(var)); MxV=max(max(var)); |
| 230 |
%fprintf('min,Max= %9.5g , %9.5g \n', mnV, MxV) |
| 231 |
%text(xtxt,ytxt,sprintf('min,Max= %9.5g , %9.5g', mnV, MxV)) |
| 232 |
|
| 233 |
title(titplot); |
| 234 |
end |
| 235 |
|
Parent Directory
| 
Revision Log
| 
Revision Graph