I would like to simulate two kinds of shock in the same DSGE model. The first shock aims at generating uncertainty between four kind of policies (therefore a discrete random variable) and the second shock is simply from a normal distribution. I know how to simulate a single non-normal shock: Drawing from non-normal distributions in model simulation
However, I am finding troubles understanding how to make these two kinds of shocks coexist into the dynare simult() function. Should I treat them separately until the random generating process and then stack them into the vector DynareResults.exo_simul?