init_timing Subroutine

    subroutine init_timing()
        != Start global timer for timing the total calculation time.

        integer :: i

        if (time_at_all) then
            global_time = MPI_WTIME()
            global_timing_on = .true.

            if (.not. allocated(timers)) allocate (timers(ntimer))

            do i = 1, itimer
                ! Have already done one run if itimer>0.  Clear existing timing info.
                timers(i)%ncalls = 0
                timers(i)%time_cpu = 0.0_dp
                timers(i)%sum_time_cpu = 0.0_dp
                timers(i)%timing_on = .false.
            end do
        end if

    end subroutine init_timing