%---------------------- Database Setup ---------------------------
%
% Read in the databases

% SPEX line list
atoms("spex");
plasma("spex");

% new database from CXC (is missing things though)
%%plasma("aped");

%---------------------- Model and Groups ---------------------------
% Set a nominal emmision model a la Capella with 2 (or more) temps
% edit_model;
% - - -
% v0.4 :
%# id    Temp         Dens         Abund        EM          redshift     Nh
%   1    6.3e+06   1.0e+12   1.0e+00   8.0e+45  0.0e+00   1.0e+18
%   2    5.0e+06   1.0e+12   1.0e+00   1.0e+45  0.0e+00   1.0e+18
%
% - - - 
% v0.5 :
%# id    Temp         Dens         Abund        Norm        Vturb       redshift     Nh     
%#       (K)          (cm^-3)                               (km/s)                   (cm^-2)
%   1    6.3000e+06   1.0000e+12   1.0000e+00   8.0000e+45  0.0         0.0000e+00   1.0000e+18
%   2    5.0000e+06   1.0000e+12   1.0000e+00   1.0000e+45  0.0         0.0000e+00   1.0000e+18
%

load_model("capella_manyT.model");
variable mod_minlam, mod_maxlam, mod_binlo, mod_binhi, mod_em;
mod_minlam = 1.0;
mod_maxlam = 181.0;
% define it on 0.002 A bins
(mod_binlo, mod_binhi) = linear_grid(mod_minlam , mod_maxlam,
	(mod_maxlam-mod_minlam)*500);
use_thermal_profile;
%%use_delta_profile;
mod_em = model_spectrum(mod_binlo, mod_binhi);

% Setup groups that I'm interested in

% Define groups
variable N_all, O_all, Ne_all, Mg_all, Si_all, S_all, Ar_all, Ca_all, Fe_all;
variable Fe_16, Fe_17, Fe_18, Fe_19;
variable Fe_22, Fe_23, Fe_24, Fe_25; 
N_all=1; O_all=2; Ne_all=3; Mg_all=4; Si_all=5; S_all=6; Fe_all=7;
Ar_all=8; Ca_all=9;
Fe_16=28; Fe_17=29; Fe_18=30; Fe_19 =31;
Fe_22=34; Fe_23=35; Fe_24=36; Fe_25=37;

% and their colors (hey, the guy is picky!)
variable cN_all, cO_all, cNe_all, cMg_all, cSi_all, cS_all, cFe_all;
variable cAr_all, cCa_all;
variable cFe_16, cFe_17, cFe_18, cFe_19;
variable cFe_22, cFe_23, cFe_24, cFe_25;
cN_all=grey; cO_all=pinkish; cNe_all=orange; cMg_all=yellow; cSi_all=white;
cS_all = purple; cAr_all=orange; cCa_all=yellow;
cFe_all=blue;
cFe_16=purple; cFe_17=ltblue; cFe_18=mediblue, cFe_19 = purpblue;
cFe_22=purple; cFe_23=ltblue; cFe_24=mediblue, cFe_25 = purpblue;


% Fill them
% limit to HETGS wavelength range
define_group(N_all, brightest(10, where( el_ion(N) and wl(1.,35.) )));
define_group(O_all, brightest(10, where( el_ion(O) and wl(1.,35.) )));
define_group(Mg_all, brightest(5, where( el_ion(Mg) and wl(1.,35.) )));
define_group(Ne_all, brightest(10, where( el_ion(Ne) and wl(1.,35.) )));
% Show these only in the short wavelength end
define_group(Si_all, brightest(4, where( el_ion(Si) and wl(1.,10.) )));
define_group(S_all, brightest(1, where( el_ion(S) and wl(1.,10.) )));
define_group(Ca_all, brightest(1, where( el_ion(Ca) and wl(1.,10.) )));
define_group(Ar_all, brightest(1, where( el_ion(Ar) and wl(1.,10.) )));
%
define_group(Fe_all, brightest(50, where( el_ion(Fe) and wl(1.,25.) )));
%
define_group(Fe_16, brightest(5, where( el_ion(Fe, 16) and wl(1.,25.) )));
define_group(Fe_17, brightest(5, where( el_ion(Fe, 17) and wl(1.,25.) )));
define_group(Fe_18, brightest(5, where( el_ion(Fe, 18) and wl(1.,25.) )));
define_group(Fe_19, brightest(5, where( el_ion(Fe, 19) and wl(1.,25.) )));
% show the ~6 keV lines of these...
define_group(Fe_22, brightest(2, where( el_ion(Fe, 22) and wl(1.,3.) )));
define_group(Fe_23, brightest(2, where( el_ion(Fe, 23) and wl(1.,3.) )));
define_group(Fe_24, brightest(2, where( el_ion(Fe, 24) and wl(1.,3.) )));
define_group(Fe_25, brightest(2, where( el_ion(Fe, 25) and wl(1.,3.) )));

% Name the groups
name_group(N_all,"N_all");
name_group(O_all,"O_all");
name_group(Ne_all,"Ne_all");
name_group(Mg_all,"Mg_all");
name_group(Si_all,"Si_all");
name_group(S_all,"S_all");
name_group(Ar_all,"Ar_all");
name_group(Ca_all,"Ca_all");

name_group(Fe_all,"Fe_all");
name_group(Fe_16,"Fe_16");
name_group(Fe_17,"Fe_17");
name_group(Fe_18,"Fe_18");
name_group(Fe_19,"Fe_19");
name_group(Fe_22,"Fe_22");
name_group(Fe_23,"Fe_23");
name_group(Fe_24,"Fe_24");
name_group(Fe_25,"Fe_25");

message("");
message(" Line groups defined:");
message("");
list_group;
message("");