[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: Negative Welch-Stetson Index



Doug Welch wrote:
> 
> On Thu, 13 Jun 2002, jg wrote:
> 
> > I ought to have commented that wsv3.pl does not output for _all_ the
> > input stars, it not listing objects for which the index = 0 for
> > instance.  Thus of the 94,500+ input objects only 33,467 of them had a
> > WS index spat out at the other end.
> 
> John,
> 
> You are correct that stars need to meet a certain set of criteria
> to make into the output. However, the index=0.000 case will make it
> through, as far as I can see. Why do you think otherwise?
> 
> Doug

Hi Doug

A section cut and paste from the wsv3.pl output run in perl on a redhat
linux 7.0 machine and then sorted on the WS index value in Excel under
windblows (so we could always blame excel ;) )  WS index is last figure.

[Damn, it's come out tab delmited, hope that's readable].

> 89987	56	56	155.8739	6.9896	11.36	0.02	10.73	0.02	-0.001
> 90061	56	56	156.0314	5.9898	11.19	0.02	10.79	0.02	-0.001
> 90329	32	32	156.5334	6.82	12.93	0.09	12.03	0.04	-0.001
> 90682	32	32	157.3337	5.4992	12.93	0.1	11.92	0.03	-0.001
> 92748	32	32	179.9254	6.8279	11.76	0.04	11.16	0.02	-0.001
> 3552	51	51	78.9004		5.5784	12.08	0.04	11.51	0.03	0.001
> 17063	92	92	93.9504		5.4796	11.41	0.02	11.11	0.02	0.001
> 26889	76	76	96.4346		8.6043	10.97	0.02	8.97	0.01	0.001
> 28787	38	38	97.1103		8.6351	12.2	0.06	11.77	0.04	0.001
> 32262	38	38	98.4623		8.1558	12.24	0.04	11.26	0.02	0.001
> 69858	56	56	112.1426	6.1565	13.15	0.11	12.55	0.1	0.001
> 80056	32	32	117.6222	9.3306	9.63	0.01	8.67	0.01	0.001

AhAh! Don't know perl at all but the following from the wsv3.pl on the
CD suggests an output of 0.000 is not going to happen :-

>   if ($ws >  0.000 && $ws < 0.001) { $ws =  0.001; }
>   if ($ws > -0.001 && $ws < 0.000) { $ws = -0.001; }

(immediately followed by the printf output statement)

So the 0.000 ones are being reported, but not noted as such?  Needs
0.000 to +/- 0.00049999 conditionals to get $ws = 0?

Cheers

John G.