This paper presents the algorithm for determining the Lemaitre-Tolman (LT) model that best fits given datasets for maximum stellar ages, and SNIa luminosities, both as functions of redshift. It then applies it to current cosmological data. Special attention must be given to the handling of the origin, and the region of the maximum diameter distances. As with a previous combination of datasets (galaxy number counts and luminosity distances versus redshift), there are relationships that must hold at the region of the maximum diameter distance, which are unlikely to be obeyed exactly by real data. We show how to make corrections that enable a self-consistent solution to be found. We address the questions of the best way to approximate discrete data with smooth functions, and how to estimate the uncertainties of the output - the 3 free functions that determine a specific LT metric. While current data does not permit any confidence in our results, we show that the method works well, and reasonable LT models do fit with or without a cosmological constant.

We consider the spherically symmetric metric with a comoving perfect fluid and non-zero pressure -- the Lemaitre metric -- and present it in the form of a calculational algorithm. We use it to review the definition of mass, and to look at the apparent horizon relations on the observer's past null cone. We show that the introduction of pressure makes it difficult to separate the mass from other physical parameters in an invariant way. Under the usual mass definition, the apparent horizon relation, that relates the diameter distance to the cosmic mass, remains the same as in the Lemaitre-Tolman case.