Index: /branches/eam_branches/ipp-20130904/psModules/src/camera/pmFPAfile.c =================================================================== --- /branches/eam_branches/ipp-20130904/psModules/src/camera/pmFPAfile.c (revision 36161) +++ /branches/eam_branches/ipp-20130904/psModules/src/camera/pmFPAfile.c (revision 36162) @@ -491,4 +491,7 @@ if (!strcasecmp(type, "CMF")) { return PM_FPA_FILE_CMF; + } + if (!strcasecmp(type, "CFF")) { + return PM_FPA_FILE_CFF; } if (!strcasecmp(type, "WCS")) { Index: /branches/eam_branches/ipp-20130904/psModules/src/camera/pmFPAfileDefine.c =================================================================== --- /branches/eam_branches/ipp-20130904/psModules/src/camera/pmFPAfileDefine.c (revision 36161) +++ /branches/eam_branches/ipp-20130904/psModules/src/camera/pmFPAfileDefine.c (revision 36162) @@ -103,7 +103,13 @@ type = psMetadataLookupStr(&status, data, "FILE.TYPE"); + if (!type) { + psError(PM_ERR_CONFIG, true, "FILE.TYPE is not defined for %s\n", name); + psFree(file); + return NULL; + } + file->type = pmFPAfileTypeFromString(type); if (file->type == PM_FPA_FILE_NONE) { - psError(PM_ERR_CONFIG, true, "FILE.TYPE is not defined for %s\n", name); + psError(PM_ERR_CONFIG, true, "FILE.TYPE %s is not registered in pmFPAfile.c:pmFPAfileTypeFromString\n", type); psFree(file); return NULL; Index: /branches/eam_branches/ipp-20130904/psModules/src/camera/pmFPAfileIO.c =================================================================== --- /branches/eam_branches/ipp-20130904/psModules/src/camera/pmFPAfileIO.c (revision 36161) +++ /branches/eam_branches/ipp-20130904/psModules/src/camera/pmFPAfileIO.c (revision 36162) @@ -318,4 +318,5 @@ case PM_FPA_FILE_CMP: case PM_FPA_FILE_CMF: + case PM_FPA_FILE_CFF: case PM_FPA_FILE_WCS: case PM_FPA_FILE_PSF: @@ -514,4 +515,8 @@ case PM_FPA_FILE_WCS: psError(PS_ERR_IO, true, "cannot write type WCS (%s)", file->name); + return false; + + case PM_FPA_FILE_CFF: + psError(PS_ERR_IO, true, "cannot write type CFF (%s)", file->name); return false; @@ -568,4 +573,5 @@ case PM_FPA_FILE_PATTERN: case PM_FPA_FILE_CMF: + case PM_FPA_FILE_CFF: case PM_FPA_FILE_WCS: case PM_FPA_FILE_PSF: @@ -644,4 +650,5 @@ case PM_FPA_FILE_CMP: case PM_FPA_FILE_CMF: + case PM_FPA_FILE_CFF: case PM_FPA_FILE_WCS: case PM_FPA_FILE_PSF: @@ -805,4 +812,5 @@ case PM_FPA_FILE_PATTERN: case PM_FPA_FILE_CMF: + case PM_FPA_FILE_CFF: case PM_FPA_FILE_WCS: case PM_FPA_FILE_PSF: @@ -1017,4 +1025,5 @@ case PM_FPA_FILE_CMP: case PM_FPA_FILE_WCS: + case PM_FPA_FILE_CFF: case PM_FPA_FILE_JPEG: case PM_FPA_FILE_KAPA: Index: /branches/eam_branches/ipp-20130904/psModules/src/objects/pmSourceIO.c =================================================================== --- /branches/eam_branches/ipp-20130904/psModules/src/objects/pmSourceIO.c (revision 36161) +++ /branches/eam_branches/ipp-20130904/psModules/src/objects/pmSourceIO.c (revision 36162) @@ -69,4 +69,5 @@ psString *xfitname, // Extension name for extended fitted measurements psString *xradname, // Extension name for radial apertures + psString *xgalname, // Extension name for galaxy shapes const pmFPAfile *file, // File of interest const pmFPAview *view // View to level of interest @@ -140,4 +141,14 @@ } *xradname = pmFPAfileNameFromRule (rule, file, view); + } + + // EXTNAME for radial apertures + if (xgalname) { + const char *rule = psMetadataLookupStr(&status, menu, "CMF.XGAL"); + if (!rule) { + psError(PS_ERR_UNKNOWN, true, "missing entry for CMF.XGAL in EXTNAME.RULES in camera.config"); + return false; + } + *xgalname = pmFPAfileNameFromRule (rule, file, view); } @@ -362,4 +373,7 @@ status &= pmSourcesWrite_##TYPE##_XRAD (file->fits, readout, sources, file->header, xradname, recipe); \ } \ + if (xgalname) { \ + status &= pmSourcesWrite_##TYPE##_XGAL (file->fits, sources, xgalname, recipe); \ + } \ } @@ -468,4 +482,5 @@ bool XFIT_OUTPUT = psMetadataLookupBool(&status, recipe, "EXTENDED_SOURCE_FITS"); bool XRAD_OUTPUT = psMetadataLookupBool(&status, recipe, "RADIAL_APERTURES"); + bool XGAL_OUTPUT = psMetadataLookupBool(&status, recipe, "GALAXY_SHAPE_ANALYSIS"); // define the EXTNAME values for the different data segments: @@ -476,8 +491,10 @@ psString xfitname = NULL; psString xradname = NULL; + psString xgalname = NULL; if (!pmSourceIOextnames(&headname, &dataname, &deteffname, XSRC_OUTPUT ? &xsrcname : NULL, XFIT_OUTPUT ? &xfitname : NULL, XRAD_OUTPUT ? &xradname : NULL, + XGAL_OUTPUT ? &xgalname : NULL, file, view)) { return false; @@ -563,4 +580,7 @@ psMetadataAddStr (outhead, PS_LIST_TAIL, "XRADNAME", PS_META_REPLACE, "name of XRAD table extension", xradname); } + if (xgalname) { + psMetadataAddStr (outhead, PS_LIST_TAIL, "XGALNAME", PS_META_REPLACE, "name of XGAL table extension", xgalname); + } // these are case-sensitive since the EXTYPE is case-sensitive @@ -609,4 +629,5 @@ psFree (xfitname); psFree (xradname); + psFree (xgalname); psFree (deteffname); @@ -620,4 +641,5 @@ psFree (xfitname); psFree (xradname); + psFree (xgalname); psFree (deteffname); return false; @@ -922,4 +944,5 @@ psString xfitname = NULL; psString xradname = NULL; + psString xgalname = NULL; psMetadata *tableHeader = NULL; @@ -988,4 +1011,5 @@ bool XFIT_OUTPUT = psMetadataLookupBool(&status, recipe, "EXTENDED_SOURCE_FITS"); bool XRAD_OUTPUT = psMetadataLookupBool(&status, recipe, "RADIAL_APERTURES"); + bool XGAL_OUTPUT = false; // psMetadataLookupBool(&status, recipe, "GALAXY_SHAPES"); if (!pmSourceIOextnames(&headname, &dataname, &deteffname, @@ -993,4 +1017,5 @@ XFIT_OUTPUT ? &xfitname : NULL, XRAD_OUTPUT ? &xradname : NULL, + XGAL_OUTPUT ? &xgalname : NULL, file, view)) { return false; Index: /branches/eam_branches/ipp-20130904/psModules/src/objects/pmSourceIO.h =================================================================== --- /branches/eam_branches/ipp-20130904/psModules/src/objects/pmSourceIO.h (revision 36161) +++ /branches/eam_branches/ipp-20130904/psModules/src/objects/pmSourceIO.h (revision 36162) @@ -21,4 +21,5 @@ bool pmSourcesWrite_##TYPE##_XFIT(psFits *fits, pmReadout *readout, psArray *sources, psMetadata *imageHeader, char *extname); \ bool pmSourcesWrite_##TYPE##_XRAD(psFits *fits, pmReadout *readout, psArray *sources, psMetadata *imageHeader, char *extname, psMetadata *recipe); \ + bool pmSourcesWrite_##TYPE##_XGAL(psFits *fits, psArray *sources, char *extname, psMetadata *recipe); \ psArray *pmSourcesRead_##TYPE (psFits *fits, psMetadata *header); \ bool pmSourcesRead_##TYPE##_XSRC (psFits *fits, pmReadout *readout, psMetadata *header, psMetadata *tableHeader, psArray *sources, long *index); \ Index: /branches/eam_branches/ipp-20130904/psModules/src/objects/pmSourceIO_CMF.c.in =================================================================== --- /branches/eam_branches/ipp-20130904/psModules/src/objects/pmSourceIO_CMF.c.in (revision 36161) +++ /branches/eam_branches/ipp-20130904/psModules/src/objects/pmSourceIO_CMF.c.in (revision 36162) @@ -1344,14 +1344,14 @@ } -// XXX this layout is still the same as PS1_DEV_1 -bool pmSourcesWrite_CMF_@CMFMODE@_GAL (psFits *fits, pmReadout *readout, psArray *sources, psMetadata *imageHeader, char *extname) +// XXX where should I record the number of columns?? +bool pmSourcesWrite_CMF_@CMFMODE@_XGAL (psFits *fits, psArray *sources, char *extname, psMetadata *recipe) { - - psArray *table; - psMetadata *row; - psF32 *PAR, *dPAR; - psF32 xPos, yPos; - psF32 xErr, yErr; - char name[64]; + bool status = false; + + // perform full non-linear fits / extended source analysis? + if (!psMetadataLookupBool (&status, recipe, "GALAXY_SHAPES")) { + psLogMsg ("psphot", PS_LOG_INFO, "galaxy shapes were not measured, skipping\n"); + return true; + } // create a header to hold the output data @@ -1364,8 +1364,8 @@ sources = psArraySort (sources, pmSourceSortByFlux); - // we are writing one row per model; we need to write out same number of columns for each row: find the max Nparams - int nParamMax = 0; + psArray *table = psArrayAllocEmpty (sources->n); + for (int i = 0; i < sources->n; i++) { - // this is the source associated with this image + pmSource *thisSource = sources->data[i]; @@ -1373,156 +1373,42 @@ pmSource *source = thisSource->parent ? thisSource->parent : thisSource; + // if we did not fit the galaxy model, modelFits will be NULL if (source->modelFits == NULL) continue; - for (int j = 0; j < source->modelFits->n; j++) { - pmModel *model = source->modelFits->data[j]; - assert (model); - nParamMax = PS_MAX (nParamMax, model->params->n); - } - } - - @>PS1_DV2@ pmChip *chip = readout->parent->parent; - - table = psArrayAllocEmpty (sources->n); - - // we write out all sources, regardless of quality. the source flags tell us the state - for (int i = 0; i < sources->n; i++) { - - pmSource *thisSource = sources->data[i]; - - // this is the "real" version of this source - pmSource *source = thisSource->parent ? thisSource->parent : thisSource; - - // XXX if no model fits are saved, write out modelEXT? - if (source->modelFits == NULL) continue; - - // We have multiple sources : need to flag the one used to subtract the light (the 'best' model) - for (int j = 0; j < source->modelFits->n; j++) { - - // choose the convolved EXT model, if available, otherwise the simple one - pmModel *model = source->modelFits->data[j]; - assert (model); - - // pmSourceExtFitPars *extPars = source->extFitPars->data[j]; - // assert (extPars); - - // skip models which were not actually fitted - if (model->flags & PM_MODEL_STATUS_BADARGS) continue; - - PAR = model->params->data.F32; - dPAR = model->dparams->data.F32; - xPos = PAR[PM_PAR_XPOS]; - yPos = PAR[PM_PAR_YPOS]; - xErr = dPAR[PM_PAR_XPOS]; - yErr = dPAR[PM_PAR_YPOS]; - - @>PS1_DV2@ psSphere ptSky = {0.0, 0.0, 0.0, 0.0}; - @>PS1_DV2@ float posAngle = 0.0; - @>PS1_DV2@ float pltScale = 0.0; - @>PS1_DV2@ pmSourceLocalAstrometry (&ptSky, &posAngle, &pltScale, chip, xPos, yPos); - @>PS1_DV2@ double raPos = ptSky.r*PS_DEG_RAD; - @>PS1_DV2@ double decPos = ptSky.d*PS_DEG_RAD; - @>PS1_DV2@ posAngle *= PS_DEG_RAD; - @>PS1_DV2@ pltScale *= PS_DEG_RAD*3600.0; - - float kronFlux = source->moments ? source->moments->KronFlux : NAN; - float kronMag = isfinite(kronFlux) ? -2.5*log10(kronFlux) : NAN; - - row = psMetadataAlloc (); - - // XXX we are not writing out the mode (flags) or the type (psf, ext, etc) - // the psMetadataAdd entry and the double quotes are used by grep to select the output fields for automatic documentation - // This set of psMetadataAdd Entries marks the "----" "Start of the XFIT segment" - psMetadataAddU32 (row, PS_LIST_TAIL, "IPP_IDET", 0, "IPP detection identifier index", source->seq); - psMetadataAddF32 (row, PS_LIST_TAIL, "X_EXT", 0, "EXT model x coordinate", xPos); - psMetadataAddF32 (row, PS_LIST_TAIL, "Y_EXT", 0, "EXT model y coordinate", yPos); - psMetadataAddF32 (row, PS_LIST_TAIL, "X_EXT_SIG", 0, "Sigma in EXT x coordinate", xErr); - psMetadataAddF32 (row, PS_LIST_TAIL, "Y_EXT_SIG", 0, "Sigma in EXT y coordinate", yErr); - @>PS1_DV2@ psMetadataAddF32 (row, PS_LIST_TAIL, "RA_EXT", 0, "EXT model ra coordinate", raPos); - @>PS1_DV2@ psMetadataAddF32 (row, PS_LIST_TAIL, "DEC_EXT", 0, "EXT model dec coordinate", decPos); - @>PS1_DV2@ float instFlux = isfinite(model->mag) ? pow(10.0, -0.4*model->mag) : NAN; - @>PS1_DV2@ psMetadataAddF32 (row, PS_LIST_TAIL, "EXT_INST_FLUX", 0, "EXT fit instrumental counts", instFlux); - - psMetadataAddF32 (row, PS_LIST_TAIL, "EXT_INST_MAG", 0, "EXT fit instrumental magnitude", model->mag); - psMetadataAddF32 (row, PS_LIST_TAIL, "EXT_INST_MAG_SIG", 0, "Sigma of PSF instrumental magnitude", model->magErr); - - @>PS1_DV2@ float calMag = isfinite(magOffset) ? model->mag + magOffset : NAN; - @>PS1_DV2@ psMetadataAdd (row, PS_LIST_TAIL, "EXT_CAL_MAG", PS_DATA_F32, "EXT Magnitude using supplied calibration", calMag); - @>PS1_DV2@ psMetadataAdd (row, PS_LIST_TAIL, "EXT_CHISQ", PS_DATA_F32, "EXT Magnitude using supplied calibration", model->chisq); - @>PS1_DV2@ psMetadataAdd (row, PS_LIST_TAIL, "EXT_NDOF", PS_DATA_S32, "EXT Magnitude using supplied calibration", model->nDOF); - - @>PS1_DV2@ psMetadataAddF32 (row, PS_LIST_TAIL, "POSANGLE", 0, "position angle at source (degrees)", posAngle); - @>PS1_DV2@ psMetadataAddF32 (row, PS_LIST_TAIL, "PLTSCALE", 0, "plate scale at source (arcsec/pixel)", pltScale); - - // psMetadataAddF32 (row, PS_LIST_TAIL, "MOMENTS_XX", 0, "second moment in x", extPars->Mxx); - // psMetadataAddF32 (row, PS_LIST_TAIL, "MOMENTS_XY", 0, "second moment in x,y", extPars->Mxy); - // psMetadataAddF32 (row, PS_LIST_TAIL, "MOMENTS_YY", 0, "second moment in y", extPars->Myy); - // psMetadataAddF32 (row, PS_LIST_TAIL, "MOMENTS_R1", 0, "first radial moment", extPars->Mrf); - // psMetadataAddF32 (row, PS_LIST_TAIL, "MOMENTS_RH", 0, "half radial moment", extPars->Mrh); - - psMetadataAddF32 (row, PS_LIST_TAIL, "PSF_INST_MAG", 0, "PSF fit instrumental magnitude", source->psfMag); - psMetadataAddF32 (row, PS_LIST_TAIL, "AP_MAG", 0, "PSF-sized aperture magnitude", source->apMag); - psMetadataAddF32 (row, PS_LIST_TAIL, "KRON_MAG", 0, "Kron Mag", kronMag); - - @PS1_DV0,PS1_DV1@ psMetadataAddF32 (row, PS_LIST_TAIL, "NPARAMS", 0, "number of model parameters", model->params->n); - @ALL,!PS1_DV0,!PS1_DV1@ psMetadataAddS32 (row, PS_LIST_TAIL, "NPARAMS", 0, "number of model parameters", model->params->n); - psMetadataAddStr (row, PS_LIST_TAIL, "MODEL_TYPE", 0, "name of model", pmModelClassGetName (model->type)); - - // XXX these should be major and minor, not 'x' and 'y' - if (model->type == pmModelClassGetType("PS_MODEL_TRAIL")) { - psMetadataAddF32 (row, PS_LIST_TAIL, "EXT_WIDTH_MAJ", 0, "EXT width (major axis), length for trail", PAR[PM_PAR_LENGTH]); - psMetadataAddF32 (row, PS_LIST_TAIL, "EXT_WIDTH_MIN", 0, "EXT width (minor axis), sigma for trail", PAR[PM_PAR_SIGMA]); // this is not fitted - psMetadataAddF32 (row, PS_LIST_TAIL, "EXT_THETA", 0, "EXT orientation angle", PAR[PM_PAR_THETA]); - psMetadataAddF32 (row, PS_LIST_TAIL, "EXT_WIDTH_MAJ_ERR",0, "EXT width error (major axis)", dPAR[PM_PAR_LENGTH]); - psMetadataAddF32 (row, PS_LIST_TAIL, "EXT_WIDTH_MIN_ERR",0, "EXT width error (minor axis)", NAN); // this is not fitted, so error is NAN - psMetadataAddF32 (row, PS_LIST_TAIL, "EXT_THETA_ERR", 0, "EXT orientation angle (error)", dPAR[PM_PAR_THETA]); - } else { - if (!isfinite(PAR[PM_PAR_SXX]) || !isfinite(PAR[PM_PAR_SYY]) || !isfinite(PAR[PM_PAR_SXY])) { - psMetadataAddF32 (row, PS_LIST_TAIL, "EXT_WIDTH_MAJ", 0, "EXT width (SXX, isnan)", PAR[PM_PAR_SXX]); - psMetadataAddF32 (row, PS_LIST_TAIL, "EXT_WIDTH_MIN", 0, "EXT width (SYY, isnan)", PAR[PM_PAR_SYY]); - psMetadataAddF32 (row, PS_LIST_TAIL, "EXT_THETA", 0, "EXT angle (SXY, isnan)", PAR[PM_PAR_SXY]); - psMetadataAddF32 (row, PS_LIST_TAIL, "EXT_WIDTH_MAJ_ERR", 0, "EXT width err (SXX, isnan)", dPAR[PM_PAR_SXX]); - psMetadataAddF32 (row, PS_LIST_TAIL, "EXT_WIDTH_MIN_ERR", 0, "EXT width err (SYY, isnan)", dPAR[PM_PAR_SYY]); - psMetadataAddF32 (row, PS_LIST_TAIL, "EXT_THETA_ERR", 0, "EXT angle err (SXY, isnan)", dPAR[PM_PAR_SXY]); - } else { - psEllipseAxes axes = pmPSF_ModelToAxes (PAR, model->type); - psMetadataAddF32 (row, PS_LIST_TAIL, "EXT_WIDTH_MAJ", 0, "EXT width (major axis), length for trail", axes.major); - psMetadataAddF32 (row, PS_LIST_TAIL, "EXT_WIDTH_MIN", 0, "EXT width (minor axis), sigma for trail", axes.minor); - psMetadataAddF32 (row, PS_LIST_TAIL, "EXT_THETA", 0, "EXT orientation angle", axes.theta); - psMetadataAddF32 (row, PS_LIST_TAIL, "EXT_WIDTH_MAJ_ERR",0, "EXT width error (major axis)", dPAR[PM_PAR_SXX]); - psMetadataAddF32 (row, PS_LIST_TAIL, "EXT_WIDTH_MIN_ERR",0, "EXT width error (minor axis)", dPAR[PM_PAR_SYY]); - psMetadataAddF32 (row, PS_LIST_TAIL, "EXT_THETA_ERR", 0, "EXT orientation angle (error)", dPAR[PM_PAR_SXY]); - } - } - - // write out the other generic parameters - for (int k = 0; k < nParamMax; k++) { - if (k == PM_PAR_I0) continue; - if (k == PM_PAR_SKY) continue; - if (k == PM_PAR_XPOS) continue; - if (k == PM_PAR_YPOS) continue; - if (k == PM_PAR_SXX) continue; - if (k == PM_PAR_SXY) continue; - if (k == PM_PAR_SYY) continue; - - snprintf (name, 64, "EXT_PAR_%02d", k); - - if (k < model->params->n) { - psMetadataAddF32 (row, PS_LIST_TAIL, name, 0, "", model->params->data.F32[k]); - } else { - psMetadataAddF32 (row, PS_LIST_TAIL, name, 0, "", NAN); - } - } - - // optionally, write out the covariance matrix values - // XXX do I need to pad this to match the biggest covar matrix? - if (model->covar) { - for (int iy = 0; iy < model->covar->numCols; iy++) { - for (int ix = iy; ix < model->covar->numCols; ix++) { - snprintf (name, 64, "EXT_COVAR_%02d_%02d", iy, ix); - psMetadataAddF32 (row, PS_LIST_TAIL, name, 0, "", model->covar->data.F32[iy][ix]); - - } - } - } + + // if we did not fit the galaxy model, galaxyFits will also be NULL + if (source->galaxyFits == NULL) continue; + + pmModel *model = source->modelFits->data[0]; + if (!model) return false; + + // X,Y coordinates are stored with the model parameters + psF32 *PAR = model->params->data.F32; + + psMetadata *row = psMetadataAlloc (); + + // we write out the x,y positions so people can link to the psf either way (position or ID) + psMetadataAddU32 (row, PS_LIST_TAIL, "IPP_IDET", 0, "IPP detection identifier index", source->seq); + psMetadataAddF32 (row, PS_LIST_TAIL, "X_FIT", 0, "model x coordinate", PAR[PM_PAR_XPOS]); + psMetadataAddF32 (row, PS_LIST_TAIL, "Y_FIT", 0, "model y coordinate", PAR[PM_PAR_YPOS]); + psMetadataAddF32 (row, PS_LIST_TAIL, "NPIX", 0, "number of pixels for fits", source->galaxyFits->nPix); + + psVector *Io = source->galaxyFits->Io; + psVector *dIo = source->galaxyFits->dIo; + psVector *chisq = source->galaxyFits->chisq; + + // we write a series of columns for each of the elements (maybe I should save them as vectors instead?) + for (int j = 0; j < Io->n; j++) { + + char field[16]; + + snprintf (field, 16, "CHI_%02d", j); + psMetadataAddF32 (row, PS_LIST_TAIL, field, 0, "chisq of model fit", chisq->data.F32[j]); + + snprintf (field, 16, "I0_%02d", j); + psMetadataAddF32 (row, PS_LIST_TAIL, field, 0, "chisq of model fit", Io->data.F32[j]); + + snprintf (field, 16, "D_I0_%02d", j); + psMetadataAddF32 (row, PS_LIST_TAIL, field, 0, "chisq of model fit", dIo->data.F32[j]); + psArrayAdd (table, 100, row); psFree (row); @@ -1542,7 +1428,7 @@ } - psTrace ("pmFPAfile", 5, "writing ext data %s\n", extname); + psTrace ("pmFPAfile", 5, "writing galaxy data %s\n", extname); if (!psFitsWriteTable (fits, outhead, table, extname)) { - psError(psErrorCodeLast(), false, "writing ext data %s\n", extname); + psError(psErrorCodeLast(), false, "writing galaxy data %s\n", extname); psFree (outhead); psFree(table); Index: /branches/eam_branches/ipp-20130904/psModules/src/objects/pmSourceIO_PS1_CAL_0.c =================================================================== --- /branches/eam_branches/ipp-20130904/psModules/src/objects/pmSourceIO_PS1_CAL_0.c (revision 36161) +++ /branches/eam_branches/ipp-20130904/psModules/src/objects/pmSourceIO_PS1_CAL_0.c (revision 36162) @@ -713,2 +713,7 @@ return true; } + +bool pmSourcesWrite_PS1_CAL_0_XGAL (psFits *fits, psArray *sources, char *extname, psMetadata *recipe) +{ + return true; +} Index: /branches/eam_branches/ipp-20130904/psModules/src/objects/pmSourceIO_PS1_DEV_0.c =================================================================== --- /branches/eam_branches/ipp-20130904/psModules/src/objects/pmSourceIO_PS1_DEV_0.c (revision 36161) +++ /branches/eam_branches/ipp-20130904/psModules/src/objects/pmSourceIO_PS1_DEV_0.c (revision 36162) @@ -255,2 +255,7 @@ return true; } + +bool pmSourcesWrite_PS1_DEV_0_XGAL(psFits *fits, psArray *sources, char *extname, psMetadata *recipe) +{ + return true; +} Index: /branches/eam_branches/ipp-20130904/psModules/src/objects/pmSourceIO_PS1_DEV_1.c =================================================================== --- /branches/eam_branches/ipp-20130904/psModules/src/objects/pmSourceIO_PS1_DEV_1.c (revision 36161) +++ /branches/eam_branches/ipp-20130904/psModules/src/objects/pmSourceIO_PS1_DEV_1.c (revision 36162) @@ -595,2 +595,7 @@ return true; } + +bool pmSourcesWrite_PS1_DEV_1_XGAL(psFits *fits, psArray *sources, char *extname, psMetadata *recipe) +{ + return true; +} Index: /branches/eam_branches/ipp-20130904/psModules/src/objects/pmSourceIO_SMPDATA.c =================================================================== --- /branches/eam_branches/ipp-20130904/psModules/src/objects/pmSourceIO_SMPDATA.c (revision 36161) +++ /branches/eam_branches/ipp-20130904/psModules/src/objects/pmSourceIO_SMPDATA.c (revision 36162) @@ -225,2 +225,7 @@ return true; } + +bool pmSourcesWrite_SMPDATA_XGAL(psFits *fits, psArray *sources, char *extname, psMetadata *recipe) +{ + return true; +}