unwise_coadd.py starting: args: ['/n/fink1/ameisner/unwise-coadds/unwise_coadd.py', '--tile', '1792m197', '--band', '2', '--outdir', '/n/fink2/ameisner/unwise-coadds/tr_neo2/e000', '--recover_warped', '--no_irsa_dl', '--reference_dir', '/n/fink2/ameisner/unwise-coadds/neo2_ref_w2', '--epoch', '0', '--warp_all', '--no_output_masks'] PBS_ARRAYID: 0 opt: {'threads1': None, 'reference_dir': '/n/fink2/ameisner/unwise-coadds/neo2_ref_w2', 'plotprefix': None, 'force': False, 'no_sanity_check': False, 'height': 0, 'hi_lo_rej': False, 'medfilt': None, 'minmax': None, 'size': 2048, 'plots2': False, 'cube1': False, 'bgmatch': False, 'pixscale': 2.75, 'int_gz': False, 'nframes': 0, 'width': 0, 'epoch': 0, 'ra': None, 'dsky': False, 'preprocess': False, 'try_download': False, 'before': 100000.0, 'do_rebin': True, 'md5': False, 'dataset': 'sequels', 'cube': False, 'after': None, 'tile': '1792m197', 'threads': None, 'allmd5': False, 'rchi_fraction': 0.01, 'plots': False, 'frame0': 0, 'wishlist': False, 'maxmem': 0, 'compare_moon_all': False, 'center': False, 'band': [2], 'output_masks': False, 'use_zp_meta': False, 'outdir': '/n/fink2/ameisner/unwise-coadds/tr_neo2/e000', 'pdf': False, 'warp_all': True, 'dec': None, 'outdir2': None, 'recover_warped': True} ['__cmp__', '__doc__', '__init__', '__module__', '__repr__', '__str__', '_update', '_update_careful', '_update_loose', 'after', 'allmd5', 'band', 'before', 'bgmatch', 'center', 'compare_moon_all', 'cube', 'cube1', 'dataset', 'dec', 'do_rebin', 'dsky', 'ensure_value', 'epoch', 'force', 'frame0', 'height', 'hi_lo_rej', 'int_gz', 'maxmem', 'md5', 'medfilt', 'minmax', 'nframes', 'no_sanity_check', 'outdir', 'outdir2', 'output_masks', 'pdf', 'pixscale', 'plotprefix', 'plots', 'plots2', 'preprocess', 'ra', 'rchi_fraction', 'read_file', 'read_module', 'recover_warped', 'reference_dir', 'size', 'threads', 'threads1', 'tile', 'try_download', 'use_zp_meta', 'warp_all', 'width', 'wishlist'] Running on host: pan.rc.fas.harvard.edu Running as user: ameisner MKL_NUM_THREADS: 1 Tile RA,Dec 179.2 -19.7 Looking for file 1792m197-atlas.fits Reading /n/fink1/ameisner/unwise-coadds/etc_p9/wise_allsky_4band_p3as_cdd.fits Read 18240 Atlas tiles cosdec: 0.941464661246 Cut to 1 Atlas tiles near RA,Dec box Wrote 1792m197-atlas.fits Read 7904566 WISE L1b frames from /n/fink1/ameisner/unwise-coadds/etc_p9/WISE-index-L1b_w2.fits figuring out which frames are nearby took 4.33437490463 seconds Cut to 1240 WISE frames near RA,Dec box Reading /n/fink1/ameisner/unwise-coadds/etc_p9/WISE-l1b-metadata-4band.fits Read 1491686 from /n/fink1/ameisner/unwise-coadds/etc_p9/WISE-l1b-metadata-4band.fits Cut to 857 near RA,Dec box Matched 217 WISE-index-L1b scan_id: |S6 frame_num: >i8 WISE-metadata scan_id: |S6 frame_num: >i8 Cut to 204 matching scan/frame Band 1 : 0 Band 2 : 204 Band 3 : 0 Band 4 : 0 Reading /n/fink1/ameisner/unwise-coadds/etc_p9/WISE-l1b-metadata-3band.fits Read 392879 from /n/fink1/ameisner/unwise-coadds/etc_p9/WISE-l1b-metadata-3band.fits Cut to 0 near RA,Dec box Reading /n/fink1/ameisner/unwise-coadds/etc_p9/WISE-l1b-metadata-2band.fits Read 901271 from /n/fink1/ameisner/unwise-coadds/etc_p9/WISE-l1b-metadata-2band.fits Cut to 1016 near RA,Dec box Matched 295 WISE-index-L1b scan_id: |S6 frame_num: >i8 WISE-metadata scan_id: |S6 frame_num: >i8 Cut to 271 matching scan/frame Band 1 : 0 Band 2 : 271 Reading /n/fink1/ameisner/unwise-coadds/etc_p9/WISE-l1b-metadata-neowiser.fits Read 2497867 from /n/fink1/ameisner/unwise-coadds/etc_p9/WISE-l1b-metadata-neowiser.fits Cut to 1534 near RA,Dec box Matched 423 WISE-index-L1b scan_id: |S6 frame_num: >i8 WISE-metadata scan_id: |S6 frame_num: >i2 Cut to 380 matching scan/frame Band 1 : 0 Band 2 : 380 Reading /n/fink1/ameisner/unwise-coadds/etc_p9/WISE-l1b-metadata-neowiser2.fits Read 2620896 from /n/fink1/ameisner/unwise-coadds/etc_p9/WISE-l1b-metadata-neowiser2.fits Cut to 1656 near RA,Dec box Matched 416 WISE-index-L1b scan_id: |S6 frame_num: >i8 WISE-metadata scan_id: |S6 frame_num: >i4 Cut to 385 matching scan/frame Band 1 : 0 Band 2 : 385 1240 of 1240 matched to metadata tables Doing coadd tile 1792m197 band 2 Coadd tile 1792m197 RA,Dec 179.200353 -19.687778 Band 2 "git describe" version info: v3.5-15-gdebad9c will attempt to recover frames using per-quadrant polynomial warps ... mkdir /n/fink2/ameisner/unwise-coadds/tr_neo2/e000/179/1792m197 Coadd IWC polygon: [[ 7.89479167e-01 3.18085660e-15] [ 7.89479167e-01 7.89479167e-01] [ 0.00000000e+00 7.89479167e-01] [ -7.89479167e-01 7.89479167e-01] [ -7.89479167e-01 -6.36171320e-15] [ -7.89479167e-01 -7.89479167e-01] [ -9.94017687e-17 -7.89479167e-01] [ 7.89479167e-01 -7.89479167e-01]] Found 1240 WISE frames in range and in band W2 cut to 972 in RA,Dec box Defined epoch breaks [55284.81127749, 55465.032709755003, 56109.436509464998, 56746.185776525002, 56928.313990170005, 57107.709382100002] Found 6 epoch breaks Cut to 1 within epoch Cut to 1 frames before MJD 100000.0 Cut out qual_frame = 0; 1 remaining Cut out planets != 0; 1 remaining 0 of 1 frames are not moon_masked /n/sw/numpy-1.6.2_python-2.7.1/lib/python2.7/site-packages/numpy/core/fromnumeric.py:2374: RuntimeWarning: invalid value encountered in divide return mean(axis, dtype, out) /n/sw/numpy-1.6.2_python-2.7.1/lib/python2.7/site-packages/numpy/core/fromnumeric.py:2374: RuntimeWarning: invalid value encountered in double_scalars return mean(axis, dtype, out) Median nan MAD nan 1 of 1 moon-masked frames have large pixel variance Cut to 1 on moon Frames: 0 00716a 59 MJD 55203.2237211 Reading L1b top-level directory locations from /n/fink1/ameisner/unwise-coadds/etc_p9/l1b_dirs.yml L1b top-level directories are: {'2band': '/n/fink2/ameisner/staging/postcryo/2band_p1bm_frm', '3band': '/n/fink2/ameisner/staging/cryo_3band/3band_p1bm_frm', '4band': '/n/fink2/ameisner/staging/allsky/4band_p1bm_frm', 'missing': 'merge_p1bm_frm', 'neo1': '/n/fink2/ameisner/staging/neowiser/p1bm_frm', 'neo2': '/n/fink2/ameisner/staging/neowiser2/p1bm_frm'} 1 of 1 scan 00716a frame 59 band 2 intfn /n/fink2/ameisner/staging/allsky/4band_p1bm_frm/6a/00716a/059/00716a059-w2-int-1b.fits XTENSION: '(null)' Image extent: [ 871 1015 165 567] Coadd extent: [ 0 161 0 351] Quadrant 2 does not intersect target Quadrant 3 does not intersect target Total pixels in coadd space: 57024.0 Cut to 1 frames intersecting target Up to coadd_wise: Wall: 0.14 s, CPU: 0.06 s, VmPeak: 845 MB, VmSize: 376 MB, VmRSS: 59 MB, VmData: 50 MB, maxrss: 0.53974 MB Creating reference image from files: /n/fink2/ameisner/unwise-coadds/neo2_ref_w2/179/1792m197/unwise-1792m197-w2-img-u.fits /n/fink2/ameisner/unwise-coadds/neo2_ref_w2/179/1792m197/unwise-1792m197-w2-std-u.fits /n/fink2/ameisner/unwise-coadds/neo2_ref_w2/179/1792m197/unwise-1792m197-w2-n-u.fits Reading zero point polynomial coefficients for phase : 3band, W2 Reading zero point polynomial coefficients for phase : 2band, W2 Reading zero point polynomial coefficients for phase : neowiser, W2 Reading zero point polynomial coefficients for phase : neowiser2, W2 Warping all quadrants relative to reference image Creating per-quadrant FirstRoundImage objects Creating per-quadrant FirstRoundImage objects Coadd round 1, image 1 of 1 intfn /n/fink2/ameisner/staging/allsky/4band_p1bm_frm/6a/00716a/059/00716a059-w2-int-1b.fits uncfn /n/fink2/ameisner/staging/allsky/4band_p1bm_frm/6a/00716a/059/00716a059-w2-unc-1b.fits.gz maskfn /n/fink2/ameisner/staging/allsky/4band_p1bm_frm/6a/00716a/059/00716a059-w2-msk-1b.fits.gz Zeropoint: 19.596 -> scale 14.5077478909 sig1: 6.7958984375 Pixels to patch: 676 Estimated sky: 163.059796251 Image median: 166.651290894 Image median w/ noise: 167.121612549 Pixels in range: 23072 Wall: 1.20 s, CPU: 1.11 s, VmPeak: 845 MB, VmSize: 571 MB, VmRSS: 254 MB, VmData: 244 MB, maxrss: 0.53974 MB not rebinning !! time to fit warp = 0.000213861465454 seconds, number of pixels = 92 total time in do_one_warp 0.000695943832397 seconds, number of pixels = 92 Subtracted polynomial warp from quadrant deleting x, y coordinates for quadrant 1 Skipping quadrant 2 Skipping quadrant 3 time to fit warp = 0.00931096076965 seconds, number of pixels = 22707 total time in do_one_warp 0.027302980423 seconds, number of pixels = 22707 deleting x, y coordinates for quadrant 4 deleting x, y coordinates for quadrant -1 Accumulating first-round coadds... Wall: 0.00 s, CPU: 0.00 s, VmPeak: 845 MB, VmSize: 571 MB, VmRSS: 254 MB, VmData: 244 MB, maxrss: 0.53974 MB Before garbage collection: Wall: 0.00 s, CPU: 0.01 s, VmPeak: 845 MB, VmSize: 635 MB, VmRSS: 318 MB, VmData: 308 MB, maxrss: 0.53974 MB After garbage collection: Wall: 0.02 s, CPU: 0.03 s, VmPeak: 845 MB, VmSize: 635 MB, VmRSS: 318 MB, VmData: 308 MB, maxrss: 0.53974 MB Coadd round 2, image 1 of 1 Wall: 0.00 s, CPU: 0.00 s, VmPeak: 845 MB, VmSize: 827 MB, VmRSS: 510 MB, VmData: 500 MB, maxrss: 0.53974 MB Before garbage collection: Wall: 0.00 s, CPU: 0.00 s, VmPeak: 845 MB, VmSize: 827 MB, VmRSS: 510 MB, VmData: 500 MB, maxrss: 0.53974 MB After garbage collection: Wall: 0.02 s, CPU: 0.02 s, VmPeak: 845 MB, VmSize: 827 MB, VmRSS: 510 MB, VmData: 500 MB, maxrss: 0.53974 MB /n/pan/www/ameisner/software/astrometry.net-0.64/astrometry/util/miscutils.py:56: RuntimeWarning: invalid value encountered in double_scalars mx = -X[1] / (2. * X[2]) Estimated coadd sky: nan coadd_wise: Wall: 3.69 s, CPU: 3.04 s, VmPeak: 959 MB, VmSize: 639 MB, VmRSS: 322 MB, VmData: 313 MB, maxrss: 0.657388 MB Updating metadata based on quadrant SecondRoundImage objects. Traceback (most recent call last): File "/n/fink1/ameisner/unwise-coadds/unwise_coadd.py", line 3135, in sys.exit(main()) File "/n/fink1/ameisner/unwise-coadds/unwise_coadd.py", line 3129, in main output_masks=opt.output_masks): File "/n/fink1/ameisner/unwise-coadds/unwise_coadd.py", line 1229, in one_coadd fitsio.write(ofn, coim.astype(np.float32), header=hdr, clobber=True, extname='coadded image, outliers patched') File "/n/pan/www/ameisner/software/lib/python2.7/site-packages/fitsio/fitslib.py", line 215, in write **keys) File "/n/pan/www/ameisner/software/lib/python2.7/site-packages/fitsio/fitslib.py", line 443, in write header=header) File "/n/pan/www/ameisner/software/lib/python2.7/site-packages/fitsio/fitslib.py", line 497, in write_image self[-1].write_keys(header) File "/n/pan/www/ameisner/software/lib/python2.7/site-packages/fitsio/fitslib.py", line 1281, in write_keys self.write_key(name,value,comment=comment) File "/n/pan/www/ameisner/software/lib/python2.7/site-packages/fitsio/fitslib.py", line 1215, in write_key str(comment)) IOError: FITSIO status = 402: bad float to string conversion Error in ffd2e: double value is a NaN or INDEF