Hello.
How can I do this in octave?
%-------------------------------------------
function [] = plotting()
clc; more off;
delimiterIn = ' ';
headerlinesIn = 1;
files_mesh = 'MESH_Coordinate*.dat';
files_solution = 'SOLUTION_*.dat';
%------------------------------%
% Importing the MESH DATA
%------------------------------%
import_mesh = dir(files_mesh);
nfiles_mesh = length(import_mesh);
for i=1:nfiles_mesh
if (i==1)
filename = import_mesh(i).name;
mydata_mesh = importdata(filename,delimiterIn,headerlinesIn);
x = mydata_mesh.data(:,1);
xIni = min(x);
xEnd = max(x);
end
if (i==2)
filename = import_mesh(i).name;
mydata_mesh = importdata(filename,delimiterIn,headerlinesIn);
y = mydata_mesh.data(:,1);
yIni = min(y);
yEnd = max(y);
end
end
[X,Y] = meshgrid(x,y);
%------------------------------%
% Importing the SOLUTION DATA
%------------------------------%
import_solution = dir(files_solution);
nfiles_solution = length(import_solution);
v = VideoWriter('Video_filename.avi');
open(v)
for i=1:nfiles_solution
filename = import_solution(i).name;
mydata_solution = importdata(filename,delimiterIn,headerlinesIn);
u = mydata_solution.data(:,1);
uMin = 0;
uMax = 5;
frame = getframe(gcf);
writeVideo(v,frame);
aFile = sprintf('%s',filename);
fprintf('Plotting file %s\n',aFile);
U = reshape(u,[size(X,2),size(Y,1)]);
U = transpose(U);
mesh(X,Y,U);
xlim([xIni,xEnd]);
ylim([yIni,yEnd]);
zlim([uMin,uMax])
pbaspect([1 1 1]);
daspect([1 1 1]);
pause(0.1);
end
close(v)
end
%------------------------------------------
Thanks
--
Cordialmente,
| Bernardo José Luján Esquivia Estudiante de Doctorado Tel: 3106027532 Fundación Universidad del Norte |