function [phi_L_ini] = g_position(xl,yl,rt,rg)

%rl=rt+800*10^3;
%tp_leo=km/7; %vitesse leo 7km/s, tp est le tps parcouru par leo pour la valeur de km

% position du GPS sur son orbite position calculee pour RO possible
[xg,yg] = f_tangeant(xl,yl,rt,rg);

% on fait le produit scalaire N^C
prod_scal = -xl*(xg-xl)-yl*(yg-yl);
norme1 = sqrt(xl^2+yl^2);
norme2 = sqrt((xg-xl)^2+(yg-yl)^2);

phi_L_ini = acos(prod_scal/(norme1*norme2));