Fractional time stepping and adjoint based gradient computation in an inverse problem for a fractionally damped wave equation.