_RL xx_ij,xx_ip1j,xx_ijp1,xx_ip1jp1
xx_ij=xG(i,j,bi,bj)
if (xx_ip1j.GT.xx_ij+180) xx_ip1j=xx_ip1j-360
if (xx_ip1j.LT.xx_ij-180) xx_ip1j=xx_ip1j+360
if (xx_ijp1.GT.xx_ij+180) xx_ijp1=xx_ijp1-360
if (xx_ijp1.LT.xx_ij-180) xx_ijp1=xx_ijp1+360
if (xx_ip1jp1.GT.xx_ij+180) xx_ip1jp1=xx_ip1jp1-360
if (xx_ip1jp1.LT.xx_ij-180) xx_ip1jp1=xx_ip1jp1+360
xx = xx_ij + ddx*( xx_ip1j - xx_ij )
& + ddy*( xx_ijp1 - xx_ij )
& + ddx*ddy*( xx_ip1jp1 - xx_ip1j - xx_ijp1 + xx_ij )