Asymmetric impulse responses

I am looking to progam a model which has non linear dynamics. In particular, i wish to explore the differences in response to positive and negative shocks to the same variable. What is the best way to do this in dynare? (second or third oder approximation in stoch_simul, or are there better alternatives?)

Please see the RBC_state_dependent_GIRF.mod on my homepage.