//typedef unsigned char Boolean ;
#define   True    1
#define   False   0

/********* type definitions *********/

typedef char  i1;  /* data type allowing for the use of  8 bits. */
typedef short i2;  /* data type allowing for the use of 16 bits. */
typedef long  i4;  /* data type allowing for the use of 32 bits. */

typedef unsigned char  u1;          /* unsigned  8 bits */
typedef unsigned short u2;          /* unsigned 16 bits */
typedef unsigned long  u4;         /* unsigned 32 bits */
typedef float    r4;
typedef double   r8;

/* emphric varition */
long   TLW;
long   HOW;
char   URA;
int    prn;
short  wn;				 /*  GPS week number */
long   tow;                          /*  sec of GPS week */
double  tgd;                          /*  grup delay (s)*/
long   aodc;                         /*  clock date issue */
long   toc;                          /*  clock data reference time(sec) */
double af2,                          /*  clock correction(sec/sec2) */
       af1,                          /*  clock correction(sec/sec) */
       af0;                          /*  clock correction(sec) */
double alph0,                        /*  ionospheric parameter(second) */
       alph1,                        /*  ionospheric parameter(sec/semi-circles) */
       alph2,                        /*  ionospheric parameter(s/semi-circles2)  */
       alph3,                        /*  ionospheric parameter(s/semi-circles3)  */
       belta0,                       /*  ionospheric parameter(s)     */
       belta1,                       /*  ionospheric parameter(s/semi-circles)  */
       belta2,                       /*  ionospheric parameter(s/semi-circles2) */
       belta3;                       /*  ionospheric parameter(s/semi-circles3) */
long   aode,aode3;                   /*  orbit date issue  in subframes 2 and 3 */
double  deltan;                       /*  mean anomoly correction
					     semi-circles per sec*pi = rads */
double m0,                           /*  mean anomoly at ref time
				     semi-circles*pi = rads */
       e,                            /*  eccentricity */
       roota;                        /*  sqr root a ( metres 1/2 ) */
long   toe;                          /*  ref time (sec) */
double cic,                          /*  harmonic correction term(rads) */
       crc,                          /*  harmonic correction term(metres) */
       cis,                          /*  harmonic correction term(rads) */
       crs,                          /*  harmonic correction term(metres) */
       cuc,                          /*  harmonic correction term(rads) */
       cus;                          /*  harmonic correction term(rads) */
double omega0,                       /*  longitude  of asccending node
					     semi-circles*pi = rads */
       omega,                        /*  argument of perigee
					     semi-circles*pi */
       i0;                           /*  inclination angle at ref time
					     semi-circles*pi */
double omegadot,                     /*  rate of right ascension
					     semi-circles/sec*pi */
       idot;                         /*  rate of inclination
					     semi-circles/sec*pi */
int     accuracy,                     /* User range accuracy */
	health,                       /* Satellite health */
	fit;                          /* Curve fit interval */
char   prnum;			 /* sv prn number */
char   res;                           /* Reserved character */
unsigned checksum;                     /* check sum */



#define PI          3.1415926535898      /* GPS value of Pi */
#define WE          7.2921151467E-5      /* WGS84 earths rotation rate (rads/sec) */

#define SQRMU       1.99649818432174E7   /* Square root of MU.               */
#define FCONST     -4.442807633E-10      /* FCONST is used in relativistic clock
					    correction.  Defined on pg. 73 of
					    ICD-GPS-200.                     */

#define TWOM5       0.03125000000000     /* 2^-5  */
#define TWOM11      (TWOM5  /  64.0)
#define TWOM12      (TWOM11  /  2.0)
#define TWOM19      (TWOM11 / 256.0)
#define TWOM20      (TWOM19 /   2.0)
#define TWOM21      (TWOM20 /   2.0)
#define TWOM23      (TWOM21 /   4.0)
#define TWOM24      (TWOM23 /   2.0)
#define TWOM27      (TWOM24 /   8.0)
#define TWOM29      (TWOM23 /  64.0)
#define TWOM30      (TWOM29 /   2.0)
#define TWOM31      (TWOM30 /   2.0)
#define TWOM33      (TWOM31 /   4.0)
#define TWOM38      (TWOM33 /  32.0)
#define TWOM43      (TWOM38 /  32.0)
#define TWOM50      (TWOM43 / 128.0)
#define TWOM55      (TWOM50 /  32.0)
#define TWOM3125    (0.25 * TWOM31)

#define TWOP11      2048.0               /* 2^11  */
#define TWOP12      (TWOP11 * 2.0)
#define TWOP14      (TWOP12 * 4.0)
#define TWOP16      (TWOP14 * 4.0)

#define TWOP12_INTEGER   4096            /* 2^12  */
#define TWOP16_INTEGER  65536            /* 2^16  */

#define PITWOM19    (PI * TWOM19)
#define PITWOM23    (PI * TWOM23)
#define PITWOM31    (PI * TWOM31)
#define PITWOM38    (PI * TWOM38)
#define PITWOM43    (PI * TWOM43)