Feature: Support nspin=4 PW partial-charge and wavefunction .cube output (Useful Information to know how to use out_pchg, out_wfc_norm, and out_wfc_re_im functions) - #7841