You can do this using the simult_-function. See https://github.com/JohannesPfeifer/DSGE_mod/blob/master/RBC_news_shock_model/RBC_news_shock_model.mod for an example.
But note that at first order, the IRF to a joint shock is simply the sum of the IRFs to the individual shocks due to linearity.