DIPOLE and Joint DOS calculation

The latest version of code in the repository (172) now produces equal results for dosjnt.ftn and excite3-orig.f after correcting the smearing factor (it should be set to what ever is used in excite), the expression is:
FUNC=BACKW*VFAC*EXP(-((ERG+EDIFF)/0.01)**2)
with VFAC set to 1.
Also the dipole moment has been checked with what is calculated in jntdos.ftn DOS against apotnl.ftn
The electronic dipole moments are equal. However, the result given in the DIPOLE file (written in apotnl.ftn) has already taken nuclear contributions into account.