23 |
if test1;%the expand faces |
if test1;%the expand faces |
24 |
for iFace=1:mygrid.nFaces; |
for iFace=1:mygrid.nFaces; |
25 |
ii=size(v0{iFace}); ii(1:2)=widefacesSize(iFace,:); |
ii=size(v0{iFace}); ii(1:2)=widefacesSize(iFace,:); |
26 |
tmp1=NaN*zeros(ii); |
tmp1=mygrid.missVal*zeros(ii); |
27 |
if iFace==1|iFace==2;%fill the Northern part of faces 1 and 2 |
if iFace==1|iFace==2;%fill the Northern part of faces 1 and 2 |
28 |
jj=[widefacesSize(iFace,2)-v0facesSize(iFace,2)+1:widefacesSize(iFace,2)]; |
jj=[widefacesSize(iFace,2)-v0facesSize(iFace,2)+1:widefacesSize(iFace,2)]; |
29 |
else; |
else; |
30 |
jj=[1:v0facesSize(iFace,2)]; |
jj=[1:v0facesSize(iFace,2)]; |
31 |
end; |
end; |
32 |
ii=[1:v0facesSize(iFace,1)]; |
ii=[1:v0facesSize(iFace,1)]; |
33 |
tmp1(ii,jj,:)=v0{iFace}; |
tmp1(ii,jj,:,:)=v0{iFace}; |
34 |
v0{iFace}=tmp1; |
v0{iFace}=tmp1; |
35 |
end; |
end; |
36 |
else;%then truncate faces |
else;%then truncate faces |
41 |
jj=[1:mygrid.facesSize(iFace,2)]; |
jj=[1:mygrid.facesSize(iFace,2)]; |
42 |
end; |
end; |
43 |
ii=[1:mygrid.facesSize(iFace,1)]; |
ii=[1:mygrid.facesSize(iFace,1)]; |
44 |
v0{iFace}=v0{iFace}(ii,jj,:); |
v0{iFace}=v0{iFace}(ii,jj,:,:); |
45 |
end; |
end; |
46 |
end; |
end; |
47 |
|
|