========================================= The start-up time of CASA may vary depending on whether the shared libraries are cached or not. ========================================= ]0;IPython: 2018S1SEQUOIACommissioning/92998IPython 5.1.0 -- An enhanced Interactive Python. CASA 5.8.0-109 -- Common Astronomy Software Applications Found an existing telemetry logfile: /home/teuben/.casa/casastats-580-109-1982ac1bc21613821e-20220204-015000.log Telemetry initialized. Telemetry will send anonymized usage statistics to NRAO. You can disable telemetry by adding the following line to your ~/.casarc file: EnableTelemetry: False Adding au casaVersion = 5.8.0-109 imported casatasks and tools using taskinit * The import of casa items did not complete. You may need to update the version of analysisUtils that you are using. See https://casaguides.nrao.edu/index.php/Analysis_Utilities au: 1.4901, casa: 5.8.0 QAC: Root /home/teuben/.casa/QAC QAC: Load contrib/tp2vis.py tp2vis for CASA5 [18-jun-2020 PJT] QAC: Skip distribute/tp2vis.py QAC: Skip tp2vis/tp2vis.py QAC: Load src/qac.py Loading QAC 27-oct-2020 QAC: Load src/ssc.py QAC: Load src/plot.py qac: version 27-oct-2020 qac_root: /home/teuben/.casa/QAC casa:5.8.0 data:/an/terra_astromake/opt/casa/casa-release-5.8.0-109.el7/data ('QAC: ', None) --> CrashReporter initialized. WARNING: No ASAP; ASAPLineFinder cannot function. ('ADMIT1: Version ', '21-dec-2017', 'loglevel ', 15) ('Removing previous results from ', 'Ori-KL_92998.nf.admit') ('Skipping ', 'admit1.apar') ('Found parameter file to execfile', 'Ori-KL_92998.nf.fits.apar') ('Skipping ', '') ('_ADMIT_ROOT=', '/lma1/teuben/lmtoy/admit') ADMIT basedir = /lma1/teuben/lmtoy/examples/2018S1SEQUOIACommissioning/92998/Ori-KL_92998.nf.admit ADMIT root = /lma1/teuben/lmtoy/admit ADMIT version = 1.0.8.5 ADMIT flowcount = 1 stale = 0 ('Starting a new ADMIT using', 'Ori-KL_92998.nf.fits') ('Found parameter file to cp:', 'Ori-KL_92998.nf.fits.apar') INFO : AT.py : Setting {'basename': 'x', 'file': 'Ori-KL_92998.nf.fits'} for Ingest_AT ('PJT0 _taskid', -1) ('PJT0 _taskid', -1) ('PJT0 _taskid', -1) ('PJT0 _taskid', -1) INFO : AT.py : Setting 'ppp' = True for CubeStats_AT ('PJT0 _taskid', -1) ('PJT0 _taskid', -1) INFO : AT.py : Setting 'numsigma' = 4.0 for CubeSum_AT INFO : AT.py : Setting 'sigma' = 99.0 for CubeSum_AT ('PJT0 _taskid', -1) ('PJT0 _taskid', -1) INFO : AT.py : Setting 'numsigma' = 4.0 for SFind2D_AT INFO : Admit.py : ADMIT run() called [flowcount 1] INFO : INFO : INFO : Executing Ingest_AT - '' (V1.2.13) INFO : INFO : INFO : Run using the following settings: INFO : box : [] INFO : file : Ori-KL_92998.nf.fits INFO : smooth : [] INFO : basename : x INFO : mask : True INFO : pb : INFO : edge : [] INFO : restfreq : -1.0 INFO : vlsr : -999999.9 INFO : usepb : True INFO : TIMING : Ingest ADMIT [ 2.12931500e+00 1.64394047e+09] TIMING : Ingest BEGIN [ 0. 0.] INFO : Ingest_AT.py : OBJECT: Ori-KL SHAPE: [ 33 33 515] ALMA OBJECT = Ori-KL ALMA DATE-OBS = 2021-12-17T02:07:27 INFO : Ingest_AT.py : basename=x INFO : utils.py : OSTYPE: linux TIMING : Ingest start [ 1.42768000e-01 1.48476839e-01 1.41080078e+03 2.17679688e+02] TIMING : Ingest importfits [ 9.45070000e-02 1.43806219e-01 1.41067578e+03 2.19011719e+02] WARNING : Ingest_AT.py : Adding dummy STOKES-I axis TIMING : Ingest adddegaxes [ 5.00350000e-02 7.31279850e-02 1.41099609e+03 2.19937500e+02] TIMING : Ingest summary-0 [ 1.30180000e-02 1.30798817e-02 1.41099609e+03 2.19937500e+02] TIMING : Ingest mask [ 2.04990000e-02 7.67390728e-02 1.41099609e+03 2.20277344e+02] TIMING : Ingest summary-1 [ 1.33340000e-02 2.43330002e-02 1.41099609e+03 2.20277344e+02] TIMING : Ingest statistics [ 2.22190000e-02 2.22809315e-02 1.41632422e+03 2.26339844e+02] 2022-02-04 02:07:50 WARN importfits::image::commonbeam (file ../../tools/images/image_cmpt.cc, line 753) This image only has one beam, so just returning that INFO : Ingest_AT.py : COMMONBEAM[3] {'major': {'value': 17.2500012, 'unit': 'arcsec'}, 'pa': {'value': 0.0, 'unit': 'deg'}, 'minor': {'value': 17.2500012, 'unit': 'arcsec'}} INFO : Ingest_AT.py : BASICS: [shape] npts min max: [ 33 33 515 1] 142140 -1.533444 91.802063 INFO : Ingest_AT.py : S/N (all data): 139.013621 INFO : Ingest_AT.py : GOOD PIXELS: 142140/560835 (25.344353% good or 74.655647% bad) WARNING : Ingest_AT.py : MASKS: ['mask0'] REGRESSION : CUBE: -1.53344 91.8021 0.660382 33 33 515 74.655647 INFO : Ingest_AT.py : TELESCOPE: LMT INFO : Ingest_AT.py : OBJECT: Ori-KL INFO : Ingest_AT.py : REFFREQTYPE: LSRK INFO : Ingest_AT.py : RA Axis 1: 0.000000 -7.499999 16.000000 INFO : Ingest_AT.py : DEC Axis 2: 0.000000 7.499999 16.000000 INFO : Ingest_AT.py : VLSRv = 0.000000 (from source catalog) INFO : Ingest_AT.py : VLSRz = 0.000000 +/- 0.000000 1 values: [ 0.] INFO : Ingest_AT.py : Freq Orig Axis 3: 86.3455 -0.000390625 0 INFO : Ingest_AT.py : Cube Orig Axis 3: type=Frequency velocity increment=1.357859 km/s @ fc=86.245076 fw=-0.201172 GHz INFO : Ingest_AT.py : RESTFREQ: 86.2434 86.2434 -1 INFO : Ingest_AT.py : VLSRc= -5.678937 VLSRf= 0.000000 VLSRv= 0.000000 VLSRz= 0.000000 WIDTH= 699.297385 INFO : Ingest_AT.py : VLSR = -5.678937 errs = 0.000000 0.000000 0.000000 width = 1.357859 TIMING : Ingest done [ 5.31228000e-01 5.36417961e-01 1.41932422e+03 2.29839844e+02] TIMING : Ingest END [ 0.899641 1.05032897] INFO : AT.py : BDP_OUT[0] = SpwCube_BDP x.im INFO : INFO : INFO : Executing CubeStats_AT - '' (V1.2.3) INFO : INFO : INFO : Run using the following settings: INFO : ppp : True INFO : psample : -1 INFO : robust : [] INFO : maxvrms : 2.0 INFO : TIMING : CubeStats ADMIT [ 3.08925800e+00 1.64394047e+09] TIMING : CubeStats BEGIN [ 0. 0.] TIMING : CubeStats imval [ 5.06090000e-02 5.08370399e-02 1.41807422e+03 2.28980469e+02] TIMING : CubeStats start [ 1.23860000e-02 1.24309063e-02 1.41807422e+03 2.28980469e+02] TIMING : CubeStats imstat0 [ 4.35420000e-02 4.37400341e-02 1.41485156e+03 2.26453125e+02] TIMING : CubeStats imstat1 [ 6.20860000e-02 6.23519421e-02 1.41698828e+03 2.28535156e+02] INFO : CubeStats_AT.py : sigma varies from 0.037103 to 0.427935; 515/515 channels ok WARNING : CubeStats_AT.py : sigma varies too much, going to clip to 0.0742055 (11.5338 > 2) INFO : CubeStats_AT.py : Computing MaxPos for PeakPointPlot /an/terra_astromake/opt/casa/casa-release-5.8.0-109.el7/lib/python2.7/site-packages/numpy/ma/core.py:4144: UserWarning: Warning: converting a masked element to nan. warnings.warn("Warning: converting a masked element to nan.") TIMING : CubeStats ppp [ 5.74970000e-01 5.76937199e-01 1.41698828e+03 2.28777344e+02] INFO : CubeStats_AT.py : CubeMax: 91.802063 @ [ 17 16 258 0] INFO : CubeStats_AT.py : CubeMin: -1.533444 @ [ 16 19 258 0] INFO : CubeStats_AT.py : CubeRMS: 0.052124 INFO : CubeStats_AT.py : RMS Sanity check 12.669329 WARNING : CubeStats_AT.py : RMS sanity check = 12.669329. Either bad sidelobes, lotsa signal, or both REGRESSION : CST: 0.052124 12.669329 INFO : CubeStats_AT.py : mean,rms,S/N=0.021580 0.052124 1761.209291 INFO : CubeStats_AT.py : RMS BAD VARIATION RATIO: 1.513547 1.059747 /an/terra_astromake/opt/casa/casa-release-5.8.0-109.el7/lib/python2.7/site-packages/matplotlib/collections.py:446: FutureWarning: elementwise comparison failed; returning scalar instead, but in the future will perform elementwise comparison if self._edgecolors == 'face': TIMING : CubeStats plotting [ 2.281325 2.45253181 1436.45703125 246.74609375] TIMING : CubeStats done [ 1.35950000e-02 1.36311054e-02 1.43645703e+03 2.46746094e+02] TIMING : CubeStats summary [ 1.21730000e-02 1.22239590e-02 1.43645703e+03 2.46746094e+02] TIMING : CubeStats END [ 3.062858 3.23689198] INFO : AT.py : BDP_OUT[0] = CubeStats_BDP x.cst INFO : INFO : INFO : Executing CubeSum_AT - '' (V1.2.2) INFO : INFO : INFO : Run using the following settings: INFO : linesum : True INFO : numsigma : 4.0 INFO : zoom : 1 INFO : pad : 5 INFO : sigma : 99.0 INFO : TIMING : CubeSum ADMIT [ 6.21195000e+00 1.64394047e+09] TIMING : CubeSum BEGIN [ 0. 0.] INFO : CubeSum_AT.py : Using constant sigma = 0.052124 TIMING : CubeSum start [ 5.19270000e-02 5.20749092e-02 1.43645703e+03 2.46746094e+02] TIMING : CubeSum immoments [ 4.82190000e-02 1.18519068e-01 1.43859375e+03 2.49082031e+02] 2022-02-04 02:07:54 WARN immoments::image::open (file ../../tools/images/image_cmpt.cc, line 4125) Another image is already open, closing first TIMING : CubeSum statistics [ 1.91640000e-02 1.92389488e-02 1.43859375e+03 2.49082031e+02] INFO : CubeSum_AT.py : Total flux: 304426.920425 (sum=5412.035873) REGRESSION : CSM: [304426.92042471049, 5412.0358727350831] TIMING : CubeSum implot [ 1.41955000e-01 5.82964301e+00 1.52316016e+03 2.61703125e+02] TIMING : CubeSum getdata [ 2.06930000e-02 2.07490921e-02 1.52316016e+03 2.61703125e+02] TIMING : CubeSum done [ 2.39121000e-01 2.57791996e-01 1.50848828e+03 2.47191406e+02] TIMING : CubeSum END [ 0.534583 6.31155896] INFO : AT.py : BDP_OUT[0] = Moment_BDP x.csm INFO : INFO : INFO : Executing SFind2D_AT - 'csm' (V1.2.2) INFO : INFO : INFO : Run using the following settings: INFO : nmax : 30 INFO : numsigma : 4.0 INFO : snmax : 35.0 INFO : region : INFO : zoom : 1 INFO : robust : ['hin', 1.5] INFO : sigma : -1.0 INFO : TIMING : SFind2D ADMIT [ 6.80639600e+00 1.64394048e+09] TIMING : SFind2D BEGIN [ 0. 0.] ('BDP', 0, ) ('BDP', 1, ) ('BDP', 2, ) No PB given TIMING : SFind2D start [ 3.32570000e-02 3.33549976e-02 1.50848828e+03 2.47191406e+02] ('ROBUST:', {'fence': 1.5, 'algorithm': 'hinges-fences'}) TIMING : SFind2D imstat [ 4.03950000e-02 4.05960083e-02 1.50848828e+03 2.47191406e+02] INFO : AT.py : Setting 'sigma' = 6.45128333512 for SFind2D_AT INFO : SFind2D_AT.py : sigma, dmin, dmax, snmax, cutoff 6.45128 -13.2911 1074.73 35 0.0285714 WARNING : SFind2D_AT.py : Temporarely patching your K.km/s units to Jy/beam for ia.findsources() 2022-02-04 02:08:00 WARN imstat::ImageRegion::fromRecord Fit did not converge, resorting to default shape estimates 2022-02-04 02:08:00 WARN imstat::ImageRegion::fromRecord Fit did not converge, resorting to default shape estimates 2022-02-04 02:08:00 WARN imstat::ImageRegion::fromRecord Fit did not converge, resorting to default shape estimates 2022-02-04 02:08:00 WARN imstat::ImageRegion::fromRecord Fit did not converge, resorting to default shape estimates 2022-02-04 02:08:01 WARN imstat::ImageRegion::fromRecord Fit did not converge, resorting to default shape estimates 2022-02-04 02:08:01 WARN imstat::ImageRegion::fromRecord There is no Stokes coordinate in the CoordinateSystem - assuming Stokes I TIMING : SFind2D findsources [ 5.77460000e-01 5.86024046e-01 1.50848828e+03 2.47191406e+02] INFO : SFind2D_AT.py : Right Ascen. Declination X(pix) Y(pix) Peak Flux Major Minor PA SNR INFO : SFind2D_AT.py : K.km/s Jy arcsec arcsec deg INFO : SFind2D_AT.py : 23:59:59.957 -00.00.00.42 16.08 15.94 0 1.01e+03 15.452 10.193 79.6 0.0 INFO : SFind2D_AT.py : 00:00:00.003 -00.00.53.13 15.99 8.92 6.55 100 37.500 37.500 90.0 1.0 INFO : SFind2D_AT.py : 00:00:02.189 -00.00.39.20 11.62 10.77 13.5 65.8 37.500 37.500 180.0 2.1 INFO : SFind2D_AT.py : 23:59:57.762 -00.00.35.49 20.48 11.27 12.9 60.7 37.500 37.500 0.0 2.0 INFO : SFind2D_AT.py : 00:00:01.808 -00.00.56.19 12.38 8.51 8.74 43.4 37.500 37.500 90.0 1.4 INFO : SFind2D_AT.py : 00:00:03.759 -00.00.04.21 8.48 15.44 8.5 24.6 37.500 37.500 0.0 1.3 INFO : SFind2D_AT.py : Wrote ds9.reg TIMING : SFind2D table [ 2.61904000e-01 2.66238928e-01 1.50848828e+03 2.47781250e+02] REGRESSION : CONTFLUX: 6 1308.17 INFO : SFind2D_AT.py : Fitted Gaussian size; NOT deconvolved source size. INFO : SFind2D_AT.py : Restoring Beam: Major axis: 17.3 arcsec , Minor axis: 17.3 arcsec , PA: 0.0 deg TIMING : SFind2D done [ 2.73647000e-01 2.92111158e-01 1.50937109e+03 2.48664062e+02] TIMING : SFind2D END [ 1.200223 1.23193717] INFO : AT.py : BDP_OUT[0] = SourceList_BDP x-csm.sl TIMING : ADMITrun END [ 6.26659 12.45491886] INFO : Admit.py : ADMIT run() finished [flowcount 1] [cpu 6.26659 12.4549 ] ('N sources in CSM:', 6) ('PJT0 _taskid', -1) ('PJT0 _taskid', -1) ('PJT0 _taskid', -1) ('PJT0 _taskid', -1) INFO : AT.py : Setting 'csub' = [0, 0] for LineSegment_AT INFO : Admit.py : ADMIT run() called [flowcount 1] INFO : INFO : INFO : Executing CubeSpectrum_AT - '' (V1.2.5) INFO : INFO : INFO : Run using the following settings: INFO : xaxis : INFO : sources : [0] INFO : pos : [] INFO : TIMING : CubeSpectrum ADMIT [ 8.21122700e+00 1.64394048e+09] TIMING : CubeSpectrum BEGIN [ 0. 0.] INFO : CubeSpectrum_AT.py : CubeStats::maxpos,val=[17, 16, 258],91.802063 TIMING : CubeSpectrum CubeStats-pos [ 4.84770000e-02 4.86199856e-02 1.50937109e+03 2.48667969e+02] INFO : CubeSpectrum_AT.py : CubeSum::maxpos,val=[17, 16],1074.725830 TIMING : CubeSpectrum Moment-pos [ 4.19460000e-02 4.20689583e-02 1.50937109e+03 2.48667969e+02] INFO : CubeSpectrum_AT.py : SourceList::maxpos,val=('23h59m59.957s', '-00d00m00.42s'),0.000000 TIMING : CubeSpectrum SourceList-pos [ 3.57130000e-02 3.58030796e-02 1.50937109e+03 2.48667969e+02] TIMING : CubeSpectrum open [ 1.21520000e-02 1.21870041e-02 1.50937109e+03 2.48667969e+02] TIMING : CubeSpectrum imval [ 4.18170000e-02 4.20169830e-02 1.50937109e+03 2.49019531e+02] TIMING : CubeSpectrum imhead [ 4.11730000e-02 4.13160324e-02 1.50937109e+03 2.49019531e+02] TIMING : CubeSpectrum imval [ 2.37556000e-01 2.56230831e-01 1.50957422e+03 2.49222656e+02] REGRESSION : CSP: [91.80206298828125, 0.0] INFO : CubeSpectrum_AT.py : Writing 2 testCubeSpectrum tables TIMING : CubeSpectrum done [ 2.54522000e-01 2.73374081e-01 1.50870312e+03 2.48351562e+02] TIMING : CubeSpectrum summary [ 1.23600000e-02 1.23898983e-02 1.50870312e+03 2.48351562e+02] TIMING : CubeSpectrum END [ 0.737854 0.77618098] INFO : AT.py : BDP_OUT[0] = CubeSpectrum_BDP x.csp INFO : INFO : INFO : Executing LineSegment_AT - '' (V1.2.3) INFO : INFO : INFO : Run using the following settings: INFO : numsigma : 5.0 INFO : minchan : 4 INFO : edgechannels : 0 INFO : smooth : [] INFO : recalcnoise : False INFO : maxgap : 3 INFO : csub : [0, 0] INFO : iterate : True INFO : segment : ADMIT INFO : TIMING : LineSegment ADMIT [ 9.00917500e+00 1.64394048e+09] TIMING : LineSegment BEGIN [ 0. 0.] TIMING : LineSegment start [ 2.43630000e-02 2.44271755e-02 1.50870312e+03 2.48351562e+02] PJT-robust unexpected array with no data PJT-robust unexpected array with no data ERROR : Admit.py : Project run() failed; (, TypeError("ufunc 'bitwise_or' not supported for the input types, and the inputs could not be safely coerced to any supported types according to the casting rule ''safe''",), ) : saving state... Traceback (most recent call last): File "/an/terra_astromake/opt/casa/casa-release-5.8.0-109.el7/lib/python2.7/init_welcome.py", line 33, in execfile(__candidates[0]) File "/lma1/teuben/lmtoy/admit/admit/test/admit1.py", line 438, in a.run() File "/lma1/teuben/lmtoy/admit/admit/Admit.py", line 797, in run self.fm.run() File "/lma1/teuben/lmtoy/admit/admit/FlowManager.py", line 367, in run task.execute(args) File "/lma1/teuben/lmtoy/admit/admit/AT.py", line 1795, in execute self.run() File "/lma1/teuben/lmtoy/admit/admit/at/LineSegment_AT.py", line 243, in run self.getkey("recalcnoise"), basicsegment) File "/lma1/teuben/lmtoy/admit/admit/util/specutil.py", line 183, in getspectrum sep, cut, noise, mean = sfinder.find() File "/lma1/teuben/lmtoy/admit/admit/util/segmentfinder/SegmentFinder.py", line 137, in find seg, cut, noi, mn = sfinder.find() File "/lma1/teuben/lmtoy/admit/admit/util/segmentfinder/ADMITSegmentFinder.py", line 334, in find dm = ma.masked_outside(dr, dmin, dmax) File "/an/terra_astromake/opt/casa/casa-release-5.8.0-109.el7/lib/python2.7/site-packages/numpy/ma/core.py", line 2117, in masked_outside condition = (xf < v1) | (xf > v2) TypeError: ufunc 'bitwise_or' not supported for the input types, and the inputs could not be safely coerced to any supported types according to the casting rule ''safe''