velylineoutmovie.m

%movie for flash lineoutsets

figure(5)
i=0;
var1='vely';
var2='tion';
xoffset=0;
x1=xoffset+FLASH.([var1,num2str(i,'%04i')])(:,1);
if (strcmp(var1, 'tele') || strcmp(var1,'tion') || strcmp(var1,'trad'))
    y1=FLASH.([var1,num2str(i,'%04i')])(:,2)./11.6e3;
else
    y1=FLASH.([var1,num2str(i,'%04i')])(:,2)/1e5; %change cm/s to km/s
end

x2=xoffset+FLASH.([var2,num2str(i,'%04i')])(:,1);

if (strcmp(var2, 'tele') || strcmp(var2,'tion') || strcmp(var2,'trad'))
    y2=FLASH.([var2,num2str(i,'%04i')])(:,2)./11.6e3;
else
    y2=FLASH.([var2,num2str(i,'%04i')])(:,2);
end

hline1=plot(x1,y1,'LineWidth',1);
%surf(Z)
axis tight manual
ax1 = gca; %get current axes
ax1.YScale='linear';
ax1_pos = ax1.Position; % position of first axes
ax1.YLim=[-20 80];
ax1.NextPlot = 'replaceChildren'; 
%after 2014b, it's OK, to use a dot to set properties. 
ax1.XAxisLocation='bottom';
ax1.YAxisLocation='left';
title('Axial velocity at 00.00 ns');
xlabel('Position (cm)');
%ylabel('Temperature (eV)')
ylabel('Axial velocity (km/s)')
%legend('Electron Temperautre');
loops = 250;
F(loops) = struct('cdata',[],'colormap',[]);
for j = 1:loops
x1=xoffset+FLASH.([var1,num2str(j,'%04i')])(:,1);
y1=FLASH.([var1,num2str(j,'%04i')])(:,2)./1e5; %change unit
%x2=xoffset+FLASH.([var2,num2str(j,'%04i')])(:,1);
%y2=FLASH.([var2,num2str(j,'%04i')])(:,2)./11.6e3;
plot(ax1,x1,y1,  'LineWidth',1);
title(['Axial velocity at ',num2str(j*0.2,'%02.2f'),' ns']);
%legend('Electron Temperautre', 'Ion Temperature');

  %  drawnow
    F(j) = getframe(gcf);%use gcf to get whole figure frame. 
end
movie(F,1)
%movie(F,2)
%myVideo=VideoWrite('try.avi');
%Undefined function 'VideoWrite' for input arguments of type 'char'.
 
myVideo=VideoWriter(['CRFshock',var1,'.avi']);
myVideo.FrameRate = 10;  % Default 30
myVideo.Quality = 75;    % Default 75
open(myVideo);
writeVideo(myVideo, F);
close(myVideo)