I want to find the position in a 3D numpy array where a 2D array first exceeds the value in the 3D array. E.g.
import numpy as np
array_3D = np.random.rand(10, 3, 3)
array_2D = np.random.rand(3, 3)
np.argmax(array_2D > array_3D)
I get the result 0, but I want a 2D array i.e. essentially argmax results for each grid cell.