localtimeOffset
p->rJD += localtimeOffset(p);
c1 = localtimeOffset(p);
p->rJD += c1 - localtimeOffset(p);