This code facilitates running transient simulations of hysteresis by regulating the rate of change of the forcing as a function of the rate of change of the result.

make clean
make test env=[airaki,manto,...]
./test_hyster.x > tmp

In R, plot the results here:

source("plot_test.r")