Thank you @mhscott for the encouragement
I managed to solve the problem by using the "parameter" and "updateParameter" commands.
Inspiration for this was one of your blog posts, and I did end up reading the C++ source to find out the arguments to the parameter command.
Here is an excerpt for future reference:
Code: Select all
...
ops.timeSeries('Constant', id_timeseries)
ops.pattern('Plain', id_pattern, id_timeseries)
ops.load(n, *[0.001, 0.0, 0.0])
ops.parameter(n, 'pattern', id_pattern, "loadAtNode", n, '1')
...
for step in range(n_steps):
...
f = slamming_force(Vrel, depth, Rc, chord)
ops.updateParameter(id_pattern, f + 0.1)