This subroutine creates the mesh data structure.
The geometry of the domain is taken from the machine module.
Number of streamlines