debugging equilibrium by offering another gradient computation. In fact, the equil gradient computation through the lpk parabolic fitting is not always the best. An option has been added in the equilibrium options list . Default is using the old deviation grad computation.
In addition, the option for lpk (in equil and modal) has been added to the options list for direct access (instead of reaching them through lpk which was not the best)