Abstract: The McMillan-Ginzburg-Landau (MGL) model for charge density waves (CDW) is employed in a systematic phenomenological study of the different phases that have been probed in recent experiments involving transition metal dichalcogenides. We implemented an efficient imaginary time evolution method to solve the MGL equations, which enabled us to investigate the role of different coupling parameters on the CDW patterns and to perform calculations with different energy functionals that lead to several experimentally observed singularities in the CDW phase profiles. In particular, by choosing the appropriate energy functionals, we were able to obtain phases that go beyond the well-known periodic phase slips (discommensurations), exhibiting also topological defects (i.e., vortex-antivortex pairs), domain walls where the CDW order parameter is suppressed, and even CDW with broken rotational symmetry. Finally, we briefly discuss the effect of these different CDW phases on the profile and critical temperature of the competing superconducting state.