How can I translate this code:

rose_number = np.copy(rosehrs)

for a in range(0,1597):
    for b in range(0,2145):
        for c in range(0,64):
            rose_number[a][b][c] = round((rosehrs[a][b][c]/totalhrs[a][b])*100,2)

Into a numpy function? Maybe using np.divide is an option with where but I’m not sure of the operations needed.