returning clock bias covariance in pntpos

src/pntpos.c

@@ -470,7 +470,11 @@ static int estpos(const obsd_t *obs, int n, const double *rs, const double *dts,
             for (j=0;j<3;j++) sol->qr[j]=(float)Q[j+j*NX];
             sol->qr[3]=(float)Q[1];    /* cov xy */
             sol->qr[4]=(float)Q[2+NX]; /* cov yz */
-            sol->qr[5]=(float)Q[2];    /* cov zx */
+            sol->qr[5]=(float)Q[2];    /* cov xz */
+            for (j=3;j<7;j++) sol->qdtr[j-3]=(float)Q[j+j*NX];
+            sol->qdtr[4]=(float)Q[3];     /* cov xdt*/
+            sol->qdtr[5]=(float)Q[3+NX];  /* cov ydt*/
+            sol->qdtr[6]=(float)Q[3+2*NX];/* cov zdt*/
             sol->ns=(unsigned char)ns;
             sol->age=sol->ratio=0.0;
             

src/rtklib.h

@@ -962,6 +962,8 @@ typedef struct {        /* solution type */
                         /* {c_ee,c_nn,c_uu,c_en,c_nu,c_ue} */
     float  qv[6];       /* velocity variance/covariance (m^2/s^2) */
     double dtr[6];      /* receiver clock bias to time systems (s) */
+    float  qdtr[7];     /* receiver clock bias to time systems variance/covariance*/
+                        /* {c_dtdt, c_GPS_GLO, c_GPS_GAL, c_GPS_CMP, c_xdt, c_ydt, czdt}*/
     unsigned char type; /* type (0:xyz-ecef,1:enu-baseline) */
     unsigned char stat; /* solution status (SOLQ_???) */
     unsigned char ns;   /* number of valid satellites */