window(screenPLT);

% Set a small range around the center wavelength...
xrange(11.8,12.2);

% Show the spatial-only, velocity-only and combined profiles
% for all 4 HETG spectra 

set_par(3,0.0); set_par(4,1.0);
()=eval_counts;
title(demo_shape+" : spatial(red), spectral(blue), combined(black)");
plot_model_flux(megm,1);
oplot_model_flux(megm,2);
set_par(3,1.0); set_par(4,0.0);
()=eval_counts;
oplot_model_flux(megm,4);
set_par(3,1.0); set_par(4,1.0);
()=eval_counts;
oplot_model_flux(megm,1);

set_par(3,0.0); set_par(4,1.0);
()=eval_counts;
title(demo_shape+" : spatial(red), spectral(blue), combined(black)");
plot_model_flux(megp,1);
oplot_model_flux(megp,2);
set_par(3,1.0); set_par(4,0.0);
()=eval_counts;
oplot_model_flux(megp,4);
set_par(3,1.0); set_par(4,1.0);
()=eval_counts;
oplot_model_flux(megp,1);

set_par(3,0.0); set_par(4,1.0);
()=eval_counts;
title(demo_shape+" : spatial(red), spectral(blue), combined(black)");
plot_model_flux(hegm,1);
oplot_model_flux(hegm,2);
set_par(3,1.0); set_par(4,0.0);
()=eval_counts;
oplot_model_flux(hegm,4);
set_par(3,1.0); set_par(4,1.0);
()=eval_counts;
oplot_model_flux(hegm,1);

set_par(3,0.0); set_par(4,1.0);
()=eval_counts;
title(demo_shape+" : spatial(red), spectral(blue), combined(black)");
plot_model_flux(hegp,1);
oplot_model_flux(hegp,2);
set_par(3,1.0); set_par(4,0.0);
()=eval_counts;
oplot_model_flux(hegp,4);
set_par(3,1.0); set_par(4,1.0);
()=eval_counts;
oplot_model_flux(hegp,1);