def METHOD_NAME(do):
    # Locate bump
    tStart = 1
    t = do.grid.t[tStart:]
    p = np.zeros(t.shape)
    pCut = -19
    fPar = do.eqsys.f_hot[-1,0,-1,:pCut]
    df = np.diff(fPar)
    maxIdx = np.argmax(df)
    #plt.semilogy(do.grid.hottail.p[:pCut-1], np.abs(df))
    #plt.show()
    bumpLocIdx = np.argwhere(df[maxIdx:]<0)[0][0]
    return do.grid.hottail.p[maxIdx+bumpLocIdx]