An easy way to calculate P and S seismic travel times based on 1D velocity model. It's based on PyKonal
conda env create -f ptt.yml
conda activate ptt
git clone https://github.com/malcolmw/pykonal mypath/to/pykonal
cd mypath/to/pykonal
pip install .
See the example.