(1.0 - tmp1 / kw1 * np.tanh(2.0 * kw1 * h)) +