https://docs.google.com/document/d/1ysHpXGzfprZuVg0V-CMAQqFnVbBuWaN48931F4Zn8JE/edit?pli=1#heading=h.fyltvlf8wds2 impulse response function fortran