Recreating sox spectrograms in python

[MJWeldy]

Good morning,
I am working on updating a data pipeline for an acoustic classification project developed in 2016. In 2016, we decided to build our spectrogram generation process using sox. I am now trying to move the project into some APIs which are trying to minimized dependencies. However, to accomplish this I am trying to develop python code to recreate the sox spectrograms in python. I read through spectrogram.c and it appears that there is a gain processes being applied, but I am having trouble putting my finger on the specifics.
I would really appreciate it if someone more familiar with the codebase could give me the FFT and gain specifics needed to recreate the spectrograms in python.
Thanks for your time,
Matt