highFac(i,j)= 0. _d 0
highFac(i,j)= cH3
highFac(i,j)= cM2
rC4km = highFac(i,j)*(rCenter+ABS(rCenter))*mskM
rC4kp = highFac(i,j)*(rCenter-ABS(rCenter))*mskP
_RL highFac(1-OLx:sNx+OLx,1-OLy:sNy+OLy)