Fixed the use of temperature in photon statistics, subroutine total_rates.
The code was still using the temper variable which has become undefined.
This can lead to divbyzero errors.
SHould use get_temperature_point. This has now been added to this
subroutine.
Removed the use of the variable temper from the output routines.

This module now uses the new temperature structures.
The quantity sum_xh0_int was not being calculated correctly for the non ALLFRAC
case. It should not be 1.0 - sum_xh1_int, but real(mesh(1)*mesh(2)*mesh(3)) - sum_xh1_int.
The global average fractions are not normalized during this calculation hence
this factor instead of 1.0
This is now fixed.
Fixed bug in handling restarts from iteration dumps. On Aug 16, 2013
the mpi_bcast of the phih_grid array was erroneously removed. This
has now been restored.

Additionally the two quantities used in checking convergence of the
global ionization fractions (prev_sum_xh1_int, prev_sum_xh0_int)
are being mpi_bcast.
Merge branch 'master' of github.com:garrelt/C2-Ray3Dm into factorization

This may be wrong...

Conflicts:
	output.F90
Not sure about these... I think they were unchanged
Modified to use the correct factorized modules of the former material
module.
Adjusting to the new factorized material modules.