function [xl_out,yl_out] = f_orbite_leo(xl_in,yl_in,tp,sens)

ae=6378.137*10^3;
rl=ae+800*10^3;

teta=atan(yl_in/xl_in);

if ( ((xl_in<0) & (yl_in>0)) | ((xl_in<0) & (yl_in<0)) )
   teta=teta+pi;
end

delta_teta=(tp*2*pi)/(1.7897*3600);

if (sens == 1)
   xl_out=rl*cos(teta+delta_teta);
   yl_out=rl*sin(teta+delta_teta);
else
   xl_out=rl*cos(teta-delta_teta);
   yl_out=rl*sin(teta-delta_teta);
   
end