The user is now allowed to define how many body forces they want to compute, whereas one body force could be assembled by 1 or more surface ids in the mesh. The user then gives each body force a unique name, which will be dumped as such in the stat file, e.g. force_uniquebodyforcename%1 is the x component of the body force which has the name uniquebodyforcename.