Re: problem using model_spectrum and brightest

Hi David:

> 
>     juan> but when I type 
> 
>     juan> g=brightest(10,where(wl(10,12)));
> 
> 
> You might not have evaluated the model, or evaluated
> it over the
> region of interest.  You can try:
> 
> 
>   notice(all_data);
>   eval_counts;
>   g=brightest(10,where(wl(10,12)));
> 
> 
> or try:
> 
>   y = eval_fun( 1.0, 40.0 ) ; 
>   g=brightest(10,where(wl(10,12)));
> 
> 
> If that doesn't work, I'll need more details on the
> state at which
> this doesn't work.

It didn't work. I give you more details:

I define a model:

# id   Temp	Density    Abund    Norm    Vturb   
redshift
#      (K)	(cm^-3) 		    (km/s)
  1    3.e6	1.e-3	   1.0      0.0008635055     10.0   
0.0
#
# abundances for component 1:
#     Si = 0.8	Na = 0.9   Fe = 0.3
#
  2    20.e6	1.e-3	   1.0      0.001866956    10.0   
0.0


then I loaded a dataset, and evaluate the model as you
suggested:

>load_model("my_model.txt")
>d=get_data_counts(4);
>flux=model_spectrum(d.bin_lo ,d.bin_hi);
>eval_counts; %or y=eval_fun(12.0,17.5);
 Parameters[Variable] = 26[21]
            Data bins = 6604
           Chi-square = 4925.318
   Reduced chi-square = 0.7481875
0

>g=brightest(5, where (wl(12.0,17.5)));

and still I'm getting
>all lines have flux = 0


I'm fitting a function to the data before doing this,
it can be that the problem?

Thanks
juan


> 
>     juan> I'm following the example on page 23 of
> the manual to
>     juan> identify the brightest lines in a model
> spectrum which
>     juan> I define as:
> 
>     juan> # id    Temp         Dens         Abund   
>     Norm   
>     juan> Vturb       redshift
>     juan> #       (K)          (cm^-3)              
>            
>     juan> (km/s)              
>     juan> 1	3.14e6	     1.e-3	  1.0	    1.0		  
> 200.0	0.0
>     juan> 2	79.7e6	     1.e-3	  1.0	    0.02960599  
>   200.0    
>     juan> 0.0
>     juan> 3	7.56e6	     1.e-3        1.0      
> 0.003366989   
>     juan> 200.0        0.0
>     juan> 4	33.4e7	     1.e-3        1.0      
> 0.02999968    
>     juan> 200.0        0.0
> 
>     juan> but when I type 
> 
>     juan> g=brightest(10,where(wl(10,12)));
> 
>     juan> I got
> 
>     juan> "all lines have flux = 0"
> 
>     juan> I'm reading just the pha2 files as in the
> example.
> 
>     juan> Thanks
>     juan> Juan
> 
>     juan>
>
########################################################
>     juan> GERARDO JUAN MANUEL LUNA
>     juan> Eberly College of Sciences
>     juan> Penn State University
>     juan> 525 Davey Lab
>     juan> University Park, PA 16802, USA
>     juan>
>
########################################################
> 
> 
>     juan> Tarjeta de crédito Yahoo! de Banco
> Supervielle.
>     juan> Solicitá tu nueva Tarjeta de crédito. De
> tu PC directo a tu casa. www.tuprimeratarjeta.com.ar
> 
>     juan> ----
>     juan> You received this message because you are
>     juan> subscribed to the isis-users list.
>     juan> To unsubscribe, send a message to
>     juan> isis-users-request_at_email.domain.hidden>     juan> with the first line of the message as:
>     juan> unsubscribe
> 
> 
> 
> ----
> 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.hidden> with the first line of the message as:
> unsubscribe
> 
> 


########################################################
GERARDO JUAN MANUEL LUNA
Eberly College of Sciences
Penn State University
525 Davey Lab
University Park, PA 16802, USA
########################################################


      Yahoo! Encuentros.

Ahora encontrar pareja es mucho más fácil, probá el nuevo Yahoo! Encuentros http://yahoo.cupidovirtual.com/servlet/NewRegistration
----
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