_RL ftmp
ftmp = 2. _d 0 * floc(-1) - floc(+0)
floc(-2) = ftmp
& + mloc(-2)*(fbar(ir-2)-ftmp)
ftmp = 2. _d 0 * floc(+0) - floc(-1)
floc(+1) = ftmp
& + mloc(+1)*(fbar(ir+1)-ftmp)
& + mloc(-2)*(fbar(ix-2)-ftmp)
& + mloc(+1)*(fbar(ix+1)-ftmp)
& + mloc(-2)*(fbar(iy-2)-ftmp)
& + mloc(+1)*(fbar(iy+1)-ftmp)
ftmp = 2. _d 0 * floc(-2) - floc(-1)
floc(-3) = ftmp
& + mloc(-3)*(fbar(ir-3)-ftmp)
ftmp = 2. _d 0 * floc(+1) - floc(+0)
floc(+2) = ftmp
& + mloc(+2)*(fbar(ir+2)-ftmp)
& + mloc(-3)*(fbar(ix-3)-ftmp)
& + mloc(+2)*(fbar(ix+2)-ftmp)
& + mloc(-3)*(fbar(iy-3)-ftmp)
& + mloc(+2)*(fbar(iy+2)-ftmp)