========================================= The start-up time of CASA may vary depending on whether the shared libraries are cached or not. ========================================= ]0;IPython: 2018-S1-MU-46/85824IPython 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 ', 'Region_J-K_85824.nfs.admit') ('Skipping ', 'admit1.apar') ('Found parameter file to execfile', 'Region_J-K_85824.nfs.fits.apar') ('Skipping ', '') ('_ADMIT_ROOT=', '/lma1/teuben/lmtoy/admit') ADMIT basedir = /lma1/teuben/lmtoy/examples/2018-S1-MU-46/85824/Region_J-K_85824.nfs.admit ADMIT root = /lma1/teuben/lmtoy/admit ADMIT version = 1.0.8.5 ADMIT flowcount = 1 stale = 0 ('Starting a new ADMIT using', 'Region_J-K_85824.nfs.fits') ('Found parameter file to cp:', 'Region_J-K_85824.nfs.fits.apar') INFO : AT.py : Setting {'basename': 'x', 'file': 'Region_J-K_85824.nfs.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 : Region_J-K_85824.nfs.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.06440800e+00 1.64394311e+09] TIMING : Ingest BEGIN [ 0. 0.] INFO : Ingest_AT.py : OBJECT: Region_J-K SHAPE: [ 91 91 689] ALMA OBJECT = Region_J-K ALMA DATE-OBS = 2019-11-01T03:09:48 INFO : Ingest_AT.py : basename=x INFO : utils.py : OSTYPE: linux TIMING : Ingest start [ 1.96195000e-01 2.03284979e-01 1.41021875e+03 2.17093750e+02] TIMING : Ingest importfits [ 1.34912000e-01 2.28847980e-01 1.41095703e+03 2.18320312e+02] WARNING : Ingest_AT.py : Adding dummy STOKES-I axis TIMING : Ingest adddegaxes [ 6.40040000e-02 2.53926039e-01 1.41556641e+03 2.23820312e+02] TIMING : Ingest summary-0 [ 1.42170000e-02 1.43020153e-02 1.41556641e+03 2.23820312e+02] TIMING : Ingest mask [ 4.32780000e-02 8.77859592e-02 1.41556641e+03 2.24167969e+02] TIMING : Ingest summary-1 [ 1.36850000e-02 1.78360939e-02 1.41556641e+03 2.24167969e+02] TIMING : Ingest statistics [ 6.62870000e-02 6.64389133e-02 1.41670312e+03 2.26398438e+02] 2022-02-04 02:51:53 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': 12.6500004, 'unit': 'arcsec'}, 'pa': {'value': 0.0, 'unit': 'deg'}, 'minor': {'value': 12.6500004, 'unit': 'arcsec'}} INFO : Ingest_AT.py : BASICS: [shape] npts min max: [ 91 91 689 1] 4869163 -0.090834 0.777210 INFO : Ingest_AT.py : S/N (all data): 31.780266 INFO : Ingest_AT.py : GOOD PIXELS: 4869163/5705609 (85.339935% good or 14.660065% bad) WARNING : Ingest_AT.py : MASKS: ['mask0'] REGRESSION : CUBE: -0.0908337 0.77721 0.0244557 91 91 689 14.660065 INFO : Ingest_AT.py : TELESCOPE: LMT INFO : Ingest_AT.py : OBJECT: Region_J-K INFO : Ingest_AT.py : REFFREQTYPE: LSRK INFO : Ingest_AT.py : RA Axis 1: 10.648750 -5.500001 45.000000 INFO : Ingest_AT.py : DEC Axis 2: 41.513700 5.500001 45.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: 115.25 0.000390625 0 INFO : Ingest_AT.py : Cube Orig Axis 3: type=Frequency velocity increment=-1.015921 km/s @ fc=115.384822 fw=0.269141 GHz INFO : Ingest_AT.py : RESTFREQ: 115.271 115.271 -1 INFO : Ingest_AT.py : VLSRc= -295.492064 VLSRf= 0.000000 VLSRv= 0.000000 VLSRz= 0.000000 WIDTH= -699.969569 INFO : Ingest_AT.py : VLSR = -295.492064 errs = 0.000000 0.000000 0.000000 width = -1.015921 TIMING : Ingest done [ 5.76561000e-01 5.84861994e-01 1.41945312e+03 2.29859375e+02] TIMING : Ingest END [ 1.121248 1.46942401] 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.24615500e+00 1.64394311e+09] TIMING : CubeStats BEGIN [ 0. 0.] TIMING : CubeStats imval [ 5.20410000e-02 5.28540611e-02 1.41870312e+03 2.29714844e+02] TIMING : CubeStats start [ 1.23660000e-02 1.24058723e-02 1.41870312e+03 2.29714844e+02] TIMING : CubeStats imstat0 [ 2.41431000e-01 2.42677212e-01 1.44747656e+03 2.58777344e+02] TIMING : CubeStats imstat1 [ 3.00041000e-01 3.00788879e-01 1.47939453e+03 2.90753906e+02] INFO : CubeStats_AT.py : sigma varies from 0.007625 to 0.027186; 689/689 channels ok WARNING : CubeStats_AT.py : sigma varies too much, going to clip to 0.0152508 (3.5652 > 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 [ 8.35340000e-01 8.38148117e-01 1.47939453e+03 2.90878906e+02] INFO : CubeStats_AT.py : CubeMax: 0.777210 @ [ 54 34 324 0] INFO : CubeStats_AT.py : CubeMin: -0.090834 @ [ 41 71 456 0] INFO : CubeStats_AT.py : CubeRMS: 0.014854 INFO : CubeStats_AT.py : RMS Sanity check 1.646395 WARNING : CubeStats_AT.py : RMS sanity check = 1.646395. Either bad sidelobes, lotsa signal, or both REGRESSION : CST: 0.014854 1.646395 INFO : CubeStats_AT.py : mean,rms,S/N=0.001907 0.014854 52.322875 INFO : CubeStats_AT.py : RMS BAD VARIATION RATIO: 2.099883 2.411432 /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.296366 2.49018884 1488.8984375 298.8125 ] TIMING : CubeStats done [ 1.35630000e-02 1.35970116e-02 1.48889844e+03 2.98812500e+02] TIMING : CubeStats summary [ 1.21940000e-02 1.22389793e-02 1.48889844e+03 2.98812500e+02] TIMING : CubeStats END [ 3.775476 3.97507095] 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 [ 7.08146700e+00 1.64394312e+09] TIMING : CubeSum BEGIN [ 0. 0.] INFO : CubeSum_AT.py : Using constant sigma = 0.014854 TIMING : CubeSum start [ 5.17730000e-02 5.24561405e-02 1.48889844e+03 2.98812500e+02] TIMING : CubeSum immoments [ 9.46210000e-02 1.67165041e-01 1.48889844e+03 2.99117188e+02] 2022-02-04 02:51:57 WARN immoments::image::open (file ../../tools/images/image_cmpt.cc, line 4125) Another image is already open, closing first TIMING : CubeSum statistics [ 1.99250000e-02 2.00028419e-02 1.48889844e+03 2.99117188e+02] INFO : CubeSum_AT.py : Total flux: 161555.511120 (sum=5340.676500) REGRESSION : CSM: [161555.51112046328, 5340.6764999175757] TIMING : CubeSum implot [ 1.43419000e-01 6.09378099e+00 1.56090234e+03 2.99160156e+02] TIMING : CubeSum getdata [ 2.13840000e-02 2.14509964e-02 1.56090234e+03 2.99160156e+02] TIMING : CubeSum done [ 2.56087000e-01 2.74775028e-01 1.56090234e+03 2.99222656e+02] TIMING : CubeSum END [ 0.601047 6.64351201] 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 [ 7.74229800e+00 1.64394312e+09] TIMING : SFind2D BEGIN [ 0. 0.] ('BDP', 0, ) ('BDP', 1, ) ('BDP', 2, ) No PB given TIMING : SFind2D start [ 3.36930000e-02 3.38041782e-02 1.56090234e+03 2.99222656e+02] ('ROBUST:', {'fence': 1.5, 'algorithm': 'hinges-fences'}) TIMING : SFind2D imstat [ 4.11240000e-02 4.13570404e-02 1.56090234e+03 2.99222656e+02] INFO : AT.py : Setting 'sigma' = 0.650839954277 for SFind2D_AT INFO : SFind2D_AT.py : sigma, dmin, dmax, snmax, cutoff 0.65084 -1.02803 13.035 35 0.199721 WARNING : SFind2D_AT.py : Temporarely patching your K.km/s units to Jy/beam for ia.findsources() 2022-02-04 02:52:04 WARN imstat::ImageRegion::fromRecord Fit did not converge, resorting to default shape estimates 2022-02-04 02:52:04 WARN imstat::ImageRegion::fromRecord There is no Stokes coordinate in the CoordinateSystem - assuming Stokes I TIMING : SFind2D findsources [ 2.46877000e-01 2.89003849e-01 1.56090234e+03 2.99222656e+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 : 00:42:38.777 +41.32.02.54 38.72 58.31 12.7 271 93.628 37.771 69.7 19.6 INFO : SFind2D_AT.py : 00:42:31.350 +41.29.47.75 53.89 33.81 10.7 175 61.125 38.719 42.8 16.4 INFO : SFind2D_AT.py : 00:42:43.190 +41.33.00.38 29.71 68.83 10.2 244 72.679 52.176 8.7 15.6 INFO : SFind2D_AT.py : 00:42:29.004 +41.29.03.70 58.68 25.80 8.42 162 74.418 42.744 19.9 12.9 INFO : SFind2D_AT.py : 00:42:25.028 +41.28.17.44 66.81 17.39 4.98 82.1 67.775 35.575 53.8 7.7 INFO : SFind2D_AT.py : 00:42:46.584 +41.32.40.34 22.78 65.19 2.68 12.5 27.500 27.500 90.0 4.1 INFO : SFind2D_AT.py : Wrote ds9.reg TIMING : SFind2D table [ 2.56037000e-01 2.60931969e-01 1.56090234e+03 2.99886719e+02] REGRESSION : CONTFLUX: 6 945.936 INFO : SFind2D_AT.py : Fitted Gaussian size; NOT deconvolved source size. INFO : SFind2D_AT.py : Restoring Beam: Major axis: 12.7 arcsec , Minor axis: 12.7 arcsec , PA: 0.0 deg TIMING : SFind2D done [ 2.81663000e-01 3.00103188e-01 1.56090234e+03 2.99937500e+02] TIMING : SFind2D END [ 0.873125 0.93897009] INFO : AT.py : BDP_OUT[0] = SourceList_BDP x-csm.sl TIMING : ADMITrun END [ 6.951805 13.66486883] INFO : Admit.py : ADMIT run() finished [flowcount 1] [cpu 6.95181 13.6649 ] ('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.82947000e+00 1.64394313e+09] TIMING : CubeSpectrum BEGIN [ 0. 0.] INFO : CubeSpectrum_AT.py : CubeStats::maxpos,val=[54, 34, 324],0.777210 TIMING : CubeSpectrum CubeStats-pos [ 4.79490000e-02 4.80649471e-02 1.56090234e+03 2.99937500e+02] INFO : CubeSpectrum_AT.py : CubeSum::maxpos,val=[39, 58],13.035013 TIMING : CubeSpectrum Moment-pos [ 4.18400000e-02 4.19740677e-02 1.56090234e+03 2.99937500e+02] INFO : CubeSpectrum_AT.py : SourceList::maxpos,val=('00h42m38.777s', '+41d32m02.54s'),12.747886 TIMING : CubeSpectrum SourceList-pos [ 3.58390000e-02 3.59530449e-02 1.56090234e+03 2.99937500e+02] TIMING : CubeSpectrum open [ 1.21210000e-02 1.21550560e-02 1.56090234e+03 2.99937500e+02] TIMING : CubeSpectrum imval [ 4.34690000e-02 4.36608791e-02 1.56090234e+03 3.00296875e+02] TIMING : CubeSpectrum imhead [ 4.16970000e-02 4.18510437e-02 1.56090234e+03 3.00296875e+02] TIMING : CubeSpectrum imval [ 2.31340000e-01 2.50246048e-01 1.56090234e+03 3.00304688e+02] TIMING : CubeSpectrum imval [ 2.65660000e-01 2.85055876e-01 1.56090234e+03 3.00367188e+02] REGRESSION : CSP: [0.5920255184173584, 0.77720963954925537, 0.5920255184173584] INFO : CubeSpectrum_AT.py : Writing 3 testCubeSpectrum tables TIMING : CubeSpectrum done [ 2.82629000e-01 3.01213026e-01 1.56090234e+03 3.00367188e+02] TIMING : CubeSpectrum summary [ 1.25210000e-02 1.25501156e-02 1.56090234e+03 3.00367188e+02] TIMING : CubeSpectrum END [ 1.027408 1.08510399] 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.91784400e+00 1.64394313e+09] TIMING : LineSegment BEGIN [ 0. 0.] TIMING : LineSegment start [ 2.42920000e-02 2.43589878e-02 1.56090234e+03 3.00367188e+02] INFO : LineSegment_AT.py : Attempting Continuum Subtraction for Input Spectra INFO : LineSegment_AT.py : Attempting Continuum Subtraction for Input CubeStats Spectra TIMING : LineSegment getspectrum [ 36.182786 36.27138591 1560.90234375 300.921875 ] INFO : LineSegment_AT.py : Detecting segments in CubeSpectrum based data WARNING : specutil.py : 0 [295, 334] WARNING : specutil.py : 1 [448, 452] WARNING : specutil.py : 2 [505, 511] WARNING : specutil.py : 3 [585, 590] WARNING : specutil.py : 0 [90, 94] WARNING : specutil.py : 1 [160, 163] WARNING : specutil.py : 2 [220, 226] WARNING : specutil.py : 3 [258, 261] WARNING : specutil.py : 4 [299, 302] WARNING : specutil.py : 5 [309, 339] WARNING : specutil.py : 6 [407, 411] WARNING : specutil.py : 7 [420, 424] WARNING : specutil.py : 8 [463, 479] WARNING : specutil.py : 9 [484, 499] WARNING : specutil.py : 10 [611, 622] WARNING : specutil.py : 11 [635, 640] WARNING : specutil.py : 12 [677, 685] WARNING : specutil.py : 0 [295, 334] WARNING : specutil.py : 1 [448, 452] WARNING : specutil.py : 2 [505, 511] WARNING : specutil.py : 3 [585, 590] INFO : LineSegment_AT.py : Detecting segments in CubeStats based data WARNING : specutil.py : 0 [152, 155] WARNING : specutil.py : 1 [281, 354] WARNING : specutil.py : 2 [506, 509] WARNING : specutil.py : 3 [567, 571] WARNING : specutil.py : 4 [597, 600] WARNING : specutil.py : 0 [649, 652] TIMING : LineSegment segment finder [ 5.168889 5.1820631 1560.90234375 300.921875 ] INFO : LineSegment_AT.py : Segment Coverage 153 / 689 = 0.222061 REGRESSION : LINESEG: [[115.3098216830261, 115.31099355810548], [115.3602123114395, 115.3887279383711], [115.44810294239309, 115.44927481747247], [115.47193106900717, 115.47349356911302], [115.483649819801, 115.48482169488037], [115.50396232117693, 115.50513419625632], [115.36568106180994, 115.38091543784191], [115.42544669085838, 115.42700919096423], [115.44771231736664, 115.4500560675254], [115.47896231948346, 115.48091544461577], [115.28560293138555, 115.28716543149139], [115.3129466832378, 115.31411855831718], [115.33638418482541, 115.33872793498418], [115.35122793583091, 115.3523998109103], [115.36724356191579, 115.36841543699516], [115.37114981218039, 115.3828685629742], [115.40943106477351, 115.41099356487935], [115.4145091901175, 115.41607169022335], [115.43130606625529, 115.43755606667867], [115.43950919181096, 115.44536856720786], [115.48911857017143, 115.49341544546249], [115.49849357080649, 115.50044669593878], [115.51489982191782, 115.5180248221295], [115.36568106180994, 115.38091543784191], [115.42544669085838, 115.42700919096423], [115.44771231736664, 115.4500560675254], [115.47896231948346, 115.48091544461577]] TIMING : LineSegment done [ 2.287808 2.41865277 1581.31640625 321.2734375 ] TIMING : LineSegment END [ 43.676097 43.90881586] INFO : AT.py : BDP_OUT[0] = LineSegment_BDP x.lseg TIMING : ADMITrun END [ 51.942277 59.00060201] INFO : Admit.py : ADMIT run() finished [flowcount 1] [cpu 51.9423 59.0006 ] Found 19 segments in LineSegment No ContinuumSub needed INFO : Admit.py : ADMIT run() called [flowcount 1] TIMING : ADMITrun END [ 52.043501 59.16218996] INFO : Admit.py : ADMIT run() finished [flowcount 1] [cpu 52.0435 59.1622 ] ('OBJECT = ', 'Region_J-K') ('PJT0 _taskid', -1) ('PJT0 _taskid', -1) INFO : AT.py : Setting 'csub' = [0, 0] for LineID_AT INFO : AT.py : Setting 'references' = etc/tier1_lines.list for LineID_AT INFO : Admit.py : ADMIT run() called [flowcount 1] INFO : INFO : INFO : Executing LineID_AT - '' (V1.2.7) INFO : INFO : INFO : Run using the following settings: INFO : numsigma : 5.0 INFO : force : [] INFO : tier1width : 0.0 INFO : online : False INFO : recomblevel : shallow INFO : references : etc/tier1_lines.list INFO : csub : [0, 0] INFO : identifylines : True INFO : iterate : True INFO : segment : ADMIT INFO : allowexotics : False INFO : minchan : 4 INFO : pattern : AUTO INFO : edgechannels : 0 INFO : smooth : [] INFO : recalcnoise : False INFO : vlsr : -999999.99 INFO : maxgap : 3 INFO : reject : [] INFO : method : {'PeakFinder': {'thresh': 0.0}} INFO : mode : ONE INFO : TIMING : LineID ADMIT [ 5.39463480e+01 1.64394317e+09] TIMING : LineID BEGIN [ 0. 0.] INFO : LineID_AT.py : Set vlsr = -295.49 for line identification. INFO : LineID_AT.py : Identifylines = True INFO : LineID_AT.py : Using vlsr = -295.492 INFO : LineID_AT.py : Attempting Continuum Subtraction for Input Spectra TIMING : LineID getspectrum-cubespecs [ 21.920838 21.97549295 1581.31640625 321.33984375] INFO : LineID_AT.py : Attempting Continuum Subtraction for Input CubeStats Spectra TIMING : LineID getspectrum-cubestats [ 14.593051 14.62913609 1581.31640625 321.33984375] TIMING : LineID segment finder [ 1.42520000e-02 1.42920017e-02 1.58131641e+03 3.21339844e+02] INFO : LineID_AT.py : Detecting segments in CubeSpectrum based data WARNING : specutil.py : 0 [295, 334] WARNING : specutil.py : 1 [448, 452] WARNING : specutil.py : 2 [505, 511] WARNING : specutil.py : 3 [585, 590] WARNING : specutil.py : 0 [90, 94] WARNING : specutil.py : 1 [160, 163] WARNING : specutil.py : 2 [220, 226] WARNING : specutil.py : 3 [258, 261] WARNING : specutil.py : 4 [299, 302] WARNING : specutil.py : 5 [309, 339] WARNING : specutil.py : 6 [407, 411] WARNING : specutil.py : 7 [420, 424] WARNING : specutil.py : 8 [463, 479] WARNING : specutil.py : 9 [484, 499] WARNING : specutil.py : 10 [611, 622] WARNING : specutil.py : 11 [635, 640] WARNING : specutil.py : 12 [677, 685] WARNING : specutil.py : 0 [295, 334] WARNING : specutil.py : 1 [448, 452] WARNING : specutil.py : 2 [505, 511] WARNING : specutil.py : 3 [585, 590] INFO : LineID_AT.py : Detecting segments in CubeStats based data WARNING : specutil.py : 0 [152, 155] WARNING : specutil.py : 1 [281, 354] WARNING : specutil.py : 2 [506, 509] WARNING : specutil.py : 3 [567, 571] WARNING : specutil.py : 4 [597, 600] WARNING : specutil.py : 0 [649, 652] INFO : LineID_AT.py : Searching for spectral peaks with method: PeakFinder INFO : LineID_AT.py : Too many peaks in CubeSpectrum for pattern finding to be useful, turning it off.[2] /an/terra_astromake/opt/casa/casa-release-5.8.0-109.el7/lib/python2.7/site-packages/scipy/optimize/minpack.py:715: OptimizeWarning: Covariance of the parameters could not be estimated category=OptimizeWarning) INFO : LineID_AT.py : Found line: CO 1-0 @ 115.2712GHz, channels 90 - 640 REGRESSION : LINEID: CO 115.27120 90 640 INFO : LineID_AT.py : Found line: U_115.3906 @ 115.390617557GHz, channels 649 - 652 REGRESSION : LINEID: U_115.3906 115.39062 649 652 INFO : LineID_AT.py : Found line: OS18O 4(1,4)-3(0,3) @ 115.40201GHz, channels 677 - 685 REGRESSION : LINEID: OS18O 115.40201 677 685 INFO : LineID_AT.py : Line Coverage 564 / 689 = 0.818578 TIMING : LineID done [ 19.464701 19.71278286 1666.3828125 407.85546875] TIMING : LineID END [ 56.006929 56.34611607] INFO : AT.py : BDP_OUT[0] = LineList_BDP x.ll TIMING : ADMITrun END [ 108.392966 115.91578197] INFO : Admit.py : ADMIT run() finished [flowcount 1] [cpu 108.393 115.916 ] Found 3 lines in LineID ('PJT0 _taskid', -1) ('PJT0 _taskid', -1) INFO : Admit.py : ADMIT run() called [flowcount 1] INFO : INFO : INFO : Executing LineCube_AT - '' (V1.2.2) INFO : INFO : INFO : Run using the following settings: INFO : equalize : False INFO : pad : 5 INFO : fpad : -1.0 INFO : TIMING : LineCube ADMIT [ 1.10253763e+02 1.64394323e+09] TIMING : LineCube BEGIN [ 0. 0.] TIMING : LineCube start [ 8.59380000e-02 8.68251324e-02 1.61694531e+03 3.58570312e+02] TIMING : LineCube pad [ 1.35170000e-02 1.35548115e-02 1.61694531e+03 3.58570312e+02] TIMING : LineCube trans-x.CO_115.27120 [ 8.64340000e-02 1.94797993e-01 1.61694531e+03 3.58691406e+02] TIMING : LineCube trans-x.U_115.3906 [ 7.05350000e-02 1.55294180e-01 1.61694531e+03 3.58691406e+02] WARNING : LineCube_AT.py : pad=5 too large, end=690 resetting to 688 TIMING : LineCube trans-x.OS18O_115.40201 [ 9.45640000e-02 3.44420910e-01 1.61694531e+03 3.58691406e+02] REGRESSION : LC: [85, 645, 644, 657, 672, 688] TIMING : LineCube done [ 2.51750000e-02 2.52399445e-02 1.61694531e+03 3.58691406e+02] TIMING : LineCube END [ 0.388563 0.83256698] INFO : AT.py : BDP_OUT[0] = LineCube_BDP x.CO_115.27120/lc.im INFO : AT.py : BDP_OUT[1] = LineCube_BDP x.U_115.3906/lc.im INFO : AT.py : BDP_OUT[2] = LineCube_BDP x.OS18O_115.40201/lc.im TIMING : ADMITrun END [ 109.013201 117.04926395] INFO : Admit.py : ADMIT run() finished [flowcount 1] [cpu 109.013 117.049 ] Found 3 lines in LineCube ('LineDir:', 0, 'x.CO_115.27120/lc.im') ('PJT0 _taskid', -1) ('PJT0 _taskid', -1) ('MOMENT_AT:', 8) INFO : AT.py : Setting 'moments' = [0, 1, 2] for Moment_AT INFO : AT.py : Setting 'mom0clip' = 2.0 for Moment_AT ('PJT0 _taskid', -1) ('PJT0 _taskid', -1) ('LineDir:', 1, 'x.U_115.3906/lc.im') ('PJT0 _taskid', -1) ('PJT0 _taskid', -1) ('MOMENT_AT:', 10) INFO : AT.py : Setting 'moments' = [0, 1, 2] for Moment_AT INFO : AT.py : Setting 'mom0clip' = 2.0 for Moment_AT ('PJT0 _taskid', -1) ('PJT0 _taskid', -1) ('LineDir:', 2, 'x.OS18O_115.40201/lc.im') ('PJT0 _taskid', -1) ('PJT0 _taskid', -1) ('MOMENT_AT:', 12) INFO : AT.py : Setting 'moments' = [0, 1, 2] for Moment_AT INFO : AT.py : Setting 'mom0clip' = 2.0 for Moment_AT ('PJT0 _taskid', -1) ('PJT0 _taskid', -1) ('PJT0 _taskid', -1) ('PJT0 _taskid', -1) INFO : AT.py : Setting 'numsigma' = [3.0] for Moment_AT INFO : AT.py : Setting 'mom0clip' = 2.0 for Moment_AT INFO : Admit.py : ADMIT run() called [flowcount 1] INFO : INFO : INFO : Executing Moment_AT - '@3' (V1.2.2) INFO : INFO : INFO : Run using the following settings: INFO : numsigma : [3.0] INFO : moments : [0] INFO : mom0clip : 2.0 INFO : zoom : 1 INFO : chans : INFO : sigma : -1.0 INFO : TIMING : Moment ADMIT [ 1.11057668e+02 1.64394323e+09] TIMING : Moment BEGIN [ 0. 0.] TIMING : Moment open [ 4.88210000e-02 4.89518642e-02 1.61694531e+03 3.58691406e+02] TIMING : Moment immoments-0 [ 1.08209000e-01 1.72924995e-01 1.63201172e+03 3.73632812e+02] TIMING : Moment mom0clip [ 6.85860000e-02 1.34312153e-01 1.63201172e+03 3.73882812e+02] TIMING : Moment implot [ 9.77720000e-02 3.47199917e-01 1.63201172e+03 3.73882812e+02] TIMING : Moment getdata [ 2.30140000e-02 2.30801105e-02 1.63201172e+03 3.73882812e+02] TIMING : Moment ren+mask_0 [ 3.46704000e-01 3.51675034e-01 1.63448047e+03 3.76351562e+02] TIMING : Moment done [ 1.50610000e-02 1.51000023e-02 1.63448047e+03 3.76351562e+02] TIMING : Moment END [ 0.72048 1.1056149] INFO : AT.py : BDP_OUT[0] = Moment_BDP x-@3.mom_0 INFO : INFO : INFO : Executing Moment_AT - '' (V1.2.2) INFO : INFO : INFO : Run using the following settings: INFO : numsigma : [2.0] INFO : moments : [0, 1, 2] INFO : mom0clip : 2.0 INFO : zoom : 1 INFO : chans : INFO : sigma : -1.0 INFO : TIMING : Moment ADMIT [ 1.11838958e+02 1.64394323e+09] TIMING : Moment BEGIN [ 0. 0.] TIMING : Moment open [ 4.84420000e-02 4.85610962e-02 1.63448047e+03 3.76351562e+02] TIMING : Moment immoments-all [ 1.12502000e-01 4.02621031e-01 1.63448047e+03 3.76351562e+02] TIMING : Moment mom0clip [ 6.90810000e-02 1.19409084e-01 1.63448047e+03 3.76351562e+02] TIMING : Moment implot [ 9.45730000e-02 3.64938974e-01 1.63448047e+03 3.76351562e+02] TIMING : Moment getdata [ 2.04390000e-02 2.11958885e-02 1.63448047e+03 3.76351562e+02] TIMING : Moment ren+mask_0 [ 2.68730000e-01 2.73782015e-01 1.63398047e+03 3.75851562e+02] 2022-02-04 02:53:51 WARN ImageRegrid::regrid You are regridding an image whose beam is not well sampled by the pixel size. Total flux can be lost when regridding such images, especially when the new pixel size is larger than the old pixel size. It is recommended to check the total flux of your input and output image, and if necessary rebin the input to have smaller pixels. TIMING : Moment makemask [ 1.49617000e-01 5.64512968e-01 1.63398047e+03 3.75914062e+02] TIMING : Moment implot [ 7.70160000e-02 3.20559025e-01 1.63398047e+03 3.75914062e+02] TIMING : Moment getdata [ 2.02170000e-02 2.03440189e-02 1.63398047e+03 3.75914062e+02] TIMING : Moment ren+mask_1 [ 2.80080000e-01 2.85331011e-01 1.63398047e+03 3.76433594e+02] 2022-02-04 02:53:52 WARN ImageRegrid::regrid You are regridding an image whose beam is not well sampled by the pixel size. Total flux can be lost when regridding such images, especially when the new pixel size is larger than the old pixel size. It is recommended to check the total flux of your input and output image, and if necessary rebin the input to have smaller pixels. TIMING : Moment makemask [ 1.51437000e-01 5.08565903e-01 1.63398047e+03 3.76433594e+02] TIMING : Moment implot [ 7.62270000e-02 3.20397139e-01 1.63398047e+03 3.76433594e+02] TIMING : Moment getdata [ 2.00430000e-02 2.01170444e-02 1.63398047e+03 3.76433594e+02] TIMING : Moment ren+mask_2 [ 2.64442000e-01 2.68919945e-01 1.63373047e+03 3.76183594e+02] REGRESSION : MOM0FLUX: x.CO_115.27120 9295.53 3953.22 0.0330578 -269.271 -269.271 22.6362 TIMING : Moment mom0flux [ 5.21690000e-02 5.23190498e-02 1.63373047e+03 3.76183594e+02] TIMING : Moment flux-spectrum [ 4.83813000e-01 4.89567995e-01 1.63348047e+03 3.75933594e+02] TIMING : Moment done [ 1.48260000e-02 1.48689747e-02 1.63348047e+03 3.75933594e+02] TIMING : Moment END [ 2.216096 4.10850501] INFO : AT.py : BDP_OUT[0] = Moment_BDP x.CO_115.27120/lc.mom_0 INFO : AT.py : BDP_OUT[1] = Moment_BDP x.CO_115.27120/lc.mom_1 INFO : AT.py : BDP_OUT[2] = Moment_BDP x.CO_115.27120/lc.mom_2 INFO : INFO : INFO : Executing Moment_AT - '@1' (V1.2.2) INFO : INFO : INFO : Run using the following settings: INFO : numsigma : [2.0] INFO : moments : [0, 1, 2] INFO : mom0clip : 2.0 INFO : zoom : 1 INFO : chans : INFO : sigma : -1.0 INFO : TIMING : Moment ADMIT [ 1.14164199e+02 1.64394323e+09] TIMING : Moment BEGIN [ 0. 0.] TIMING : Moment open [ 4.85410000e-02 4.86640930e-02 1.63348047e+03 3.75933594e+02] TIMING : Moment immoments-all [ 6.86090000e-02 1.96924925e-01 1.63348047e+03 3.75933594e+02] TIMING : Moment mom0clip [ 6.72210000e-02 1.25940084e-01 1.63348047e+03 3.75933594e+02] TIMING : Moment implot [ 9.89590000e-02 3.41333866e-01 1.63348047e+03 3.75933594e+02] TIMING : Moment getdata [ 2.06950000e-02 2.08010674e-02 1.63348047e+03 3.75933594e+02] TIMING : Moment ren+mask_0 [ 2.65552000e-01 2.69582987e-01 1.63348047e+03 3.75933594e+02] 2022-02-04 02:53:55 WARN ImageRegrid::regrid You are regridding an image whose beam is not well sampled by the pixel size. Total flux can be lost when regridding such images, especially when the new pixel size is larger than the old pixel size. It is recommended to check the total flux of your input and output image, and if necessary rebin the input to have smaller pixels. TIMING : Moment makemask [ 1.50033000e-01 5.03468037e-01 1.63348047e+03 3.75933594e+02] TIMING : Moment implot [ 7.58990000e-02 3.25937033e-01 1.63348047e+03 3.75933594e+02] TIMING : Moment getdata [ 1.99420000e-02 2.00109482e-02 1.63348047e+03 3.75933594e+02] TIMING : Moment ren+mask_1 [ 2.58314000e-01 2.63098001e-01 1.63348047e+03 3.75933594e+02] 2022-02-04 02:53:56 WARN ImageRegrid::regrid You are regridding an image whose beam is not well sampled by the pixel size. Total flux can be lost when regridding such images, especially when the new pixel size is larger than the old pixel size. It is recommended to check the total flux of your input and output image, and if necessary rebin the input to have smaller pixels. TIMING : Moment makemask [ 1.52850000e-01 5.66740036e-01 1.63348047e+03 3.75933594e+02] TIMING : Moment implot [ 7.72190000e-02 3.14831972e-01 1.63348047e+03 3.75933594e+02] TIMING : Moment getdata [ 2.02450000e-02 2.03139782e-02 1.63348047e+03 3.75933594e+02] TIMING : Moment ren+mask_2 [ 2.55877000e-01 2.60810137e-01 1.63348047e+03 3.75933594e+02] REGRESSION : MOM0FLUX: x.U_115.3906 -139.626 6.07861 0.0330578 -296.241 -296.241 1.24897 TIMING : Moment mom0flux [ 5.08340000e-02 5.13868332e-02 1.63348047e+03 3.75933594e+02] TIMING : Moment flux-spectrum [ 2.55375000e-01 2.59870052e-01 1.63348047e+03 3.75933594e+02] TIMING : Moment done [ 1.49370000e-02 1.49779320e-02 1.63348047e+03 3.75933594e+02] TIMING : Moment END [ 1.913478 3.61712003] INFO : AT.py : BDP_OUT[0] = Moment_BDP x.U_115.3906/lc-@1.mom_0 INFO : AT.py : BDP_OUT[1] = Moment_BDP x.U_115.3906/lc-@1.mom_1 INFO : AT.py : BDP_OUT[2] = Moment_BDP x.U_115.3906/lc-@1.mom_2 INFO : INFO : INFO : Executing Moment_AT - '@2' (V1.2.2) INFO : INFO : INFO : Run using the following settings: INFO : numsigma : [2.0] INFO : moments : [0, 1, 2] INFO : mom0clip : 2.0 INFO : zoom : 1 INFO : chans : INFO : sigma : -1.0 INFO : TIMING : Moment ADMIT [ 1.16186703e+02 1.64394324e+09] TIMING : Moment BEGIN [ 0. 0.] TIMING : Moment open [ 4.88170000e-02 4.89299297e-02 1.63348047e+03 3.75933594e+02] TIMING : Moment immoments-all [ 6.96070000e-02 3.68674040e-01 1.63348047e+03 3.75933594e+02] TIMING : Moment mom0clip [ 6.73620000e-02 1.42038107e-01 1.63348047e+03 3.75933594e+02] TIMING : Moment implot [ 9.72900000e-02 3.64435911e-01 1.63348047e+03 3.75933594e+02] TIMING : Moment getdata [ 2.30390000e-02 2.31199265e-02 1.63348047e+03 3.75933594e+02] TIMING : Moment ren+mask_0 [ 2.64874000e-01 2.69795179e-01 1.63348047e+03 3.75933594e+02] 2022-02-04 02:53:59 WARN ImageRegrid::regrid You are regridding an image whose beam is not well sampled by the pixel size. Total flux can be lost when regridding such images, especially when the new pixel size is larger than the old pixel size. It is recommended to check the total flux of your input and output image, and if necessary rebin the input to have smaller pixels. TIMING : Moment makemask [ 1.49500000e-01 4.42434788e-01 1.63348047e+03 3.75933594e+02] TIMING : Moment implot [ 7.55040000e-02 3.22351217e-01 1.63348047e+03 3.75933594e+02] TIMING : Moment getdata [ 1.98830000e-02 1.99508667e-02 1.63348047e+03 3.75933594e+02] TIMING : Moment ren+mask_1 [ 2.51497000e-01 2.56025076e-01 1.63348047e+03 3.75933594e+02] 2022-02-04 02:54:00 WARN ImageRegrid::regrid You are regridding an image whose beam is not well sampled by the pixel size. Total flux can be lost when regridding such images, especially when the new pixel size is larger than the old pixel size. It is recommended to check the total flux of your input and output image, and if necessary rebin the input to have smaller pixels. TIMING : Moment makemask [ 1.49894000e-01 4.23143864e-01 1.63348047e+03 3.75933594e+02] TIMING : Moment implot [ 7.55290000e-02 3.20048094e-01 1.63348047e+03 3.75933594e+02] TIMING : Moment getdata [ 1.99490000e-02 2.00231075e-02 1.63348047e+03 3.75933594e+02] TIMING : Moment ren+mask_2 [ 2.47012000e-01 2.51562834e-01 1.63323047e+03 3.75683594e+02] REGRESSION : MOM0FLUX: x.OS18O_115.40201 -417.188 0.662747 0.0330578 -294.2 -294.2 0.0691986 TIMING : Moment mom0flux [ 5.12550000e-02 5.14061451e-02 1.63323047e+03 3.75683594e+02] TIMING : Moment flux-spectrum [ 2.61646000e-01 2.65436888e-01 1.63323047e+03 3.75683594e+02] TIMING : Moment done [ 1.47030000e-02 1.47490501e-02 1.63323047e+03 3.75683594e+02] TIMING : Moment END [ 1.899709 3.61655307] INFO : AT.py : BDP_OUT[0] = Moment_BDP x.OS18O_115.40201/lc-@2.mom_0 INFO : AT.py : BDP_OUT[1] = Moment_BDP x.OS18O_115.40201/lc-@2.mom_1 INFO : AT.py : BDP_OUT[2] = Moment_BDP x.OS18O_115.40201/lc-@2.mom_2 INFO : INFO : INFO : Executing CubeSpectrum_AT - '@1' (V1.2.5) INFO : INFO : INFO : Run using the following settings: INFO : xaxis : INFO : sources : [0] INFO : pos : [] INFO : TIMING : CubeSpectrum ADMIT [ 1.18195666e+02 1.64394324e+09] TIMING : CubeSpectrum BEGIN [ 0. 0.] INFO : CubeSpectrum_AT.py : CubeSum::maxpos,val=[39, 58],14.251700 TIMING : CubeSpectrum Moment-pos [ 5.41200000e-02 5.42709827e-02 1.63323047e+03 3.75683594e+02] TIMING : CubeSpectrum open [ 1.22960000e-02 1.23331547e-02 1.63323047e+03 3.75683594e+02] TIMING : CubeSpectrum imval [ 4.12270000e-02 4.14218903e-02 1.63323047e+03 3.75683594e+02] TIMING : CubeSpectrum imhead [ 4.34110000e-02 4.35700417e-02 1.63323047e+03 3.75683594e+02] REGRESSION : CSP: [0.5920255184173584] INFO : CubeSpectrum_AT.py : Writing 1 testCubeSpectrum tables TIMING : CubeSpectrum done [ 2.62542000e-01 2.67728090e-01 1.63323047e+03 3.75683594e+02] TIMING : CubeSpectrum summary [ 1.25090000e-02 1.25408173e-02 1.63323047e+03 3.75683594e+02] TIMING : CubeSpectrum END [ 0.438435 0.44425297] INFO : AT.py : BDP_OUT[0] = CubeSpectrum_BDP x.CO_115.27120/lc-@1.csp INFO : INFO : INFO : Executing CubeSpectrum_AT - '@2' (V1.2.5) INFO : INFO : INFO : Run using the following settings: INFO : xaxis : INFO : sources : [0] INFO : pos : [] INFO : TIMING : CubeSpectrum ADMIT [ 1.18694559e+02 1.64394324e+09] TIMING : CubeSpectrum BEGIN [ 0. 0.] INFO : CubeSpectrum_AT.py : CubeSum::maxpos,val=[23, 21],0.381633 TIMING : CubeSpectrum Moment-pos [ 5.39680000e-02 5.41110039e-02 1.63323047e+03 3.75683594e+02] TIMING : CubeSpectrum open [ 1.22510000e-02 1.23019218e-02 1.63323047e+03 3.75683594e+02] TIMING : CubeSpectrum imval [ 3.61100000e-02 3.62949371e-02 1.63323047e+03 3.75683594e+02] TIMING : CubeSpectrum imhead [ 4.21940000e-02 4.23550606e-02 1.63323047e+03 3.75683594e+02] REGRESSION : CSP: [0.051480114459991455] INFO : CubeSpectrum_AT.py : Writing 1 testCubeSpectrum tables TIMING : CubeSpectrum done [ 2.61370000e-01 2.65883923e-01 1.63323047e+03 3.75683594e+02] TIMING : CubeSpectrum summary [ 1.27420000e-02 1.27780437e-02 1.63323047e+03 3.75683594e+02] TIMING : CubeSpectrum END [ 0.430986 0.43611097] INFO : AT.py : BDP_OUT[0] = CubeSpectrum_BDP x.U_115.3906/lc-@2.csp INFO : INFO : INFO : Executing CubeSpectrum_AT - '@3' (V1.2.5) INFO : INFO : INFO : Run using the following settings: INFO : xaxis : INFO : sources : [0] INFO : pos : [] INFO : TIMING : CubeSpectrum ADMIT [ 1.19186508e+02 1.64394324e+09] TIMING : CubeSpectrum BEGIN [ 0. 0.] INFO : CubeSpectrum_AT.py : CubeSum::maxpos,val=[57, 51],0.331697 TIMING : CubeSpectrum Moment-pos [ 5.47540000e-02 5.49008846e-02 1.63323047e+03 3.75683594e+02] TIMING : CubeSpectrum open [ 1.24040000e-02 1.24409199e-02 1.63323047e+03 3.75683594e+02] TIMING : CubeSpectrum imval [ 3.66090000e-02 3.67891788e-02 1.63323047e+03 3.75683594e+02] TIMING : CubeSpectrum imhead [ 4.23450000e-02 4.25000191e-02 1.63323047e+03 3.75683594e+02] REGRESSION : CSP: [0.052873451262712479] INFO : CubeSpectrum_AT.py : Writing 1 testCubeSpectrum tables TIMING : CubeSpectrum done [ 2.70212000e-01 2.74788857e-01 1.63323047e+03 3.75683594e+02] TIMING : CubeSpectrum summary [ 1.24630000e-02 1.24950409e-02 1.63323047e+03 3.75683594e+02] TIMING : CubeSpectrum END [ 0.441055 0.44622087] INFO : AT.py : BDP_OUT[0] = CubeSpectrum_BDP x.OS18O_115.40201/lc-@3.csp TIMING : ADMITrun END [ 117.97832 131.8125298] INFO : Admit.py : ADMIT run() finished [flowcount 1] [cpu 117.978 131.813 ]