The paper will present the framework of fully automated two/three dimensional ice accretion simulation package NSMB3D-ICE, with emphasis on the remeshing step. The NSMB3D-ICE Navier-Stokes code, coupled to an Eulerian droplet module and iterative Messinger thermodynamic model, can perform multi time-steps ice accretion simulations via an automated multi-block elliptic/parabolic grid generation code (NSGRID3D). Attention is paid to the efficiency and robustness of the numerical procedure especially for complex 3D glaze ice simulation. The new automated multi time-step icing code NSMB3D-ICE/NSGRID3D is validated and verified using several icing case studies such as the GLC305-3D rime and glaze ice cases.