========================================= The start-up time of CASA may vary depending on whether the shared libraries are cached or not. ========================================= ]0;IPython: 2022S1SEQUOIACommissioning/94052IPython 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-20220308-184842.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 ', 'R-Cas_94052.nf.admit') ('Skipping ', 'admit1.apar') ('Found parameter file to execfile', 'R-Cas_94052.nf.fits.apar') ('Skipping ', '') ('_ADMIT_ROOT=', '/lma1/teuben/lmtoy/admit') ADMIT basedir = /lma1/teuben/lmtoy/examples/2022S1SEQUOIACommissioning/94052/R-Cas_94052.nf.admit ADMIT root = /lma1/teuben/lmtoy/admit ADMIT version = 1.0.8.5 ADMIT flowcount = 1 stale = 0 ('Starting a new ADMIT using', 'R-Cas_94052.nf.fits') ('Found parameter file to cp:', 'R-Cas_94052.nf.fits.apar') INFO : AT.py : Setting {'basename': 'x', 'file': 'R-Cas_94052.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 : R-Cas_94052.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.29381100e+00 1.64728993e+09] TIMING : Ingest BEGIN [ 0. 0.] INFO : Ingest_AT.py : OBJECT: R-Cas SHAPE: [ 97 97 2062] ALMA OBJECT = R-Cas ALMA DATE-OBS = 2022-02-22T16:01:05 INFO : Ingest_AT.py : basename=x INFO : utils.py : OSTYPE: linux TIMING : Ingest start [ 3.50806000e-01 3.57424974e-01 1.41041016e+03 2.11894531e+02] TIMING : Ingest importfits [ 2.56568000e-01 3.31405163e-01 1.41111328e+03 2.13718750e+02] WARNING : Ingest_AT.py : Adding dummy STOKES-I axis TIMING : Ingest adddegaxes [ 1.27527000e-01 1.51062965e-01 1.42086328e+03 2.24218750e+02] TIMING : Ingest summary-0 [ 1.42140000e-02 1.43418312e-02 1.42086328e+03 2.24218750e+02] TIMING : Ingest mask [ 1.26568000e-01 1.37189150e-01 1.42173828e+03 2.25714844e+02] TIMING : Ingest summary-1 [ 1.46440000e-02 1.89168453e-02 1.42173828e+03 2.25714844e+02] TIMING : Ingest statistics [ 1.92037000e-01 1.92487001e-01 1.42396094e+03 2.28300781e+02] 2022-03-14 20:32:11 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: [ 97 97 2062 1] 11487402 -1.135214 68.046791 INFO : Ingest_AT.py : S/N (all data): 420.188317 INFO : Ingest_AT.py : GOOD PIXELS: 11487402/19401358 (59.209268% good or 40.790732% bad) WARNING : Ingest_AT.py : MASKS: ['mask0'] REGRESSION : CUBE: -1.13521 68.0468 0.161944 97 97 2062 40.790732 INFO : Ingest_AT.py : TELESCOPE: LMT INFO : Ingest_AT.py : OBJECT: R-Cas INFO : Ingest_AT.py : REFFREQTYPE: LSRK INFO : Ingest_AT.py : RA Axis 1: 359.603600 -7.499999 48.000000 INFO : Ingest_AT.py : DEC Axis 2: 51.388800 7.499999 48.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.3372 -9.76562e-05 0 INFO : Ingest_AT.py : Cube Orig Axis 3: type=Frequency velocity increment=0.339465 km/s @ fc=86.236577 fw=-0.201367 GHz INFO : Ingest_AT.py : RESTFREQ: 86.2434 86.2434 -1 INFO : Ingest_AT.py : VLSRc= 23.864173 VLSRf= 0.000000 VLSRv= 0.000000 VLSRz= 0.000000 WIDTH= 699.976211 INFO : Ingest_AT.py : VLSR = 23.864173 errs = 0.000000 0.000000 0.000000 width = 0.339465 TIMING : Ingest done [ 7.05157000e-01 7.11146116e-01 1.42696094e+03 2.31515625e+02] TIMING : Ingest END [ 1.799516 1.92600393] 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 [ 4.15377900e+00 1.64728993e+09] TIMING : CubeStats BEGIN [ 0. 0.] TIMING : CubeStats imval [ 6.19720000e-02 6.22251034e-02 1.42571094e+03 2.30601562e+02] TIMING : CubeStats start [ 1.25590000e-02 1.26099586e-02 1.42571094e+03 2.30601562e+02] TIMING : CubeStats imstat0 [ 7.98159000e-01 8.00117970e-01 1.45086328e+03 2.56207031e+02] TIMING : CubeStats imstat1 [ 8.66749000e-01 8.68767977e-01 1.46425391e+03 2.69812500e+02] INFO : CubeStats_AT.py : sigma varies from 0.114606 to 0.145667; 2062/2062 channels ok INFO : CubeStats_AT.py : Computing MaxPos for PeakPointPlot TIMING : CubeStats ppp [ 2.412759 2.42111111 1464.25390625 269.8125 ] INFO : CubeStats_AT.py : CubeMax: 68.046791 @ [ 47 48 1035 0] INFO : CubeStats_AT.py : CubeMin: -1.135214 @ [ 18 86 138 0] INFO : CubeStats_AT.py : CubeRMS: 0.127923 INFO : CubeStats_AT.py : RMS Sanity check 1.265944 REGRESSION : CST: 0.127923 1.265944 INFO : CubeStats_AT.py : mean,rms,S/N=0.000678 0.127923 531.934677 INFO : CubeStats_AT.py : RMS BAD VARIATION RATIO: 2.759680 2.834592 /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.306641 2.48172593 1474.0078125 277.66796875] TIMING : CubeStats done [ 1.38100000e-02 1.38480663e-02 1.47400781e+03 2.77667969e+02] TIMING : CubeStats summary [ 1.24020000e-02 1.24619007e-02 1.47400781e+03 2.77667969e+02] TIMING : CubeStats END [ 6.49731 6.68516493] 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 [ 1.07108170e+01 1.64728994e+09] TIMING : CubeSum BEGIN [ 0. 0.] INFO : CubeSum_AT.py : Using constant sigma = 0.127923 TIMING : CubeSum start [ 5.25370000e-02 5.26828766e-02 1.47400781e+03 2.77667969e+02] TIMING : CubeSum immoments [ 2.14527000e-01 2.68045187e-01 1.47400781e+03 2.78046875e+02] 2022-03-14 20:32:19 WARN immoments::image::open (file ../../tools/images/image_cmpt.cc, line 4125) Another image is already open, closing first TIMING : CubeSum statistics [ 2.00100000e-02 2.01218128e-02 1.47400781e+03 2.78046875e+02] INFO : CubeSum_AT.py : Total flux: 135162.912382 (sum=2402.896989) REGRESSION : CSM: [135162.91238209314, 2402.8969890530752] TIMING : CubeSum implot [ 1.47274000e-01 6.06929898e+00 1.54601172e+03 2.78074219e+02] TIMING : CubeSum getdata [ 2.12100000e-02 2.12881565e-02 1.54601172e+03 2.78074219e+02] TIMING : CubeSum done [ 2.32812000e-01 2.55270958e-01 1.54601172e+03 2.78187500e+02] TIMING : CubeSum END [ 0.701869 6.70024586] 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 [ 1.14728340e+01 1.64728995e+09] TIMING : SFind2D BEGIN [ 0. 0.] ('BDP', 0, ) ('BDP', 1, ) ('BDP', 2, ) No PB given TIMING : SFind2D start [ 3.34800000e-02 3.35850716e-02 1.54601172e+03 2.78187500e+02] ('ROBUST:', {'fence': 1.5, 'algorithm': 'hinges-fences'}) TIMING : SFind2D imstat [ 4.09240000e-02 4.11748886e-02 1.54601172e+03 2.78187500e+02] INFO : AT.py : Setting 'sigma' = 0.219800873397 for SFind2D_AT INFO : SFind2D_AT.py : sigma, dmin, dmax, snmax, cutoff 0.219801 -1.73031 199.734 35 0.0285714 WARNING : SFind2D_AT.py : Temporarely patching your K.km/s units to Jy/beam for ia.findsources() 2022-03-14 20:32:26 WARN imstat::ImageRegion::fromRecord There is no Stokes coordinate in the CoordinateSystem - assuming Stokes I TIMING : SFind2D findsources [ 9.90340000e-02 1.08254194e-01 1.54601172e+03 2.78187500e+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:58:25.513 +51.23.18.71 47.19 47.87 146 283 22.997 18.432 86.9 664.1 INFO : SFind2D_AT.py : 23:58:25.511 +51.23.18.70 47.19 47.87 146 282 22.848 18.365 86.5 664.1 INFO : SFind2D_AT.py : 23:58:25.549 +51.23.18.71 47.14 47.87 146 270 21.809 18.049 85.7 664.1 INFO : SFind2D_AT.py : Wrote ds9.reg TIMING : SFind2D table [ 1.82020000e-01 1.86734915e-01 1.54601172e+03 2.78187500e+02] REGRESSION : CONTFLUX: 3 835.388 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 WARNING : SFind2D_AT.py : LogScaling applied /lma1/teuben/lmtoy/admit/admit/at/SFind2D_AT.py:429: RuntimeWarning: invalid value encountered in log10 data = np.where(data<0,-np.log10(1-data),+np.log10(1+data)) TIMING : SFind2D done [ 3.00506000e-01 3.20414066e-01 1.54601172e+03 2.78953125e+02] TIMING : SFind2D END [ 0.6694 0.70363617] INFO : AT.py : BDP_OUT[0] = SourceList_BDP x-csm.sl TIMING : ADMITrun END [ 10.332075 16.74209404] INFO : Admit.py : ADMIT run() finished [flowcount 1] [cpu 10.3321 16.7421 ] ('N sources in CSM:', 3) ('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 [ 1.24281620e+01 1.64728995e+09] TIMING : CubeSpectrum BEGIN [ 0. 0.] INFO : CubeSpectrum_AT.py : CubeStats::maxpos,val=[47, 48, 1035],68.046791 TIMING : CubeSpectrum CubeStats-pos [ 4.80920000e-02 4.82120514e-02 1.54601172e+03 2.78953125e+02] INFO : CubeSpectrum_AT.py : CubeSum::maxpos,val=[47, 48],199.734421 TIMING : CubeSpectrum Moment-pos [ 4.21800000e-02 4.23099995e-02 1.54601172e+03 2.78953125e+02] INFO : CubeSpectrum_AT.py : SourceList::maxpos,val=('23h58m25.513s', '+51d23m18.71s'),145.963531 TIMING : CubeSpectrum SourceList-pos [ 3.59070000e-02 3.60019207e-02 1.54601172e+03 2.78953125e+02] TIMING : CubeSpectrum open [ 1.21380000e-02 1.21750832e-02 1.54601172e+03 2.78953125e+02] TIMING : CubeSpectrum imval [ 5.39060000e-02 5.41660786e-02 1.54601172e+03 2.79386719e+02] TIMING : CubeSpectrum imhead [ 4.13080000e-02 4.14707661e-02 1.54601172e+03 2.79386719e+02] TIMING : CubeSpectrum imval [ 2.37807000e-01 2.56999016e-01 1.54576172e+03 2.79183594e+02] REGRESSION : CSP: [68.046791076660156, 68.046791076660156] INFO : CubeSpectrum_AT.py : Writing 2 testCubeSpectrum tables TIMING : CubeSpectrum done [ 2.84368000e-01 3.04475069e-01 1.54576172e+03 2.79183594e+02] TIMING : CubeSpectrum summary [ 1.23280000e-02 1.23620033e-02 1.54576172e+03 2.79183594e+02] TIMING : CubeSpectrum END [ 0.78015 0.82032204] 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.32682020e+01 1.64728995e+09] TIMING : LineSegment BEGIN [ 0. 0.] TIMING : LineSegment start [ 2.42990000e-02 2.43721008e-02 1.54576172e+03 2.79183594e+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 [ 16.474974 16.5152998 1545.76171875 279.56640625] INFO : LineSegment_AT.py : Detecting segments in CubeSpectrum based data WARNING : specutil.py : 0 [1023, 1057] WARNING : specutil.py : 0 [1023, 1057] INFO : LineSegment_AT.py : Detecting segments in CubeStats based data WARNING : specutil.py : 0 [1025, 1056] WARNING : specutil.py : 0 [278, 279] TIMING : LineSegment segment finder [ 3.098927 3.10740519 1545.76171875 279.56640625] INFO : LineSegment_AT.py : Segment Coverage 35 / 2062 = 0.0169738 REGRESSION : LINESEG: [[86.234086588342478, 86.237113931775951], [86.309965486658825, 86.310063142898628], [86.23398893210269, 86.237309244255542], [86.23398893210269, 86.237309244255542]] TIMING : LineSegment done [ 2.20365 2.31569386 1561.74609375 295.40625 ] TIMING : LineSegment END [ 21.814254 21.9752059] INFO : AT.py : BDP_OUT[0] = LineSegment_BDP x.lseg TIMING : ADMITrun END [ 33.250381 39.91922712] INFO : Admit.py : ADMIT run() finished [flowcount 1] [cpu 33.2504 39.9192 ] Found 2 segments in LineSegment No ContinuumSub needed INFO : Admit.py : ADMIT run() called [flowcount 1] TIMING : ADMITrun END [ 33.35203 40.08479095] INFO : Admit.py : ADMIT run() finished [flowcount 1] [cpu 33.352 40.0848 ] ('OBJECT = ', 'R-Cas') ('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 [ 3.54708560e+01 1.64728997e+09] TIMING : LineID BEGIN [ 0. 0.] INFO : LineID_AT.py : Set vlsr = 23.86 for line identification. INFO : LineID_AT.py : Identifylines = True INFO : LineID_AT.py : Using vlsr = 23.8642 INFO : LineID_AT.py : Attempting Continuum Subtraction for Input Spectra TIMING : LineID getspectrum-cubespecs [ 5.722472 5.73703694 1561.74609375 295.43359375] INFO : LineID_AT.py : Attempting Continuum Subtraction for Input CubeStats Spectra TIMING : LineID getspectrum-cubestats [ 10.900076 10.93450284 1561.74609375 295.43359375] TIMING : LineID segment finder [ 1.45880000e-02 1.46470070e-02 1.56174609e+03 2.95433594e+02] INFO : LineID_AT.py : Detecting segments in CubeSpectrum based data WARNING : specutil.py : 0 [1023, 1057] WARNING : specutil.py : 0 [1023, 1057] INFO : LineID_AT.py : Detecting segments in CubeStats based data WARNING : specutil.py : 0 [1025, 1056] WARNING : specutil.py : 0 [278, 279] INFO : LineID_AT.py : Searching for spectral peaks with method: PeakFinder INFO : LineID_AT.py : Found line: SiO 2-1 @ 86.24337GHz, channels 1023 - 1057 REGRESSION : LINEID: SiO 86.24337 1023 1057 INFO : LineID_AT.py : Found line: U_86.3168 @ 86.3168365066GHz, channels 278 - 279 REGRESSION : LINEID: U_86.3168 86.31684 278 279 INFO : LineID_AT.py : Line Coverage 37 / 2062 = 0.0179437 TIMING : LineID done [ 14.771008 14.99448609 1670.1640625 404.80078125] TIMING : LineID END [ 31.422148 31.69470978] INFO : AT.py : BDP_OUT[0] = LineList_BDP x.ll TIMING : ADMITrun END [ 65.288032 72.37793612] INFO : Admit.py : ADMIT run() finished [flowcount 1] [cpu 65.288 72.3779 ] 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 [ 6.73616420e+01 1.64729000e+09] TIMING : LineCube BEGIN [ 0. 0.] TIMING : LineCube start [ 2.07999000e-01 2.08599091e-01 1.62693750e+03 3.61667969e+02] TIMING : LineCube pad [ 1.40720000e-02 1.41220093e-02 1.62693750e+03 3.61667969e+02] TIMING : LineCube trans-x.SiO_86.24337 [ 8.45520000e-02 2.01071978e-01 1.62693750e+03 3.61785156e+02] TIMING : LineCube trans-x.U_86.3168 [ 7.23140000e-02 1.25529051e-01 1.62693750e+03 3.61785156e+02] REGRESSION : LC: [1018, 1062, 273, 284] TIMING : LineCube done [ 2.48770000e-02 2.49409676e-02 1.62693750e+03 3.61785156e+02] TIMING : LineCube END [ 0.416142 0.58663011] INFO : AT.py : BDP_OUT[0] = LineCube_BDP x.SiO_86.24337/lc.im INFO : AT.py : BDP_OUT[1] = LineCube_BDP x.U_86.3168/lc.im TIMING : ADMITrun END [ 65.906153 73.22687697] INFO : Admit.py : ADMIT run() finished [flowcount 1] [cpu 65.9062 73.2269 ] Found 2 lines in LineCube ('LineDir:', 0, 'x.SiO_86.24337/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_86.3168/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 [ 6.81273810e+01 1.64729000e+09] TIMING : Moment BEGIN [ 0. 0.] TIMING : Moment open [ 4.92010000e-02 4.93259430e-02 1.62693750e+03 3.61785156e+02] TIMING : Moment immoments-0 [ 2.23221000e-01 2.66345024e-01 1.63884766e+03 3.73785156e+02] TIMING : Moment mom0clip [ 6.84350000e-02 1.14454985e-01 1.63884766e+03 3.73886719e+02] TIMING : Moment implot [ 9.58890000e-02 3.50378036e-01 1.63884766e+03 3.73886719e+02] TIMING : Moment getdata [ 2.08340000e-02 2.09059715e-02 1.63884766e+03 3.73886719e+02] TIMING : Moment ren+mask_0 [ 3.32866000e-01 3.38657856e-01 1.63996484e+03 3.74910156e+02] TIMING : Moment done [ 1.50220000e-02 1.50661469e-02 1.63996484e+03 3.74910156e+02] TIMING : Moment END [ 0.817766 1.16749382] 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.90062390e+01 1.64729000e+09] TIMING : Moment BEGIN [ 0. 0.] TIMING : Moment open [ 4.81800000e-02 4.83169556e-02 1.63996484e+03 3.74910156e+02] TIMING : Moment immoments-all [ 7.15300000e-02 1.88380003e-01 1.63996484e+03 3.74910156e+02] TIMING : Moment mom0clip [ 6.70790000e-02 9.87770557e-02 1.63996484e+03 3.74910156e+02] TIMING : Moment implot [ 9.20980000e-02 3.43288898e-01 1.63996484e+03 3.74910156e+02] TIMING : Moment getdata [ 2.04440000e-02 2.05230713e-02 1.63996484e+03 3.74910156e+02] TIMING : Moment ren+mask_0 [ 2.57596000e-01 2.63177156e-01 1.63996484e+03 3.74910156e+02] 2022-03-14 20:33:25 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.48327000e-01 3.51189852e-01 1.63996484e+03 3.75355469e+02] TIMING : Moment implot [ 7.88440000e-02 3.12722921e-01 1.63996484e+03 3.75761719e+02] TIMING : Moment getdata [ 2.01040000e-02 2.01790333e-02 1.63996484e+03 3.75761719e+02] TIMING : Moment ren+mask_1 [ 2.58110000e-01 2.65694141e-01 1.63996484e+03 3.75761719e+02] 2022-03-14 20:33:26 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.48447000e-01 3.73209000e-01 1.63996484e+03 3.75761719e+02] TIMING : Moment implot [ 7.69350000e-02 3.01884890e-01 1.63996484e+03 3.75761719e+02] TIMING : Moment getdata [ 2.01360000e-02 2.02100277e-02 1.63996484e+03 3.75761719e+02] TIMING : Moment ren+mask_2 [ 2.60062000e-01 2.65140057e-01 1.63971484e+03 3.75511719e+02] REGRESSION : MOM0FLUX: x.SiO_86.24337 2754.24 1997.14 0.0177778 25.2249 25.2249 0.0420509 TIMING : Moment mom0flux [ 5.13100000e-02 5.14769554e-02 1.63971484e+03 3.75511719e+02] TIMING : Moment flux-spectrum [ 2.73127000e-01 2.78117895e-01 1.63971484e+03 3.75511719e+02] TIMING : Moment done [ 1.49190000e-02 1.49610043e-02 1.63971484e+03 3.75511719e+02] TIMING : Moment END [ 1.919584 3.22963595] INFO : AT.py : BDP_OUT[0] = Moment_BDP x.SiO_86.24337/lc.mom_0 INFO : AT.py : BDP_OUT[1] = Moment_BDP x.SiO_86.24337/lc.mom_1 INFO : AT.py : BDP_OUT[2] = Moment_BDP x.SiO_86.24337/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 [ 7.10351450e+01 1.64729001e+09] TIMING : Moment BEGIN [ 0. 0.] TIMING : Moment open [ 4.81770000e-02 4.82859612e-02 1.63971484e+03 3.75511719e+02] TIMING : Moment immoments-all [ 6.85820000e-02 2.10186005e-01 1.63971484e+03 3.75511719e+02] TIMING : Moment mom0clip [ 6.73200000e-02 9.15620327e-02 1.63971484e+03 3.75511719e+02] TIMING : Moment implot [ 9.10770000e-02 3.33803892e-01 1.63971484e+03 3.75511719e+02] TIMING : Moment getdata [ 2.02310000e-02 2.03061104e-02 1.63971484e+03 3.75511719e+02] TIMING : Moment ren+mask_0 [ 2.60966000e-01 2.65969992e-01 1.63971484e+03 3.75511719e+02] 2022-03-14 20:33:29 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.48219000e-01 3.46889019e-01 1.63971484e+03 3.75511719e+02] TIMING : Moment implot [ 7.66170000e-02 2.95742035e-01 1.63971484e+03 3.75511719e+02] TIMING : Moment getdata [ 2.01250000e-02 2.01890469e-02 1.63971484e+03 3.75511719e+02] TIMING : Moment ren+mask_1 [ 2.52367000e-01 2.57299900e-01 1.63971484e+03 3.75511719e+02] 2022-03-14 20:33:30 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.48112000e-01 3.53309870e-01 1.63971484e+03 3.75511719e+02] TIMING : Moment implot [ 7.69150000e-02 2.96504021e-01 1.63971484e+03 3.75511719e+02] TIMING : Moment getdata [ 1.99500000e-02 2.00161934e-02 1.63971484e+03 3.75511719e+02] TIMING : Moment ren+mask_2 [ 2.57663000e-01 2.62599945e-01 1.63971484e+03 3.75511719e+02] REGRESSION : MOM0FLUX: x.U_86.3168 20.3201 10.0121 0.0177778 23.5037 23.5037 0.731677 TIMING : Moment mom0flux [ 5.14470000e-02 5.16078472e-02 1.63971484e+03 3.75511719e+02] TIMING : Moment flux-spectrum [ 2.59770000e-01 2.64991045e-01 1.63971484e+03 3.75511719e+02] TIMING : Moment done [ 1.47820000e-02 1.48320198e-02 1.63971484e+03 3.75511719e+02] TIMING : Moment END [ 1.894639 3.16647601] INFO : AT.py : BDP_OUT[0] = Moment_BDP x.U_86.3168/lc-@1.mom_0 INFO : AT.py : BDP_OUT[1] = Moment_BDP x.U_86.3168/lc-@1.mom_1 INFO : AT.py : BDP_OUT[2] = Moment_BDP x.U_86.3168/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 [ 7.30378460e+01 1.64729001e+09] TIMING : CubeSpectrum BEGIN [ 0. 0.] INFO : CubeSpectrum_AT.py : CubeSum::maxpos,val=[47, 48],200.128494 TIMING : CubeSpectrum Moment-pos [ 5.43800000e-02 5.45411110e-02 1.63971484e+03 3.75511719e+02] TIMING : CubeSpectrum open [ 1.23730000e-02 1.24149323e-02 1.63971484e+03 3.75511719e+02] TIMING : CubeSpectrum imval [ 3.61610000e-02 3.63478661e-02 1.63971484e+03 3.75511719e+02] TIMING : CubeSpectrum imhead [ 4.21450000e-02 4.23080921e-02 1.63971484e+03 3.75511719e+02] REGRESSION : CSP: [68.046791076660156] INFO : CubeSpectrum_AT.py : Writing 1 testCubeSpectrum tables TIMING : CubeSpectrum done [ 2.59525000e-01 2.64437914e-01 1.63971484e+03 3.75562500e+02] TIMING : CubeSpectrum summary [ 1.25680000e-02 1.26039982e-02 1.63971484e+03 3.75562500e+02] TIMING : CubeSpectrum END [ 0.429475 0.4350121] INFO : AT.py : BDP_OUT[0] = CubeSpectrum_BDP x.SiO_86.24337/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 [ 7.35278530e+01 1.64729001e+09] TIMING : CubeSpectrum BEGIN [ 0. 0.] INFO : CubeSpectrum_AT.py : CubeSum::maxpos,val=[14, 34],0.380805 TIMING : CubeSpectrum Moment-pos [ 5.43710000e-02 5.45320511e-02 1.63971484e+03 3.75562500e+02] TIMING : CubeSpectrum open [ 1.22900000e-02 1.23260021e-02 1.63971484e+03 3.75562500e+02] TIMING : CubeSpectrum imval [ 3.58200000e-02 3.60050201e-02 1.63971484e+03 3.75562500e+02] TIMING : CubeSpectrum imhead [ 4.20460000e-02 4.22089100e-02 1.63971484e+03 3.75562500e+02] REGRESSION : CSP: [0.44393777847290039] INFO : CubeSpectrum_AT.py : Writing 1 testCubeSpectrum tables TIMING : CubeSpectrum done [ 2.61013000e-01 2.66053200e-01 1.63971484e+03 3.75562500e+02] TIMING : CubeSpectrum summary [ 1.24560000e-02 1.24919415e-02 1.63971484e+03 3.75562500e+02] TIMING : CubeSpectrum END [ 0.43025 0.43590498] INFO : AT.py : BDP_OUT[0] = CubeSpectrum_BDP x.U_86.3168/lc-@2.csp TIMING : ADMITrun END [ 72.076174 82.43299794] INFO : Admit.py : ADMIT run() finished [flowcount 1] [cpu 72.0762 82.433 ]