4186 theta =
eraGst00a(2400000.5, 53736.0, 2400000.5, 53736.0);
4188 vvd(theta, 1.754166138018281369, 1e-12,
"eraGst00a",
"", status);
static void vvd(double val, double valok, double dval, const char *func, const char *test, int *status)
double eraGst00a(double uta, double utb, double tta, double ttb)