How to change the electrolyte_exchange_current_density function?

[fangyuer]

Hello, everyone! I'd like to know if I can replace the Function Parameter like below.

def graphite_LGM50_electrolyte_exchange_current_density_Chen2020(c_e, c_s_surf, c_s_max, T):

As we all know, when we calculate the exchange current density, we need to set the Electrochemical Reaction Constant. However, I don't know the precise value of them. So I want to use the parameter estimation method to get them. But I don't know how to set the parameter "m_ref" as the input of the exchange current density function, like below.

m_ref = 6.48e-7 # (A/m2)(m3/mol)**1.5 - includes ref concentrations

Whether I can set the Electrochemical Reaction Constant as the parameter of the model, such as the Negative/Positive current collector thickness. It then, make the Electrochemical Reaction Constant be used as the input parameter of the exchange current density function.

Please teach me how to realize it! Thank you!