4136 theta =
eraGmst06(2400000.5, 53736.0, 2400000.5, 53736.0);
4138 vvd(theta, 1.754174971870091203, 1e-12,
"eraGmst06",
"", status);
static void vvd(double val, double valok, double dval, const char *func, const char *test, int *status)
double eraGmst06(double uta, double utb, double tta, double ttb)