store the output of all the steps in a array of array
change wave representation of PtM1, Pt, PtP1
encapsulate the stuff into functions
add receivers to record the wave field
full RTM implementation, but can not pass valgrind test