Problem while changing exchange current density

[akariwrynn]

Hello everyone,
I am trying to observe the temperature effect on li-ion batteries behavior and trying to change the current density of the exchange to do it. Since it is a function, I tried to change like this:

from pybamm import *
import numpy as np
from matplotlib import pyplot as plt
import scipy.signal
import math

amplitude=1e-3
period=10
points_per_period=128
frequency_array=np.logspace(-3,2,60)

mychemistry=pybamm.parameter_sets.Chen2020
model=pybamm.lithium_ion.DFN()
params=pybamm.ParameterValues(mychemistry)
solver=pybamm.CasadiSolver()



def pos_electrolyte_exchange(c_e, c_s_surf,T):

    m_ref = 3.42e-6  # (A/m2)(mol/m3)**1.5 - includes ref concentrations
    E_r = 17800
    arrhenius = exp(E_r / constants.R * (1 / 298.15 - 1 / 300))

    c_p_max = params["Maximum concentration in positive electrode [mol.m-3]"]

    return (m_ref * arrhenius * c_e ** 0.5 * c_s_surf ** 0.5 * (c_p_max - c_s_surf) ** 0.5)

params["Positive electrode exchange-current density [A.m-2]"] = pos_electrolyte_exchange

sim = pybamm.Simulation(model, parameter_values=params,solver=solver)
sim.solve(t_eval=(0,3600))

but getting errors like :

pos_electrolyte_exchange() takes 3 positional arguments but 4 were given and other previous 

and this

FunctionParameter(-0x5d8b0a551ed8dfd, Positive electrode exchange-current density [A.m-2], children=['maximum(Positive electrode porosity times concentration [mol.m-3] / Positive electrode porosity, 1e-08)', 'maximum(minimum(boundary value(Positive particle concentration [mol.m-3]), 0.99999999 * Maximum concentration in positive electrode [mol.m-3]), 1e-08 * Maximum concentration in positive electrode [mol.m-3])', 'Maximum concentration in positive electrode [mol.m-3]', 'broadcast(Ambient temperature [K])'], domains={'primary': ['positive electrode'], 'secondary': ['current collector']})

any suggestions on how to solve this?
Also, I want to ask that is it possible to change the Temperature term in all events for example when I change Ambient temperature parameter does it affect both diffusion and kinetic related events? Or I should change the diffusion part like I am trying to change exchange current density?

Thanks,

[rtimms]

The function should take c_e, c_s_surf, c_s_max, T as arguments, i.e.

def pos_electrolyte_exchange(c_e, c_s_surf, c_s_max, T):
    ...

see

PyBaMM/pybamm/input/parameters/lithium_ion/Chen2020_composite.py

Line 6 in cf09d5e

def graphite_LGM50_electrolyte_exchange_current_density_Chen2020(