========================================= The start-up time of CASA may vary depending on whether the shared libraries are cached or not. ========================================= CASA Version 4.2.2 (r30986) Compiled on: Thu 2014/09/04 17:59:42 UTC casaviewer: could not start Xvfb ___________________________________________________________________ For help use the following commands: tasklist - Task list organized by category taskhelp - One line summary of available tasks help taskname - Full help for task toolhelp - One line summary of available tools help par.parametername - Full help for parameter name ___________________________________________________________________ will execute script /home/teuben/ADMIT/ALMAvis/astute/bin/test1 Activating auto-logging. Current session state plus future input saved. Filename : ipython-20140918-173801.log Mode : backup Output logging : False Raw input log : False Timestamping : False State : active *** Loading ATNF ASAP Package... *** ... ASAP (4.2.0a rev#30794) import complete *** ################################################## Major interface changes to SINGLE DISH tasks have been taken place in CASA 4.2.2 release The interface of the following tasks are modified: sdbaseline, sdcal, sdcal2, sdfit, sdflag, sdgrid, sdimaging, sdmath, sdplot, sdreduce, sdsave, and sdstat. Additionally, a new task called sdaverage is available. Task sdsmooth has been incorporated in the new task and removed. The tasks with old interfaces are available with name {taskname}old. They will be kept until CASA 4.3 release and removed from later releases. Users are advised to update existing scripts. ################################################## ['/home/teuben/ADMIT/ALMAvis/astute/bin/test1', '999'] Trying mode 999 TRY1a BDP(FILE,fits,None) AT(FILE,[0],[1]) AT_file.init AT_file.run AT_DATE=FILE=Thu Sep 18 13:38:01 EDT 2014 BDP.pdump(fits.p) BDP(FILE,cim,None) AT_flow.init AT(FLOW,[1],[1]) BDP cim depends on BDP_FILE(fits): special flow keys: ('self', 'keyvals', 'kv') special flow vals: None AT_flow.run(FLOW) AT_DATE=FLOW=Thu Sep 18 13:38:01 EDT 2014 work_flow: 1 -> 1 BDP.pdump(cim.p) BDP.pdump(fits.p) BDP(SUMMARY,summary,None) AT(SUMMARY,[1],[1]) BDP summary depends on BDP_FILE(cim): AT_DATE=SUMMARY=Thu Sep 18 13:38:01 EDT 2014 ============================================================================ AXIS ['Right Ascension' 'Declination' 'Frequency' 'Stokes'] UNIT ['rad' 'rad' 'Hz' ''] SHAP [640 640 246 1] INCR [ -1.93925472e-06 1.93925472e-06 7.59541583e+06 1.00000000e+00] BMAJ {'value': 4.401422977446, 'unit': 'arcsec'} COOR [ 2.07486613e-01 -4.41362243e-01 1.13779323e+11 1.00000000e+00] ============================================================================ BDP.pdump(summary.p) BDP.pdump(cim.p) BDP(CUBESTATS,cubestats,None) AT(CUBESTATS,[1],[1]) BDP cubestats depends on BDP_FILE(cim): AT_DATE=CUBESTATS=Thu Sep 18 13:38:02 EDT 2014 casa::imhead(cim) *** reffreqtype in get before *** reffreqtype in get LSRK casa::imstat(cim) Cube Stats: number of points: 17119379 mean, sigma: 0.273594 1.36351 (6.57325) mJy/beam min: -16.2805467844 mJy/beam @: [329 336 193 0] max: 659.840106964 mJy/beam @: [327 313 194 0] Plane Stats: N mean std (all) N mean std (robust) mean :: (246, 0.27362257102568238, 0.74791136119946566, 156, 0.0080285586976625214, 0.016183423031758848) madm :: (246, 1.3777742695207216, 0.17143201773732017, 175, 1.3074084298153008, 0.062798216056489761) max :: (246, 44.276835518033522, 101.25720383457779, 168, 7.4753716831938144, 0.80379082265788104) sigma :: (246, 3.814451094711198, 5.301189993368979, 178, 1.9461221242259847, 0.10515835170011939) s/n :: (246, 26.467562079832884, 53.579195721194772, 162, 5.6764872009853438, 0.47796643970575098) Freq range : 111.918 113.779 GHz Mean range : -0.0775488 3.98419 mJy/beam Noise range : 1.19297 2.12259 mJy/beam Peak range : 5.93869 659.84 mJy/beam S/N range : 4.48543 310.866 mJy/beam SIGNAL CONFIDENCE FACTOR: 638.51585978 SIGNAL/NOISE: 504.693171557 BDP.pdump(cubestats.p) BDP.pdump(cim.p) Masking outside -6.61525 7.16249 17119379/100761600 (16.99 %) data are non-zero GaussPlot(0.273623,1.37777,583311) BDP(CUBESPECTRUM,cubespectrum,None) AT(CUBESPECTRUM,[2],[1]) BDP cubespectrum depends on BDP_FILE(cim): BDP cubespectrum depends on BDP_CUBESTATS(cubestats): AT_DATE=CUBESPECTRUM=Thu Sep 18 13:38:25 EDT 2014 IMVALS: [327, 313, 0, 0] [327, 313, 245, 0] *** reffreqtype in get before *** reffreqtype in get LSRK Freq range : 111.918 113.779 GHz Data range : -5.24892 659.84 mJy/beam BDP.pdump(cubespectrum.p) BDP.pdump(cim.p) BDP.pdump(cubestats.p) BDP(FILE,linelist,None) AT(LINELIST,[2],[1]) BDP linelist depends on BDP_FILE(cim): BDP linelist depends on BDP_CUBESTATS(cubestats): AT_DATE=LINELIST=Thu Sep 18 13:38:27 EDT 2014 Median: 0.770049050286 Q1,2,3: 0.73404132456 0.770744472539 0.935958951962 f1,f3: 0.532123697157 1.13787657937 ROBUST: 0.773196091105 0.754880325624 0.0753404509951 SEGMENTS: f=1 nsigma=2 ngap=2 nline=3 data range: 0.651804180037 2.49257272125 cutoff: 0.923876993096 Found 3 lines Line in channels 41 - 52 @ 112.268 GHz +/- 0.0235659 GHz log(S/N) = 1.32607 Line in channels 141 - 165 @ 113.077 GHz +/- 0.0476785 GHz log(S/N) = 2.23727 Line in channels 185 - 207 @ 113.405 GHz +/- 0.0430925 GHz log(S/N) = 2.49257 GaussPlot(0.773196,0.0753405,39) BDP.pdump(linelist.p) BDP.pdump(cim.p) BDP.pdump(cubestats.p) nothing to print here yet BDP(FILE,linecube,None) AT(LINECUBE,[2],[1]) BDP linecube depends on BDP_FILE(cim): BDP linecube depends on BDP_FILE(linelist): AT_DATE=LINECUBE=Thu Sep 18 13:38:30 EDT 2014 Found 3 lines Cutting a cube linecube.U-112.268 chans=41~52 @ 112.268 GHz Cutting a cube linecube.U-113.077 chans=141~165 @ 113.077 GHz BDP(FILE,linecube.U-113.077,None) Cutting a cube linecube.U-113.405 chans=185~207 @ 113.405 GHz BDP(FILE,linecube.U-113.405,None) BDP.pdump(linecube.U-112.268.p) BDP.pdump(linecube.U-113.077.p) BDP.pdump(linecube.U-113.405.p) BDP.pdump(cim.p) BDP.pdump(linelist.p) nothing to print here yet Found 3 linecubes BDP(FILE,linecube.U-112.268.mom0,None) BDP(FILE,linecube.U-112.268.mom1,None) BDP(FILE,linecube.U-112.268.mom2,None) AT(MOMENTS,[1],[3]) BDP linecube.U-112.268.mom0 depends on BDP_FILE(linecube.U-112.268): BDP linecube.U-112.268.mom1 depends on BDP_FILE(linecube.U-112.268): BDP linecube.U-112.268.mom2 depends on BDP_FILE(linecube.U-112.268): AT_DATE=MOMENTS=Thu Sep 18 13:38:31 EDT 2014 BDP.pdump(linecube.U-112.268.mom0.p) BDP.pdump(linecube.U-112.268.mom1.p) BDP.pdump(linecube.U-112.268.mom2.p) BDP.pdump(linecube.U-112.268.p) nothing to print here yet BDP(FILE,linecube.U-113.077.mom0,None) BDP(FILE,linecube.U-113.077.mom1,None) BDP(FILE,linecube.U-113.077.mom2,None) AT(MOMENTS,[1],[3]) BDP linecube.U-113.077.mom0 depends on BDP_FILE(linecube.U-113.077): BDP linecube.U-113.077.mom1 depends on BDP_FILE(linecube.U-113.077): BDP linecube.U-113.077.mom2 depends on BDP_FILE(linecube.U-113.077): AT_DATE=MOMENTS=Thu Sep 18 13:38:39 EDT 2014 BDP.pdump(linecube.U-113.077.mom0.p) BDP.pdump(linecube.U-113.077.mom1.p) BDP.pdump(linecube.U-113.077.mom2.p) BDP.pdump(linecube.U-113.077.p) nothing to print here yet BDP(FILE,linecube.U-113.405.mom0,None) BDP(FILE,linecube.U-113.405.mom1,None) BDP(FILE,linecube.U-113.405.mom2,None) AT(MOMENTS,[1],[3]) BDP linecube.U-113.405.mom0 depends on BDP_FILE(linecube.U-113.405): BDP linecube.U-113.405.mom1 depends on BDP_FILE(linecube.U-113.405): BDP linecube.U-113.405.mom2 depends on BDP_FILE(linecube.U-113.405): AT_DATE=MOMENTS=Thu Sep 18 13:38:48 EDT 2014 BDP.pdump(linecube.U-113.405.mom0.p) BDP.pdump(linecube.U-113.405.mom1.p) BDP.pdump(linecube.U-113.405.mom2.p) BDP.pdump(linecube.U-113.405.p) nothing to print here yet ADMIT: pickle saving TRY1a.p END mcrit=6 leaving casapy...