Index: trunk/Ohana/src/addstar/src/ReadStarsSDSS.c
===================================================================
--- trunk/Ohana/src/addstar/src/ReadStarsSDSS.c	(revision 28938)
+++ trunk/Ohana/src/addstar/src/ReadStarsSDSS.c	(revision 28939)
@@ -11,4 +11,8 @@
   NAME = (TYPE *) gfits_get_bintable_column_data (table.header, &table, #NAME, type, &Nrow, &Ncol); \
   assert (NAME); assert (!strcmp (type, #TYPE)); assert (Nrow == Nstars); assert (Ncol == NFILTER);
+
+# define GET_COLUMN_5_NOASSERT(NAME,TYPE) \
+  TYPE *NAME; \
+  NAME = (TYPE *) gfits_get_bintable_column_data (table.header, &table, #NAME, type, &Nrow, &Ncol);
 
 # define GET_COLUMN_1(NAME,TYPE) \
@@ -123,5 +127,9 @@
   GET_COLUMN_5 (flags2, int);
 
+#ifdef notyet
   GET_COLUMN_5 (prob_psf, float);
+#else
+  GET_COLUMN_5_NOASSERT (prob_psf, float);
+#endif
 
   GET_COLUMN_1 (ra, double);
@@ -160,8 +168,12 @@
       stars[N].measure.FWx       = ShortPixels(seeing[j]); // reported in arcsec?
       stars[N].measure.FWy       = ShortPixels(seeing[j]);
-      stars[N].measure.psfChisq  = prob_psf[N]; // XXX not really the correct value...
+      if (prob_psf) {
+          stars[N].measure.psfChisq  = prob_psf[N]; // XXX not really the correct value...
+      } else {
+          stars[N].measure.psfChisq  = NAN;
+      }
       stars[N].measure.detID     = N;
       stars[N].measure.t         = tzero[j] + clockRate*rowc[N]; // time since row 0
-      stars[N].measure.dt        = 53.907456; // is this 2048*clockRate ?
+      stars[N].measure.dt        = 4.32912209; // 2.5 * log(53.907456) the sdss exposure time // old comment is 53907456 is this 2048*clockRate ?
 
       SetSDSSFlags (&stars[N], flags[N], flags2[N]);