Implementation of the callback function

[jiwei08]

Is it able to implement the callback function on u[t,x,y,z] when running op.apply(), like here? I would not like to use save=nt because the size of 4-d arrray u[t,x,y,z] is too large. Any help would be appreciated!