Changed the assignment to phi%heat at the end of heat_lookuptable to
be cumulative.
Removed spurious writing of average density to standard output.
Added debugging output to be activated with the -DMPILOG option
for the preprocessor
Module to report memory usage
Added the report_memory module file to the list of files to be compiled.
Made proper series for test_kyl programs.
Edited number of redshifts and time step, for testing.
Removed unnecessary executables from the repository
Output the grid of heating rates (phiheat)

In evolve needed to add initialization to zero and phiheat to
the list of public entities of the module
For debugging changed to "short" version of 3 time steps of 1 Myr.
Serious bug: using undefined variable NFlux when calling the heat_lookuptable.
This should have been NormFLux(nsrc) and NormFluxPL(nsrc).

Removed unused variable NFlux and replaced it in the call to heat_lookuptable
with the correct variables.
Hunting down a bug which manifested itself as a different result between
OMP and no OMP version of the code.

It turned out be due to the scaling factors in the table lookup functions.
These are now calculated in photoion_rates and passed to the lookup
functions.

In the process many write statements were added and commmented out.