code: def cell_thermal_mass(temperature, conductivity):
      
  
      alpha = 0.03  # Thermal anomaly amplitude.
      beta = 1.0 / 7  # Thermal anomaly time constant (1/beta).
  
      sample_interval = 1 / 15.0
      a = 2 * alpha / (sample_interval * beta + 2)
      b = 1 - (2 * a / alpha)
      dCodT = 0.1 * (1 + 0.006 * [temperature - 20])
      dT = np.diff(temperature)
      ctm = -1.0 * b * conductivity + a * (dCodT) * dT  # [S/m]
      return ctm
docs: Sample interval is measured in seconds.
      Temperature in degrees.
      CTM is calculated in S/m.