Add routine to calculate and plot transient reflectivity. Replaces !179 (closed) and closes #49 (closed)