An easy to use immersed boundary method in 2D, with full implementations in MATLAB and Python that contains over 75 built-in examples, including multiple options for fiber-structure models and advection-diffusion, Boussinesq approximations, and/or artificial forcing.
The main function in IBM_Driver.py only accepts six parameters. Why is the main function in line 160 of pyIB2d called with seven parameters?