3D version of initial.f90 (IC of mhd wave)