Mercurial > octave-nkf

comparison scripts/plot/slice.m @ 10549:95c3e38098bf

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Untabify .m scripts
author Rik <code@nomad.inbox5.com>
date Fri, 23 Apr 2010 11:28:50 -0700
parents 923c7cb7f13f
children be55736a0783
comparison
equal deleted inserted replaced
10548:479536c5bb10 10549:95c3e38098bf
133 hs = zeros(ns,1); 133 hs = zeros(ns,1);
134 [ny, nx, nz] = size (v); 134 [ny, nx, nz] = size (v);
135 if (length(sz) > 0) 135 if (length(sz) > 0)
136 for i = 1:length(sz) 136 for i = 1:length(sz)
137 [xi, yi, zi] = meshgrid (squeeze (x(1,:,1)), 137 [xi, yi, zi] = meshgrid (squeeze (x(1,:,1)),
138 squeeze (y(:,1,1)), sz(i)); 138 squeeze (y(:,1,1)), sz(i));
139 vz = squeeze (interp3 (x, y, z, v, xi, yi, zi, method)); 139 vz = squeeze (interp3 (x, y, z, v, xi, yi, zi, method));
140 tmp(sidx++) = surface (xi, yi, sz(i) * ones (size (yi)), vz); 140 tmp(sidx++) = surface (xi, yi, sz(i) * ones (size (yi)), vz);
141 endfor 141 endfor
142 endif 142 endif
143 143
144 if (length (sy) > 0) 144 if (length (sy) > 0)
145 for i = length(sy):-1:1 145 for i = length(sy):-1:1
146 [xi, yi, zi] = meshgrid (squeeze (x(1,:,1)), sy(i), squeeze (z(1,1,:))); 146 [xi, yi, zi] = meshgrid (squeeze (x(1,:,1)), sy(i), squeeze (z(1,1,:)));
147 vy = squeeze (interp3 (x, y, z, v, xi, yi, zi, method)); 147 vy = squeeze (interp3 (x, y, z, v, xi, yi, zi, method));
148 tmp(sidx++) = surface (squeeze (xi), 148 tmp(sidx++) = surface (squeeze (xi),
149 squeeze (sy(i) * ones (size (zi))), 149 squeeze (sy(i) * ones (size (zi))),
150 squeeze (zi), vy); 150 squeeze (zi), vy);
151 endfor 151 endfor
152 endif 152 endif
153 153
154 if (length (sx) > 0) 154 if (length (sx) > 0)
155 for i = length(sx):-1:1 155 for i = length(sx):-1:1
156 [xi, yi, zi] = meshgrid (sx(i), squeeze (y(:,1,1)), squeeze (z(1,1,:))); 156 [xi, yi, zi] = meshgrid (sx(i), squeeze (y(:,1,1)), squeeze (z(1,1,:)));
157 vx = squeeze (interp3 (x, y, z, v, xi, yi, zi, method)); 157 vx = squeeze (interp3 (x, y, z, v, xi, yi, zi, method));
158 tmp(sidx++) = surface (squeeze (sx(i) * ones (size (zi))), 158 tmp(sidx++) = surface (squeeze (sx(i) * ones (size (zi))),
159 squeeze (yi), squeeze(zi), vx); 159 squeeze (yi), squeeze(zi), vx);
160 endfor 160 endfor
161 endif 161 endif
162 else 162 else
163 vi = interp3 (x, y, z, v, sx, sy, sz); 163 vi = interp3 (x, y, z, v, sx, sy, sz);
164 tmp = surface (sx, sy, sz, vi); 164 tmp = surface (sx, sy, sz, vi);
165 endif 165 endif
166 166
167 if (! ishold ()) 167 if (! ishold ())
168 set (ax, "view", [-37.5, 30.0], "box", "off", "xgrid", "on", 168 set (ax, "view", [-37.5, 30.0], "box", "off", "xgrid", "on",
169 "ygrid", "on", "zgrid", "on"); 169 "ygrid", "on", "zgrid", "on");
170 endif 170 endif
171 171
172 if (nargout > 0) 172 if (nargout > 0)
173 h = tmp; 173 h = tmp;
174 endif 174 endif