4111 theta =
eraGmst00(2400000.5, 53736.0, 2400000.5, 53736.0);
4113 vvd(theta, 1.754174972210740592, 1e-12,
"eraGmst00",
"", status);
double eraGmst00(double uta, double utb, double tta, double ttb)
static void vvd(double val, double valok, double dval, const char *func, const char *test, int *status)