changed compute_temperature.F90