Update volume phase function

I have a Blender environment setup with a Mix Shader combining Volume Absorption and Volume Scatter nodes to render a volume in the Cycles render engine. I am trying to change the volume Henyey Greenstein phase function for experiments. In the source code, this corresponds to the single_peaked_henyey_greenstein() function in volume.h: https://github.com/blender/blender/blob/blender-v4.0-release/intern/cycles/kernel/closure/volume.h

I need more than the anisotropy parameter g for the experiments so possibly add more arguments to the function.

What is the best approach for this? Is this possible to change the phase within Blender or should I modify the source code? Ideally I would also be able to change the additional parameters in the shader editor.