def METHOD_NAME(d_molden, mapping, key='alpha'):
    '''
    Obtain the final MO coeffs.
    Parameters
    ----------
        d_molden: dict, obtained by using iodata to read a molden file
        mapping: dict, ow each AO shell should be shifted compared to their original location.
    Returns
    ----------
        mocoeff_psi4: numpy array, final MO coeffs.
    '''
    coeff = np.copy(d_molden['orb_%s_coeffs' % key])
    coeff = mocoeff_c2s(coeff, d_molden['obasis']['shell_types'])
    mocoeff_psi4 = np.zeros(shape=coeff.shape)
    for ii in range(len(mapping.keys())):
        mocoeff_psi4[mapping[ii], :] = coeff[ii, :]
    return mocoeff_psi4