propagate a predefined python surface over time
Here is the blend file with the modified script. Basically, you have to:


create a few shape keys, with incremented angle given to sin function
uncheck the Relative box for basis shapekey
make the interpolation type linear
set the value of evaluation time for the end frame one less than the maximum evaluation time