Source code and Mac-OS binaries
This release reflects the changes made for Bell, Mamajek & Naylor.
1) The grid now needs membership probabilities for all the stars.
These can be given by having a colour labelled "Prior" in the input
file, which allows you to have different probabilities for each star.
If such a column is not present then grid will prompt for a value to
apply to all stars.
For details of how to use this feature see Section 4.1 and Appendix C
of Bell, Mamajek & Naylor.
2) The output catalogues grid writes now have the prior probability
used written in, but also the calculated probability of membership from
the stars CMD position, called "Post" (short for posterior probability).
3) The programme to calculate Pr(tau^2), called tau^2 has been
completely re-written with a new algorithm, though the prompts
remain the same.
The old version is still available as old_tau2.
4) Some of the file names have changed, which I hope makes them more
logical (see manual).
5) There is a program to create the cluster format input file from a
simple ascii file, called ascii2cluster.
6) The mass function used by monte is now that of Dabringhausen J., Hilker M., Kroupa P., 2008, MNRAS, 386, 864 (it was a Salpeter function).
Note that the binaries in this version require libgfortran.3.dylib whereas recent updates of gcc (8.X) which come with High Sierra have libgfortran.5.dylib, so you may need to have a copy of the older version. The obvious alternative is to simply recompile the code.
Source code and Mac-OS binaries
This release is primarily to update the available isochrones to include those from Bell et al (2014), and update Example 2 to use them rather than the older Jeffries Pleiades tuning.1) Many of the sub-routines have been re-grouped. 2) A bug which made it impossible to use user isochrones in monte cured. 3) A bug whereby all the singe stars were written out twice in the ascii file which gives the first 10,000 stars from monte cured. 4) The Bolometric luminosity of the Sun now picked up from the BC file headers (if its there) and set to Eric Mamajek's value.
Source code and Mac-OS binariesFor this release large parts of the code which create the isochrones have been re-written, which has resulted in a simpler structure, within which it is easy to add new models. As a result there is a far richer range of options available, yet the code is lightly shorter than the previous release. Old hands, used to previous versions should note the following points.
1) The absolute values of tau^2 have changed. Of course this does not make any practical difference as only relative values used, but there was a conceptual "bug" in which when normalising the isochrone image I did not divide by the pixel size. 2) Monte is structured so that the binary fraction can be a function of primary mass. Now the binary fraction for stars greater than 14.4Mo is 1.5 times the normal binary fraction, and the mass ratio for such stars is distributed as described in Naylor (2009; IAU paper). 3) The extinction vector for the Tycho colours has been improved. The difference between the new values and those in Mayne & Naylor (2008) is small (delta Av = 0.05*E(B-V)) so does not affect the results in Mayne & Naylor or Naylor (2009), but would change things if you had large extinctions. 4) Some more FITS keywords have been added to the output of monte so that both DS9 and the latest version of GAIA give you the co-ordinates in colour-magnitude space. 5) The scripts have been updated and uncer changed slightly to make the code compile on 64-bit machines.
cmdfit.tar.gz_21-08-09cmdfit_mac.tar.gz_21-08-09 Bug found by Phill Cargile fixed, which caused a crash on certain Bessell atmospheres.
cmdfit.tar.gz_27-06-09cmdfit_mac.tar.gz_27-06-09 (1) This release is of the version of the software used in Naylor (2009).
cmdfit.tar.gz_29-10-08Documentation fix for example 2.
To be on the safe side, I switched off the extrapolation of the colour-Teff relationships.
This was a major update in techniques, which were presented in
(1) I use a new normalisation, which retains the mass function in the
grids, but then normalises the grid to sum to one between the brightest and
This avoids the problem of the normalisation going infinite if the isochrone
(2) To generate the uncertainties I now use the uncer method (replacing bootstrap), for which there is a brief description. Its faster and more accurate.
(3) The calculation of the expected value of tau^2 has been improved, and is now done in a program called tau2, which replaces tau. The old estimate became increasingly poor as the number of free parameters increased.
(1) A correction to Example 2.
The extra uncertainty to be added to the colour was given as 0.03, it
should be 0.042.
Both the web-page and the sample files have been updated to reflect this
(1) V vs V-I for the Geneva isochrones.
(2) V vs B-V for the Siess and Baraffe isochrones.
(3) The option to use Bessell et al (1998) colour and bolometric corrections.
(4) Grid now imposes the upper limit on tau-squared to be the tau-squared of the best fitting data point (i.e. that with the lowest tau-squared) plus some value. We do this so that if a different normalisation for tau-squared is used, which effectively adds a constant to all the tau-squared values, the answer remains the same. With the chi-squared-like normalisation the lowest tau-squared value was so close to zero this has little effect on the fits.
(1) The Padova option was using the z=0.004 tracks from the Cioni et
al (2006) papers, not as stated on the web pages Girardi et al (2002).
The code has been changed so the web page is correct, and the web
page now explicitly states its the z=0.019 models.
(1) Corrected a problem whereby if the mass-range you requested outran
that available in the isochrones for model numbers greater than 20
(i.e. those which come ready calibrated in colour and magnitude), you
got some strange low-level effects in the 2D isochrones.
(2) Updated Example 1 so the numbers are correct for the new code (the differences are within the uncertainties).
(3) Added .in files to Example 1, which contain the precise inputs needed to run the programs.
(1) Corrected bundle so it includes files which are soft links on my
(2) Compiled all software statically, to improve portability.