Fix unit bug when upgrading numpy to > 1.21

[Adnan-Ali-Ahmad]

Hi,

I've upgraded my numpy version to '1.23.3' and this is causing some weird behavior when dividing Arrays/Vectors.
See the example below:

data["hydro"]["velocity"] / data["amr"]["position"]
Out:
'' Min: 1.279e-14 Max: 1.712e-11 [cm / s] (340488,), {x,y,z}

The units are not being properly converted. The reason this happens is because np.divide.__name__ no longer returns 'true_divide', but now returns 'divide'.

Fixing this bug simply requires replacing the string in APPLY_OP_TO_UNIT.

[Adnan-Ali-Ahmad]

On second thought, it's better to also keep "true_divide" in APPLY_OP_TO_UNIT for compatibility with older versions of numpy.

[nvaytet]

On second thought, it's better to also keep "true_divide" in APPLY_OP_TO_UNIT for compatibility with older versions of numpy.

Yes, that's what I was about to comment.