update:
added driver for 2band calculations at fixed density.
small bug fixed for parallel calculations, i.e. use ED_MPI_ variables
in DMFT_ED.f90