おそらくやりたいことは、n, mに変数rightのインデックスの連番をつくりたいのかなぁと予想し、下記に修正コードを書きました。
right(isnan(right))=0;
s=surf(right);
zlim([inR axR]);
% view([-30.6 57.0])
s.EdgeColor = 'none';
for n=1:size(right,1) % rightの1次元目のインデックスの連番を作成
for m=1:size(right,2)% rightの2次元目のインデックスの連番を作成
if Data(n,m,face(n,m))==0
Data(n,m,face(n,m))=1;
end
end
end
