#define I96 96#define I50 50#define I36 36#define I10 10#define I5 5#define I3 3#define FL1 1575.42e6   /* MHz */#define FL2 1227.60e6   /* MHz */#define C 2.99792458e8  /* light velocity */#define Lambda1 (C/FL1) /* m */#define Lambda2 (C/FL2) /* m */#define PI 3.14159265358979324e0int getline();               /* rinex.c */void r_rinexh();void r_rinexd();void getdat();void getdat_o();void t_hms();void id_i12j12();void odat_tc();void sdtdat();double line_interpo();double line2_interpo();double line2_integrate();void id_obs();void cycleslips_test();void gps_dat();void r_con();               /* controlf.c */int indexx();void w_info();void DOP();void wr_condat();void igs();                 /* igs.c + broadcast.c + orbit.c> */void broadcast();int getorb();int refsat();void igsdat();double lagrange();void r_igs_h();void r_igs_d();void r_eph();void eph_check();  void DdE_e();void broadcast_orbit();double Kepler_equation();void borbc();int findisat();void rotation();           /* tranf.c */void TRANF();              double zenith1();int chol_inv();            /* ad_core.c */double ls();void eq_dd_s();            /* eqdd_s.c */void edit_norm_dd();       /* normdd_s.c */void ambifix();            /* ambifix.c */void NINi();void min2m0();void getNk();void NisNk();double dfix();void earthrot();           /* singlep.c */double distance();void transrot();void transrot1();void singlep();double twoi2f();void PC12();#define Isat 26            /* header of main */#define Isato 15#define Ista 6#define I_LA 20#define Iepoch 96#define Iepoches 10#define Islips 2000#define Iswitch 6#define Iline 130#define Iunknown 150       /* maximum unknowns number */#define N Iunknown#define M (10*N) #define Nlow (N*(N+2)/2)   /* if N counted from 0, then N += 1 */#define IR (Isato*Ista)#define Rlow (IR*(IR+2)/2)int O_rinexh=1,O_rinexd=2,O_getdat=2,O_getdato=1,O_odat_tc=1; /*Debug*/int O_sdtdat=1,O_cycle=2,O_gpsdat=6,O_rcon=1;int O_igsdat=1,O_lagrange=1,O_lagrangef=1,O_rigs=1,O_getorb=1;int O_r_eph=1,O_broadcast=2,O_borbc=1,O_eph_check=1;int O_broadcast_orbit=1,O_refsat=5;int O_tro=1,O_rot=1,O_rela=1,O_Rotats=1,O_SL_eph=1,O_tide=2;int O_gcslsn=0,O_gcsls=0,O_cls=2,O_rls=0,O_ls=1,O_chol=3,O_gaussj=3;int O_eq=7,O_norm=3,O_model=0,O_main=6,O_norm_cp=5,O_teun=2,O_kalman=2;int O_screen=1,O_singlep=7,O_specula=4,O_ambifix=4,O_ambifib=4;int C1,C2,L1,L2,P1,P2,D1,D2,Ob_num,lastn_num,Islip=0,Ref_sv=0,Inv_n=1;/*common*/int Doppler_sign=-1,Kstation=0,Code=3,Igs_v=1,Run_n=1,Ncancel_n=0,Coord_gl=3,    N_n=50,CP=0,K2=0,Nfix=0,Condin=0,CPDV=0,M_cyctest=0,Combine_code=0,    Velocity=0,Ambiguity=0,Ifilter=0,Ifilter1=0,Robust=0,TideCoef=0;int Ncancel_n1=0,Ncancel_n2=0,Super=0;int Frequency_n=1,Tsta,PC1,PC2,PD1,PD2,ReceivC=1,DopplerI=0,Schue=0,Iearthrot=2;int Ixyz=4,Icoord,Iluisa=1,Istate=0,Flysta=0,Igood=0,Irfila=0,Iotf=0;int Itidesign=1,Itropsign=1,Iclocksign1=1,Iclocksign2=1,Iclocksign3=1;int Iddfly=8,Imulstatic=0,Ispecula=0,Idisn=0,Idis1=0;double Cycleslip=1.1,Wc=0.1,Wp=100.,Wd=0.,Distance=4.566,Condic[6];double Factor1=1.,Factor2=-1.,Dt_real,truedis[Ista];double M0best,Ratio,M0round,M0teun;double N00f[Iunknown];long N00n[Iunknown],N00_n[Iunknown],N00id[Iunknown],N00tb[Iunknown],     N00te[Iunknown],N00_n1[Iunknown],N00nc[Iunknown];long N0xid[Iunknown],N0xtb[Iunknown],N0xte[Iunknown],     N0x[Iunknown],Fi0x[Iunknown],N0vv[Iunknown];double Nk[Iunknown];double TropCoef[Ista],X_0[Iunknown];