xmf = 0. _d 0
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
Cheapmask(i,j,bi,bj) = xmf*cheapaml_taurelax
_RL local,localt,xmf