/* tle_pos() accuracy --------------------------------------------------------*/
static void utest3(void)
{
const char *file1="../data/tle/brdc3050.12*";
const char *file2="../data/tle/TLE_GNSS_20121101.txt";
const char *file3="../data/tle/igs17127.erp";
const double ep[6]={2012,10,31,0,0,0};
nav_t nav={0};
erp_t erp={0};
tle_t tle={0};
gtime_t time;
char sat[32];
double rs1[6],rs2[6],ds[6],dts[2],var;
int i,j,k,stat,svh;
readrnx(file1,0,"",NULL,&nav,NULL);
assert(nav.n>0);
stat=readerp(file3,&erp);
assert(stat);
stat=tle_read(file2,&tle);
assert(stat);
for (i=0;i<MAXSAT;i++) {
satno2id(i+1,sat);
fprintf(OUT,"SAT=%s\n",sat);
for (j=0;j<96;j++) {
time=timeadd(epoch2time(ep),900.0*j);
if (!satpos(time,time,i+1,EPHOPT_BRDC,&nav,rs1,dts,&var,&svh)) continue;
if (satsys(i+1,NULL)==SYS_QZS) svh&=0xFE;
if (svh) continue;
stat=tle_pos(time,sat,"","",&tle,&erp,rs2);
assert(stat);
for (k=0;k<3;k++) ds[k]=rs2[k]-rs1[k];
fprintf(OUT,"%6.0f %11.3f %11.3f %11.3f %11.3f\n",900.0*j,
ds[0]/1e3,ds[1]/1e3,ds[2]/1e3,norm(ds,3)/1e3);
assert(norm(ds,3)/1e3<300.0);
}
fprintf(OUT,"\n");
}
fprintf(OUT,"%s utest3 : OK\n",__FILE__);
}
/* open procssing session ----------------------------------------------------*/
static int openses(const prcopt_t *popt, const solopt_t *sopt,
const filopt_t *fopt, nav_t *nav, pcvs_t *pcvs, pcvs_t *pcvr)
{
char *ext;
trace(3,"openses :\n");
/* read satellite antenna parameters */
if (*fopt->satantp&&!(readpcv(fopt->satantp,pcvs))) {
showmsg("error : no sat ant pcv in %s",fopt->satantp);
trace(1,"sat antenna pcv read error: %s\n",fopt->satantp);
return 0;
}
/* read receiver antenna parameters */
if (*fopt->rcvantp&&!(readpcv(fopt->rcvantp,pcvr))) {
showmsg("error : no rec ant pcv in %s",fopt->rcvantp);
trace(1,"rec antenna pcv read error: %s\n",fopt->rcvantp);
return 0;
}
/* read dcb parameters */
if (*fopt->dcb) {
readdcb(fopt->dcb,nav);
}
/* read ionosphere data file */
if (*fopt->iono&&(ext=strrchr(fopt->iono,'.'))) {
if (strlen(ext)==4&&(ext[3]=='i'||ext[3]=='I')) {
readtec(fopt->iono,nav,0);
}
#ifdef EXTSTEC
else if (!strcmp(ext,".stec")||!strcmp(ext,".STEC")) {
stec_read(fopt->iono,nav);
}
#endif
}
/* open geoid data */
if (sopt->geoid>0&&*fopt->geoid) {
if (!opengeoid(sopt->geoid,fopt->geoid)) {
showmsg("error : no geoid data %s",fopt->geoid);
trace(2,"no geoid data %s\n",fopt->geoid);
}
}
/* read erp data */
if (*fopt->eop) {
if (!readerp(fopt->eop,&nav->erp)) {
showmsg("error : no erp data %s",fopt->eop);
trace(2,"no erp data %s\n",fopt->eop);
}
}
return 1;
}
