def NMPB(obs, mod, paxis, axis=None):
    """ Normalized Mean Peak Bias (%)"""
    modmax = mod.max(axis=paxis)
    obsmax = obs.max(axis=paxis)
    return (modmax - obsmax).mean(axis=axis) / obsmax.mean(axis=axis) * 100.0
