tp = findtp(data1[iy], min_h)