Monte Carlo Simulation after DSGE Estimation

Hello,
I am using the Smets-Wouters 2007 model for a project where I start by using the estimated model to simulate 1000 datasets for a Monte Carlo simulation. I will then fit VARs and LPs to each of these datasets and look at the IRFs and bias and variance compared to the original model to see which method is best. I know stoch_simul can simulate datasets but also solves models (which I don’t think I need), but can it be used after estimation() so that it uses the estimated model for the simulation? Are there any additional options that may be needed? Or should I be doing a different (manual) procedure following the estimation() to get the simulated datasets? I’ve seen simult- as another option. Just want to make sure the simulations are correct–hope this is an easy question.
Thank you.

See

You may want to use the simul_replic-option to obtain the datasets. See e.g. DSGE_mod/Hansen_1985/Hansen_1985.mod at master · JohannesPfeifer/DSGE_mod · GitHub