Fixed a bug (angle[] is only for 2D simulation).
- When Pe is smaller than Pe_switch, some parameters are scaled.
- Introducing zero_shear for the relaxation process to generate initial configuration.
  Otherwise, it can generate configuration for sheared periodic box, which we cannot use for initial configuration.