LMTOY>> seq_pipeline: 11-jan-2022
LMTOY>> READING lmtoy_functions 24-dec-2021 via /lma1/teuben/lmtoy/examples/seq_pipeline.sh
Working directory ./2018S1SEQUOIACommissioning/82695
LMTOY>> Hang on, creating a bootstrap lmtoy_82695.rc from path=/n/lma1/lmt/data_lmt/
# GLOBS: /n/lma1/lmt/data_lmt//ifproc/ifproc*082695*.nc
#
# version=11-jan-2022
# ifproc="/n/lma1/lmt/data_lmt//ifproc/ifproc_2019-03-21_082695_01_0000.nc"
date_obs="2019-03-21T03:51:07"
# skytime=676.413 sec
# inttime=597.72 sec
obsnum="82695"
calobsnum="82694"
obspgm="Map"
ProjectId="2018S1SEQUOIACommissioning"
# SkyOff=900 0
# HPBW=16 arcsec
# XYLength=400 400 arcsec
# XYRamp=0 0 arcsec
# XYoff=0 0 arcsec
# XYstep=1 0.5
numbands=1
vlsr=-20 # km/s
skyfreq=115.27 # GHz
restfreq=115.271 # Ghz
src="IRC+10216"
resolution=11 # arcsec
cell=5.5 # arcsec
x_extent=400 # arcsec
y_extent=400 # arcsec
instrument="SEQ"
#
LMTOY>> this is your startup lmtoy_82695.rc file:
# seq_pipeline: 11-jan-2022
# DATE: 2022-01-11T21:59:23.788815224
# obsnum=82695
# Using ifproc=/n/lma1/lmt/data_lmt//ifproc/ifproc_2019-03-21_082695_01_0000.nc
# path=/n/lma1/lmt/data_lmt/
# GLOBS: /n/lma1/lmt/data_lmt//ifproc/ifproc*082695*.nc
#
# version=11-jan-2022
# ifproc="/n/lma1/lmt/data_lmt//ifproc/ifproc_2019-03-21_082695_01_0000.nc"
date_obs="2019-03-21T03:51:07"
# skytime=676.413 sec
# inttime=597.72 sec
obsnum="82695"
calobsnum="82694"
obspgm="Map"
ProjectId="2018S1SEQUOIACommissioning"
# SkyOff=900 0
# HPBW=16 arcsec
# XYLength=400 400 arcsec
# XYRamp=0 0 arcsec
# XYoff=0 0 arcsec
# XYstep=1 0.5
numbands=1
vlsr=-20 # km/s
skyfreq=115.27 # GHz
restfreq=115.271 # Ghz
src="IRC+10216"
resolution=11 # arcsec
cell=5.5 # arcsec
x_extent=400 # arcsec
y_extent=400 # arcsec
instrument="SEQ"
#
# based on vlsr=-20, dv=100, dw=250
b_order=0
b_regions=[[-370,-120],[80,330]]
l_regions=[[-120,80]]
slice=[-370,330]
v_range=-120,80
pix_list=0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15
rmax=3
otf_a=1.1
otf_b=4.75
otf_c=2
sample=-1
otf_cal=0
edge=0
LMTOY>> Sleeping for 5 seconds, you can abort, edit lmtoy_82695.rc, then continuing
LMTOY>> process_otf_map2 in 2 seconds
Traceback (most recent call last):
File "/lma1/teuben/lmtoy/SpectralLineReduction/lmtslr/ifproc/ifproc.py", line 581, in __init__
self.chop = self.nc.variables['Data.Msip1mm.BeamChopperActPos'][:]
KeyError: 'Data.Msip1mm.BeamChopperActPos'
Traceback (most recent call last):
File "/lma1/teuben/lmtoy/SpectralLineReduction/lmtslr/ifproc/ifproc.py", line 712, in __init__
self.chop = self.nc.variables['Data.Msip1mm.BeamChopperActPos'][:]
KeyError: 'Data.Msip1mm.BeamChopperActPos'
/lma1/teuben/lmtoy/anaconda3/lib/python3.8/site-packages/numpy/core/fromnumeric.py:748: UserWarning: Warning: 'partition' will ignore the 'mask' of the MaskedArray.
a.partition(kth, axis=axis, kind=kind, order=order)
/lma1/teuben/lmtoy/SpectralLineReduction/lmtslr/spec/spec.py:547: RuntimeWarning: invalid value encountered in true_divide
self.tsys_spectrum = 280 * sky_spectrum / (hot_spectrum -
/lma1/teuben/lmtoy/SpectralLineReduction/lmtslr/spec/spec.py:547: RuntimeWarning: divide by zero encountered in true_divide
self.tsys_spectrum = 280 * sky_spectrum / (hot_spectrum -
/lma1/teuben/lmtoy/SpectralLineReduction/lmtslr/spec/spec.py:392: RuntimeWarning: invalid value encountered in true_divide
spectra.append((self.raw_spec[i,:] - ref) / ref)
/lma1/teuben/lmtoy/SpectralLineReduction/lmtslr/spec/spec.py:392: RuntimeWarning: divide by zero encountered in true_divide
spectra.append((self.raw_spec[i,:] - ref) / ref)
/lma1/teuben/lmtoy/SpectralLineReduction/lmtslr/spec/spec.py:434: RuntimeWarning: invalid value encountered in multiply
self.reduced_spectra = np.array(spectra) * tsys_spectrum
{'--b_order': '0',
'--b_regions': '[[-370,-120],[80,330]]',
'--bank': '0',
'--eliminate_list': '0',
'--help': False,
'--l_regions': '[[-120,80]]',
'--obsnum': '82695',
'--output': 'IRC+10216_82695.nc',
'--path': '/n/lma1/lmt/data_lmt/',
'--pix_list': '0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15',
'--restfreq': None,
'--sample': None,
'--save_tsys': False,
'--slice': '[-370,330]',
'--stype': '2',
'--tsys': '250.0',
'--use_cal': True,
'--use_otf_cal': False,
'--x_axis': 'VLSR'}
program process_otf_map options
data_path = /n/lma1/lmt/data_lmt/
output = IRC+10216_82695.nc
use_otf_cal = False
path = /n/lma1/lmt/data_lmt/
slice = [-370, 330]
b_order = 0
l_regions = [[-120, 80]]
output_file_name = IRC+10216_82695.nc
bank = 0
b_regions = [[-370, -120], [80, 330]]
pix_list = [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15]
use_cal = True
tsys = 250.0
stype = 2
obsnum = 82695
x_axis = VLSR
eliminate_list = [0]
2019-03-21T03:51:07 begin /n/lma1/lmt/data_lmt/ifproc/ifproc_2019-03-21_082695_01_0000.nc
2019-03-21T04:02:24 end /n/lma1/lmt/data_lmt/ifproc/ifproc_2019-03-21_082695_01_0000.nc
BufPos 1 37.8 sec 82695
BufPos 0 277.8 sec 82695
BufPos 1 38.3 sec 82695
BufPos 0 276.8 sec 82695
BufPos 1 37.2 sec 82695
BufPos 0 8.5 sec 82695
before read npix
from pixels npix = 16
from xlen npix = 16
TRACKING Sequoia PIXEL 10
Map Parameters: Az Continuous
HPBW= 16.0 XLength= 400.0 YLength= 400.0 XStep= 1.00 YStep= 0.50
/n/lma1/lmt/data_lmt/ifproc/ifproc_2019-03-21_082695_01_0000.nc does not have bs parameters
82695 is a Map observation
get chop
'Data.Msip1mm.BeamChopperActPos'
no chop
2019-03-21T03:50:45 begin /n/lma1/lmt/data_lmt/ifproc/ifproc_2019-03-21_082694_00_0001.nc
2019-03-21T03:50:52 end /n/lma1/lmt/data_lmt/ifproc/ifproc_2019-03-21_082694_00_0001.nc
BufPos 3 4.7 sec 82694
BufPos 2 2.6 sec 82694
before read npix
from pixels npix = 16
from xlen npix = 16
TRACKING Sequoia PIXEL 10
/n/lma1/lmt/data_lmt/ifproc/ifproc_2019-03-21_082694_00_0001.nc does not have map parameters
/n/lma1/lmt/data_lmt/ifproc/ifproc_2019-03-21_082694_00_0001.nc does not have bs parameters
82694 is a Cal observation
get chop cal
'Data.Msip1mm.BeamChopperActPos'
no chop cal
read_roach /n/lma1/lmt/data_lmt/spectrometer/roach0/roach0_82695_1_0_IRC+10216_2019-03-21_075102.nc nspec,nchan=25740,2048
r:0 inp:0 pix:0 time_offset:-0.030000
r:0 inp:1 pix:1 time_offset:-0.030000
r:0 inp:2 pix:2 time_offset:-0.030000
r:0 inp:3 pix:3 time_offset:-0.030000
read_roach /n/lma1/lmt/data_lmt/spectrometer/roach1/roach1_82695_1_0_IRC+10216_2019-03-21_075102.nc nspec,nchan=25752,2048
r:1 inp:0 pix:4 time_offset:-0.030000
r:1 inp:1 pix:5 time_offset:-0.030000
r:1 inp:2 pix:6 time_offset:-0.030000
r:1 inp:3 pix:7 time_offset:-0.030000
read_roach /n/lma1/lmt/data_lmt/spectrometer/roach2/roach2_82695_1_0_IRC+10216_2019-03-21_075102.nc nspec,nchan=25764,2048
r:2 inp:0 pix:8 time_offset:-0.030000
r:2 inp:1 pix:9 time_offset:-0.030000
r:2 inp:2 pix:10 time_offset:-0.030000
r:2 inp:3 pix:11 time_offset:-0.030000
read_roach /n/lma1/lmt/data_lmt/spectrometer/roach3/roach3_82695_1_0_IRC+10216_2019-03-21_075102.nc nspec,nchan=25772,2048
r:3 inp:0 pix:12 time_offset:-0.030000
r:3 inp:1 pix:13 time_offset:-0.030000
r:3 inp:2 pix:14 time_offset:-0.030000
r:3 inp:3 pix:15 time_offset:-0.030000
RESTFREQ 115.271 GHz
read_roach /n/lma1/lmt/data_lmt/spectrometer/roach0/roach0_82694_0_1_IRC+10216_2019-03-21_075043.nc nspec,nchan=232,2048
r:0 inp:0 pix:0 time_offset:-0.030000
r:0 inp:1 pix:1 time_offset:-0.030000
r:0 inp:2 pix:2 time_offset:-0.030000
r:0 inp:3 pix:3 time_offset:-0.030000
read_roach /n/lma1/lmt/data_lmt/spectrometer/roach1/roach1_82694_0_1_IRC+10216_2019-03-21_075044.nc nspec,nchan=236,2048
r:1 inp:0 pix:4 time_offset:-0.030000
r:1 inp:1 pix:5 time_offset:-0.030000
r:1 inp:2 pix:6 time_offset:-0.030000
r:1 inp:3 pix:7 time_offset:-0.030000
read_roach /n/lma1/lmt/data_lmt/spectrometer/roach2/roach2_82694_0_1_IRC+10216_2019-03-21_075044.nc nspec,nchan=240,2048
r:2 inp:0 pix:8 time_offset:-0.030000
r:2 inp:1 pix:9 time_offset:-0.030000
r:2 inp:2 pix:10 time_offset:-0.030000
r:2 inp:3 pix:11 time_offset:-0.030000
read_roach /n/lma1/lmt/data_lmt/spectrometer/roach3/roach3_82694_0_1_IRC+10216_2019-03-21_075044.nc nspec,nchan=240,2048
r:3 inp:0 pix:12 time_offset:-0.030000
r:3 inp:1 pix:13 time_offset:-0.030000
r:3 inp:2 pix:14 time_offset:-0.030000
r:3 inp:3 pix:15 time_offset:-0.030000
TSYS[0] = 456.186 +/- 23.7841 (2048 channels)
TSYS[1] = 397.114 +/- 19.8351 (2048 channels)
TSYS[2] = 417.307 +/- 20.7008 (2048 channels)
TSYS[3] = 396.862 +/- 21.1605 (2048 channels)
TSYS[4] = 418.743 +/- 28.4395 (2048 channels)
TSYS[5] = 416.628 +/- 23.4353 (2048 channels)
TSYS[6] = 463.078 +/- 30.4436 (2048 channels)
TSYS[7] = 402.422 +/- 18.2658 (2048 channels)
TSYS[8] = 400.786 +/- 23.0008 (2048 channels)
TSYS[9] = 414.749 +/- 23.7305 (2048 channels)
TSYS[10] = 412.256 +/- 24.7739 (2048 channels)
TSYS[11] = 411.703 +/- 23.4238 (2048 channels)
TSYS[12] = 379.153 +/- 18.3571 (2048 channels)
TSYS[13] = 434.848 +/- 19.9951 (2048 channels)
TSYS[14] = 425.174 +/- 21.8982 (2048 channels)
TSYS[15] = 431.067 +/- 24.4929 (2048 channels)
Spectral Band velocity range: -1059.8 1019.8 km/s
Total Number of OTF Spectra = 87904
NSPEC: 87904
NCHAN: 689
FAST_NC: True
Looping over 16 pixel list [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15]:
Processing 1 CAL's for Tsys
Pix Nspec Mean Std MAD_std Min Max RMS_max Warnings
TSYS[0] slice: 456.35 (10.4741) minmax: 436.516 479.47
0 5493 -0.023 2.411 2.403 -13.803 14.863 2.402 4.138
TSYS[1] slice: 396.612 (7.10052) minmax: 382.134 410.055
1 5493 0.017 2.106 2.085 -20.409 22.884 2.090 7.890 *M 1.4
TSYS[2] slice: 416.489 (7.99118) minmax: 401.413 430.718
2 5493 0.014 2.201 2.186 -18.366 14.722 2.188 5.375
TSYS[3] slice: 396.372 (8.00738) minmax: 380.648 411.661
3 5493 0.017 2.098 2.082 -14.634 13.057 2.086 4.798
TSYS[4] slice: 418.074 (9.62367) minmax: 397.511 436.476
4 5494 0.010 2.204 2.191 -12.977 19.232 2.197 4.785
TSYS[5] slice: 420.527 (9.00789) minmax: 400.314 437.001
5 5494 0.135 2.557 2.485 -14.513 14.968 2.596 5.826 *M 1.7
TSYS[6] slice: 468.918 (11.0601) minmax: 440.629 487.227
6 5494 0.014 2.473 2.459 -20.274 18.224 2.454 6.339
TSYS[7] slice: 407.377 (8.04546) minmax: 386.422 423.662
7 5494 0.011 2.151 2.140 -11.682 14.409 2.137 5.596
TSYS[8] slice: 399.611 (7.85109) minmax: 384.44 416.12
8 5494 0.008 2.111 2.098 -11.506 14.556 2.103 4.018
TSYS[9] slice: 413.179 (10.146) minmax: 395.5 433.922
9 5494 0.009 2.186 2.169 -20.877 17.842 2.176 7.431 *M 1.3
TSYS[10] slice: 411.849 (11.3658) minmax: 391.192 434.103
10 5494 0.016 2.177 2.165 -16.854 17.293 2.166 5.850
TSYS[11] slice: 413.388 (10.607) minmax: 393.19 432.642
11 5494 0.017 2.181 2.169 -15.025 15.348 2.170 5.686
TSYS[12] slice: 378.766 (6.6079) minmax: 363.434 391.583
12 5495 0.008 1.999 1.989 -10.458 10.277 1.992 3.097
TSYS[13] slice: 434.674 (7.64801) minmax: 418.149 450.811
13 5495 0.008 2.296 2.283 -18.188 19.510 2.287 7.244 *M 1.3
TSYS[14] slice: 424.291 (7.86416) minmax: 406.965 439.786
14 5495 0.007 2.244 2.227 -22.656 15.765 2.231 6.597 *M 1.2
TSYS[15] slice: 429.424 (9.02467) minmax: 411.365 448.269
15 5495 0.016 2.265 2.254 -13.246 14.302 2.257 5.101
CPU TIME: 44.9526 sec
CPU TIME: 45.0729 sec
Warnings:
*P ratio: ratio of std/mad too high for data
*M ratio: ratio of std/mad too high for RMS
X-range: -234.593 263.861 Y-range: -259.949 -5.11764 arcsec
Written IRC+10216_82695.nc
LMTOY>> view_spec_file
{'--help': False,
'--input': 'IRC+10216_82695.nc',
'--pix_list': '0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15',
'--plot_range': '-1,3',
'--plots': 'IRC+10216_82695_specviews',
'--rms_cut': '-4',
'--show_pixel': None,
'--skip_tsys': False}
SLR Plots initialized with 1 method=IRC+10216_82695_specviews
Max sequence=5495
Plots saving in IRC+10216_82695_specviews.1.png
Plots saving in IRC+10216_82695_specviews.2.png
Plots saving in IRC+10216_82695_specviews.3.png
Plots saving in IRC+10216_82695_specviews.4.png
Plots saving in IRC+10216_82695_specviews.5.png
NEW TSYS PIC
SHAPE (1, 16, 689)
PIXEL: 0 [ 0 1 2 ... 5490 5491 5492]
PIXEL: 1 [ 5493 5494 5495 ... 10983 10984 10985]
PIXEL: 2 [10986 10987 10988 ... 16476 16477 16478]
PIXEL: 3 [16479 16480 16481 ... 21969 21970 21971]
PIXEL: 4 [21972 21973 21974 ... 27463 27464 27465]
PIXEL: 5 [27466 27467 27468 ... 32957 32958 32959]
PIXEL: 6 [32960 32961 32962 ... 38451 38452 38453]
PIXEL: 7 [38454 38455 38456 ... 43945 43946 43947]
PIXEL: 8 [43948 43949 43950 ... 49439 49440 49441]
PIXEL: 9 [49442 49443 49444 ... 54933 54934 54935]
PIXEL: 10 [54936 54937 54938 ... 60427 60428 60429]
PIXEL: 11 [60430 60431 60432 ... 65921 65922 65923]
PIXEL: 12 [65924 65925 65926 ... 71416 71417 71418]
PIXEL: 13 [71419 71420 71421 ... 76911 76912 76913]
PIXEL: 14 [76914 76915 76916 ... 82406 82407 82408]
PIXEL: 15 [82409 82410 82411 ... 87901 87902 87903]
Plots saving in IRC+10216_82695_specviews.6.png
/lma1/teuben/lmtoy/anaconda3/lib/python3.8/site-packages/numpy/core/fromnumeric.py:748: UserWarning: Warning: 'partition' will ignore the 'mask' of the MaskedArray.
a.partition(kth, axis=axis, kind=kind, order=order)
SLR Plots initialized with 3 method=IRC+10216_82695_specpoint,png,1
MAD: 0 2.3936436 0.15301155746673512 3.0056898474969453
MAD: 1 2.0809395 0.14702731925744209 2.669048808356062
MAD: 2 2.1806033 0.14895466849021152 2.7764219397231753
MAD: 3 2.0754359 0.1516729290312472 2.682127592971302
MAD: 4 2.1894493 0.149882023045401 2.7889774024843383
MAD: 5 2.4491549 0.21628640626756226 3.3143004788910497
MAD: 6 2.4435034 0.16004651472131143 3.0836894387070233
MAD: 7 2.1321795 0.1483244138029987 2.72547715388448
MAD: 8 2.094389 0.1549254744627608 2.7140908596430355
MAD: 9 2.1656165 0.15144846928566102 2.771410389441228
MAD: 10 2.1529245 0.150253883915254 2.7539400733197072
MAD: 11 2.1606004 0.1523936745765231 2.7701751221189586
MAD: 12 1.9751909 0.14852024167552585 2.569271844616532
MAD: 13 2.282539 0.1530705886124247 2.894821245288322
MAD: 14 2.2230742 0.15326146776615152 2.836120068833771
MAD: 15 2.2471201 0.15546753591135357 2.8689902856284464
Plots saving in IRC+10216_82695_specpoint.1.png
/lma1/teuben/lmtoy/anaconda3/lib/python3.8/site-packages/numpy/core/fromnumeric.py:748: UserWarning: Warning: 'partition' will ignore the 'mask' of the MaskedArray.
a.partition(kth, axis=axis, kind=kind, order=order)
/lma1/teuben/lmtoy/anaconda3/lib/python3.8/site-packages/numpy/ma/core.py:5243: RuntimeWarning: Mean of empty slice.
result = super(MaskedArray, self).mean(axis=axis,
/lma1/teuben/lmtoy/anaconda3/lib/python3.8/site-packages/numpy/core/_methods.py:162: RuntimeWarning: invalid value encountered in true_divide
ret = um.true_divide(
SLR Plots initialized with 3 method=IRC+10216_82695_specpoint,png,2
MAD: 0 2.3936436 0.15301155746673512 3.0056898474969453
MAD: 1 2.0809395 0.14702731925744209 2.669048808356062
MAD: 2 2.1806033 0.14895466849021152 2.7764219397231753
MAD: 3 2.0754359 0.1516729290312472 2.682127592971302
MAD: 4 2.1894493 0.149882023045401 2.7889774024843383
MAD: 5 2.4491549 0.21628640626756226 3.3143004788910497
MAD: 6 2.4435034 0.16004651472131143 3.0836894387070233
MAD: 7 2.1321795 0.1483244138029987 2.72547715388448
MAD: 8 2.094389 0.1549254744627608 2.7140908596430355
MAD: 9 2.1656165 0.15144846928566102 2.771410389441228
MAD: 10 2.1529245 0.150253883915254 2.7539400733197072
MAD: 11 2.1606004 0.1523936745765231 2.7701751221189586
MAD: 12 1.9751909 0.14852024167552585 2.569271844616532
MAD: 13 2.282539 0.1530705886124247 2.894821245288322
MAD: 14 2.2230742 0.15326146776615152 2.836120068833771
MAD: 15 2.2471201 0.15546753591135357 2.8689902856284464
Plots saving in IRC+10216_82695_specpoint.2.png
LMTOY>> make_spec_fits (waterfall)
pix [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15]
Min nsamp 5493
Binning [1, 1]
FITS file will be 5493 x 689 x 16
Written waterfall cube to IRC+10216_82695.wf.fits
pix [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15]
Min nsamp 5493
Binning [10, 1]
FITS file will be 549 x 689 x 16
New shape: (16, 689, 549)
Written waterfall cube to IRC+10216_82695.wf10.fits
LMTOY>> grid_data
{'--cell': '5.5',
'--edge': '0',
'--help': False,
'--input': 'IRC+10216_82695.nc',
'--n_samples': '256',
'--noise_sigma': '1',
'--otf_a': '1.1',
'--otf_b': '4.75',
'--otf_c': '2',
'--otf_select': '1',
'--output': 'IRC+10216_82695.fits',
'--pix_list': '0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15',
'--program_path': 'spec_driver_fits',
'--resolution': '11',
'--rmax': '3',
'--rms_cut': '-4',
'--sample': '-1',
'--weight': 'IRC+10216_82695.wt.fits',
'--x_extent': '400',
'--y_extent': '400'}
STDOUT ****************************
spec_driver_fits LMTSLR 26-sep-2021
CMD: spec_driver_fits -i IRC+10216_82695.nc -o IRC+10216_82695.fits -w IRC+10216_82695.wt.fits -l 11 -c 5.5 -u [0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15] -z -4 -s 1 -x 400 -y 400 -f 1 -r 3 -n 256 -0 1.1 -1 4.75 -2 2 -b -1 -e 0
file list=IRC+10216_82695.nc
DEBUG:5.5
sample_list=-1
Warning: sample list is not a multiple of 3
nsamples=1
PJT: OTF->fuzzy_edge = 0
Processing 1 SpecFiles:
Opening SpecFile file IRC+10216_82695.nc
file: IRC+10216_82695.nc nspec= 87904 nchan= 689
SpecFile version 6-mar-2021
X-range: -234.593 263.861 Y-range: -259.949 -5.11764 arcsec
X-ramp: -151.193 180.461 Y-ramp: -176.549 -88.5176 arcsec
MapSize: 498.454 x 254.831 arcsec
82695
IRC+10216
DATE-OBS 2019-03-21T03:51:07
CF.n_cells= 7 delta=0.129412 cell=5.5 resol=11 oft_select=1
C-AXIS 2 starts at 329.985
C-AXIS 1 starts at -401.5
C-AXIS 0 starts at -401.5
P-AXIS 1 starts at -401.5
P-AXIS 0 starts at -401.5
P-AXIS 1 starts at -401.5
P-AXIS 0 starts at -401.5
WARNING fuzzy_edge=0
pixel: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
used?: 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
rms_stats: n=87904 rms_cutoff=-4
Cube Completed, 87806/87904 Spectra accepted = 0.999
Weighting Completed, fuzzy_edge=0
Weight at 0.000000 0.000000 is 0.087984
write to IRC+10216_82695.fits
write weights to IRC+10216_82695.wt.fits
STDERR ****************************
Exit Code: 0
LMTOY>> Some NEMO post-processing
### nemo Debug Info: FITS file: Image size 147 147 689
### Warning [fitsccd]: Setting map minimum from 0 to -678.918
### Warning [fitsccd]: Setting map maximum from 0 to 622.681
### nemo Debug Info: There were 12054055 blank values in the image (80.961636 %)
### nemo Debug Info: FITS file: Image size 147 147
### Warning [fitsccd]: Setting map minimum from 0 to -0.925052
### Warning [fitsccd]: Setting map maximum from 0 to 18.5612
### nemo Debug Info: There were no blank values set in the image
LMTOY>> STATS IRC+10216_82695.ccd centerbox robust
### Warning [ccdsub]: Sampling will be done
### Warning [ccdsub]: Attempting to fix the WCS
### nemo Debug Info: WCS Corner: -0.055 0.055 329985
Number of points : 1828606
Min and Max : -9.089821 11.094045
Mean and dispersion : 0.020324 0.407110
Skewness and kurtosis : 0.618963 15.067163
Sum and Sum*Dx*Dy*Dz* : 37163.936800 88.125614
Points per beam (map) : 1
N Robust : 1773062
Mean Robust : 0.013104
Sigma Robust : 0.329061
Median Robust : 0.012060
1944358/3772964 out-of-range points discarded
LMTOY>> STATS IRC+10216_82695.wt.ccd centerbox robust
### Warning [ccdsub]: Sampling will be done
### Warning [ccdsub]: Attempting to fix the WCS
### nemo Debug Info: WCS Corner: -0.055 0.055 0
Number of points : 3076
Min and Max : -0.255010 18.561224
Mean and dispersion : 7.483670 4.472125
Skewness and kurtosis : -0.328111 -1.097391
Sum and Sum*Dx*Dy* : 23019.768131 -0.053731
Points per beam (map) : 1
N Robust : 3076
Mean Robust : 7.483670
Sigma Robust : 4.472125
Median Robust : 8.330002
2400/5476 out-of-range points discarded
### nemo Debug Info: Keeping 147*147*689 cube
### Warning [ccdfits]: Replacing header from IRC+10216_82695.wt.fits
SLABS: [[-370,-120],[80,330]] == -370,-120,80,330
### Warning [ccdslice]: axis=1 not fully supported yet
### nemo Debug Info: Axis: 689 0 329985 -1015.92
### nemo Debug Info: Axis min/max: 329985 -368969
### nemo Debug Info: slab -370 -120 -> 689 442
### nemo Debug Info: slab 80 330 -> 246 0
### nemo Debug Info: Total planes: 495
### Warning [ccdslice]: Code not converted to fix reference pixel value
### Warning [ccdfits]: Replacing header from IRC+10216_82695.wt.fits
### Warning [ccdfits]: Replacing header from IRC+10216_82695.wt.fits
### Warning [ccdfits]: Replacing header from IRC+10216_82695.fits
[0;31mQAC_STATS: IRC+10216_82695-full 0.0119256 3.41516 -678.918 622.681 33803.6 0.0291235[0m
### Warning [ccdsub]: Sampling will be done
### Warning [ccdsub]: Attempting to fix the WCS
### nemo Debug Info: WCS Corner: -0.055 0.055 329985
[0;31mQAC_STATS: IRC+10216_82695-cent 0.0203236 0.40711 -9.08982 11.094 37163.9 0.069258[0m
### nemo Debug Info: FITS file: Image size 147 147 689
### Warning [fitsccd]: Setting map minimum from 0 to -0.366753
### Warning [fitsccd]: Setting map maximum from 0 to 1.24726
### nemo Debug Info: There were no blank values set in the image
spectab : # 73 73 689 min/mean/sig/max 0 0 0 0 0 0 -nan
specstab: # 73 73 689 min/mean/sig/max -0.00898264 0.00332505 0.0184213 0.106165 -6.83785e-07 0.00208153 0.0799003
Traceback (most recent call last):
File "/lma1/teuben/lmtoy/examples/fitsplot.py", line 21, in
f = aplpy.FITSFigure(fitsfile)
File "/lma1/teuben/lmtoy/anaconda3/lib/python3.8/site-packages/aplpy/decorators.py", line 28, in wrapper
return f(*args, **kwargs)
File "/lma1/teuben/lmtoy/anaconda3/lib/python3.8/site-packages/aplpy/core.py", line 232, in __init__
self._get_hdu(data, hdu, north, convention=convention,
File "/lma1/teuben/lmtoy/anaconda3/lib/python3.8/site-packages/aplpy/core.py", line 309, in _get_hdu
raise IOError("File not found: " + filename)
OSError: File not found: IRC+10216_82695.mom0.fits
Traceback (most recent call last):
File "/lma1/teuben/lmtoy/examples/fitsplot.py", line 21, in
f = aplpy.FITSFigure(fitsfile)
File "/lma1/teuben/lmtoy/anaconda3/lib/python3.8/site-packages/aplpy/decorators.py", line 28, in wrapper
return f(*args, **kwargs)
File "/lma1/teuben/lmtoy/anaconda3/lib/python3.8/site-packages/aplpy/core.py", line 232, in __init__
self._get_hdu(data, hdu, north, convention=convention,
File "/lma1/teuben/lmtoy/anaconda3/lib/python3.8/site-packages/aplpy/core.py", line 309, in _get_hdu
raise IOError("File not found: " + filename)
OSError: File not found: IRC+10216_82695.mom1.fits
Traceback (most recent call last):
File "/lma1/teuben/lmtoy/examples/fitsplot.py", line 21, in
f = aplpy.FITSFigure(fitsfile)
File "/lma1/teuben/lmtoy/anaconda3/lib/python3.8/site-packages/aplpy/decorators.py", line 28, in wrapper
return f(*args, **kwargs)
File "/lma1/teuben/lmtoy/anaconda3/lib/python3.8/site-packages/aplpy/core.py", line 232, in __init__
self._get_hdu(data, hdu, north, convention=convention,
File "/lma1/teuben/lmtoy/anaconda3/lib/python3.8/site-packages/aplpy/core.py", line 309, in _get_hdu
raise IOError("File not found: " + filename)
OSError: File not found: IRC+10216_82695.mom2.fits
INFO: Auto-setting vmin to -1.647e+00 [aplpy.core]
INFO: Auto-setting vmax to 1.523e+01 [aplpy.core]
INFO: Auto-setting vmin to -1.723e+00 [aplpy.core]
INFO: Auto-setting vmax to 1.550e+01 [aplpy.core]
Writing IRC+10216_82695.wt.png
INFO: Auto-setting vmin to -1.765e+00 [aplpy.core]
INFO: Auto-setting vmax to 1.959e+01 [aplpy.core]
INFO: Auto-setting vmin to -1.819e+00 [aplpy.core]
INFO: Auto-setting vmax to 2.019e+01 [aplpy.core]
Writing IRC+10216_82695.wt2.png
INFO: Auto-setting vmin to -1.709e+02 [aplpy.core]
INFO: Auto-setting vmax to 1.897e+03 [aplpy.core]
INFO: Auto-setting vmin to -3.011e+02 [aplpy.core]
INFO: Auto-setting vmax to 3.342e+03 [aplpy.core]
Writing IRC+10216_82695.wt3.png
Traceback (most recent call last):
File "/lma1/teuben/lmtoy/examples/fitsplot.py", line 21, in
f = aplpy.FITSFigure(fitsfile)
File "/lma1/teuben/lmtoy/anaconda3/lib/python3.8/site-packages/aplpy/decorators.py", line 28, in wrapper
return f(*args, **kwargs)
File "/lma1/teuben/lmtoy/anaconda3/lib/python3.8/site-packages/aplpy/core.py", line 232, in __init__
self._get_hdu(data, hdu, north, convention=convention,
File "/lma1/teuben/lmtoy/anaconda3/lib/python3.8/site-packages/aplpy/core.py", line 309, in _get_hdu
raise IOError("File not found: " + filename)
OSError: File not found: IRC+10216_82695.wtn.fits
INFO: Auto-setting vmin to -1.075e-01 [aplpy.core]
INFO: Auto-setting vmax to 1.194e+00 [aplpy.core]
INFO: Auto-setting vmin to -1.072e-01 [aplpy.core]
INFO: Auto-setting vmax to 1.190e+00 [aplpy.core]
Writing IRC+10216_82695.wtr.png
WARNING: FITSFixedWarning: RADECSYS= 'FK5 ' /
the RADECSYS keyword is deprecated, use RADESYSa. [astropy.wcs.wcs]
INFO: Setting slices=[0] [aplpy.core]
INFO: Auto-setting vmin to -1.439e+04 [aplpy.core]
INFO: Auto-setting vmax to 2.977e+04 [aplpy.core]
INFO: Auto-setting vmin to -1.492e+04 [aplpy.core]
INFO: Auto-setting vmax to 3.046e+04 [aplpy.core]
Writing IRC+10216_82695.mom0.png
LMTOY>> Created IRC+10216_82695.nf.fits and IRC+10216_82695.nfs.fits
LMTOY>> ADMIT post-processing
LMTOY>> lmtoy_admit: 31-oct-2021
/lma1/teuben/lmtoy/examples/lmtoy_admit.sh: line 33: export: `IRC+10216_82695.nfs.fits': not a valid identifier
No PDIR directory used, all work in the current directory
ADMIT = /lma1/teuben/lmtoy/admit
version = 1.0.8.5
CASA_PATH = /n/astromake/opt/casa/58 linux admit lma.astro.umd.edu
CASA_ROOT = /n/astromake/opt/casa/58
prefix = /an/terra_astromake/opt/casa/casa-release-5.8.0-109.el7
version = 5.8.0-rel-109
revision = 1
Trying CASA6, no worries if this fails
Traceback (most recent call last):
File "", line 1, in
ImportError: No module named casatools
LMTOY>> ADMIT post-processing
Removing suspect directory TempLattice*_*
Checking casapy log files Found 1
Checking ipython log files Found 1
Checking casa log files Found 1
Cleaning 1 casapy*.log files
Cleaning 1 ipython-*.log files
Cleaning 1 casa-*.log files
### (1/1) : processing IRC+10216_82695.nfs.fits
IRC+10216_82695.nfs.fits [147, 147, 689] IRC+10216 0.0 0.0 0.0003299848 0.055
/lma1/teuben/lmtoy/admit/admit/test/admit1.py --basename x IRC+10216_82695.nfs.fits
194.170u 8.136s 4:43.70 71.3% 0+0k 416464+936360io 944pf+0w
Logfile: IRC+10216_82695.nfs.fits.log
LMTOY>> lmtoy_admit: 31-oct-2021
/lma1/teuben/lmtoy/examples/lmtoy_admit.sh: line 33: export: `IRC+10216_82695.nf.fits': not a valid identifier
No PDIR directory used, all work in the current directory
ADMIT = /lma1/teuben/lmtoy/admit
version = 1.0.8.5
CASA_PATH = /n/astromake/opt/casa/58 linux admit lma.astro.umd.edu
CASA_ROOT = /n/astromake/opt/casa/58
prefix = /an/terra_astromake/opt/casa/casa-release-5.8.0-109.el7
version = 5.8.0-rel-109
revision = 1
Trying CASA6, no worries if this fails
Traceback (most recent call last):
File "", line 1, in
ImportError: No module named casatools
LMTOY>> ADMIT post-processing
Removing suspect directory TempLattice*_*
Checking casapy log files Found 1
Checking ipython log files Found 1
Checking casa log files Found 1
Cleaning 1 casapy*.log files
Cleaning 1 ipython-*.log files
Cleaning 1 casa-*.log files
### (1/1) : processing IRC+10216_82695.nf.fits
IRC+10216_82695.nf.fits [147, 147, 689] IRC+10216 0.0 0.0 0.0003299848 0.055
/lma1/teuben/lmtoy/admit/admit/test/admit1.py --basename x IRC+10216_82695.nf.fits
64.187u 2.360s 3:00.35 36.8% 0+0k 1032+293712io 2pf+0w
Logfile: IRC+10216_82695.nf.fits.log
LMTOY>> Created IRC+10216_82695.fits and IRC+10216_82695.wt.fits
LMTOY>> Parameter file used: lmtoy_82695.rc
LMTOY>> Making summary index.html:
Making index.html for obsnum=82695
Writing index.html
Writing index_pipeline.html
Writing index_admit.html
Writing index_pars.html
Writing index_log.html