Changeset 19948

Timestamp:

Oct 7, 2008, 10:06:25 AM (18 years ago)

Author:

eugene

Message:

handle chisq case for model missing

File:

• trunk/psModules/src/objects/pmSourceIO_PS1_DEV_1.c (modified) (5 diffs)

trunk/psModules/src/objects/pmSourceIO_PS1_DEV_1.c

@@ -3 +3 @@
  *  @author EAM, IfA
  *
- *  @version $Revision: 1.11 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2008-10-06 13:02:43 $
+ *  @version $Revision: 1.12 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2008-10-07 20:06:25 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -61 +61 @@
     psF32 xPos, yPos;
     psF32 xErr, yErr;
-    psF32 errMag;
+    psF32 errMag, chisq;
 
     // let's write these out in S/N order
@@ -84 +84 @@
             yErr = dPAR[PM_PAR_YPOS];
             axes = pmPSF_ModelToAxes (PAR, 20.0);
+            chisq = model->chisq;
 
             // need to determine the PSF photometry error: source->errMag is the error on the 'best' model mag.
@@ -95 +96 @@
             axes.minor = NAN;
             axes.theta = NAN;
+            chisq = NAN;
             errMag = NAN;
         }
@@ -114 +116 @@
         psMetadataAdd (row, PS_LIST_TAIL, "SKY_SIGMA",        PS_DATA_F32, "Sigma of sky level",                         source->skyErr);
 
-        psMetadataAdd (row, PS_LIST_TAIL, "PSF_CHISQ",        PS_DATA_F32,  "Chisq of PSF-fit",                          model->chisq);
+        psMetadataAdd (row, PS_LIST_TAIL, "PSF_CHISQ",        PS_DATA_F32,  "Chisq of PSF-fit",                          chisq);
         psMetadataAdd (row, PS_LIST_TAIL, "CR_NSIGMA",        PS_DATA_F32,  "Nsigma deviations from PSF to CF",          source->crNsigma);
         psMetadataAdd (row, PS_LIST_TAIL, "EXT_NSIGMA",       PS_DATA_F32,  "Nsigma deviations from PSF to EXT",         source->extNsigma);