def METHOD_NAME(d_molden, mapping, key='alpha'):
    coeff = np.copy(d_molden['orb_%s_coeffs' % key])
    coeff = mocoeff_c2s(coeff, d_molden['obasis']['shell_types'])
    mooeff_psi4 = np.zeros(shape=coeff.shape)
    for ii in range(len(mapping.keys())):
        #         print(ii, mapping[ii])
        mooeff_psi4[mapping[ii], :] = coeff[ii, :]
    return mooeff_psi4