ComputeGIRFs function

Compute GIRFs for all models