diff --git a/source/stellar_astrophysics.tracks.file.F90 b/source/stellar_astrophysics.tracks.file.F90 index 3033da3d7b..6fdbd0bb78 100644 --- a/source/stellar_astrophysics.tracks.file.F90 +++ b/source/stellar_astrophysics.tracks.file.F90 @@ -445,7 +445,7 @@ subroutine fileInterpolationCompute(self,initialMass,metallicity,age,interpolati interpolationFactorsMass(iMetallicity,:)=[0.0d0,1.0d0] massOutOfRange=.true. else - call self%interpolatorMass(iMetallicity)%linearFactors(initialMass,interpolationIndicesMass(iMetallicity,1),interpolationFactorsMass(iMetallicity,:)) + call self%interpolatorMass(jMetallicity)%linearFactors(initialMass,interpolationIndicesMass(iMetallicity,1),interpolationFactorsMass(iMetallicity,:)) interpolationIndicesMass(iMetallicity,2)=interpolationIndicesMass(iMetallicity,1)+1 end if ! Loop over masses.