spinning HD example
ext_corepos[Z] now holds core polarization