code: def linterp(self, setx, sety, x):
          
          if math.isnan(sety[0]) or math.isnan(setx[0]):
              return np.nan
          #if math.isnan(sety[0]):
          #    sety[0] = 0.
          #if math.isnan(sety[1]):
          #    sety[1] = 0.
          return sety[0] + (x - setx[0]) * ( (sety[1]-sety[0]) / (setx[1]-setx[0]) )
docs: Linear interp of model data values between time steps