ald.py now can run displacement jobs. Virtual crystal added.

ald.py successfully runs with fulldisp=disp. Virtual function averages harmonic and cubic force constants.
Next step is to have ald code run displacement jobs using qepy.