[T, H] = np.meshgrid(Tp, Hm0)