_traceback = 1;
plasma(aped);
%require("he_modifier.sl");

create_aped_fun( "Aped_1", default_plasma_state ) ;

define ionpop_modifier(pars,state,last_ionpop){
   variable n = _isis_max_proton_number+1;
   variable new_ionpop = Float_Type[n,n];
   %new_ionpop[*,1]=pars[0];
   %new_ionpop[*,2]=pars[1];
   %new_ionpop[*,-2]=pars[0];
   %new_ionpop[*,-3]=pars[1];
   variable i;
   _for i(2,n-1,1){
      % H-like:
      new_ionpop[i,i-1]=pars[0];
      % He-like:
      new_ionpop[i,i-2]=pars[1];
      new_ionpop[i,*]/=(pars[0]+pars[1]);
   }
   return new_ionpop;
}

create_aped_ionpop_modifier("ionpop",&ionpop_modifier,["xh","xhe"]);

%fit_fun( "Aped_1( 1, ionpop(1), He_triplets(1) )" );
fit_fun( "Aped_1( 1, ionpop(1) )" );
%fit_fun("Aped_1(1)");

%set_par("Aped_1(1).density",1e11,1,0,0);
set_par( "ionpop(1).xh",0.5,1,0,0);
set_par( "ionpop(1).xhe",0.5,1,0,0);

variable atom=10;
variable id_r = where( trans( atom, atom-1,     7, 1 ) ) ;
variable line_inf = line_info( id_r[0] ) ;
variable w = line_inf.lambda ;
eval_fun( w-0.05, w+0.05 );