Update SS316L.py to use the energy_density_total_enthalpy model and the new implementation of prepare_interpolation_arrays

Add bounds check to prevent out-of-bounds access. Ensure idx_E_neq can't exceed size-2 for safe interpolation.

Modify loop condition to ensure left and right indices remain adjacent. Change index updates to maintain proper ordering. Use back() method for consistent array access.

Update InterpolationArrayContainer to conditionally include appropriate methods, add unified interpolate() menthod

Refactor prepare_interpolation_arrays to support both binary search and double lookup methods with automatic selection.