CODE UPDATED: removed full diagonalization channel from the code.
old code is saved in a suitable branch.

small udpate: lambda_SOC=0.0 modification in the BHZ code.
BUG[solved]: small bug in the printing routine 'print_imp_gf' inside the module
	     ED_GREENS_FUNCTION.f90.
	     The files for the impurity self energies (matsubara and real
	     frequencies) were opened for any value of the input flag
	     'ed_verbose' and were closed only for ed_verbose<3.
	     This caused the increase of the values of  the corresponding
	     unit variable with the number of DMFT-loop.
	     When the number of loop was too large or the solver was used
	     in a real space setup the values of the unit variables
	     reached very fast the upper bound of 900.
	     The problem was not present in the superconducting case.

FIXED 	   : Bug fixed, now the impurity self-energies are printed for
	     ed_verbose<=3.