 First read in the best fitting model FITS file.
Give input file ("?" gives help)>  Reading file best_model.fit
 No gradient file found, using natural normalisation.
 > Give the file of data points (or type end for chunk of image).
 The available colours are
 2 B-V       
 3 U-B       
 Give the number of the colour do you want?
 Will add V to B-V to create B
 > Give the number of free parameters.
 Doing star  2
 Modal tau is  22.52535  for star  1  out of  16
 Doing star  4
 Modal tau is  45.52646  for star  2  out of  16
 Doing star  5
 Modal tau is  69.120766  for star  3  out of  16
 Doing star  6
 Modal tau is  92.25028  for star  4  out of  16
 Doing star  9
 Modal tau is  115.207634  for star  5  out of  16
 Doing star  11
 Modal tau is  137.86983  for star  6  out of  16
 Doing star  12
 Modal tau is  161.2609  for star  7  out of  16
 Doing star  13
 Modal tau is  184.36424  for star  8  out of  16
 Doing star  14
 Modal tau is  207.28773  for star  9  out of  16
 Doing star  15
 Modal tau is  230.5095  for star  10  out of  16
 Doing star  16
 Modal tau is  252.72806  for star  11  out of  16
 Doing star  17
 Modal tau is  276.76636  for star  12  out of  16
 Doing star  18
 Modal tau is  299.81647  for star  13  out of  16
 Doing star  24
 Modal tau is  323.8882  for star  14  out of  16
 Doing star  31
 Modal tau is  347.52393  for star  15  out of  16
 Doing star  32
 Modal tau is  372.18787  for star  16  out of  16
 The sigma of the distribution is  5.5106964 which should
 be approximately the square root of twice the number of free 
 parameters  15 i.e.  5.477226
  
 Reading best-fit tau^2 value from grid.fit
 Reading file grid.fit
 Pr( 381.7701 ) is  0.029761888  for  1  free parameters.
