What really should happen is a morphing transition between the sinusoid and the linear hard-clipper-like shape:
You can still horizontally scale the sinusoid up until -1.499 dB
Then switch to the morphing function.
I worked through many solutions after "The Pan Law Discussion".
This really has to be done piecewise and I'll post a fairly simple function that accomplishes what is needed. It works just like a variable knee soft clipper such as GClip. Can be simplified further.
https://www.desmos.com/calculator/xj1ed3hdf6