========================================= The start-up time of CASA may vary depending on whether the shared libraries are cached or not. ========================================= ]0;IPython: 2022VLBI1mmCommissioning/93231IPython 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 ', 'W-Hya_93231_0.nf.admit') ('Skipping ', 'admit1.apar') ('Found parameter file to execfile', 'W-Hya_93231_0.nf.fits.apar') ('Skipping ', '') ('_ADMIT_ROOT=', '/lma1/teuben/lmtoy/admit') ADMIT basedir = /lma1/teuben/lmtoy/examples/2022VLBI1mmCommissioning/93231/W-Hya_93231_0.nf.admit ADMIT root = /lma1/teuben/lmtoy/admit ADMIT version = 1.0.8.5 ADMIT flowcount = 1 stale = 0 ('Starting a new ADMIT using', 'W-Hya_93231_0.nf.fits') ('Found parameter file to cp:', 'W-Hya_93231_0.nf.fits.apar') INFO : AT.py : Setting {'basename': 'x', 'file': 'W-Hya_93231_0.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 : W-Hya_93231_0.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.23999600e+00 1.64394099e+09] TIMING : Ingest BEGIN [ 0. 0.] INFO : Ingest_AT.py : OBJECT: W-Hya SHAPE: [ 81 81 1379] ALMA OBJECT = W-Hya ALMA DATE-OBS = 2022-01-24T07:23:34 INFO : Ingest_AT.py : basename=x INFO : utils.py : OSTYPE: linux TIMING : Ingest start [ 2.33618000e-01 2.40031004e-01 1.41113672e+03 2.17289062e+02] TIMING : Ingest importfits [ 1.65620000e-01 2.29859114e-01 1.41113672e+03 2.18714844e+02] WARNING : Ingest_AT.py : Adding dummy STOKES-I axis TIMING : Ingest adddegaxes [ 7.12360000e-02 1.00524902e-01 1.41113672e+03 2.19554688e+02] TIMING : Ingest summary-0 [ 1.42640000e-02 1.43530369e-02 1.41113672e+03 2.19554688e+02] TIMING : Ingest mask [ 5.81530000e-02 7.27150440e-02 1.41113672e+03 2.19554688e+02] TIMING : Ingest summary-1 [ 1.35440000e-02 3.23739052e-02 1.41113672e+03 2.19554688e+02] TIMING : Ingest statistics [ 4.44200000e-02 4.45520878e-02 1.41113672e+03 2.20578125e+02] 2022-02-04 02:16:34 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': 6.9000012, 'unit': 'arcsec'}, 'pa': {'value': 0.0, 'unit': 'deg'}, 'minor': {'value': 6.9000012, 'unit': 'arcsec'}} INFO : Ingest_AT.py : BASICS: [shape] npts min max: [ 81 81 1379 1] 1071483 -1.489022 2.783515 INFO : Ingest_AT.py : S/N (all data): 17.002727 INFO : Ingest_AT.py : GOOD PIXELS: 1071483/9047619 (11.842707% good or 88.157293% bad) WARNING : Ingest_AT.py : MASKS: ['mask0'] REGRESSION : CUBE: -1.48902 2.78352 0.16371 81 81 1379 88.157293 INFO : Ingest_AT.py : TELESCOPE: LMT INFO : Ingest_AT.py : OBJECT: W-Hya INFO : Ingest_AT.py : REFFREQTYPE: LSRK INFO : Ingest_AT.py : RA Axis 1: 0.000000 -3.000000 40.000000 INFO : Ingest_AT.py : DEC Axis 2: 0.000000 3.000000 40.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: 230.776 -0.000390625 0 INFO : Ingest_AT.py : Cube Orig Axis 3: type=Frequency velocity increment=0.507970 km/s @ fc=230.506667 fw=-0.538672 GHz INFO : Ingest_AT.py : RESTFREQ: 230.538 230.538 -1 INFO : Ingest_AT.py : VLSRc= 40.746068 VLSRf= 0.000000 VLSRv= 0.000000 VLSRz= 0.000000 WIDTH= 700.490906 INFO : Ingest_AT.py : VLSR = 40.746068 errs = 0.000000 0.000000 0.000000 width = 0.507970 TIMING : Ingest done [ 5.68291000e-01 5.94938040e-01 1.41512500e+03 2.24914062e+02] TIMING : Ingest END [ 1.18129 1.34152293] 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.48241800e+00 1.64394099e+09] TIMING : CubeStats BEGIN [ 0. 0.] TIMING : CubeStats imval [ 5.75800000e-02 5.78219891e-02 1.41387500e+03 2.24191406e+02] TIMING : CubeStats start [ 1.26910000e-02 1.27320290e-02 1.41387500e+03 2.24191406e+02] TIMING : CubeStats imstat0 [ 1.43430000e-01 1.43865108e-01 1.41373828e+03 2.24636719e+02] TIMING : CubeStats imstat1 [ 2.94117000e-01 2.94954062e-01 1.45823047e+03 2.69312500e+02] INFO : CubeStats_AT.py : sigma varies from 0.071626 to 0.255730; 1379/1379 channels ok WARNING : CubeStats_AT.py : sigma varies too much, going to clip to 0.143252 (3.57035 > 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 [ 1.585759 1.59159493 1458.23046875 269.3125 ] INFO : CubeStats_AT.py : CubeMax: 2.783515 @ [ 41 49 688 0] INFO : CubeStats_AT.py : CubeMin: -1.489022 @ [38 72 4 0] INFO : CubeStats_AT.py : CubeRMS: 0.154706 INFO : CubeStats_AT.py : RMS Sanity check 1.058200 REGRESSION : CST: 0.154706 1.058200 INFO : CubeStats_AT.py : mean,rms,S/N=0.001940 0.154706 17.992281 INFO : CubeStats_AT.py : RMS BAD VARIATION RATIO: 0.820002 0.529815 /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.328357 2.49864507 1467.984375 277.30078125] TIMING : CubeStats done [ 1.35260000e-02 1.35648251e-02 1.46798438e+03 2.77300781e+02] TIMING : CubeStats summary [ 1.22210000e-02 1.22680664e-02 1.46798438e+03 2.77300781e+02] TIMING : CubeStats END [ 4.459853 4.63765812] 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 [ 8.00216600e+00 1.64394100e+09] TIMING : CubeSum BEGIN [ 0. 0.] INFO : CubeSum_AT.py : Using constant sigma = 0.154706 TIMING : CubeSum start [ 5.23270000e-02 5.24990559e-02 1.46798438e+03 2.77300781e+02] TIMING : CubeSum immoments [ 1.20252000e-01 2.83983946e-01 1.46798438e+03 2.77746094e+02] 2022-02-04 02:16:40 WARN immoments::image::open (file ../../tools/images/image_cmpt.cc, line 4125) Another image is already open, closing first TIMING : CubeSum statistics [ 2.01620000e-02 2.02438831e-02 1.46798438e+03 2.77746094e+02] INFO : CubeSum_AT.py : Total flux: -107.373363 (sum=-11.930372) REGRESSION : CSM: [-107.37336293943326, -11.930371751077473] TIMING : CubeSum implot [ 1.69893000e-01 6.34973216e+00 1.53998828e+03 2.77785156e+02] TIMING : CubeSum getdata [ 3.16770000e-02 3.17778587e-02 1.53998828e+03 2.77785156e+02] TIMING : CubeSum done [ 2.88714000e-01 3.08145046e-01 1.53998828e+03 2.77851562e+02] TIMING : CubeSum END [ 0.697129 7.06054497] 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 [ 8.75973800e+00 1.64394101e+09] TIMING : SFind2D BEGIN [ 0. 0.] ('BDP', 0, ) ('BDP', 1, ) ('BDP', 2, ) No PB given TIMING : SFind2D start [ 3.41020000e-02 3.42772007e-02 1.53998828e+03 2.77851562e+02] ('ROBUST:', {'fence': 1.5, 'algorithm': 'hinges-fences'}) TIMING : SFind2D imstat [ 4.12690000e-02 4.15008068e-02 1.53998828e+03 2.77851562e+02] INFO : AT.py : Setting 'sigma' = 0.65903369134 for SFind2D_AT INFO : SFind2D_AT.py : sigma, dmin, dmax, snmax, cutoff 0.659034 -39.7263 36.2996 35 0.0726216 WARNING : SFind2D_AT.py : Temporarely patching your K.km/s units to Jy/beam for ia.findsources() 2022-02-04 02:16:47 WARN imstat::ImageRegion::fromRecord No sources were found TIMING : SFind2D findsources [ 9.09360000e-02 1.22226000e-01 1.53998828e+03 2.78027344e+02] REGRESSION : CONTFLUX: 0 0 INFO : SFind2D_AT.py : Fitted Gaussian size; NOT deconvolved source size. INFO : SFind2D_AT.py : Restoring Beam: Major axis: 6.9 arcsec , Minor axis: 6.9 arcsec , PA: 0.0 deg TIMING : SFind2D done [ 2.71806000e-01 2.95084000e-01 1.53998828e+03 2.78085938e+02] TIMING : SFind2D END [ 0.451901 0.50690913] INFO : AT.py : BDP_OUT[0] = SourceList_BDP x-csm.sl TIMING : ADMITrun END [ 7.41932 14.23855996] INFO : Admit.py : ADMIT run() finished [flowcount 1] [cpu 7.41932 14.2386 ] ('N sources in CSM:', 0) ('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 [ 9.46434200e+00 1.64394101e+09] TIMING : CubeSpectrum BEGIN [ 0. 0.] INFO : CubeSpectrum_AT.py : CubeStats::maxpos,val=[41, 49, 688],2.783515 TIMING : CubeSpectrum CubeStats-pos [ 5.82130000e-02 5.83589077e-02 1.53998828e+03 2.78085938e+02] INFO : CubeSpectrum_AT.py : CubeSum::maxpos,val=[44, 64],36.299595 TIMING : CubeSpectrum Moment-pos [ 4.43270000e-02 4.44600582e-02 1.53998828e+03 2.78085938e+02] TIMING : CubeSpectrum open [ 1.21790000e-02 1.22148991e-02 1.53998828e+03 2.78085938e+02] TIMING : CubeSpectrum imval [ 5.01180000e-02 5.03771305e-02 1.53998828e+03 2.78742188e+02] TIMING : CubeSpectrum imhead [ 4.19880000e-02 4.21559811e-02 1.53998828e+03 2.78742188e+02] TIMING : CubeSpectrum imval [ 2.69610000e-01 2.88204908e-01 1.53998828e+03 2.78843750e+02] REGRESSION : CSP: [2.783515453338623, 1.4877701997756958] INFO : CubeSpectrum_AT.py : Writing 2 testCubeSpectrum tables TIMING : CubeSpectrum done [ 2.90558000e-01 3.09509039e-01 1.53998828e+03 2.78843750e+02] TIMING : CubeSpectrum summary [ 1.24120000e-02 1.24459267e-02 1.53998828e+03 2.78843750e+02] TIMING : CubeSpectrum END [ 0.791582 0.82993579] 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 [ 1.03163910e+01 1.64394101e+09] TIMING : LineSegment BEGIN [ 0. 0.] TIMING : LineSegment start [ 2.44140000e-02 2.44801044e-02 1.53998828e+03 2.78843750e+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 [ 13.002557 13.03654981 1539.98828125 279.19140625] INFO : LineSegment_AT.py : Detecting segments in CubeSpectrum based data WARNING : specutil.py : 0 [679, 699] WARNING : specutil.py : 0 [0, 60] INFO : LineSegment_AT.py : Detecting segments in CubeStats based data WARNING : specutil.py : 0 [0, 33] WARNING : specutil.py : 1 [40, 43] WARNING : specutil.py : 2 [680, 697] WARNING : specutil.py : 0 [0, 28] TIMING : LineSegment segment finder [ 3.228434 3.23689318 1539.98828125 279.19140625] INFO : LineSegment_AT.py : Segment Coverage 82 / 1379 = 0.0594634 REGRESSION : LINESEG: [[230.7629166111586, 230.77580723671574], [230.75901036098978, 230.76018223604041], [230.50354159994839, 230.51018222523538], [230.76486973624301, 230.77580723671574], [230.50276034991461, 230.51057285025229], [230.75236973570276, 230.77580723671574]] TIMING : LineSegment done [ 2.170756 2.2811749 1568.06640625 307.19140625] TIMING : LineSegment END [ 18.438568 18.591537] INFO : AT.py : BDP_OUT[0] = LineSegment_BDP x.lseg TIMING : ADMITrun END [ 26.950732 34.01343203] INFO : Admit.py : ADMIT run() finished [flowcount 1] [cpu 26.9507 34.0134 ] Found 2 segments in LineSegment No ContinuumSub needed INFO : Admit.py : ADMIT run() called [flowcount 1] TIMING : ADMITrun END [ 27.050761 34.16893196] INFO : Admit.py : ADMIT run() finished [flowcount 1] [cpu 27.0508 34.1689 ] ('OBJECT = ', 'W-Hya') ('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 [ 2.91204810e+01 1.64394103e+09] TIMING : LineID BEGIN [ 0. 0.] INFO : LineID_AT.py : Set vlsr = 40.75 for line identification. INFO : LineID_AT.py : Identifylines = True INFO : LineID_AT.py : Using vlsr = 40.7461 INFO : LineID_AT.py : Attempting Continuum Subtraction for Input Spectra TIMING : LineID getspectrum-cubespecs [ 5.847185 5.86190891 1568.06640625 307.19921875] INFO : LineID_AT.py : Attempting Continuum Subtraction for Input CubeStats Spectra TIMING : LineID getspectrum-cubestats [ 7.307291 7.32567191 1568.06640625 307.19921875] TIMING : LineID segment finder [ 1.50940000e-02 1.51340961e-02 1.56806641e+03 3.07199219e+02] INFO : LineID_AT.py : Detecting segments in CubeSpectrum based data WARNING : specutil.py : 0 [679, 699] WARNING : specutil.py : 0 [0, 60] INFO : LineID_AT.py : Detecting segments in CubeStats based data WARNING : specutil.py : 0 [0, 33] WARNING : specutil.py : 1 [40, 43] WARNING : specutil.py : 2 [680, 697] WARNING : specutil.py : 0 [0, 28] INFO : LineID_AT.py : Searching for spectral peaks with method: PeakFinder /an/terra_astromake/opt/casa/casa-release-5.8.0-109.el7/lib/python2.7/site-packages/matplotlib/text.py:1206: RuntimeWarning: divide by zero encountered in double_scalars cwd *= 1+dashpad/np.sqrt(np.dot(cwd,cwd)) /an/terra_astromake/opt/casa/casa-release-5.8.0-109.el7/lib/python2.7/site-packages/matplotlib/text.py:1206: RuntimeWarning: invalid value encountered in multiply cwd *= 1+dashpad/np.sqrt(np.dot(cwd,cwd)) INFO : LineID_AT.py : Found line: CO 2-1 @ 230.538GHz, channels 679 - 699 REGRESSION : LINEID: CO 230.53800 679 699 INFO : LineID_AT.py : Found line: H2NCH2CN 16(5,11)-16(4,12) @ 230.79742GHz, channels 0 - 60 REGRESSION : LINEID: H2NCH2CN 230.79742 0 60 INFO : LineID_AT.py : Line Coverage 82 / 1379 = 0.0594634 TIMING : LineID done [ 15.938676 16.17301202 1668.4140625 408.70703125] TIMING : LineID END [ 29.123502 29.39101887] INFO : AT.py : BDP_OUT[0] = LineList_BDP x.ll TIMING : ADMITrun END [ 56.602254 64.05316901] INFO : Admit.py : ADMIT run() finished [flowcount 1] [cpu 56.6023 64.0532 ] Found 2 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 [ 5.86275620e+01 1.64394106e+09] TIMING : LineCube BEGIN [ 0. 0.] TIMING : LineCube start [ 6.42260000e-02 6.47649765e-02 1.62313281e+03 3.63503906e+02] TIMING : LineCube pad [ 1.34180000e-02 1.34589672e-02 1.62313281e+03 3.63503906e+02] TIMING : LineCube trans-x.CO_230.53800 [ 7.46890000e-02 1.86520100e-01 1.62313281e+03 3.63566406e+02] WARNING : LineCube_AT.py : pad=5 too large, start=-5 resetting to 0 TIMING : LineCube trans-x.H2NCH2CN_230.79742 [ 9.85230000e-02 4.55888033e-01 1.62313281e+03 3.63566406e+02] REGRESSION : LC: [674, 704, 0, 65] TIMING : LineCube done [ 2.48520000e-02 2.49149799e-02 1.62313281e+03 3.63566406e+02] TIMING : LineCube END [ 0.288026 0.7579 ] INFO : AT.py : BDP_OUT[0] = LineCube_BDP x.CO_230.53800/lc.im INFO : AT.py : BDP_OUT[1] = LineCube_BDP x.H2NCH2CN_230.79742/lc.im TIMING : ADMITrun END [ 57.095015 65.09293604] INFO : Admit.py : ADMIT run() finished [flowcount 1] [cpu 57.095 65.0929 ] Found 2 lines in LineCube ('LineDir:', 0, 'x.CO_230.53800/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.H2NCH2CN_230.79742/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) ('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 - '@2' (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 [ 5.92640720e+01 1.64394106e+09] TIMING : Moment BEGIN [ 0. 0.] TIMING : Moment open [ 4.94280000e-02 4.95457649e-02 1.62313281e+03 3.63566406e+02] TIMING : Moment immoments-0 [ 1.29042000e-01 1.88088179e-01 1.63389062e+03 3.74324219e+02] TIMING : Moment mom0clip [ 6.95310000e-02 1.24019861e-01 1.63389062e+03 3.74476562e+02] TIMING : Moment implot [ 9.90300000e-02 3.38202000e-01 1.63389062e+03 3.74476562e+02] TIMING : Moment getdata [ 2.24810000e-02 2.25601196e-02 1.63389062e+03 3.74476562e+02] TIMING : Moment ren+mask_0 [ 5.00698000e-01 5.09245872e-01 1.63601172e+03 3.76597656e+02] TIMING : Moment done [ 1.64550000e-02 1.64990425e-02 1.63601172e+03 3.76597656e+02] TIMING : Moment END [ 0.899154 1.2607038] INFO : AT.py : BDP_OUT[0] = Moment_BDP x-@2.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 [ 6.02241580e+01 1.64394106e+09] TIMING : Moment BEGIN [ 0. 0.] TIMING : Moment open [ 4.85490000e-02 4.86650467e-02 1.63601172e+03 3.76597656e+02] TIMING : Moment immoments-all [ 7.08470000e-02 2.52609968e-01 1.63601172e+03 3.76597656e+02] TIMING : Moment mom0clip [ 6.96020000e-02 1.51921988e-01 1.63601172e+03 3.76597656e+02] TIMING : Moment implot [ 9.87780000e-02 3.47354889e-01 1.63601172e+03 3.76597656e+02] TIMING : Moment getdata [ 2.81490000e-02 2.82571316e-02 1.63601172e+03 3.76597656e+02] TIMING : Moment ren+mask_0 [ 3.44935000e-01 3.49473953e-01 1.63601172e+03 3.76597656e+02] 2022-02-04 02:17:41 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.52765000e-01 7.34120131e-01 1.63601172e+03 3.77238281e+02] TIMING : Moment implot [ 8.25010000e-02 3.20882797e-01 1.63601172e+03 3.77328125e+02] TIMING : Moment getdata [ 2.18460000e-02 2.19221115e-02 1.63601172e+03 3.77328125e+02] TIMING : Moment ren+mask_1 [ 3.02870000e-01 3.08063984e-01 1.63601172e+03 3.77328125e+02] 2022-02-04 02:17:42 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.57172000e-01 4.96177912e-01 1.63601172e+03 3.77328125e+02] TIMING : Moment implot [ 8.20550000e-02 3.08596134e-01 1.63601172e+03 3.77328125e+02] TIMING : Moment getdata [ 2.18230000e-02 2.19018459e-02 1.63601172e+03 3.77328125e+02] TIMING : Moment ren+mask_2 [ 3.04248000e-01 3.09207201e-01 1.63601172e+03 3.77328125e+02] REGRESSION : MOM0FLUX: x.CO_230.53800 730.995 343.548 0.111111 39.9185 39.9185 0.66235 TIMING : Moment mom0flux [ 5.23950000e-02 5.30378819e-02 1.63601172e+03 3.77328125e+02] TIMING : Moment flux-spectrum [ 3.31206000e-01 3.36648941e-01 1.63601172e+03 3.77328125e+02] TIMING : Moment done [ 1.65660000e-02 1.66170597e-02 1.63601172e+03 3.77328125e+02] TIMING : Moment END [ 2.198712 4.11793399] INFO : AT.py : BDP_OUT[0] = Moment_BDP x.CO_230.53800/lc.mom_0 INFO : AT.py : BDP_OUT[1] = Moment_BDP x.CO_230.53800/lc.mom_1 INFO : AT.py : BDP_OUT[2] = Moment_BDP x.CO_230.53800/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 [ 6.25318280e+01 1.64394106e+09] TIMING : Moment BEGIN [ 0. 0.] TIMING : Moment open [ 4.85360000e-02 4.86481190e-02 1.63601172e+03 3.77328125e+02] TIMING : Moment immoments-all [ 7.26640000e-02 2.04943895e-01 1.63601172e+03 3.77328125e+02] TIMING : Moment mom0clip [ 6.96890000e-02 1.24480009e-01 1.63601172e+03 3.77328125e+02] TIMING : Moment implot [ 9.65070000e-02 3.29553127e-01 1.63601172e+03 3.77328125e+02] TIMING : Moment getdata [ 2.18560000e-02 2.19318867e-02 1.63601172e+03 3.77328125e+02] TIMING : Moment ren+mask_0 [ 2.90303000e-01 2.95135975e-01 1.63601172e+03 3.77328125e+02] 2022-02-04 02:17:45 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.55318000e-01 5.00236034e-01 1.63601172e+03 3.77328125e+02] TIMING : Moment implot [ 8.50170000e-02 3.10328007e-01 1.63601172e+03 3.77328125e+02] TIMING : Moment getdata [ 2.14910000e-02 2.15661526e-02 1.63601172e+03 3.77328125e+02] TIMING : Moment ren+mask_1 [ 3.12408000e-01 3.16857815e-01 1.63601172e+03 3.77328125e+02] 2022-02-04 02:17:46 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.52985000e-01 4.80935097e-01 1.63601172e+03 3.77328125e+02] TIMING : Moment implot [ 8.11940000e-02 3.10870886e-01 1.63601172e+03 3.77328125e+02] TIMING : Moment getdata [ 2.18950000e-02 2.19781399e-02 1.63601172e+03 3.77328125e+02] TIMING : Moment ren+mask_2 [ 3.17532000e-01 3.22016001e-01 1.63601172e+03 3.77328125e+02] REGRESSION : MOM0FLUX: x.H2NCH2CN_230.79742 -109.28 141.981 0.111111 42.9213 42.9213 0.245618 TIMING : Moment mom0flux [ 5.25070000e-02 5.26640415e-02 1.63601172e+03 3.77328125e+02] TIMING : Moment flux-spectrum [ 3.34582000e-01 3.39637995e-01 1.63601172e+03 3.77328125e+02] TIMING : Moment done [ 1.59480000e-02 1.59959793e-02 1.63601172e+03 3.77328125e+02] TIMING : Moment END [ 2.162815 3.73021317] INFO : AT.py : BDP_OUT[0] = Moment_BDP x.H2NCH2CN_230.79742/lc-@1.mom_0 INFO : AT.py : BDP_OUT[1] = Moment_BDP x.H2NCH2CN_230.79742/lc-@1.mom_1 INFO : AT.py : BDP_OUT[2] = Moment_BDP x.H2NCH2CN_230.79742/lc-@1.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 [ 6.48233870e+01 1.64394107e+09] TIMING : CubeSpectrum BEGIN [ 0. 0.] INFO : CubeSpectrum_AT.py : CubeSum::maxpos,val=[41, 49],22.436716 TIMING : CubeSpectrum Moment-pos [ 6.15990000e-02 6.17887974e-02 1.63601172e+03 3.77328125e+02] TIMING : CubeSpectrum open [ 1.26150000e-02 1.26500130e-02 1.63601172e+03 3.77328125e+02] TIMING : CubeSpectrum imval [ 3.88090000e-02 3.90181541e-02 1.63601172e+03 3.77328125e+02] TIMING : CubeSpectrum imhead [ 4.27630000e-02 4.29389477e-02 1.63601172e+03 3.77328125e+02] REGRESSION : CSP: [2.783515453338623] INFO : CubeSpectrum_AT.py : Writing 1 testCubeSpectrum tables TIMING : CubeSpectrum done [ 3.03197000e-01 3.08434010e-01 1.63601172e+03 3.77328125e+02] TIMING : CubeSpectrum summary [ 1.25110000e-02 1.25470161e-02 1.63601172e+03 3.77328125e+02] TIMING : CubeSpectrum END [ 0.483727 0.48964286] INFO : AT.py : BDP_OUT[0] = CubeSpectrum_BDP x.CO_230.53800/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 [ 6.53675950e+01 1.64394107e+09] TIMING : CubeSpectrum BEGIN [ 0. 0.] INFO : CubeSpectrum_AT.py : CubeSum::maxpos,val=[44, 64],34.591694 TIMING : CubeSpectrum Moment-pos [ 5.52990000e-02 5.54630756e-02 1.63601172e+03 3.77328125e+02] TIMING : CubeSpectrum open [ 1.23760000e-02 1.24118328e-02 1.63601172e+03 3.77328125e+02] TIMING : CubeSpectrum imval [ 3.67110000e-02 3.69000435e-02 1.63601172e+03 3.77328125e+02] TIMING : CubeSpectrum imhead [ 4.22920000e-02 4.24549580e-02 1.63601172e+03 3.77328125e+02] REGRESSION : CSP: [1.4877701997756958] INFO : CubeSpectrum_AT.py : Writing 1 testCubeSpectrum tables TIMING : CubeSpectrum done [ 2.59209000e-01 2.63945103e-01 1.63601172e+03 3.77328125e+02] TIMING : CubeSpectrum summary [ 1.25990000e-02 1.26299858e-02 1.63601172e+03 3.77328125e+02] TIMING : CubeSpectrum END [ 0.430892 0.43624902] INFO : AT.py : BDP_OUT[0] = CubeSpectrum_BDP x.H2NCH2CN_230.79742/lc-@2.csp TIMING : ADMITrun END [ 63.9705 75.90421104] INFO : Admit.py : ADMIT run() finished [flowcount 1] [cpu 63.9705 75.9042 ]