Error with pybamm base thermal model source code

Hi,

Just thought I’d mention that changing the parameter:

Negative electrode OCP entropic change [V.K-1] to a non-zero value gives the error:

ShapeError: Cannot find shape (original error: operands could not be broadcast together with shapes (20,1) (11,1) )

when using the thermal lumped model. This stems from the reversible heat loss calculation: Q_rev_n += a_j_n * T_n * dUdT_n

Does anyone have any suggestions to why this occurs or some potential fixes to the source code ?

All the best,

Alex