subroutine output_ftlm() integer :: i real(dp) :: beta if (iProcIndex /= root) return write(stdout, '(1x,a4,18X,a11,11x,a11)') "Beta", "E_numerator", "Denominator" beta = 0.0_dp do i = 1, nbeta_ftlm + 1 write(stdout, '(es17.10,5x,es17.10,5x,es17.10)') beta, ftlm_e_num(i), ftlm_trace(i) beta = beta + delta_beta_ftlm end do end subroutine output_ftlm