Inversion of phaseless total field data using a two-step strategy based on the iterative multiscaling approach