for ( double timeJD = startTime; timeJD <= endTime; timeJD += step ) {
// Calculate coordinates and add transit to orbit
boolean relativity = false;
GaiaSimuTime gst = new GaiaSimuTime(new GaiaTime());
gst.setJD(timeJD);
GVector3d comrs = model.getStar().getAstrometry().getCoMRSPosition(gst, relativity);
GVector3d comrsRef = model.getAstrometry().getCoMRSPosition(gst, relativity);
GVector2d lpc = Conversion.fromCoMRStoLPC(comrs, comrsRef.getLongitude(), comrsRef.getLatitude());