Re: Prob.: timelag calc. using sitar package

Dear Prof. Nowak,

i tried to run isis with '-g' argument but still i get the same error 
message.
my isis ver. is 1.4.9 and sitar ver. is 0.8.0

about the sizes of diff. variables , here is check,
----------------
isis> printf("%d %d %d %d %d %d\n", length(aflo1), length(apsd1),
               length(apsd2), length(acpd), length(npsd1), length(npsd2));
36 36 36 36 1 1
---------------

one more thing , i tried to see in sitar , where the error message is
being generated. and i find that in sitar_lags() routine, in the 
begining it tries to check the length of all the arguments to be equal or 
not. it fails there, possibly for 'navg' variable.

i have attached  my lc files and script(list of commands) for your 
convenience.

regards,
nilay

>
>
>> i am trying to calculate timelag between lightcurves from two diff. energy bands for one XMM-Newton data set.
>> I give the following command at isis prompt and it quits with
>> following error message.
>> -----------------------
>> isis> (lag,dlag,g,dg) = sitar_lags(aflo1,apsd1,apsd2,acpd,npsd1,npsd2);
>> Inconsistent array sizes
>> Stack Underflow Error
>> -----------------------
>> here, i have not supplied one "optional" variable (navg - no. of averaged segments).
>> now if I supply that optional variable too, i will still get the same
>> error.
>>
>> i dont know how to get around it. Pl. help.
>
> I'll need a little more help to figure this out.  Run isis as:
>
> unix%> isis -g
>
> you'll get a little more debug info that way, like the line in the code that actually causes the problem.  Also, it would be good to know the lengths of the arrays in the input, to see that they are indeed all the same length.  Could be a bug in my code, or that something got off in the set up, but there's not enough to go on at the moment.
>
> So, send me the info, and we'll work from there.
>


----
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