One need to check, whether maxFS is not equal
zero. In this case ratio=0 and is a divider
for shearforce, which leads to NaNs.

Clean some headers and indentions.

One need also to implement a "fool proofing" from adding
several equal IGs in InteractionLoop.