code: def get_interpolated_params(m_frac, ph, params):    
      
      slope, offset = np.polyfit(np.arange(1,4), params[:3,ph,0], deg=1)
      a   = slope*m_frac + offset
      
      slope, offset = np.polyfit(np.arange(1,4), params[:3,ph,1], deg=1)
      loc   = slope*m_frac + offset
      
      slope, offset = np.polyfit(np.arange(1,4), params[:3,ph,2], deg=1)
      scale   = slope*m_frac + offset
  
      return (a, loc, scale)
docs: Get parameters describing a 1-D capactity distribution 
      for fractional number of modules.