import seaborn as sns

def boxplot(df, wells, analytes, col_wrap=3):
    data = df[df.location_id.isin(wells) & df.param_name.isin(analytes)]
    g = sns.factorplot(x="location_id", y="analysis_result", 
                       col='param_name', data=data, kind="box",
                       sharey=False, sharex=False, col_wrap=col_wrap)
    g.set_titles("{col_name}")    
    g.set(xlabel='Location ID', ylabel = 'Analysis Result')    
    plt.show()

Create plots.py