LMTOY>> seq_pipeline: 2-feb-2022
LMTOY>> READING lmtoy_functions 11-mar-2022 via /lma1/teuben/lmtoy/examples/seq_pipeline.sh
Working directory ./2018S1SEQUOIACommissioning/79448
LMTOY>> Hang on, creating a bootstrap lmtoy_79448.rc from path=/n/lma1/lmt/data_lmt/
#
# version=2-feb-2022
# ifproc="/n/lma1/lmt/data_lmt//ifproc/ifproc_2018-11-16_079448_01_0000.nc"
date_obs="2018-11-16T06:48:52"
# skytime=685.807 sec
# inttime=528.4 sec
obsnum="79448"
calobsnum="79447"
obspgm="Map"
obsgoal="Pointing"
ProjectId="2018S1SEQUOIACommissioning"
# SkyOff=900 0
# HPBW=16 arcsec
# XYLength=400 400 arcsec
# XYRamp=0 0 arcsec
# XYoff=0 0 arcsec
# XYstep=1 0.7
numbands=1
vlsr=-20 # km/s
skyfreq=115.2877722202265 # GHz
restfreq=115.271204 # 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_79448.rc file:
# seq_pipeline: 2-feb-2022
# DATE: 2022-03-27T22:37:49.873474349
# obsnum=79448
# Using ifproc=/n/lma1/lmt/data_lmt//ifproc/ifproc_2018-11-16_079448_01_0000.nc
# path=/n/lma1/lmt/data_lmt/
#
# version=2-feb-2022
# ifproc="/n/lma1/lmt/data_lmt//ifproc/ifproc_2018-11-16_079448_01_0000.nc"
date_obs="2018-11-16T06:48:52"
# skytime=685.807 sec
# inttime=528.4 sec
obsnum="79448"
calobsnum="79447"
obspgm="Map"
obsgoal="Pointing"
ProjectId="2018S1SEQUOIACommissioning"
# SkyOff=900 0
# HPBW=16 arcsec
# XYLength=400 400 arcsec
# XYRamp=0 0 arcsec
# XYoff=0 0 arcsec
# XYstep=1 0.7
numbands=1
vlsr=-20 # km/s
skyfreq=115.2877722202265 # GHz
restfreq=115.271204 # 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_79448.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 537, 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 668, 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,-300],[-280,-120],[80,330]]',
'--bank': '0',
'--eliminate_list': '0',
'--help': False,
'--l_regions': '[[-120,80]]',
'--obsnum': '79448',
'--output': 'IRC+10216_79448.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
stype = 2
pix_list = [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15]
x_axis = VLSR
b_regions = [[-370, -300], [-280, -120], [80, 330]]
path = /n/lma1/lmt/data_lmt/
data_path = /n/lma1/lmt/data_lmt/
eliminate_list = [0]
slice = [-370, 330]
output = IRC+10216_79448.nc
use_otf_cal = False
bank = 0
use_cal = True
tsys = 250.0
b_order = 0
output_file_name = IRC+10216_79448.nc
obsnum = 79448
l_regions = [[-120, 80]]
2018-11-16T06:48:52 begin /n/lma1/lmt/data_lmt/ifproc/ifproc_2018-11-16_079448_01_0000.nc
2018-11-16T07:00:18 end /n/lma1/lmt/data_lmt/ifproc/ifproc_2018-11-16_079448_01_0000.nc
BufPos 1 16.9 sec 79448
BufPos 0 150.6 sec 79448
BufPos 1 17.7 sec 79448
BufPos 0 151.6 sec 79448
BufPos 1 17.3 sec 79448
BufPos 0 150.9 sec 79448
BufPos 1 17.4 sec 79448
BufPos 0 153.3 sec 79448
BufPos 1 10.1 sec 79448
before read npix
from pixels npix = 16
from xlen npix = 16
TRACKING Sequoia PIXEL 10
Map Parameters: Ra Continuous
HPBW= 16.0 XLength= 400.0 YLength= 400.0 XStep= 1.00 YStep= 0.70
/n/lma1/lmt/data_lmt/ifproc/ifproc_2018-11-16_079448_01_0000.nc does not have bs parameters
79448 is a Map observation
get chop
'Data.Msip1mm.BeamChopperActPos'
no chop
2018-11-16T06:48:30 begin /n/lma1/lmt/data_lmt/ifproc/ifproc_2018-11-16_079447_00_0001.nc
2018-11-16T06:48:38 end /n/lma1/lmt/data_lmt/ifproc/ifproc_2018-11-16_079447_00_0001.nc
BufPos 3 5.1 sec 79447
BufPos 2 2.6 sec 79447
before read npix
from pixels npix = 16
from xlen npix = 16
TRACKING Sequoia PIXEL 10
/n/lma1/lmt/data_lmt/ifproc/ifproc_2018-11-16_079447_00_0001.nc does not have map parameters
/n/lma1/lmt/data_lmt/ifproc/ifproc_2018-11-16_079447_00_0001.nc does not have bs parameters
79447 is a Cal observation
get chop cal
'Data.Msip1mm.BeamChopperActPos'
no chop cal
read_roach /n/lma1/lmt/data_lmt/spectrometer/roach0/roach0_79448_1_0_IRC+10216_2018-11-16_114845.nc nspec,nchan=25108,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_79448_1_0_IRC+10216_2018-11-16_114845.nc nspec,nchan=25112,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_79448_1_0_IRC+10216_2018-11-16_114845.nc nspec,nchan=25140,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_79448_1_0_IRC+10216_2018-11-16_114845.nc nspec,nchan=25144,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_79447_0_1_IRC+10216_2018-11-16_114830.nc nspec,nchan=236,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_79447_0_1_IRC+10216_2018-11-16_114830.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_79447_0_1_IRC+10216_2018-11-16_114830.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_79447_0_1_IRC+10216_2018-11-16_114830.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] = 228.711 +/- 14.6717 (2048 channels)
TSYS[1] = 186.476 +/- 10.7159 (2048 channels)
TSYS[2] = 237.321 +/- 15.4345 (2048 channels)
TSYS[3] = 176.894 +/- 9.99853 (2048 channels)
TSYS[4] = 196.23 +/- 14.0989 (2048 channels)
TSYS[5] = 188.974 +/- 11.6229 (2048 channels)
TSYS[6] = 244.756 +/- 20.2181 (2048 channels)
TSYS[7] = 191.019 +/- 8.66106 (2048 channels)
TSYS[8] = 180.528 +/- 12.6098 (2048 channels)
TSYS[9] = 190.449 +/- 11.9417 (2048 channels)
TSYS[10] = 187.345 +/- 12.224 (2048 channels)
TSYS[11] = 189.03 +/- 10.584 (2048 channels)
TSYS[12] = 169.13 +/- 9.37749 (2048 channels)
TSYS[13] = 205.447 +/- 11.5503 (2048 channels)
TSYS[14] = 196.334 +/- 11.6919 (2048 channels)
TSYS[15] = 198.554 +/- 14.0957 (2048 channels)
Spectral Band velocity range: -1059.8 1019.8 km/s
Total Number of OTF Spectra = 91816
NSPEC: 91816
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: 228.221 (5.52475) minmax: 218.315 239.078
0 5739 0.029 1.253 1.219 -6.867 16.610 1.224 1.964 *M 3.7
TSYS[1] slice: 186.218 (3.82251) minmax: 178.93 193.181
1 5739 0.082 1.247 1.032 -8.401 18.525 1.075 2.903 *P 1.2 *M 11.7
TSYS[2] slice: 236.981 (5.89269) minmax: 225.707 247.79
2 5739 -0.020 1.457 1.297 -11.996 42.242 1.366 3.786 *M 11.7
TSYS[3] slice: 176.677 (3.8792) minmax: 169.197 183.713
3 5739 0.001 1.001 0.956 -5.648 16.210 0.964 1.920 *M 5.6
TSYS[4] slice: 197.745 (6.24488) minmax: 185.273 207.345
4 5738 0.023 1.113 1.060 -6.573 23.085 1.081 2.182 *M 5.9
TSYS[5] slice: 191.114 (4.29401) minmax: 181.14 197.63
5 5738 0.054 1.087 1.034 -6.455 16.790 1.026 1.790 *M 4.2
TSYS[6] slice: 247.739 (6.13411) minmax: 232.428 257.318
6 5738 -0.027 1.427 1.336 -9.379 37.046 1.357 2.708 *M 6.4
TSYS[7] slice: 191.52 (3.30968) minmax: 185.798 198.447
7 5738 0.003 1.062 1.027 -6.529 15.740 1.034 2.084 *M 5.1
TSYS[8] slice: 179.962 (5.18863) minmax: 171.029 189.979
8 5739 0.029 1.014 0.966 -5.757 16.611 0.976 1.906 *M 4.9
TSYS[9] slice: 190.915 (5.61381) minmax: 181.127 200.316
9 5739 0.047 1.138 1.042 -8.260 17.819 1.059 2.476 *M 7.6
TSYS[10] slice: 188.455 (5.97086) minmax: 177.665 198.259
10 5739 0.014 1.083 1.019 -7.303 17.356 1.044 2.370 *M 7.6
TSYS[11] slice: 190.503 (3.27864) minmax: 183.323 196.734
11 5739 0.009 1.069 1.025 -6.606 17.375 1.041 2.273 *M 6.2
TSYS[12] slice: 169.003 (3.48854) minmax: 162.366 175.419
12 5738 0.029 0.933 0.901 -4.684 16.786 0.899 1.458 *M 3.0
TSYS[13] slice: 205.027 (3.89592) minmax: 196.987 212.257
13 5738 0.036 1.208 1.109 -7.209 38.975 1.159 2.832 *M 9.4
TSYS[14] slice: 195.495 (3.92732) minmax: 187.8 203.465
14 5738 0.016 1.147 1.063 -8.136 16.101 1.091 2.747 *M 8.4
TSYS[15] slice: 197.349 (5.43977) minmax: 187.91 207.655
15 5738 0.003 1.109 1.063 -6.845 16.836 1.074 2.069 *M 5.3
CPU TIME: 48.2921 sec
CPU TIME: 48.4461 sec
Warnings:
*P ratio: ratio of std/mad too high for data
*M ratio: ratio of std/mad too high for RMS
X-range: -244.574 269.081 Y-range: -256.818 267.553 arcsec
Written IRC+10216_79448.nc
LMTOY>> view_spec_file
{'--help': False,
'--input': 'IRC+10216_79448.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_79448_specviews',
'--rms_cut': '-4',
'--show_pixel': None,
'--skip_tsys': False}
SLR Plots initialized with 1 method=IRC+10216_79448_specviews
Max sequence=5739
Plots saving in IRC+10216_79448_specviews.1.png
Plots saving in IRC+10216_79448_specviews.2.png
Plots saving in IRC+10216_79448_specviews.3.png
Plots saving in IRC+10216_79448_specviews.4.png
Plots saving in IRC+10216_79448_specviews.5.png
NEW TSYS PIC
SHAPE (1, 16, 689)
PIXEL: 0 [ 0 1 2 ... 5736 5737 5738]
PIXEL: 1 [ 5739 5740 5741 ... 11475 11476 11477]
PIXEL: 2 [11478 11479 11480 ... 17214 17215 17216]
PIXEL: 3 [17217 17218 17219 ... 22953 22954 22955]
PIXEL: 4 [22956 22957 22958 ... 28691 28692 28693]
PIXEL: 5 [28694 28695 28696 ... 34429 34430 34431]
PIXEL: 6 [34432 34433 34434 ... 40167 40168 40169]
PIXEL: 7 [40170 40171 40172 ... 45905 45906 45907]
PIXEL: 8 [45908 45909 45910 ... 51644 51645 51646]
PIXEL: 9 [51647 51648 51649 ... 57383 57384 57385]
PIXEL: 10 [57386 57387 57388 ... 63122 63123 63124]
PIXEL: 11 [63125 63126 63127 ... 68861 68862 68863]
PIXEL: 12 [68864 68865 68866 ... 74599 74600 74601]
PIXEL: 13 [74602 74603 74604 ... 80337 80338 80339]
PIXEL: 14 [80340 80341 80342 ... 86075 86076 86077]
PIXEL: 15 [86078 86079 86080 ... 91813 91814 91815]
Plots saving in IRC+10216_79448_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_79448_specpoint,png,1
PJT -1 [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15] IRC+10216_79448_specpoint,png,1 False
MAD: 0 1.1841571 0.04290821983813239 1.3557900124549465
MAD: 1 0.96640736 0.03445218495805559 1.1042160984786151
MAD: 2 1.2296911 0.043396729079707326 1.4032780641230895
MAD: 3 0.9175047 0.032779518003545734 1.0486227998546065
MAD: 4 1.0242816 0.03721825365750154 1.1731546356093152
MAD: 5 0.990745 0.03560620850829566 1.1331698420249734
MAD: 6 1.2787111 0.04678076903934387 1.465834156708523
MAD: 7 0.9937161 0.03558362997876918 1.1360506406349864
MAD: 8 0.9360806 0.03387892890702508 1.0715962906174192
MAD: 9 0.9910272 0.03575846998138424 1.1340610563056943
MAD: 10 0.9773009 0.03449902104670152 1.1152969665262837
MAD: 11 0.98971033 0.035736465856718504 1.1326561943900089
MAD: 12 0.8786584 0.03213136839567102 1.0071838874697079
MAD: 13 1.0653353 0.03747311267571828 1.215227724445549
MAD: 14 1.0152559 0.036810382021378325 1.162497456125064
MAD: 15 1.0261832 0.036539130372135506 1.1723397690547652
Plots saving in IRC+10216_79448_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)
SLR Plots initialized with 3 method=IRC+10216_79448_specpoint,png,2
PJT -1 [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15] IRC+10216_79448_specpoint,png,2 False
MAD: 0 1.1841571 0.04290821983813239 1.3557900124549465
MAD: 1 0.96640736 0.03445218495805559 1.1042160984786151
MAD: 2 1.2296911 0.043396729079707326 1.4032780641230895
MAD: 3 0.9175047 0.032779518003545734 1.0486227998546065
MAD: 4 1.0242816 0.03721825365750154 1.1731546356093152
MAD: 5 0.990745 0.03560620850829566 1.1331698420249734
MAD: 6 1.2787111 0.04678076903934387 1.465834156708523
MAD: 7 0.9937161 0.03558362997876918 1.1360506406349864
MAD: 8 0.9360806 0.03387892890702508 1.0715962906174192
MAD: 9 0.9910272 0.03575846998138424 1.1340610563056943
MAD: 10 0.9773009 0.03449902104670152 1.1152969665262837
MAD: 11 0.98971033 0.035736465856718504 1.1326561943900089
MAD: 12 0.8786584 0.03213136839567102 1.0071838874697079
MAD: 13 1.0653353 0.03747311267571828 1.215227724445549
MAD: 14 1.0152559 0.036810382021378325 1.162497456125064
MAD: 15 1.0261832 0.036539130372135506 1.1723397690547652
Plots saving in IRC+10216_79448_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 5738
Binning [1, 1]
FITS file will be 5738 x 689 x 16
Written waterfall cube to IRC+10216_79448.wf.fits
pix [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15]
Min nsamp 5738
Binning [10, 1]
FITS file will be 573 x 689 x 16
New shape: (16, 689, 573)
Written waterfall cube to IRC+10216_79448.wf10.fits
LMTOY>> grid_data
{'--cell': '6.25',
'--edge': '0',
'--help': False,
'--input': 'IRC+10216_79448.nc',
'--n_samples': '256',
'--noise_sigma': '1',
'--otf_a': '1.1',
'--otf_b': '4.75',
'--otf_c': '2',
'--otf_select': '1',
'--output': 'IRC+10216_79448.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': '12.5',
'--rmax': '3',
'--rms_cut': '-4',
'--sample': '-1',
'--weight': 'IRC+10216_79448.wt.fits',
'--x_extent': '300',
'--y_extent': '300'}
STDOUT ****************************
spec_driver_fits LMTSLR 4-mar-2022
CMD: spec_driver_fits -i IRC+10216_79448.nc -o IRC+10216_79448.fits -w IRC+10216_79448.wt.fits -l 12.5 -c 6.25 -u [0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15] -z -4 -s 1 -x 300 -y 300 -f 1 -r 3 -n 256 -0 1.1 -1 4.75 -2 2 -b -1 -e 0
file list=IRC+10216_79448.nc
DEBUG:6.25
sample_list=-1
Warning: sample list is not a multiple of 3
nsamples=1
Processing 1 SpecFiles:
Opening SpecFile file IRC+10216_79448.nc
file: IRC+10216_79448.nc nspec= 91816 nchan= 689
SpecFile version 6-mar-2021
X-range: -244.574 269.081 Y-range: -256.818 267.553 arcsec
X-ramp: -161.174 185.681 Y-ramp: -173.418 184.153 arcsec
MapSize: 513.655 x 524.372 arcsec
79448
IRC+10216
DATE-OBS 2018-11-16T06:48:52
CF.n_cells= 7 delta=0.147059 cell=6.25 resol=12.5 oft_select=1
C-AXIS 2 starts at 329.985
C-AXIS 1 starts at -300
C-AXIS 0 starts at -300
P-AXIS 1 starts at -300
P-AXIS 0 starts at -300
P-AXIS 1 starts at -300
P-AXIS 0 starts at -300
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=91816 rms_cutoff=-4
Cube Completed, 85003/91816 Spectra accepted = 0.926
Weighting Completed, fuzzy_edge=0
Weight at 0.000000 0.000000 is 21.385635
write to IRC+10216_79448.fits
write weights to IRC+10216_79448.wt.fits
STDERR ****************************
Exit Code: 0
LMTOY>> Some NEMO post-processing
### nemo Debug Info: FITS file: Image size 97 97 689
### Warning [fitsccd]: Setting map minimum from 0 to -55.9662
### Warning [fitsccd]: Setting map maximum from 0 to 45.8299
### nemo Debug Info: There were 1814826 blank values in the image (27.994473 %)
### nemo Debug Info: FITS file: Image size 97 97
### Warning [fitsccd]: Setting map minimum from 0 to -2.97536
### Warning [fitsccd]: Setting map maximum from 0 to 62.107
### nemo Debug Info: There were no blank values set in the image
LMTOY>> STATS IRC+10216_79448.ccd centerbox robust
### Warning [ccdsub]: Sampling will be done
### Warning [ccdsub]: Attempting to fix the WCS
### nemo Debug Info: WCS Corner: 146.949 13.3187 329985
Number of points : 1722500
Min and Max : -0.947408 15.263989
Mean and dispersion : 0.039418 0.373722
Skewness and kurtosis : 10.928029 248.488978
Sum and Sum*Dx*Dy*Dz* : 67897.132908 207.905684
Points per beam (map) : 1
N Robust : 1672804
Mean Robust : 0.002406
Sigma Robust : 0.183928
Median Robust : 0.001407
0/1722500 out-of-range points discarded
LMTOY>> STATS IRC+10216_79448.wt.ccd centerbox robust
### Warning [ccdsub]: Sampling will be done
### Warning [ccdsub]: Attempting to fix the WCS
### nemo Debug Info: WCS Corner: 146.949 13.3187 0
Number of points : 2500
Min and Max : 18.639803 59.736450
Mean and dispersion : 24.342390 3.779290
Skewness and kurtosis : 5.772878 39.494164
Sum and Sum*Dx*Dy* : 60855.973917 -0.183425
Points per beam (map) : 1
N Robust : 2310
Mean Robust : 23.735237
Sigma Robust : 1.051138
Median Robust : 23.753420
0/2500 out-of-range points discarded
### nemo Debug Info: Keeping 97*97*689 cube
### Warning [ccdfits]: Replacing header from IRC+10216_79448.wt.fits
SLABS: [[-370,-300],[-280,-120],[80,330]] == -370,-300,-280,-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 -300 -> 689 620
### nemo Debug Info: slab -280 -120 -> 600 442
### nemo Debug Info: slab 80 330 -> 246 0
### nemo Debug Info: Total planes: 476
### Warning [ccdslice]: Code not converted to fix reference pixel value
### Warning [ccdfits]: Replacing header from IRC+10216_79448.wt.fits
### Warning [ccdfits]: Replacing header from IRC+10216_79448.wt.fits
### Warning [ccdfits]: Replacing header from IRC+10216_79448.fits
[0;31mQAC_STATS: IRC+10216_79448-full 0.0154308 0.436052 -55.9662 45.8299 72030.6 0.0721877[0m
### Warning [ccdsub]: Sampling will be done
### Warning [ccdsub]: Attempting to fix the WCS
### nemo Debug Info: WCS Corner: 146.949 13.3187 329985
[0;31mQAC_STATS: IRC+10216_79448-cent 0.0394178 0.373722 -0.947408 15.264 67897.1 0.213951[0m
### nemo Debug Info: FITS file: Image size 97 97 689
### Warning [fitsccd]: Setting map minimum from 0 to -0.0714179
### Warning [fitsccd]: Setting map maximum from 0 to 6.07828
### nemo Debug Info: There were no blank values set in the image
spectab : # 48 48 689 min/mean/sig/max -0.601834 0.544644 2.67454 15.264 0.000488855 0.473258 0.125122
specstab: # 48 48 689 min/mean/sig/max -0.0560382 0.219458 1.06654 6.07828 -5.92235e-07 0.111252 0.0737588
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_79448.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_79448.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_79448.mom2.fits
INFO: Auto-setting vmin to -6.893e+00 [aplpy.core]
INFO: Auto-setting vmax to 5.566e+01 [aplpy.core]
INFO: Auto-setting vmin to -6.893e+00 [aplpy.core]
INFO: Auto-setting vmax to 5.566e+01 [aplpy.core]
Writing IRC+10216_79448.wt.png
INFO: Auto-setting vmin to -6.202e+00 [aplpy.core]
INFO: Auto-setting vmax to 6.884e+01 [aplpy.core]
INFO: Auto-setting vmin to -6.202e+00 [aplpy.core]
INFO: Auto-setting vmax to 6.884e+01 [aplpy.core]
Writing IRC+10216_79448.wt2.png
INFO: Auto-setting vmin to -8.536e+04 [aplpy.core]
INFO: Auto-setting vmax to 9.475e+05 [aplpy.core]
INFO: Auto-setting vmin to -8.536e+04 [aplpy.core]
INFO: Auto-setting vmax to 9.475e+05 [aplpy.core]
Writing IRC+10216_79448.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_79448.wtn.fits
INFO: Auto-setting vmin to -1.158e-01 [aplpy.core]
INFO: Auto-setting vmax to 1.286e+00 [aplpy.core]
INFO: Auto-setting vmin to -1.158e-01 [aplpy.core]
INFO: Auto-setting vmax to 1.286e+00 [aplpy.core]
Writing IRC+10216_79448.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.808e+04 [aplpy.core]
INFO: Auto-setting vmax to 1.296e+05 [aplpy.core]
INFO: Auto-setting vmin to -1.808e+04 [aplpy.core]
INFO: Auto-setting vmax to 1.296e+05 [aplpy.core]
Writing IRC+10216_79448.mom0.png
LMTOY>> Created IRC+10216_79448.nf.fits and IRC+10216_79448.nfs.fits
LMTOY>> ADMIT post-processing
LMTOY>> lmtoy_admit: 31-oct-2021
/lma1/teuben/lmtoy/examples/lmtoy_admit.sh: line 33: export: `IRC+10216_79448.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
### (1/1) : processing IRC+10216_79448.nfs.fits
IRC+10216_79448.nfs.fits [97, 97, 689] IRC+10216 146.9892 13.27877 0.0003299848 0.024
/lma1/teuben/lmtoy/admit/admit/test/admit1.py --basename x IRC+10216_79448.nfs.fits
161.450u 6.859s 5:35.05 50.2% 0+0k 163168+538040io 124pf+0w
Logfile: IRC+10216_79448.nfs.fits.log
LMTOY>> lmtoy_admit: 31-oct-2021
/lma1/teuben/lmtoy/examples/lmtoy_admit.sh: line 33: export: `IRC+10216_79448.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
### (1/1) : processing IRC+10216_79448.nf.fits
IRC+10216_79448.nf.fits [97, 97, 689] IRC+10216 146.9892 13.27877 0.0003299848 0.024
/lma1/teuben/lmtoy/admit/admit/test/admit1.py --basename x IRC+10216_79448.nf.fits
61.094u 2.870s 2:26.79 43.5% 0+0k 552+183280io 0pf+0w
Logfile: IRC+10216_79448.nf.fits.log
LMTOY>> Created IRC+10216_79448.fits and IRC+10216_79448.wt.fits
LMTOY>> Parameter file used: lmtoy_79448.rc
LMTOY>> Making summary index.html:
Making index.html for obsnum=79448
IRC+10216_79448
IRC+10216_79448_specviews
IRC+10216_79448_specpoint
Writing index.html
Writing index_pipeline.html
Writing index_admit.html
Writing index_pars.html
Writing index_log.html