Re: usage of apedx model

Hi Juan, 


    juan> I'm trying to use apedx to model an spectral line. I'm
    juan> following the example on LSF_script.sl from David
    juan> Huenemoerder:


    juan> But I'm getting only continuum fitting and no lines.

    juan> If I ask isis for the value of Aped_Filter I got:

    isis> print(Aped_Filter);
    juan> {type=NULL, lines=NULL}

That looks like a bug - it should evaluate lines and continuum for
type=NULL.  I will have to dig into that.


I'd like to fast-forward from that script to a more-modern method,
since the aped plasma model functions have improved a lot since I
wrote that.  It might do what you want, but there might be an easier
way.

Tell me if this suffices:

p = default_plasma_state;

p.temperature = 10^[6.8, 7.2, 7.5]; % use 3T; adjust for your application.
p.norm = p.temperature * 0.0 + 0.001 ;  % norms for each T
p.elem = [N, O, Ne, Fe ] ; % for any non-solar mix ; add more, or fewer.
p.elem_abund = [ 1.0, 0.5, 1.5, 0.2 ] ;  % relative to solar for the above.

create_aped_fun( "Aped_3", p ) ; % create a model function (the name is arbitrary)

fit_fun( "Aped_3(1)" ) ;   % instantiate the functiona
list_par;                  % scan the parameters

% and you're ready to go.


If you want arbitrary line and/or continuum evaluation, then you need
to use some hooks - see the create_aped_fun help, or I have some
replacement functions (but largely undocumented).



-- Dave

David Huenemoerder  617-253-4283 (o); -253-8084 (f); http://space.mit.edu/home/dph
MIT Kavli Institute for Astrophysics and Space Research
70 Vassar St., NE80-6065,
Cambridge, MA  02139
[Admin. Asst.: Elaine Tirrell, 617-253-7480, egt_at_email.domain.hidden
----
You received this message because you are
subscribed to the isis-users list.
To unsubscribe, send a message to
isis-users-request_at_email.domain.hiddenwith the first line of the message as:
unsubscribe