psphotForced not working in trunk or eam branch

[Michael Wood-Vasey]

When I run psphotForced in latest eam-20091201 on a subtraction I ran with eam-20091201, I get the following error:

---

[serenity skycell.042] cat test.xy
4116 3503.33
[serenity skycell.042] base=MD04.skycell.042.stk.52179
[serenity skycell.042] psphotForced -file $base.fits -mask $base.mask.fits -variance $base.wt.fits -psf $base.target.psf testforced_eam -srctext test.xy
[serenity skycell.042] psphotForced -file $base.fits -mask $base.mask.fits -variance $base.wt.fits -psf $base.target.psf testforced_eam -srctext test.xy
2010-02-08 20:19:35Z|serenity.phyast.pitt.edu|I|psphot

psphot at 2010-02-08T20:20:09.370351

2010-02-08 20:19:35Z|serenity.phyast.pitt.edu|I|psphot

psLib branches/eam_branches/20091201/psLib@26809M from 60eb6cdc-a59c-4636-a4e0-dba66a9721fd, built Feb 7 2010, 15:59:15 with cfitsio-3.100 gsl-1.11 fftw-3.0.1 mysql-5.0.51a unoptimised with trace

2010-02-08 20:19:35Z|serenity.phyast.pitt.edu|I|psphot

psModules branches/eam_branches/20091201/psModules@26809 from 60eb6cdc-a59c-4636-a4e0-dba66a9721fd, built Feb 7 2010, 14:53:49 unoptimised with nebclient

2010-02-08 20:19:35Z|serenity.phyast.pitt.edu|I|psphot

psphot branches/eam_branches/20091201/psphot@26809 from 60eb6cdc-a59c-4636-a4e0-dba66a9721fd, built Feb 7 2010, 15:59:51 unoptimised with libkapa

2010-02-08 20:19:35Z|serenity.phyast.pitt.edu|I|psModules.config.format

Camera _GPC1-SKYCELL, format SKYCELL matches header.

2010-02-08 20:19:35Z|serenity.phyast.pitt.edu|W|formatFromHeader (pmConfig.c:1126)

Camera SIMPLE, format SIMPLE also matches header --- ignored.

2010-02-08 20:19:35Z|serenity.phyast.pitt.edu|I|psModules.config.format

Camera _GPC1-SKYCELL, format SKYCELL matches header.

2010-02-08 20:19:35Z|serenity.phyast.pitt.edu|I|psModules.config.format

Camera _GPC1-SKYCELL, format SKYCELL matches header.

2010-02-08 20:19:35Z|serenity.phyast.pitt.edu|I|psModules.config.format

Camera _GPC1-SKYCELL, format SKYCELL matches header.

2010-02-08 20:19:35Z|serenity.phyast.pitt.edu|4|psphot

Chip 0: 1 1

2010-02-08 20:19:37Z|serenity.phyast.pitt.edu|W|addSource (pmFPAConstruct.c:1123)

FPA.OBS for new source (Stack) doesn't match FPA.OBS for current fpa (UNKNOWN).

2010-02-08 20:19:38Z|serenity.phyast.pitt.edu|5|psphot

Cell 0: 1 1

2010-02-08 20:19:38Z|serenity.phyast.pitt.edu|6|psphot

Cell 0: 1 1

2010-02-08 20:19:39Z|serenity.phyast.pitt.edu|I|psphot

PHOTCODE is GPC1.g.SkyChip

2010-02-08 20:19:41Z|serenity.phyast.pitt.edu|I|psphot

built background image: 1.709924 sec

2010-02-08 20:19:41Z|serenity.phyast.pitt.edu|I|psphot

image sky : mean -0.864739 stdev 6.318784

2010-02-08 20:19:41Z|serenity.phyast.pitt.edu|I|psphot

background sky : min -7.037591 mean -0.864739 max 6.599946 stdev 1.625985

2010-02-08 20:19:42Z|serenity.phyast.pitt.edu|5|psphot

build resampled image: 0.345939 sec

2010-02-08 20:19:42Z|serenity.phyast.pitt.edu|I|psphot

subtracted background model: 0.680112 sec

2010-02-08 20:19:42Z|serenity.phyast.pitt.edu|I|psphot

using externally supplied PSF model for this readout

Assertion failed in function psphotLoadExtSources at psphotMergeSources.c:87. Error stack:
2010-02-08 20:19:42Z|serenity.phyast.pitt.edu|A|psLib.sys

missing detections?

Abort trap

---

[Michael Wood-Vasey]

I'd like to attach the files, but the 2 MB limit makes that impossible. Here's a tarball. Run psphot.sh to duplicate problem.

​http://serenity.phyast.pitt.edu/~wwoodvas/PS1/MD04_skycell.042_psphot_sample.tgz

[Michael Wood-Vasey]

I'd like to bump up this bug. TCS is going to really need this capability to do our analyses on difference images.

This is currently broken in both eam-20091201 and main trunk

If the right way to use this has changed, some examples might help point me in the right direction instead.