% file: mysim_values.sl   A subset of vh1xray.sl 3/8/13

% - - - - - - -
% - - - - - - -
% Example values for the user-set public/global variables.
% Can either edit the values here or, preferably,
% copy this section to a file, e.g., mysim_values.sl,
% and .source it to set the values after vh1xray is loaded
% (see USAGE above and output of vhx_help;).

% - - - Input and output files
% Point to the VH-1 simulation dir, the dirs above it,
% the output filenames prefix, the history file
vhx_abovedir = "/nfs/spectra/d1/Science/Vikram/Tycho_2013";
vhx_simdir = "tycho6_110216";
vhx_simpre = "tycho";
vhx_simhst = "tychohst";
% Output products go into:
vhx_outdir = "./Tycho_out";

% - - - Set start and ending time steps
% 1000 will be added to these to generate the step's filename
vhx_simbeg = NULL;
vhx_simend = NULL;  % the RS moves out of range later
% output only every nth simulation step (can be NULL to be s/w set)
vhx_simnth = NULL;

% - - - Information about the source
% Distance, systemic velocity, and foreground NH
vhx_dkpc = NULL; % ! in kpc
vhx_systvel = 0.0; % ! in km/s
vhx_nh22 = 0.0; % ! x10^22 /cm^2
% Fraction of the simulation that the source fills,
% e.g., if it fills +/-1 half_ang degrees about the equatorial plane
%       then this fraction would be:  sin(half_ang*(PI/180.0));
vhx_frac4pi = 1.0;  % 1.0 for full-4pi

% - - - CD and shocked region
% User define the CD radius at the given-CD time step
vhx_cd0step = NULL;
vhx_cd0rad =  NULL;  % ! in pc
% Temperature above-which material is considered shocked:
vhx_Tshock = 1.e6; % ! in K

% - - - Mass-cut related, incl. Clumping
% The mass-cuts are defined by these two arrays:
% (if 'mass=NULL; then they can be automatically set by vhx_make_cuts)
shl_mass=NULL, shl_icmp=NULL;
% Set mass-cut numbers for the vhx_make_cuts routine:
vhx_ncutsej=12;
vhx_ncutscsm=12;
vhx_ncutsnear=3;  % number of cuts close to the CD, 1 or more.
% Indicate shells that will be clumped in their spectral model,
% this has to be set after the mass-cuts are setup ('mass, 'icmp).
% shl_clumped[ishl] = 0 for a not clumped shell; NULL for no clumping.
shl_clumped=NULL;
% The volume fraction and density factor for the clumping
% (The values are ignored if no clumped shells are specified.)
vhx_clmpfill=0.2;  % volume filling fraction of the clumps
vhx_clmpdens=5.0;  % density factor for the clumps

% - - - Abundances, Modeling, etc.
% Reference abundance set, will be passed to XSPEC models too.
% mxs_abund_pkg implements either: "angr" or "wilm".
vhx_abundstr="angr";
% The nominal XSPEC abundance Zs:
%             H He C N O Ne Mg Si S  Ar Ca Fe Ni
% (The first abundzs must be 1=H, but the rest can
%  be just a subset of the zs, e.g. [1,8,26] for H,O,Fe only.)
vhx_abundzs= [1,2, 6,7,8,10,12,14,16,18,20,26,28];
% Define the Ejecta and CSM relative abundances:
% (an overall scale factor does not matter)
% Set them to all 1.0:
vhx_abundej= 1.0 + 0.0*vhx_abundzs;
vhx_abundcsm= 1.0 + 0.0*vhx_abundzs;
% Other model-related values; if they are NULL then the
% vhx routines will not change them from what the user/isis
% have set them to.
% The XSPEC "NEIVERS" string for NEI models:
vhx_xneiverstr = "1.1";  % "1.0" or "1.1" or "2.0" (or NULL)
% The XSPEC cross-sections version:
vhx_xxsectstr = NULL;  % NULL or "vern" or "bcmc"
% Beta-minimum value for tetp_equilib_bmin function
% (will also use this beta if tetp_equilib_bmin is not available.)
vhx_betamin = 0.06;  % nominal value is 0.06

% - - - Plot options for the Locations plots
% Plot the cummulative mass from outside-in (0)
% or offset it to be zero at the CD (1):
vhx_massrelcd = 0;
% Overplot the density of a reference step - set step number
% or NULL for no overplotting.
vhx_refstep=NULL;
% Can change the (log) Y-axis values if desired,
% (though the axis is used for density, velocity, mass,
%  and temp, so usually a range of at least 10 to
%  10^8 would be used)
vhx_pltymin=0.02;
vhx_pltymax=2.1e10;
% Can set the X-axis min value (or NULL for auto-scaling)
vhx_pltxmin=NULL;  % ! in pc
% Fix the X-axis max value (or NULL for auto-scaling)
vhx_pltxmax=NULL;  % ! in pc

% - - - Energy bands and Spectral Plot
% Arrays of E_lo and E_hi, the first one covers a wide range,
% others set as desired. (Hard-coded to lengths of 3 for now.)
vhx_fluxelos = [0.5, 0.5, 2.0];
vhx_fluxehis = [8.0, 2.0, 8.0];
% Min,Max for the Spectral plots  (erg/s/cm^2/keV)
vhx_fluxmin = 1.e-16;
vhx_fluxmax = 1.e-8; 

% - - - Lightcurve plot option
vhx_lc_xlog = 0;  % linear years (0) or log days (1)
vhx_lc_showejcsm = 0; % show (1) or not (0) separate ej,CSM contrib.s

% End of User defined parameters
% - - - - - - -
% - - - - - - -