Time-dependent multi-orbital mean-field for fragmented Bose-Einstein condensates