How to use time 't' in a current function?

[cyasing]

I tried programming a current function like this:

def controlled_current(one_c, c_rate):
    def current(t):
        return float(one_c if t <= 1000 else (-1*one_c*c_rate[int(t-1001)]))    
    return current

where the c_rate will be an array of c-rates over a hundred seconds or so. I passed it as an input to the solver via current function, but it doesn't work due to not being able to use time "t" as a number (All it returns is a Time object). What should I do? The way I input is through the param["Current function [A]"] = controlled_current(one_c, c_rate).

[cyasing]

Edited the question based on changes I made that worked for simpler trial of function.

[RuiheLi]

Hi~ Yes "t" is by default a variable in PyBaMM. I think the best reference is this notebook, where you can find a drive cycle section and custom current section.
Also, it may be very helpful if you can kindly post your full scripts here, or explain a bit more how you are trying to control the current? Thanks!

[cyasing]

Turns out that the only way to access time dependent current is through step simulation.

[cyasing]

In case a function requires to make time-based decisions.

[rtimms]

You can use < and > to define a heaviside function. E.g.

def current(t):
    return t < 100

would return a current of 1 when t is less than 100 seconds and 0 otherwise. You can combine as many as you like to make the function you want.

[rtimms]

Alternatively you could make an experiment. If you have an array of rates you could make the experiment string programmatically, e.g. construct it in a loop.

[cyasing]

I did go for the second option in the end. I'll see if the simpler expression works.