def nse(obs,mod):
    """
    Nash-Sutcliffe-Efficiency = 1 - sum((obs - mod)**2) / sum((obs - mean(obs))**2)
    """
    return 1. - ((obs-mod)**2).sum(axis=0)/((obs-obs.mean(axis=0))**2).sum(axis=0)