LMTOY>> seq_pipeline: 2-feb-2022
LMTOY>> READING lmtoy_functions 11-mar-2022 via /home/lmtslr/lmtoy/examples/seq_pipeline.sh
Working directory /home/lmtslr/work_lmt/2022S1SEQUOIACommissioning/95125
LMTOY>> Hang on, creating a bootstrap lmtoy_95125.rc from path=/home/lmtslr/data_lmt3
#
# version=2-feb-2022
# ifproc="/home/lmtslr/data_lmt3/ifproc/ifproc_2022-03-04_095125_00_0001.nc"
date_obs="2022-03-03T22:31:01"
# skytime=1547.58 sec
# inttime=1071.52 sec
obsnum="95125"
calobsnum="95124"
obspgm="Map"
obsgoal="Pointing"
ProjectId="2022S1SEQUOIACommissioning"
# SkyOff=-720 0
# HPBW=15 arcsec
# XYLength=600 600 arcsec
# XYRamp=0 0 arcsec
# XYoff=0 0 arcsec
# XYstep=1 0.5
numbands=1
vlsr=-26.5 # km/s
skyfreq=115.27516851282712 # GHz
restfreq=115.271204 # Ghz
src="IRC+10216"
resolution=11 # arcsec
cell=5.5 # arcsec
x_extent=600 # arcsec
y_extent=600 # arcsec
instrument="SEQ"
#
LMTOY>> this is your startup lmtoy_95125.rc file:
# seq_pipeline: 2-feb-2022
# DATE: 2022-03-16T07:27:36.712045180
# obsnum=95125
# Using ifproc=/home/lmtslr/data_lmt3/ifproc/ifproc_2022-03-04_095125_00_0001.nc
# path=/home/lmtslr/data_lmt3
#
# version=2-feb-2022
# ifproc="/home/lmtslr/data_lmt3/ifproc/ifproc_2022-03-04_095125_00_0001.nc"
date_obs="2022-03-03T22:31:01"
# skytime=1547.58 sec
# inttime=1071.52 sec
obsnum="95125"
calobsnum="95124"
obspgm="Map"
obsgoal="Pointing"
ProjectId="2022S1SEQUOIACommissioning"
# SkyOff=-720 0
# HPBW=15 arcsec
# XYLength=600 600 arcsec
# XYRamp=0 0 arcsec
# XYoff=0 0 arcsec
# XYstep=1 0.5
numbands=1
vlsr=-26.5 # km/s
skyfreq=115.27516851282712 # GHz
restfreq=115.271204 # Ghz
src="IRC+10216"
resolution=11 # arcsec
cell=5.5 # arcsec
x_extent=600 # arcsec
y_extent=600 # arcsec
instrument="SEQ"
#
# based on vlsr=-26.5, dv=100, dw=250
b_order=0
b_regions=[[-376.5,-126.5],[73.5,323.5]]
l_regions=[[-126.5,73.5]]
slice=[-376.5,323.5]
v_range=-126.5,73.5
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_95125.rc, then continuing
LMTOY>> process_otf_map2 in 2 seconds
{'--b_order': '0',
'--b_regions': '[[-376.5,-126.5],[73.5,323.5]]',
'--bank': '0',
'--eliminate_list': '0',
'--help': False,
'--l_regions': '[[-126.5,73.5]]',
'--obsnum': '95125',
'--output': 'IRC+10216_95125.nc',
'--path': '/home/lmtslr/data_lmt3',
'--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': '[-376.5,323.5]',
'--stype': '2',
'--tsys': '250.0',
'--use_cal': True,
'--use_otf_cal': False,
'--x_axis': 'VLSR'}
program process_otf_map options
path = /home/lmtslr/data_lmt3
output_file_name = IRC+10216_95125.nc
stype = 2
b_order = 0
b_regions = [[-376.5, -126.5], [73.5, 323.5]]
bank = 0
data_path = /home/lmtslr/data_lmt3
pix_list = [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15]
eliminate_list = [0]
l_regions = [[-126.5, 73.5]]
use_cal = True
slice = [-376.5, 323.5]
use_otf_cal = False
obsnum = 95125
x_axis = VLSR
output = IRC+10216_95125.nc
tsys = 250.0
2022-03-03T22:31:01 begin /home/lmtslr/data_lmt3/ifproc/ifproc_2022-03-04_095125_00_0001.nc
2022-03-03T22:56:49 end /home/lmtslr/data_lmt3/ifproc/ifproc_2022-03-04_095125_00_0001.nc
BufPos 1 16.7 sec 95125
BufPos 0 79.2 sec 95125
BufPos 1 30.4 sec 95125
BufPos 0 81.6 sec 95125
BufPos 1 30.9 sec 95125
BufPos 0 79.4 sec 95125
BufPos 1 29.0 sec 95125
BufPos 0 81.1 sec 95125
BufPos 1 30.9 sec 95125
BufPos 0 81.3 sec 95125
BufPos 1 30.3 sec 95125
BufPos 0 81.0 sec 95125
BufPos 1 32.9 sec 95125
BufPos 0 81.5 sec 95125
BufPos 1 34.2 sec 95125
BufPos 0 80.2 sec 95125
BufPos 1 31.0 sec 95125
BufPos 0 80.7 sec 95125
BufPos 1 32.0 sec 95125
BufPos 0 81.6 sec 95125
BufPos 1 31.8 sec 95125
BufPos 0 82.1 sec 95125
BufPos 1 30.6 sec 95125
BufPos 0 81.8 sec 95125
BufPos 1 31.9 sec 95125
BufPos 0 81.7 sec 95125
BufPos 1 33.8 sec 95125
BufPos 0 58.0 sec 95125
BufPos 1 10.1 sec 95125
before read npix
from pixels npix = 16
from xlen npix = 16
TRACKING Sequoia PIXEL 10
Map Parameters: Az Continuous
HPBW= 15.0 XLength= 600.0 YLength= 600.0 XStep= 1.00 YStep= 0.50
/home/lmtslr/data_lmt3/ifproc/ifproc_2022-03-04_095125_00_0001.nc does not have bs parameters
95125 is a Map observation
get chop
'Data.Msip1mm.BeamChopperActPos'
no chop
2022-03-03T22:30:26 begin /home/lmtslr/data_lmt3/ifproc/ifproc_2022-03-04_095124_00_0001.nc
2022-03-03T22:30:33 end /home/lmtslr/data_lmt3/ifproc/ifproc_2022-03-04_095124_00_0001.nc
BufPos 3 5.2 sec 95124
BufPos 2 2.6 sec 95124
before read npix
from pixels npix = 16
from xlen npix = 16
TRACKING Sequoia PIXEL 10
/home/lmtslr/data_lmt3/ifproc/ifproc_2022-03-04_095124_00_0001.nc does not have map parameters
/home/lmtslr/data_lmt3/ifproc/ifproc_2022-03-04_095124_00_0001.nc does not have bs parameters
95124 is a Cal observation
get chop cal
'Data.Msip1mm.BeamChopperActPos'
no chop cal
read_roach /home/lmtslr/data_lmt3/spectrometer/roach0/roach0_95125_0_1_IRC+10216_2022-03-04_043052.nc nspec,nchan=52564,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 /home/lmtslr/data_lmt3/spectrometer/roach1/roach1_95125_0_1_IRC+10216_2022-03-04_043052.nc nspec,nchan=52628,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 /home/lmtslr/data_lmt3/spectrometer/roach2/roach2_95125_0_1_IRC+10216_2022-03-04_043052.nc nspec,nchan=52716,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 /home/lmtslr/data_lmt3/spectrometer/roach3/roach3_95125_0_1_IRC+10216_2022-03-04_043052.nc nspec,nchan=52784,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 /home/lmtslr/data_lmt3/spectrometer/roach0/roach0_95124_0_1_IRC+10216_2022-03-04_043024.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 /home/lmtslr/data_lmt3/spectrometer/roach1/roach1_95124_0_1_IRC+10216_2022-03-04_043024.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 /home/lmtslr/data_lmt3/spectrometer/roach2/roach2_95124_0_1_IRC+10216_2022-03-04_043024.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 /home/lmtslr/data_lmt3/spectrometer/roach3/roach3_95124_0_1_IRC+10216_2022-03-04_043024.nc nspec,nchan=244,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] = 219.692 +/- 18.0833 (2048 channels)
TSYS[1] = 186.935 +/- 8.96398 (2048 channels)
TSYS[2] = 239.16 +/- 13.4734 (2048 channels)
TSYS[3] = 175.525 +/- 9.32802 (2048 channels)
TSYS[4] = 192.356 +/- 12.0121 (2048 channels)
TSYS[5] = 183.423 +/- 10.0581 (2048 channels)
TSYS[6] = 252.37 +/- 21.7521 (2048 channels)
TSYS[7] = 184.404 +/- 9.58293 (2048 channels)
TSYS[8] = 175.363 +/- 10.0199 (2048 channels)
TSYS[9] = 185.528 +/- 9.64461 (2048 channels)
TSYS[10] = 179.988 +/- 9.88101 (2048 channels)
TSYS[11] = 184.848 +/- 10.035 (2048 channels)
TSYS[12] = 169.075 +/- 8.96375 (2048 channels)
TSYS[13] = 206.698 +/- 9.85703 (2048 channels)
TSYS[14] = 193.61 +/- 9.59571 (2048 channels)
TSYS[15] = 195.275 +/- 11.5 (2048 channels)
Spectral Band velocity range: -1066.3 1013.3 km/s
Total Number of OTF Spectra = 162032
NSPEC: 162032
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: 220.304 (4.72759) minmax: 211.471 229.589
0 10128 0.005 1.155 1.149 -6.154 17.777 1.140 1.286
TSYS[1] slice: 186.307 (3.37689) minmax: 179.727 192.792
1 10128 0.009 0.975 0.967 -5.093 18.009 0.964 1.085
TSYS[2] slice: 238.913 (5.54625) minmax: 227.365 249.371
2 10128 0.010 1.244 1.238 -6.868 17.703 1.234 1.397
TSYS[3] slice: 175.369 (3.68155) minmax: 167.831 182.798
3 10128 0.010 0.921 0.911 -4.541 19.022 0.908 1.020
TSYS[4] slice: 192.788 (4.24721) minmax: 183.899 200.619
4 10126 0.009 1.008 1.001 -5.197 19.766 0.996 1.125
TSYS[5] slice: 183.801 (3.53085) minmax: 176.643 189.698
5 10126 0.009 0.964 0.954 -4.750 19.875 0.951 1.070
TSYS[6] slice: 252.027 (8.00354) minmax: 237.159 267.295
6 10126 0.009 1.313 1.305 -7.235 20.568 1.303 1.450
TSYS[7] slice: 184.931 (4.03555) minmax: 176.927 192.646
7 10126 0.009 0.966 0.959 -5.055 18.278 0.955 1.069
TSYS[8] slice: 175.435 (3.58124) minmax: 168.117 182.039
8 10128 0.009 0.921 0.912 -4.721 19.093 0.908 1.022
TSYS[9] slice: 185.231 (3.47574) minmax: 178.049 192.031
9 10128 0.009 0.972 0.963 -5.011 20.076 0.959 1.083
TSYS[10] slice: 179.95 (3.85466) minmax: 172.851 186.525
10 10128 0.009 0.944 0.935 -4.742 19.788 0.932 1.056
TSYS[11] slice: 184.972 (3.94241) minmax: 177.242 191.737
11 10128 0.009 0.970 0.961 -5.352 18.692 0.957 1.070
TSYS[12] slice: 168.963 (3.31848) minmax: 161.997 175.167
12 10126 0.008 0.888 0.878 -4.774 17.922 0.876 0.976
TSYS[13] slice: 206.73 (3.64753) minmax: 199.084 213.891
13 10126 0.009 1.079 1.073 -5.405 16.926 1.069 1.204
TSYS[14] slice: 193.154 (3.35868) minmax: 186.235 199.699
14 10126 0.009 1.014 1.005 -5.338 19.008 1.002 1.116 Traceback (most recent call last):
File "/home/lmtslr/lmtoy/SpectralLineReduction/lmtslr/ifproc/ifproc.py", line 537, in __init__
self.chop = self.nc.variables['Data.Msip1mm.BeamChopperActPos'][:]
KeyError: 'Data.Msip1mm.BeamChopperActPos'
Traceback (most recent call last):
File "/home/lmtslr/lmtoy/SpectralLineReduction/lmtslr/ifproc/ifproc.py", line 668, in __init__
self.chop = self.nc.variables['Data.Msip1mm.BeamChopperActPos'][:]
KeyError: 'Data.Msip1mm.BeamChopperActPos'
/home/lmtslr/lmtoy/anaconda3/lib/python3.8/site-packages/numpy/core/fromnumeric.py:753: UserWarning: Warning: 'partition' will ignore the 'mask' of the MaskedArray.
a.partition(kth, axis=axis, kind=kind, order=order)
/home/lmtslr/lmtoy/SpectralLineReduction/lmtslr/spec/spec.py:547: RuntimeWarning: invalid value encountered in true_divide
self.tsys_spectrum = 280 * sky_spectrum / (hot_spectrum -
/home/lmtslr/lmtoy/SpectralLineReduction/lmtslr/spec/spec.py:547: RuntimeWarning: divide by zero encountered in true_divide
self.tsys_spectrum = 280 * sky_spectrum / (hot_spectrum -
/home/lmtslr/lmtoy/SpectralLineReduction/lmtslr/spec/spec.py:392: RuntimeWarning: invalid value encountered in true_divide
spectra.append((self.raw_spec[i,:] - ref) / ref)
/home/lmtslr/lmtoy/SpectralLineReduction/lmtslr/spec/spec.py:392: RuntimeWarning: divide by zero encountered in true_divide
spectra.append((self.raw_spec[i,:] - ref) / ref)
/home/lmtslr/lmtoy/SpectralLineReduction/lmtslr/spec/spec.py:434: RuntimeWarning: invalid value encountered in multiply
self.reduced_spectra = np.array(spectra) * tsys_spectrum
TSYS[15] slice: 195.081 (4.30867) minmax: 187.098 203.147
15 10126 0.009 1.023 1.014 -5.243 18.806 1.010 1.133
CPU TIME: 160.908 sec
CPU TIME: 161.934 sec
Warnings:
*P ratio: ratio of std/mad too high for data
*M ratio: ratio of std/mad too high for RMS
X-range: -383.8 394.351 Y-range: -378.471 347.82 arcsec
Written IRC+10216_95125.nc
LMTOY>> view_spec_file
{'--help': False,
'--input': 'IRC+10216_95125.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_95125_specviews',
'--rms_cut': '-4',
'--show_pixel': None,
'--skip_tsys': False}
SLR Plots initialized with 1 method=IRC+10216_95125_specviews
Max sequence=10128
Plots saving in IRC+10216_95125_specviews.1.png
Plots saving in IRC+10216_95125_specviews.2.png
Plots saving in IRC+10216_95125_specviews.3.png
Plots saving in IRC+10216_95125_specviews.4.png
Plots saving in IRC+10216_95125_specviews.5.png
NEW TSYS PIC
SHAPE (1, 16, 689)
PIXEL: 0 [ 0 1 2 ... 10125 10126 10127]
PIXEL: 1 [10128 10129 10130 ... 20253 20254 20255]
PIXEL: 2 [20256 20257 20258 ... 30381 30382 30383]
PIXEL: 3 [30384 30385 30386 ... 40509 40510 40511]
PIXEL: 4 [40512 40513 40514 ... 50635 50636 50637]
PIXEL: 5 [50638 50639 50640 ... 60761 60762 60763]
PIXEL: 6 [60764 60765 60766 ... 70887 70888 70889]
PIXEL: 7 [70890 70891 70892 ... 81013 81014 81015]
PIXEL: 8 [81016 81017 81018 ... 91141 91142 91143]
PIXEL: 9 [ 91144 91145 91146 ... 101269 101270 101271]
PIXEL: 10 [101272 101273 101274 ... 111397 111398 111399]
PIXEL: 11 [111400 111401 111402 ... 121525 121526 121527]
PIXEL: 12 [121528 121529 121530 ... 131651 131652 131653]
PIXEL: 13 [131654 131655 131656 ... 141777 141778 141779]
PIXEL: 14 [141780 141781 141782 ... 151903 151904 151905]
PIXEL: 15 [151906 151907 151908 ... 162029 162030 162031]
Plots saving in IRC+10216_95125_specviews.6.png
/home/lmtslr/lmtoy/anaconda3/lib/python3.8/site-packages/numpy/core/fromnumeric.py:753: 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_95125_specpoint,png,1
MAD: 0 1.1399972 0.03625643481066694 1.2850229831238933
MAD: 1 0.96382415 0.030957947635133237 1.0876559434870052
MAD: 2 1.2336667 0.03885787024002351 1.3890981393615833
MAD: 3 0.9078587 0.02918153432577371 1.0245848666655826
MAD: 4 0.9962938 0.031938412547447635 1.1240474333776567
MAD: 5 0.95009124 0.029739590740488934 1.069049605752178
MAD: 6 1.3029687 0.0414856399229961 1.4689113001552412
MAD: 7 0.9548574 0.030149627441087293 1.0754559187647459
MAD: 8 0.9080374 0.029151311793100296 1.0246426712599257
MAD: 9 0.9588224 0.031162435765560954 1.083472172242389
MAD: 10 0.9320334 0.030249220406943833 1.0530303012368452
MAD: 11 0.9563801 0.030896795609957794 1.0799673113003048
MAD: 12 0.8755839 0.027774419156888438 0.9866815637277735
MAD: 13 1.0690732 0.033862845570924026 1.204524582509526
MAD: 14 1.0020592 0.03221461291549077 1.1309176729296633
MAD: 15 1.0101984 0.032788620111339185 1.1413528351664262
Plots saving in IRC+10216_95125_specpoint.1.png
/home/lmtslr/lmtoy/anaconda3/lib/python3.8/site-packages/numpy/core/fromnumeric.py:753: 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_95125_specpoint,png,2
MAD: 0 1.1399972 0.03625643481066694 1.2850229831238933
MAD: 1 0.96382415 0.030957947635133237 1.0876559434870052
MAD: 2 1.2336667 0.03885787024002351 1.3890981393615833
MAD: 3 0.9078587 0.02918153432577371 1.0245848666655826
MAD: 4 0.9962938 0.031938412547447635 1.1240474333776567
MAD: 5 0.95009124 0.029739590740488934 1.069049605752178
MAD: 6 1.3029687 0.0414856399229961 1.4689113001552412
MAD: 7 0.9548574 0.030149627441087293 1.0754559187647459
MAD: 8 0.9080374 0.029151311793100296 1.0246426712599257
MAD: 9 0.9588224 0.031162435765560954 1.083472172242389
MAD: 10 0.9320334 0.030249220406943833 1.0530303012368452
MAD: 11 0.9563801 0.030896795609957794 1.0799673113003048
MAD: 12 0.8755839 0.027774419156888438 0.9866815637277735
MAD: 13 1.0690732 0.033862845570924026 1.204524582509526
MAD: 14 1.0020592 0.03221461291549077 1.1309176729296633
MAD: 15 1.0101984 0.032788620111339185 1.1413528351664262
Plots saving in IRC+10216_95125_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 10126
Binning [1, 1]
FITS file will be 10126 x 689 x 16
Written waterfall cube to IRC+10216_95125.wf.fits
pix [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15]
Min nsamp 10126
Binning [10, 1]
FITS file will be 1012 x 689 x 16
New shape: (16, 689, 1012)
Written waterfall cube to IRC+10216_95125.wf10.fits
LMTOY>> grid_data
{'--cell': '5.5',
'--edge': '0',
'--help': False,
'--input': 'IRC+10216_95125.nc',
'--n_samples': '256',
'--noise_sigma': '1',
'--otf_a': '1.1',
'--otf_b': '4.75',
'--otf_c': '2',
'--otf_select': '1',
'--output': 'IRC+10216_95125.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_95125.wt.fits',
'--x_extent': '600',
'--y_extent': '600'}
STDOUT ****************************
spec_driver_fits LMTSLR 4-mar-2022
CMD: spec_driver_fits -i IRC+10216_95125.nc -o IRC+10216_95125.fits -w IRC+10216_95125.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 600 -y 600 -f 1 -r 3 -n 256 -0 1.1 -1 4.75 -2 2 -b -1 -e 0
file list=IRC+10216_95125.nc
DEBUG:5.5
sample_list=-1
Warning: sample list is not a multiple of 3
nsamples=1
Processing 1 SpecFiles:
Opening SpecFile file IRC+10216_95125.nc
file: IRC+10216_95125.nc nspec= 162032 nchan= 689
SpecFile version 6-mar-2021
X-range: -383.8 394.35 Y-range: -378.471 347.82 arcsec
X-ramp: -300.4 310.95 Y-ramp: -295.071 264.42 arcsec
MapSize: 778.15 x 726.292 arcsec
95125
IRC+10216
DATE-OBS 2022-03-03T22:31:01
CF.n_cells= 7 delta=0.129412 cell=5.5 resol=11 oft_select=1
C-AXIS 2 starts at 323.485
C-AXIS 1 starts at -599.5
C-AXIS 0 starts at -599.5
P-AXIS 1 starts at -599.5
P-AXIS 0 starts at -599.5
P-AXIS 1 starts at -599.5
P-AXIS 0 starts at -599.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=162032 rms_cutoff=-4
Cube Completed, 162020/162032 Spectra accepted = 1.000
Weighting Completed, fuzzy_edge=0
Weight at 0.000000 0.000000 is 15.261176
write to IRC+10216_95125.fits
write weights to IRC+10216_95125.wt.fits
STDERR ****************************
Exit Code: 0
LMTOY>> Some NEMO post-processing
### nemo Debug Info: FITS file: Image size 219 219 689
### Warning [fitsccd]: Setting map minimum from 0 to -92.8787
### Warning [fitsccd]: Setting map maximum from 0 to 90.2376
### nemo Debug Info: There were 20915284 blank values in the image (63.293092 %)
### nemo Debug Info: FITS file: Image size 219 219
### Warning [fitsccd]: Setting map minimum from 0 to -4.35213
### Warning [fitsccd]: Setting map maximum from 0 to 68.1624
### nemo Debug Info: There were no blank values set in the image
LMTOY>> STATS IRC+10216_95125.ccd centerbox robust
### Warning [ccdsub]: Sampling will be done
### Warning [ccdsub]: Attempting to fix the WCS
### nemo Debug Info: WCS Corner: -0.0825 0.0825 323485
Number of points : 8336900
Min and Max : -2.469591 17.220646
Mean and dispersion : 0.011364 0.285189
Skewness and kurtosis : 6.744974 221.667699
Sum and Sum*Dx*Dy*Dz* : 94737.065168 224.646869
Points per beam (map) : 1
N Robust : 8198941
Mean Robust : 0.002110
Sigma Robust : 0.222115
Median Robust : 0.001655
0/8336900 out-of-range points discarded
LMTOY>> STATS IRC+10216_95125.wt.ccd centerbox robust
### Warning [ccdsub]: Sampling will be done
### Warning [ccdsub]: Attempting to fix the WCS
### nemo Debug Info: WCS Corner: -0.0825 0.0825 0
Number of points : 12100
Min and Max : 3.236960 68.162407
Mean and dispersion : 16.911558 5.391124
Skewness and kurtosis : 2.896578 13.390643
Sum and Sum*Dx*Dy* : 204629.854847 -0.477628
Points per beam (map) : 1
N Robust : 10188
Mean Robust : 15.960218
Sigma Robust : 1.714760
Median Robust : 16.183936
0/12100 out-of-range points discarded
### nemo Debug Info: Keeping 219*219*689 cube
### Warning [ccdfits]: Replacing header from IRC+10216_95125.wt.fits
SLABS: [[-376.5,-126.5],[73.5,323.5]] == -376.5,-126.5,73.5,323.5
### Warning [ccdslice]: axis=1 not fully supported yet
### nemo Debug Info: Axis: 689 0 323485 -1015.92
### nemo Debug Info: Axis min/max: 323485 -375469
### nemo Debug Info: slab -376.5 -126.5 -> 689 442
### nemo Debug Info: slab 73.5 323.5 -> 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_95125.wt.fits
### Warning [ccdfits]: Replacing header from IRC+10216_95125.wt.fits
### Warning [ccdfits]: Replacing header from IRC+10216_95125.fits
[0;31mQAC_STATS: IRC+10216_95125-full 0.00803632 0.529516 -92.8787 90.2376 97479.3 0.0326366[0m
### Warning [ccdsub]: Sampling will be done
### Warning [ccdsub]: Attempting to fix the WCS
### nemo Debug Info: WCS Corner: -0.0825 0.0825 323485
[0;31mQAC_STATS: IRC+10216_95125-cent 0.0113636 0.285189 -2.46959 17.2206 94737.1 0.0590684[0m
### nemo Debug Info: FITS file: Image size 219 219 689
### Warning [fitsccd]: Setting map minimum from 0 to -0.0768285
### Warning [fitsccd]: Setting map maximum from 0 to 5.14993
### nemo Debug Info: There were no blank values set in the image
spectab : # 109 109 689 min/mean/sig/max -0.641189 0.576443 2.9307 17.2206 -0.00011001 0.579528 0.139826
specstab: # 109 109 689 min/mean/sig/max -0.0605518 0.183715 0.916832 5.14993 -1.0606e-05 0.098631 0.0760692
Traceback (most recent call last):
File "/home/lmtslr/lmtoy/examples/fitsplot.py", line 21, in
f = aplpy.FITSFigure(fitsfile)
File "/home/lmtslr/lmtoy/anaconda3/lib/python3.8/site-packages/aplpy/decorators.py", line 26, in wrapper
return f(*args, **kwargs)
File "/home/lmtslr/lmtoy/anaconda3/lib/python3.8/site-packages/aplpy/core.py", line 239, in __init__
self._get_hdu(data, hdu, north, convention=convention,
File "/home/lmtslr/lmtoy/anaconda3/lib/python3.8/site-packages/aplpy/core.py", line 316, in _get_hdu
raise IOError("File not found: " + filename)
OSError: File not found: IRC+10216_95125.mom0.fits
Traceback (most recent call last):
File "/home/lmtslr/lmtoy/examples/fitsplot.py", line 21, in
f = aplpy.FITSFigure(fitsfile)
File "/home/lmtslr/lmtoy/anaconda3/lib/python3.8/site-packages/aplpy/decorators.py", line 26, in wrapper
return f(*args, **kwargs)
File "/home/lmtslr/lmtoy/anaconda3/lib/python3.8/site-packages/aplpy/core.py", line 239, in __init__
self._get_hdu(data, hdu, north, convention=convention,
File "/home/lmtslr/lmtoy/anaconda3/lib/python3.8/site-packages/aplpy/core.py", line 316, in _get_hdu
raise IOError("File not found: " + filename)
OSError: File not found: IRC+10216_95125.mom1.fits
Traceback (most recent call last):
File "/home/lmtslr/lmtoy/examples/fitsplot.py", line 21, in
f = aplpy.FITSFigure(fitsfile)
File "/home/lmtslr/lmtoy/anaconda3/lib/python3.8/site-packages/aplpy/decorators.py", line 26, in wrapper
return f(*args, **kwargs)
File "/home/lmtslr/lmtoy/anaconda3/lib/python3.8/site-packages/aplpy/core.py", line 239, in __init__
self._get_hdu(data, hdu, north, convention=convention,
File "/home/lmtslr/lmtoy/anaconda3/lib/python3.8/site-packages/aplpy/core.py", line 316, in _get_hdu
raise IOError("File not found: " + filename)
OSError: File not found: IRC+10216_95125.mom2.fits
INFO: Auto-setting vmin to -5.241e+00 [aplpy.core]
INFO: Auto-setting vmax to 4.910e+01 [aplpy.core]
INFO: Auto-setting vmin to -5.421e+00 [aplpy.core]
INFO: Auto-setting vmax to 5.089e+01 [aplpy.core]
Writing IRC+10216_95125.wt.png
INFO: Auto-setting vmin to -5.434e+00 [aplpy.core]
INFO: Auto-setting vmax to 6.032e+01 [aplpy.core]
INFO: Auto-setting vmin to -5.796e+00 [aplpy.core]
INFO: Auto-setting vmax to 6.433e+01 [aplpy.core]
Writing IRC+10216_95125.wt2.png
INFO: Auto-setting vmin to -2.681e+03 [aplpy.core]
INFO: Auto-setting vmax to 2.976e+04 [aplpy.core]
INFO: Auto-setting vmin to -3.321e+03 [aplpy.core]
INFO: Auto-setting vmax to 3.686e+04 [aplpy.core]
Writing IRC+10216_95125.wt3.png
Traceback (most recent call last):
File "/home/lmtslr/lmtoy/examples/fitsplot.py", line 21, in
f = aplpy.FITSFigure(fitsfile)
File "/home/lmtslr/lmtoy/anaconda3/lib/python3.8/site-packages/aplpy/decorators.py", line 26, in wrapper
return f(*args, **kwargs)
File "/home/lmtslr/lmtoy/anaconda3/lib/python3.8/site-packages/aplpy/core.py", line 239, in __init__
self._get_hdu(data, hdu, north, convention=convention,
File "/home/lmtslr/lmtoy/anaconda3/lib/python3.8/site-packages/aplpy/core.py", line 316, in _get_hdu
raise IOError("File not found: " + filename)
OSError: File not found: IRC+10216_95125.wtn.fits
INFO: Auto-setting vmin to -1.120e-01 [aplpy.core]
INFO: Auto-setting vmax to 1.244e+00 [aplpy.core]
INFO: Auto-setting vmin to -1.127e-01 [aplpy.core]
INFO: Auto-setting vmax to 1.251e+00 [aplpy.core]
Writing IRC+10216_95125.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 -9.219e+03 [aplpy.core]
INFO: Auto-setting vmax to 4.086e+04 [aplpy.core]
INFO: Auto-setting vmin to -1.005e+04 [aplpy.core]
INFO: Auto-setting vmax to 5.175e+04 [aplpy.core]
Writing IRC+10216_95125.mom0.png
LMTOY>> Created IRC+10216_95125.nf.fits and IRC+10216_95125.nfs.fits
LMTOY>> ADMIT post-processing
LMTOY>> lmtoy_admit: 31-oct-2021
/home/lmtslr/lmtoy/examples/lmtoy_admit.sh: line 33: export: `IRC+10216_95125.nfs.fits': not a valid identifier
No PDIR directory used, all work in the current directory
/bin/csh: /home/lmtslr/lmtoy/admit/casa/lib/libtinfo.so.5: no version information available (required by /bin/csh)
ADMIT = /home/lmtslr/lmtoy/admit
version = 1.0.8.5
CASA_PATH = /home/lmtslr/lmtoy/admit/casa linux admit malt
CASA_ROOT = /home/lmtslr/lmtoy/admit/casa
prefix = /home/lmtslr/lmtoy/admit/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
/bin/csh: /home/lmtslr/lmtoy/admit/casa/lib/libtinfo.so.5: no version information available (required by /bin/csh)
### (1/1) : processing IRC+10216_95125.nfs.fits
IRC+10216_95125.nfs.fits [219, 219, 689] IRC+10216 0.0 0.0 0.0003234848 0.123
/home/lmtslr/lmtoy/admit/admit/test/admit1.py --basename x IRC+10216_95125.nfs.fits
298.642u 9.735s 18:38.15 27.5% 0+0k 38232+1260120io 276pf+0w
Logfile: IRC+10216_95125.nfs.fits.log
LMTOY>> lmtoy_admit: 31-oct-2021
/home/lmtslr/lmtoy/examples/lmtoy_admit.sh: line 33: export: `IRC+10216_95125.nf.fits': not a valid identifier
No PDIR directory used, all work in the current directory
/bin/csh: /home/lmtslr/lmtoy/admit/casa/lib/libtinfo.so.5: no version information available (required by /bin/csh)
ADMIT = /home/lmtslr/lmtoy/admit
version = 1.0.8.5
CASA_PATH = /home/lmtslr/lmtoy/admit/casa linux admit malt
CASA_ROOT = /home/lmtslr/lmtoy/admit/casa
prefix = /home/lmtslr/lmtoy/admit/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
/bin/csh: /home/lmtslr/lmtoy/admit/casa/lib/libtinfo.so.5: no version information available (required by /bin/csh)
### (1/1) : processing IRC+10216_95125.nf.fits
IRC+10216_95125.nf.fits [219, 219, 689] IRC+10216 0.0 0.0 0.0003234848 0.123
/home/lmtslr/lmtoy/admit/admit/test/admit1.py --basename x IRC+10216_95125.nf.fits
126.941u 4.697s 10:10.56 21.5% 0+0k 8+644808io 0pf+0w
Logfile: IRC+10216_95125.nf.fits.log
LMTOY>> Created IRC+10216_95125.fits and IRC+10216_95125.wt.fits
LMTOY>> Parameter file used: lmtoy_95125.rc
LMTOY>> Making summary index.html:
Making index.html for obsnum=95125
Writing index.html
Writing index_pipeline.html
Writing index_admit.html
Writing index_pars.html
Writing index_log.html