def docontours_all(mc_lognormal, mc_simple, mc_lognormal_freemach):
    docontours_multi(
        mc_lognormal,
        start=10000,
        savename=savepath +
        "mc_lognormal_withmach_multipanel_abundance%s_opr%s.pdf" %
        (abundance, opr),
        dosave=True,
        parnames=('tauoneone_mu', 'tautwotwo_mu', 'meandens', 'sigma', 'mach',
                  'b', 'tau_ratio'))
    #docontours_multi(mc_lognormal,start=10000,savename=savepath+"mc_lognormal_withmach_multipanel_deviance.pdf", dosave=True,
    #                 parnames=('tauoneone_mu','tautwotwo_mu','meandens','sigma','mach','b','tau_ratio','deviance'))
    #docontours_multi(mc_lognormal,start=10000,savename=savepath+"mc_lognormal_withmach_multipanel_scalefactors.pdf", dosave=True,
    #                 parnames=('tauoneone_mu','tautwotwo_mu','meandens','sigma','mach','b','tau_ratio','Metropolis_sigma_adaptive_scale_factor','Metropolis_b_adaptive_scale_factor','Metropolis_meandens_adaptive_scale_factor','deviance'))
    #docontours_multi(mc_lognormal,start=10000,savename=savepath+"mc_lognormal_withmach_multipanel_deviance.pdf", dosave=True,
    #                 parnames=('tauoneone_mu','tautwotwo_mu','meandens','sigma','mach','b','tau_ratio','deviance'))
    docontours_multi(
        mc_simple,
        start=10000,
        savename=savepath +
        "mc_lognormal_justtau_multipanel_abundance%s_opr%s.pdf" %
        (abundance, opr),
        dosave=True,
        parnames=('tauoneone_mu', 'tautwotwo_mu', 'meandens', 'sigma',
                  'tau_ratio'))

    docontours_multi(
        mc_lognormal_freemach,
        start=10000,
        savename=savepath +
        "mc_lognormal_freemach_multipanel_abundance%s_opr%s.pdf" %
        (abundance, opr),
        dosave=True,
        parnames=('tauoneone_mu', 'tautwotwo_mu', 'meandens', 'sigma', 'mach',
                  'b', 'tau_ratio'))
def docontours_all(mc_hopkins_freemach, mc_hopkins, mc_hopkins_simple):
    docontours_multi(
        mc_hopkins_freemach,
        start=10000,
        savename=savepath + "mc_hopkins_freemach_multipanel_abundance%s_opr%s.pdf" % (abundance, opr),
        dosave=True,
        parnames=("tauoneone_mu", "tautwotwo_mu", "meandens", "sigma", "tau_ratio", "mach", "Tval", "b"),
    )
    # docontours_multi(mc_hopkins_freemach,start=10000,savename=savepath+"mc_hopkins_freemach_multipanel_deviance.pdf", dosave=True,
    #                 parnames=('tauoneone_mu','tautwotwo_mu','meandens','sigma','tau_ratio','mach','Tval','b','deviance'))

    docontours_multi(
        mc_hopkins,
        start=10000,
        savename=savepath + "mc_hopkins_withmach_multipanel_abundance%s_opr%s.pdf" % (abundance, opr),
        dosave=True,
        parnames=("tauoneone_mu", "tautwotwo_mu", "meandens", "sigma", "tau_ratio", "Tval", "b"),
    )
    # docontours_multi(mc_hopkins,start=10000,savename=savepath+"mc_hopkins_withmach_multipanel_deviance.pdf", dosave=True,
    #                 parnames=('tauoneone_mu','tautwotwo_mu','meandens','sigma','tau_ratio','Tval','b','deviance'))
    docontours_multi(
        mc_hopkins_simple,
        start=10000,
        savename=savepath + "mc_hopkins_justtau_multipanel_abundance%s_opr%s.pdf" % (abundance, opr),
        dosave=True,
        parnames=("tauoneone_mu", "tautwotwo_mu", "meandens", "sigma", "tau_ratio", "Tval"),
    )