You cannot find it because there is no "Saturation current" in a real BJT. There will be many mode parameters in an Ebers-Moll model which you will be unable to find in a datasheet.
Also note that there is no fixed point at which a BJT suddenly enters / goes out of saturation. It's more of a gradual thing. This behavior is not in the Ebers-moll model.
It's the same with diodes, you know the formula for a diode I guess, it has this current \$Is\$. It's not in any datasheet. Also, if you would try to measure it you will find that it varies over almost everything like forward current, dopings of the diode, temperature, whatnot. What you do in practice is assume a certain value for \$Is\$ such that we get a realistic forward voltage at a certain forward current.
It is probably the same with that Saturation current in the Ebers-Moll model. When modelling a transistor we tweak it such that it gives a transistor which has a realistic behaviour.