dynare thesimyp1.mod;
irf1 = oo_.irfs;
save thesimyp1.mat irf1;

dynare thesimyp2.mod;
irf2 = oo_.irfs;
save thesimyp2.mat irf2;

load('thesimyp1.mat','irf1');
load('thesimyp2.mat','irf2');

ending_cell={'_eps_a','_eps_r','_eps_g','_eps_tau'};
    
%'Hiato do produto', 'Inflação Doméstico', 'Consumo Doméstico',
%'PIB','Inflação IPC','Termo de Troca', 'Salário Real', 'Horas de
%Trabalho','Dívida','Câmbio Nominal','Balanço Comercial','Taxa de Juros Real'
    

for ii=1:length(ending_cell)
    HOR=1:options_.irf;
    var_titles  = {'Hiato do produto', 'Inflação Doméstico', 'Consumo Doméstico', 'PIB', 'Inflação IPC', 'Termo de Troca', 'Salário Real', 'Horas de Trabalho','Dívida', 'Câmbio Nominal','Balanço Comercial','Taxa de Juros Real'};    
    figure
    for jj=1:length(var)
        subplot(3,3,jj)
        eval(['irf1.' var{1,jj},ending_cell{1,ii}]);
        eval(['irf2.' var{1,jj},ending_cell{1,ii}]);
        hold on
        plot(HOR,[eval(['irf1.' var{1,jj},ending_cell{1,ii}])],'-k',HOR,[eval(['irf2.' var{1,jj},ending_cell{1,ii}])],'--r','LineWidth',2);
        title([var_titles{1,jj}] )
    end
    
    
    legend('modelo hp', 'modelo tc', 'AutoUpdate','off')
end