If the block TDKickFunction is present in the input file, and the variable
TDDeltaUserDefined is not present in the input file, the kick function to
be applied at time zero of the time-propagation will not be a "dipole" function
(i.e. $\phi \rightarrow e^{ikx} \phi$, but a general multipole in the form $r^l Y_{lm}(r)$.
Each line has three columns: integers l and m that defines the
multipole, and a weight. Any number of lines may be given, and the kick will be the sum of those
multipoles with the given weights.
This feature allows calculation of quadrupole, octupole, etc., response functions.