#computation of the phonon dispersion
PHON_NWRITE = -2
LPHON_DISPERSION = .TRUE.
PHON_DOS=1
LPHON_READ_FORCE_CONSTANTS=.TRUE.
LPHON_POLAR=.TRUE.
PHON_DIELECTRIC = \
    3.33662190      0.00000000      0.00000000\
   -0.00000000      3.33662190     -0.00000000\
   -0.00000000     -0.00000000      3.33662190
PHON_BORN_CHARGES = \
    2.00873216      0.00000000      0.00000000\
    0.00000000      2.00873216      0.00000000\
    0.00000000      0.00000000      2.00873216\
   -2.00873216      0.00000000      0.00000000\
    0.00000000     -2.00873216      0.00000000\
    0.00000000      0.00000000     -2.00873216
