1702 double utc1, utc2, dut1, elong, phi, hm, xp, yp, phpa, tc, rh, wl,
1710 elong = -0.527800806;
1719 (void)
eraApio13(utc1, utc2, dut1, elong, phi, hm, xp, yp,
1720 phpa, tc, rh, wl, &astrom);
1722 ob1 = 2.710085107986886201;
1723 ob2 = 0.1717653435758265198;
1724 eraAtoiq(
"R", ob1, ob2, &astrom, &ri, &di);
1725 vvd(ri, 2.710121574449135955, 1e-12,
1726 "eraAtoiq",
"R/ri", status);
1727 vvd(di, 0.1729371839114567725, 1e-12,
1728 "eraAtoiq",
"R/di", status);
1730 ob1 = -0.09247619879782006106;
1731 ob2 = 0.1717653435758265198;
1732 eraAtoiq(
"H", ob1, ob2, &astrom, &ri, &di);
1733 vvd(ri, 2.710121574449135955, 1e-12,
1734 "eraAtoiq",
"H/ri", status);
1735 vvd(di, 0.1729371839114567725, 1e-12,
1736 "eraAtoiq",
"H/di", status);
1738 ob1 = 0.09233952224794989993;
1739 ob2 = 1.407758704513722461;
1740 eraAtoiq(
"A", ob1, ob2, &astrom, &ri, &di);
1741 vvd(ri, 2.710121574449135955, 1e-12,
1742 "eraAtoiq",
"A/ri", status);
1743 vvd(di, 0.1729371839114567728, 1e-12,
1744 "eraAtoiq",
"A/di", status);
static void vvd(double val, double valok, double dval, const char *func, const char *test, int *status)
void eraAtoiq(const char *type, double ob1, double ob2, eraASTROM *astrom, double *ri, double *di)
int eraApio13(double utc1, double utc2, double dut1, double elong, double phi, double hm, double xp, double yp, double phpa, double tc, double rh, double wl, eraASTROM *astrom)