#include "prototype.h"

/****** Function g_ini */

void g_ini(int *flag, double *x, double *y, double *angle_rad, double *phi, double phi_l_rad,
	   double r_leo, double r1, double x1, double y1)

  {
      
    double alpha, teta;

    *phi = asin((r_leo/r1)*sin(phi_l_rad));
    *flag=1;
 
      if(*phi == 0)
	{
	  
	  alpha=*phi-phi_l-rad;
	  teta=atan(yl/xl);

	  if(xl<0)
	    teta=teta+PI;
	  else
	    teta=teta;

	  *angle_rad=teta+alpha;
	  *x=r1*cos(*angle_rad);
	  *y=r1*sin(*angle_rad);

	}
      else
	{
	  *flag=0;
	  *phi=0;
	  *angle_rad=0;
	  *x=xl;
	  *y=yl;
	}

  }