Hi, Lahcen

The projection of d-orbitals are calculated by means postproccesing tools (projwfc.x in PP dir). You must make scf calculations and "band" calculation using pw.x. After that, you can plot the bands using projwfc.x, with the add of use kresolveddos=.true. Be careful with the k-points grid in band calculation and spin-orbit coupling (use non-collinear calculation).

An alternative is the use of thermo_pw package (see example 19 for Al). This suite use quantum-espresso and plot the orbitals. The thermo file need the inclusion of

what='scf_2d_bands',

lprojpbs=.TRUE.

I hope this brief text can be useful.

Thank you so much