A new time-integration method is presented for a dislocation-density-based CP model. An implicit method with the consistent tangent moduli is implemented. The computational procedure is validated with experiments of aluminum materials. Indentation of single crystals and simple shear of the multicrystal are performed.