% file: dia_mysim.sl     dia = "do it all" 
%
% This file is edited to contains all of the processing steps to process
% a simulation. This allows convenient reprocessing or the processing
% of new similar simulations.
% 
%   Usage:
% [linux] isis
%   . . . 
%  isis> .source dia_mysim
%  . . . [ All of the processing steps are carried out. ]
% 

% Can (manually) delete previous contents of the output dir by doing:
% isis> !   rm -I MySim_out/*

% Load the software etc
.load vh1xray
;

% Set the user-parameters for the VH-1 simulation
.source mysim_values
;


% Can change any of the nominal values that are in mysim_values.sl
% e.g., change the output directory:
vhx_outdir = "./MySim_out";


% Reads the hst file, sets things up
vhx_setup;

% Manually look at a step, useful to explore a simulation
%% vhx_read_step(90); vhx_plot_step;


% Manually determine and set the CD location at a specific time step
%
%% vhx_cd0step = 32;
%% vhx_cursor_cd0;
%
% Copy the result here
% (or put the values in mysim_values.sl and these lines are not needed)
vhx_cd0step = 32;
vhx_cd0rad = 0.0403747;


% Can now calc the CD at all time steps:
vhx_calc_cds;

% Plot RS, CD, FS:
vhx_plot_cds;
vhx_psplot_cds;


% Setup the mass cuts
vhx_make_cuts;

% If desired, turn on clumping for CSM
shl_clumped = 0*shl_icmp;  % define and reset clumped array
shl_clumped[ where(shl_icmp == 2) ] = 1;

% Look at the amount of shocked mass from CD to RS (shocked ejecta)
% and from CD to FS (shocked CSM)
vhx_plot_mass;
vhx_psplot_mass;


% Make "locations" plots of all steps and assemble to a movie
vhx_plot_allsteps;
% command to make movie is output to screen, e.g.,:
%% ! convert -delay 10 ./MySim_out/locs_*.gif ./MySim_out/mov_locs.gif


% Go through the time steps and calculate and write out kT, tau, etc
% values for the mass-cut shells (every vhx_simnth one is written.)
vhx_calc_shells;


% Make a single par file at some time step:
vhx_make_par(100);

% Make plot of the mass-cut shells' kT-tau values at a step:
vhx_plot_kTtau(100);
vhx_psplot_kTtau(100);


% Make all par files, spectra plots (and movie) and the flux-loc file:
vhx_calc_allspec;
%  make movie:
%%  ! convert -delay 10 ./MySim_out/spectra_*.gif ./MySim_out/mov_spectra.gif


% Plot the lightcurves
% (uses data in the flux-loc file: ./MySim_out/fluxlocs_vs_time_Obs.txt) 
vhx_plot_lc;
vhx_psplot_lc;


% - - - end of dia_mysim - - -