Filter shocks from actual data

Let me add that you can use the simult_ function for user-define scenarios like large shocks. See the thread at