epsinvI_spin.f90 adjusted for epsinvI calculation