impulseresponse_yshock.f90 produces the data used in irfyshock.m to make Figures 6 and 7.
This fortran file requires bea_hpy.txt which has the input for the series of shocks.
impulseresponse_wealthshock.f90 plus matlab code produces Figure 8.
Next run impulseresponse_rshock.f90, impulseresponse_Tshockperm.f90 and impulsresponse_purchasesubsidy.f90 to compute the data for irfs_multiple_shock, which makes Fig 9.
Finally, impulse_yshocknof.f90 and the corresponding matlab file make Figure 10.

Note that in all cases, the actual text file output from the fortran code is labeled as impulseresponse.txt, so you will need to rename it with the above naming conventions for import to matlab.
For example, rename the impulsresponse.txt output from impulseresponse_yshock.f90 to impulseresponse_yshock.txt for input to the matlab script.

makedistboombust.f90 creates the durable gaps at different points in time, which are used by computedistributionandhazard99vs09_truemodel.m to make Figure 12.
