#include "prototype.h"

/****** Function g_position */

void g_position(double *phi_L_ini, double xl, double yl, double ae, double r_gps);

{

  
  double norme1, norme2, prod_scal;
  double xg, yg;
 
 
  f_tangeant(&xg, &yg, xl, yl, ae, r_gps);
  prod_scal=-xl*(xg-xl)-yl*(yg-yl);
  norme1=sqrt((xl*xl)+(yl*yl));
  norme2=sqrt((xg-xl)*(xg-xl)+(yg-yl)*(yg-yl));

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

}