exp( sqrt(omega0*om*k) * (I_u - 1.0d0) * n*dt )

exp( sqrt(2.0d0*pi*omega0*k/dble(nd)) * (I_u - 1.0d0) * n)
から
exp( sqrt(omega0*om*k) * (I_u - 1.0d0) * n*dt )