% file: x1822_movie.sl
%
% Create projected images from a range of phases to
% use to make a movie :-)
%  
% Can also get light-curves from the components
% by saving the sum of the opaque-adjusted component arrays
% and plotting them.

variable phases, iphase;

phases = [0.0:1.0:#101];

% light curve vlaues:
lc_corona = Float_Type[length(phases)];
lc_disk = Float_Type[length(phases)];
lc_bulge = Float_Type[length(phases)];

_for iphase (0,length(phases)-1,1)
{
   s3d_ps.phase = phases[iphase];
   s3d_update();

   % Output images?
   if ( 1 > 0) 
     {
	% Make the images
	s3d_proj_comps;
     }

	% Save the component norms factors  
   lc_corona[iphase] = s3d_comp[1].norm*sum(s3d_opaque_mask * s3d_comp[1].arr);
   lc_disk[iphase] = s3d_comp[2].norm*sum(s3d_opaque_mask * s3d_comp[2].arr);
   lc_bulge[iphase] = s3d_comp[5].norm*sum(s3d_opaque_mask * s3d_comp[5].arr);
  	
}

% Make the plot

variable PSwin = open_plot("x1822_lc.ps/cps");
xrange(0.0,2.0);
yrange(0.,1.05*max(lc_corona+lc_disk+lc_bulge));
label("Phase", "Relative norm", "X1822 Model Light Curve(Blu) = NS-Corona(Grn) + Disk(Red) + Bulge(Orange)");
%
plot(phases, lc_corona+lc_disk+lc_bulge,4);
oplot(phases+1.0, lc_corona+lc_disk+lc_bulge,4);
oplot(phases, lc_corona, 3);
oplot(phases+1.0, lc_corona, 3);
oplot(phases, lc_disk, 2);
oplot(phases+1.0, lc_disk, 2);
oplot(phases, lc_bulge, 8);
oplot(phases+1.0, lc_bulge, 8);

close_plot(PSwin);