half_K = GM_background_K
& *( half_K
_RL half_K
GM_PsiX(i,j,k,bi,bj) = half_K*sigmaX_W
GM_PsiY(i,j,k,bi,bj) = half_K*sigmaY_W