Further updatis the output.F90 file. This file is still not finished
and should not compile.
Synchronized output.F90 with the version from the non-factorized branch.
Only change is the use commands:
  use density_module, only: ndens
  use ionfractions_module, only: xh
  use temperature_module, only: temper, temperature_grid
  use radiation_sed_parameters, only: teff,rstar,lstar,S_star

This commit missed one fix, namely the "use evolve" one.
Several minor fixes to pass compilation:

The new module for evolve data is evolve_data:
  use evolve_data, only: phih_grid

names of teff, rstar, lstar variables need "_". However, these are
not being used so they were removed from the use block.

The statement
 use temperature_module, only: temperature_states_dbl
was missing

This module now compiles.