108 |
%titall='Dyncore test-case 5 (cs-32)' ; |
%titall='Dyncore test-case 5 (cs-32)' ; |
109 |
|
|
110 |
%========================================================= |
%========================================================= |
111 |
ng=0; fxb=100; fyb=60; |
ng=0; fxb=100; fyb=60; fdx=100; fdy=40; fsc=1.; |
112 |
%fyb=160; fxb=-2600; |
%fyb=-360; fxb=-2600; fdy=60; fsc=1.5; |
113 |
|
|
114 |
for jv=1:nbV, |
for jv=1:nbV, |
115 |
%------------------- |
%------------------- |
127 |
ng=ng+1; |
ng=ng+1; |
128 |
%- reset "ng" to jv for fix fig number (independent of which one is in list_on): |
%- reset "ng" to jv for fix fig number (independent of which one is in list_on): |
129 |
%ng=jv; |
%ng=jv; |
130 |
figure(ng); set(ng,'position',[fxb+100*ng fyb+40*ng 500 700]);clf; |
figure(ng); set(ng,'position',[fxb+fdx*ng fyb+fdy*ng [500 700]*fsc]);clf; |
131 |
var=squeeze(vv1(1+kl,:,1,:,:)); dd=zeros(5,Nexp); av=zeros(5,Nexp); |
var=squeeze(vv1(1+kl,:,1,:,:)); dd=zeros(5,Nexp); av=zeros(5,Nexp); |
132 |
for n=1:Nexp, |
for n=1:Nexp, |
133 |
tmp=var(isA(n):ieA(n),:,n); [I]=find(tmp(:,1)==undef); tmp(I,:)=0; |
tmp=var(isA(n):ieA(n),:,n); [I]=find(tmp(:,1)==undef); tmp(I,:)=0; |