IF (iG.LE.xmw) xmf = 1. _d 0 - (iG-1 )*recipMW
IF (iG.GE.Nx-xmw+1) xmf = 1. _d 0 - (Nx-iG)*recipMW
xmf = 1. _d 0 - (Ny-jG)*recipMW
xmf = 1. _d 0 - (iG-1 )*recipMW *(Ny-jG)*recipMW
xmf = 1. _d 0 - (Nx-iG)*recipMW *(Ny-jG)*recipMW
xmf = 1. _d 0 - (jG-1)*recipMW
xmf = 1. _d 0 - (iG-1 )*recipMW*(jG-1)*recipMW
xmf = 1. _d 0 - (Nx-iG)*recipMW*(jG-1)*recipMW
relaxMask(i,j,bi,bj) = xmf*cheapaml_taurelax
xmf = 0. _d 0
IF (iG.LE.xmw) xmf = oneRL - (iG-1 )*tmpVar
IF (iG.GE.Nx-xmw+1) xmf = oneRL - (Nx-iG)*tmpVar
& xmf = MAX( xmf, oneRL - (jG-1 )*tmpVar )
& xmf = MAX( xmf, oneRL - (Ny-jG)*tmpVar )
relaxMask(i,j,bi,bj) = xmf
_RL xmf, tmpVar