problem with f_scaling in rhoBeCu_nist_mat()