longitud1=XC(i-1,j,bi,bj)*Coord2longitude
& -SIN( tidal_freq*myTime + 2.*longitud1 )
_RL Coord2longitude,longitud1,longitud2