баг там знатный с первого и даже пятого взгляда кажется что все OK !
поэтому** все и пользуются с багом**
extern int ionocorr(gtime_t time, const nav_t *nav, int sat, const double *pos,
const double *azel, int ionoopt, double *ion, double *var)
{
.....
/* broadcast model */
// if (ionoopt==IONOOPT_BRDC)
{
double ion2 ;
ion2 = ionmodel_sim( time,nav->ion_gps,pos,azel);
*ion=ionmodel(time,nav->ion_gps,pos,azel);
trace(2," ion_old=%2.2f ion_sim=%2.2f ",*ion, ion2 );
*ion = ion2;
*var=SQR(*ion);
// return 1;
}
/* sbas ionosphere model */
if (ionoopt==IONOOPT_SBAS)
{
sbsioncorr(time,nav,pos,azel,ion,var);
trace(2," ion_sbas=%2.2f \n",*ion);
return 1;
}
.....
}
ionmodel_sim и ionmodel реализуют один и тот же алгоритм klobuchar model
только в ionmodel_sim нет бага
вот в реализации тропосферной задержки нет багов
SDCM
2 sat=25 azel= 2.9 -> 2 trp_saas=47.16 2 trp_sbas=36.04
2 sat= 8 azel= 5.2 -> 2 trp_saas=26.44 2 trp_sbas=24.18
2 sat=13 azel= 7.5 -> 2 trp_saas=18.33 2 trp_sbas=17.69
2 sat=10 azel= 8.0 -> 2 trp_saas=17.10 2 trp_sbas=16.62
2 sat=27 azel=10.7 -> 2 trp_saas=12.78 2 trp_sbas=12.69
2 sat=15 azel=11.8 -> 2 trp_saas=11.66 2 trp_sbas=11.63
2 sat=13 azel=15.5 -> 2 trp_saas= 8.95 2 trp_sbas= 9.01
2 sat= 5 azel=15.0 -> 2 trp_saas= 9.23 2 trp_sbas= 9.28
2 sat=21 azel=16.2 -> 2 trp_saas= 8.54 2 trp_sbas= 8.61
2 sat=29 azel=16.4 -> 2 trp_saas= 8.48 2 trp_sbas= 8.55
2 sat=15 azel=25.8 -> 2 trp_saas= 5.48 2 trp_sbas= 5.57
2 sat=16 azel=35.0 -> 2 trp_saas= 4.16 2 trp_sbas= 4.23
2 sat=26 azel=36.1 -> 2 trp_saas= 4.05 2 trp_sbas= 4.13
2 sat=10 azel=36.3 -> 2 trp_saas= 4.03 2 trp_sbas= 4.10
2 sat=27 azel=36.6 -> 2 trp_saas= 4.01 2 trp_sbas= 4.08
2 sat=20 azel=38.0 -> 2 trp_saas= 3.87 2 trp_sbas= 3.94
2 sat=21 azel=40.4 -> 2 trp_saas= 3.68 2 trp_sbas= 3.75
2 sat=29 azel=42.0 -> 2 trp_saas= 3.56 2 trp_sbas= 3.63
2 sat=16 azel=47.4 -> 2 trp_saas= 3.24 2 trp_sbas= 3.30
2 sat=26 azel=53.9 -> 2 trp_saas= 2.95 2 trp_sbas= 3.00
2 sat=18 azel=63.8 -> 2 trp_saas= 2.66 2 trp_sbas= 2.71
2 sat=20 azel=67.1 -> 2 trp_saas= 2.59 2 trp_sbas= 2.64
2 sat=18 azel=86.7 -> 2 trp_saas= 2.39 2 trp_sbas= 2.43
наглядно видать что тропосферная задержка больше чем ионосферная.