%; Time-stamp: <2002-03-04 10:24:15 dph> 
%; File: Define_model.sl
%; Author: D. Huenemoerder
%; Original version: 2002.02.07
%;====================================================================
%
% purpose:  define a multi-thermal plasma model, w/ non-cosmic abundances.
% based on AR Lac model.  but get rid of cooler bump.

autoload("aped_set_abund",   "aped_fit_models.sl");
autoload("aped_multi_fit",   "aped_fit_models.sl");
autoload("aped_multi_init",  "aped_fit_models.sl");
autoload("aped_multi_setup", "aped_fit_models.sl");

define Define_model()
 {

  variable mdl_temps = [
  3.1623e+05, 3.5481e+05, 3.9811e+05, 4.4668e+05, 5.0119e+05,
  5.6234e+05, 6.3096e+05, 7.0795e+05, 7.9433e+05, 8.9125e+05,
  1.0000e+06, 1.1220e+06, 1.2589e+06, 1.4125e+06, 1.5849e+06,
  1.7783e+06, 1.9953e+06, 2.2387e+06, 2.5119e+06, 2.8184e+06,
  3.1623e+06, 3.5481e+06, 3.9811e+06, 4.4668e+06, 5.0119e+06,
  5.6234e+06, 6.3096e+06, 7.0795e+06, 7.9433e+06, 8.9125e+06,
  1.0000e+07, 1.1220e+07, 1.2589e+07, 1.4125e+07, 1.5849e+07,
  1.7783e+07, 1.9953e+07, 2.2387e+07, 2.5119e+07, 2.8184e+07,
  3.1623e+07, 3.5481e+07, 3.9811e+07, 4.4668e+07, 5.0119e+07,
  5.6234e+07, 6.3096e+07, 7.0795e+07, 7.9433e+07, 8.9125e+07,
  1.0000e+08, 1.1220e+08, 1.2589e+08, 1.4125e+08, 1.5849e+08,
  1.7783e+08, 1.9953e+08, 2.2387e+08, 2.5119e+08, 2.8184e+08,
  3.1623e+08
  ];

  variable mdl_wgts = 1.e-14 * [
	1.6490e+08, 1.9159e+08, 2.2260e+08, 2.5862e+08, 3.0048e+08, 
	3.4910e+08, 4.0560e+08, 4.7124e+08, 5.4750e+08, 6.3611e+08, 
	7.3905e+08, 8.5866e+08, 9.9762e+08, 1.1591e+09, 1.5690e+09, 
	2.2144e+09, 3.0775e+09, 4.0476e+09, 5.0489e+09, 6.0191e+09, 
	6.7495e+09, 7.0727e+09, 7.0825e+09, 7.1557e+09, 8.1084e+09, 
	1.1626e+10, 2.2204e+10, 4.7552e+10, 8.5778e+10, 1.0971e+11, 
	9.4804e+10, 5.8027e+10, 3.0932e+10, 2.3546e+10, 3.3256e+10, 
	5.9225e+10, 9.5934e+10, 1.2512e+11, 1.2942e+11, 1.0919e+11, 
	7.7390e+10, 4.8367e+10, 3.0400e+10, 2.2060e+10, 1.8345e+10, 
	1.6769e+10, 1.6448e+10, 1.6555e+10, 1.6778e+10, 1.7394e+10, 
	1.8255e+10, 1.8292e+10, 1.6109e+10, 1.1762e+10, 7.5169e+09, 
	5.5995e+09, 6.2153e+09, 2.1845e+09, 1.8802e+09, 1.6183e+09, 
	1.3929e+09
   ];


% scale factor, for unabsorbed flux of 1.34e-11 (to match some known
%                      source - empirically determined scale factor)
  variable sfact = 0.854; %  
  mdl_wgts *= sfact;

% set relative abundances; Solar == 1.0; linear scale.
%
  variable elem_abund =[ 0.5,  1.0,  1.0, 3.0,  2.0,  1.0, 1.0,  1.0, 1.0  ];
  variable elem =      [ Fe,   Mg,   N,    Ne,  O,   Si,  S,   Ar,  Ca  ];

% for speedier model evaluation, make the DEM coarser:

  variable ncomps = length(mdl_wgts);
  mdl_wgts = mdl_wgts[ [0:ncomps-1:4] ] * 4.;   % fewer bins, proportionally larger value
  mdl_temps = mdl_temps[ [0:ncomps-1:4] ];
%

  aped_multi_setup( mdl_temps, mdl_wgts );
  aped_set_abund( elem, elem_abund );

}