LMTOY>> seq_pipeline: 11-jul-2023
LMTOY>> lmtoy_functions 11-jul-2023 via /home/teuben/LMT/lmtoy/bin/seq_pipeline.sh
Working directory /home/teuben/LMT/lmtoy/work_lmt/2018S1SEQUOIACommissioning/79448
LMTOY>> creating new ./lmtoy_79448.rc
LMTOY>> this is your startup ./lmtoy_79448.rc file:
# seq_pipeline: 11-jul-2023
extent=0
dv=100
dw=800
birdies=0
numbands=0
map_coord_use=-1
pix_list=0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15
rms_cut=-4
location=0,0
resolution=12.5
cell=6.25
nppb=-1
rmax=3
otf_select=1
otf_a=1.1
otf_b=4.75
otf_c=2
noise_sigma=1
b_order=0
stype=2
sample=-1
otf_cal=0
edge=0
bank=-1
#
# version=21-jun-2023
rawnc="/home/teuben/LMT/lmtoy/data_lmt/ifproc/ifproc_2018-11-16_079448_01_0000.nc"
date_obs="2018-11-16T06:48:52"
skytime=685.807
inttime=528.4
obsnum=79448
calobsnum=79447
obspgm="Map"
map_coord="Ra"
map_motion="Continuous"
ProjectId="2018S1SEQUOIACommissioning"
ra=146.989192 # deg
dec=13.278767 # deg
# SkyOff=900 0
# HPBW=16 arcsec
# XYLength=400 400 arcsec
# XYRamp=0 0 arcsec
# XYoff=0 0 arcsec
# XYstep=1 0.7
nrows=36
rowtint=12
tint_on=432 # excluding 1.6sec (?) switchback
numbands=1
vlsr=-20 # km/s
skyfreq=115.2877722202265 # GHz
restfreq=115.271204 # Ghz
src="IRC+10216"
resolution=11 # arcsec
nppb=2 # number of points per beam
cell=5.5 # arcsec
x_extent=400 # arcsec
y_extent=400 # arcsec
instrument="SEQ"
tau=0.062
PIName="LMT"
projectTitle="Commissioning"
#
#
date="2023-07-12T12:05:18" # begin
pdir=/home/teuben/LMT/lmtoy/work_lmt/2018S1SEQUOIACommissioning/79448
obsnum=79448
dv=100
dw=800
restart=1
tap=1
# Using ifproc=/home/teuben/LMT/lmtoy/data_lmt/ifproc/ifproc_2018-11-16_079448_01_0000.nc
# data_lmt=/home/teuben/LMT/lmtoy/data_lmt
vlsr=-20 # set
restfreq=115.271204 # set
# based on vlsr=-20, dv=100, dw=800
b_regions=[[-920,-120],[80,880]]
l_regions=[[-120,80]]
slice=[-920,880]
v_range=-120,80
#xy_extent should be same
pix_list=0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15
# resolution hack?
resolution=11
cell=5.5
map_coord_use=-1
LMTOY>> Sleeping for 2 seconds, you can abort, edit ./lmtoy_79448.rc, then continuing
LMTOY>> numbands=1
LMTOY>> rc1=./lmtoy_79448.rc
LMTOY>> process_otf_map2 in 2 seconds bank=-1
WARNING: resetting restfreq not supported yet
{'-': None,
'--b_order': '0',
'--b_regions': '[[-920,-120],[80,880]]',
'--bank': '-1',
'--eliminate_list': '0',
'--help': False,
'--l_regions': '[[-120,80]]',
'--map_coord': '-1',
'--obsnum': '79448',
'--output': 'IRC+10216_79448.nc',
'--path': '/home/teuben/LMT/lmtoy/data_lmt',
'--pix_list': '0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15',
'--restfreq': '-1',
'--sample': None,
'--save_tsys': False,
'--slice': '[-920,880]',
'--stype': '2',
'--tsys': '250.0',
'--use_cal': True,
'--use_otf_cal': False,
'--x_axis': 'VLSR'}
program process_otf_map options
pix_list = [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15]
eliminate_list = [0]
output = IRC+10216_79448.nc
obsnum = 79448
bank = -1
map_coord = -1
data_path = /home/teuben/LMT/lmtoy/data_lmt
l_regions = [[-120, 80]]
output_file_name = IRC+10216_79448.nc
b_regions = [[-920, -120], [80, 880]]
path = /home/teuben/LMT/lmtoy/data_lmt
restfreq = -1
x_axis = VLSR
b_order = 0
use_otf_cal = False
use_cal = True
slice = [-920, 880]
tsys = 250.0
stype = 2
PJT map_coord -1
2018-11-16T06:48:52 obs-start /home/teuben/LMT/lmtoy/data_lmt/ifproc/ifproc_2018-11-16_079448_01_0000.nc
2018-11-16T07:00:18 obs-stop /home/teuben/LMT/lmtoy/data_lmt/ifproc/ifproc_2018-11-16_079448_01_0000.nc
delta1 685.8 sec
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
delta2 685.8 sec
Old data, assuming Header.SpecBackend.DumpTime = 0.1
before read npix
from pixels npix = 16
from xlen npix = 16
TRACKING Sequoia PIXEL 10
PJT IFProc map_coord Ra
MapCoord: Ra Pointing 0 79448 -1
Warning: pointing data, assumed Az/El
Map Parameters: Ra Continuous
HPBW= 16.0 XLength= 400.0 YLength= 400.0 XStep= 1.00 YStep= 0.70 ScanAngle= 0.00
/home/teuben/LMT/lmtoy/data_lmt/ifproc/ifproc_2018-11-16_079448_01_0000.nc does not have bs parameters
PJT: line_rest_freq [115.271204 85.640456]
79448 is a Map observation
PJT IFProcData map_coord Ra
MapCoord: Ra Pointing 0 79448 -1
Warning: pointing data, assumed Az/El
get chop
'Data.Msip1mm.BeamChopperActPos'
no chop
2018-11-16T06:48:30 obs-start /home/teuben/LMT/lmtoy/data_lmt/ifproc/ifproc_2018-11-16_079447_00_0001.nc
2018-11-16T06:48:38 obs-stop /home/teuben/LMT/lmtoy/data_lmt/ifproc/ifproc_2018-11-16_079447_00_0001.nc
delta1 7.7 sec
BufPos 3 5.1 sec 79447
BufPos 2 2.6 sec 79447
delta2 7.7 sec
Old data, assuming Header.SpecBackend.DumpTime = 0.1
before read npix
from pixels npix = 16
from xlen npix = 16
TRACKING Sequoia PIXEL 10
'Header.Map.HPBW'
/home/teuben/LMT/lmtoy/data_lmt/ifproc/ifproc_2018-11-16_079447_00_0001.nc does not have map parameters
/home/teuben/LMT/lmtoy/data_lmt/ifproc/ifproc_2018-11-16_079447_00_0001.nc does not have bs parameters
PJT: line_rest_freq [115.271204 85.640456]
79447 is a Cal observation
get chop cal
'Data.Msip1mm.BeamChopperActPos'
no chop cal
deltaR 686.1
deltaS 2509.4
deltaI 475.4
read_roach /home/teuben/LMT/lmtoy/data_lmt/spectrometer/roach4/roach4_79448_1_0_IRC+10216_2018-11-16_114845.nc nspec,nchan=25108,2048
r:4 inp:0 pix:0 time_offset:-0.030000
r:4 inp:1 pix:1 time_offset:-0.030000
r:4 inp:2 pix:2 time_offset:-0.030000
r:4 inp:3 pix:3 time_offset:-0.030000
deltaR 686.1
deltaS 2509.8
deltaI 473.1
read_roach /home/teuben/LMT/lmtoy/data_lmt/spectrometer/roach5/roach5_79448_1_0_IRC+10216_2018-11-16_114845.nc nspec,nchan=25112,2048
r:5 inp:0 pix:4 time_offset:-0.030000
r:5 inp:1 pix:5 time_offset:-0.030000
r:5 inp:2 pix:6 time_offset:-0.030000
r:5 inp:3 pix:7 time_offset:-0.030000
deltaR 686.1
deltaS 2512.6
deltaI 477.0
read_roach /home/teuben/LMT/lmtoy/data_lmt/spectrometer/roach6/roach6_79448_1_0_IRC+10216_2018-11-16_114845.nc nspec,nchan=25140,2048
r:6 inp:0 pix:8 time_offset:-0.030000
r:6 inp:1 pix:9 time_offset:-0.030000
r:6 inp:2 pix:10 time_offset:-0.030000
r:6 inp:3 pix:11 time_offset:-0.030000
deltaR 686.2
deltaS 2513.0
deltaI 476.6
read_roach /home/teuben/LMT/lmtoy/data_lmt/spectrometer/roach7/roach7_79448_1_0_IRC+10216_2018-11-16_114845.nc nspec,nchan=25144,2048
r:7 inp:0 pix:12 time_offset:-0.030000
r:7 inp:1 pix:13 time_offset:-0.030000
r:7 inp:2 pix:14 time_offset:-0.030000
r:7 inp:3 pix:15 time_offset:-0.030000
PJT frequency_offset, restfreq 0.0 -1
RESTFREQ 85.6405 GHz
deltaR 7.9
deltaS 23.6
deltaI 4.5
read_roach /home/teuben/LMT/lmtoy/data_lmt/spectrometer/roach4/roach4_79447_0_1_IRC+10216_2018-11-16_114830.nc nspec,nchan=236,2048
r:4 inp:0 pix:0 time_offset:-0.030000
r:4 inp:1 pix:1 time_offset:-0.030000
r:4 inp:2 pix:2 time_offset:-0.030000
r:4 inp:3 pix:3 time_offset:-0.030000
deltaR 7.9
deltaS 23.6
deltaI 4.5
read_roach /home/teuben/LMT/lmtoy/data_lmt/spectrometer/roach5/roach5_79447_0_1_IRC+10216_2018-11-16_114830.nc nspec,nchan=236,2048
r:5 inp:0 pix:4 time_offset:-0.030000
r:5 inp:1 pix:5 time_offset:-0.030000
r:5 inp:2 pix:6 time_offset:-0.030000
r:5 inp:3 pix:7 time_offset:-0.030000
deltaR 8.0
deltaS 24.0
deltaI 4.5
read_roach /home/teuben/LMT/lmtoy/data_lmt/spectrometer/roach6/roach6_79447_0_1_IRC+10216_2018-11-16_114830.nc nspec,nchan=240,2048
r:6 inp:0 pix:8 time_offset:-0.030000
r:6 inp:1 pix:9 time_offset:-0.030000
r:6 inp:2 pix:10 time_offset:-0.030000
r:6 inp:3 pix:11 time_offset:-0.030000
deltaR 8.0
deltaS 24.4
deltaI 4.6
read_roach /home/teuben/LMT/lmtoy/data_lmt/spectrometer/roach7/roach7_79447_0_1_IRC+10216_2018-11-16_114830.nc nspec,nchan=244,2048
r:7 inp:0 pix:12 time_offset:-0.030000
r:7 inp:1 pix:13 time_offset:-0.030000
r:7 inp:2 pix:14 time_offset:-0.030000
r:7 inp:3 pix:15 time_offset:-0.030000
PJT frequency_offset, restfreq 0.0 -1
SPEC TSYS[16] = 228.711 +/- 14.6717 (2048 channels)
SPEC TSYS[17] = 186.476 +/- 10.7159 (2048 channels)
SPEC TSYS[18] = 237.321 +/- 15.4345 (2048 channels)
SPEC TSYS[19] = 176.894 +/- 9.99853 (2048 channels)
SPEC TSYS[20] = 196.23 +/- 14.0989 (2048 channels)
SPEC TSYS[21] = 188.974 +/- 11.6229 (2048 channels)
SPEC TSYS[22] = 244.756 +/- 20.2181 (2048 channels)
SPEC TSYS[23] = 191.019 +/- 8.66106 (2048 channels)
SPEC TSYS[24] = 180.528 +/- 12.6098 (2048 channels)
SPEC TSYS[25] = 190.449 +/- 11.9417 (2048 channels)
SPEC TSYS[26] = 187.345 +/- 12.224 (2048 channels)
SPEC TSYS[27] = 189.03 +/- 10.584 (2048 channels)
SPEC TSYS[28] = 169.13 +/- 9.37749 (2048 channels)
SPEC TSYS[29] = 205.447 +/- 11.5503 (2048 channels)
SPEC TSYS[30] = 196.334 +/- 11.6919 (2048 channels)
SPEC TSYS[31] = 198.554 +/- 14.0957 (2048 channels)
Warning: Line Rest Frequency of Cal observation does not match map 85.640456 115.271204
WARNING: CAL MAY NOT BE CORRECT
Velocity Range rawdata: -1419.55 1379.55 km/s = 2048 channels
Velocity Range selected: -920 880 km/s = 1317 channels
Total Number of OTF Spectra = 91816
PJT warning: assuming RA/DEC for map_coord = 0
NSPEC: 91816
NCHAN: 1317
ICHAN: 365
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
ho x -306.9130495417308 269.6288285844834
ho y -279.79445767441996 301.8275102472652
TSYS[0] slice: 228.58 (10.2469) minmax: 211.116 261.657
PJT-pos 269.0552299376975 -18.720269716913794
0 5739 0.019 1.299 1.230 -9.842 16.646 1.265 2.624 *M 10.2
ho x -306.9130495417308 269.6288285844834
ho y -279.79445767441996 301.8275102472652
TSYS[1] slice: 186.372 (7.53724) minmax: 172.67 210.768
1 5739 0.037 1.165 1.024 -8.725 18.193 1.064 2.602 *M 14.7
ho x -306.9130495417308 269.6288285844834
ho y -279.79445767441996 301.8275102472652
TSYS[2] slice: 237.356 (11.352) minmax: 218.179 368.738
2 5739 0.010 1.476 1.299 -28.506 43.437 1.409 3.675 *M 7.8
ho x -306.9130495417308 269.6288285844834
ho y -279.79445767441996 301.8275102472652
TSYS[3] slice: 176.799 (7.17279) minmax: 163.626 191.947
3 5739 0.002 0.986 0.951 -5.797 16.231 0.962 1.735 *M 7.0
ho x Traceback (most recent call last):
File "/home/teuben/LMT/lmtoy/SpectralLineReduction/lmtslr/ifproc/ifproc.py", line 593, in __init__
self.chop = self.nc.variables['Data.Msip1mm.BeamChopperActPos'][:]
KeyError: 'Data.Msip1mm.BeamChopperActPos'
Traceback (most recent call last):
File "/home/teuben/LMT/lmtoy/SpectralLineReduction/lmtslr/ifproc/ifproc.py", line 754, in __init__
self.chop = self.nc.variables['Data.Msip1mm.BeamChopperActPos'][:]
KeyError: 'Data.Msip1mm.BeamChopperActPos'
/home/teuben/LMT/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/teuben/LMT/lmtoy/SpectralLineReduction/lmtslr/spec/spec.py:580: RuntimeWarning: invalid value encountered in true_divide
self.tsys_spectrum = 280 * sky_spectrum / (hot_spectrum -
/home/teuben/LMT/lmtoy/SpectralLineReduction/lmtslr/spec/spec.py:580: RuntimeWarning: divide by zero encountered in true_divide
self.tsys_spectrum = 280 * sky_spectrum / (hot_spectrum -
/home/teuben/LMT/lmtoy/SpectralLineReduction/lmtslr/spec/spec.py:432: RuntimeWarning: invalid value encountered in true_divide
spectra.append((self.raw_spec[i,:] - ref) / ref)
/home/teuben/LMT/lmtoy/SpectralLineReduction/lmtslr/spec/spec.py:432: RuntimeWarning: divide by zero encountered in true_divide
spectra.append((self.raw_spec[i,:] - ref) / ref)
/home/teuben/LMT/lmtoy/SpectralLineReduction/lmtslr/spec/spec.py:467: RuntimeWarning: invalid value encountered in multiply
self.reduced_spectra = np.array(spectra) * tsys_spectrum
-306.94779776555174 269.5831435468327
ho y -279.6572070705403 301.80545810616337
TSYS[4] slice: 196.994 (10.7622) minmax: 177.841 274.381
4 5738 0.021 1.175 1.066 -15.592 24.071 1.126 2.668 *M 9.7
ho x -306.94779776555174 269.5831435468327
ho y -279.6572070705403 301.80545810616337
TSYS[5] slice: 189.769 (8.25446) minmax: 172.895 203.057
5 5738 0.023 1.063 1.021 -7.404 16.859 1.030 1.925 *M 7.4
ho x -306.94779776555174 269.5831435468327
ho y -279.6572070705403 301.80545810616337
TSYS[6] slice: 246.279 (13.8697) minmax: 218.036 387.956
6 5738 0.005 1.478 1.330 -31.492 38.656 1.432 3.448 *M 6.0
ho x -306.94779776555174 269.5831435468327
ho y -279.6572070705403 301.80545810616337
TSYS[7] slice: 191.478 (5.46347) minmax: 180.326 202.561
7 5738 0.005 1.047 1.022 -6.955 15.733 1.029 1.847 *M 6.2
ho x -306.8764817179339 269.60836040281544
ho y -279.6653090011365 302.0355612673752
TSYS[8] slice: 180.613 (9.60861) minmax: 164.726 198.782
8 5739 0.021 1.041 0.974 -8.997 16.594 1.001 2.203 *M 10.3
ho x -306.8764817179339 269.60836040281544
ho y -279.6653090011365 302.0355612673752
TSYS[9] slice: 190.855 (9.23891) minmax: 175.188 207.011
9 5739 0.021 1.099 1.032 -8.553 17.461 1.057 2.334 *M 10.3
ho x -306.8764817179339 269.60836040281544
ho y -279.6653090011365 302.0355612673752
TSYS[10] slice: 187.981 (9.43341) minmax: 170.423 212.883
10 5739 0.018 1.088 1.019 -7.292 17.335 1.048 2.449 *M 12.1
ho x -306.8764817179339 269.60836040281544
ho y -279.6653090011365 302.0355612673752
TSYS[11] slice: 189.746 (7.27148) minmax: 174.056 203.415
11 5739 0.008 1.053 1.018 -6.378 17.383 1.031 2.011 *M 7.7
ho x -306.8713520166912 269.5608652952257
ho y -279.8002836378797 302.15557791425806
TSYS[12] slice: 169.171 (6.84001) minmax: 157.443 182.828
12 5738 0.015 0.950 0.908 -6.626 16.827 0.925 1.776 *M 8.1
ho x -306.8713520166912 269.5608652952257
ho y -279.8002836378797 302.15557791425806
TSYS[13] slice: 205.225 (8.24847) minmax: 190.98 312.662
13 5738 0.013 1.226 1.113 -21.913 38.321 1.190 2.750 *M 6.0
ho x -306.8713520166912 269.5608652952257
ho y -279.8002836378797 302.15557791425806
TSYS[14] slice: 195.98 (8.06578) minmax: 182.17 250.268
14 5738 0.017 1.163 1.067 -12.828 16.048 1.112 2.770 *M 12.4
ho x -306.8713520166912 269.5608652952257
ho y -279.8002836378797 302.15557791425806
TSYS[15] slice: 198.31 (10.4222) minmax: 181.76 219.002
15 5738 0.006 1.108 1.066 -6.620 16.887 1.084 2.065 *M 8.1
CPU TIME: 37.181 sec
CPU TIME: 39.6398 sec
Warnings:
*P ratio: ratio of std/mad too high for data
*M ratio: ratio of std/mad too high for RMS
X-range: -359.687 331.804 Y-range: -356.455 340.486 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,1',
'--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
PJT-seq for first (array([ 0, 5739, 11478, 17217, 22956, 28694, 34432, 40170, 45908,
51647, 57386, 63125, 68864, 74602, 80340, 86078]),)
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
Plots saving in IRC+10216_79448_specviews.6.png
LMTOY>> view_spec_point
/home/teuben/LMT/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_79448_specpoint,png,1
MAD: 0 1.1930885 0.026710356244839983 1.2999299564735005
MAD: 1 0.9733612 0.02274616737584391 1.0643458636371341
MAD: 2 1.2732093 0.05954793332431095 1.5114010667170437
MAD: 3 0.923233 0.021394460183526563 1.008810813356024
MAD: 4 1.042393 0.03217524409003462 1.1710939454916083
MAD: 5 0.98765105 0.02270176146160885 1.0784580959372252
MAD: 6 1.3191054 0.06601926685552169 1.5831824541560955
MAD: 7 0.9967439 0.023271792008421287 1.089831085498895
MAD: 8 0.9447026 0.021822082509861788 1.0319909553141053
MAD: 9 0.9965455 0.02325164365330568 1.0895520682159754
MAD: 10 0.9808667 0.02204459811591926 1.0690450630721975
MAD: 11 0.9904771 0.022930772261091747 1.0822001741578924
MAD: 12 0.88376987 0.020176545138775145 0.9644760503594829
MAD: 13 1.1019962 0.050420109867090036 1.303676622863746
MAD: 14 1.0311316 0.025526110161926487 1.133236065823182
MAD: 15 1.0363319 0.023375273253335226 1.1298329850268907
Plots saving in IRC+10216_79448_specpoint.1.png
LMTOY>> view_spec_point
/home/teuben/LMT/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_79448_specpoint,png,2
MAD: 0 1.1930885 0.026710356244839983 1.2999299564735005
MAD: 1 0.9733612 0.02274616737584391 1.0643458636371341
MAD: 2 1.2732093 0.05954793332431095 1.5114010667170437
MAD: 3 0.923233 0.021394460183526563 1.008810813356024
MAD: 4 1.042393 0.03217524409003462 1.1710939454916083
MAD: 5 0.98765105 0.02270176146160885 1.0784580959372252
MAD: 6 1.3191054 0.06601926685552169 1.5831824541560955
MAD: 7 0.9967439 0.023271792008421287 1.089831085498895
MAD: 8 0.9447026 0.021822082509861788 1.0319909553141053
MAD: 9 0.9965455 0.02325164365330568 1.0895520682159754
MAD: 10 0.9808667 0.02204459811591926 1.0690450630721975
MAD: 11 0.9904771 0.022930772261091747 1.0822001741578924
MAD: 12 0.88376987 0.020176545138775145 0.9644760503594829
MAD: 13 1.1019962 0.050420109867090036 1.303676622863746
MAD: 14 1.0311316 0.025526110161926487 1.133236065823182
MAD: 15 1.0363319 0.023375273253335226 1.1298329850268907
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 1317 x 16
bumpie: 0 1004 1.3243725541879914
bumpie: 1 653 1.014861935585865
bumpie: 2 1004 3.6159087931407607
bumpie: 3 656 1.0183182718784274
bumpie: 4 1004 2.754705513042831
bumpie: 5 655 1.010690987618855
bumpie: 6 1004 3.783742205689004
bumpie: 7 660 1.0150976453997693
bumpie: 8 655 1.0252041396008138
bumpie: 9 656 1.0134061682798352
bumpie: 10 1004 1.4728024856982722
bumpie: 11 658 1.027031299862106
bumpie: 12 658 1.0176925140325133
bumpie: 13 1004 3.662766408157238
bumpie: 14 1004 2.556188067808433
bumpie: 15 658 1.0397645932472805
Written waterfall cube to IRC+10216_79448.wf.fits
chan0: 365
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 1317 x 16
bumpie: 0 1004 1.3261891314819836
bumpie: 1 653 1.0150481531777447
bumpie: 2 1004 3.615347135543742
bumpie: 3 656 1.0182025619399173
bumpie: 4 1004 2.7571837971996787
bumpie: 5 655 1.0109092856850284
bumpie: 6 1004 3.7832893691506544
bumpie: 7 660 1.0155431460702224
bumpie: 8 655 1.02499225595796
bumpie: 9 656 1.0135513247982018
bumpie: 10 1004 1.4731829303447757
bumpie: 11 658 1.026809167515583
bumpie: 12 658 1.0172768918766761
bumpie: 13 1004 3.662352265922719
bumpie: 14 1004 2.5557053477907448
bumpie: 15 658 1.0398009563655313
New shape: (16, 1317, 573)
Written waterfall cube to IRC+10216_79448.wf10.fits
chan0: 365
mv: cannot stat 'stats_wf0.png': No such file or directory
### nemo Debug Info: Reading 2 column(s)
LMTOY>> grid_data native
{'--cell': '5.5',
'--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': '11',
'--rmax': '3',
'--rms_cut': '-4',
'--sample': '-1',
'--weight': 'IRC+10216_79448.wt.fits',
'--x_extent': '400',
'--y_extent': '400'}
STDOUT ****************************
spec_driver_fits LMTSLR 4-dec-2022
CMD: spec_driver_fits -i IRC+10216_79448.nc -o IRC+10216_79448.fits -w IRC+10216_79448.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_79448.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_79448.nc"
file: IRC+10216_79448.nc nspec= 91816 nchan= 1317 npix=16 ncal=1
SpecFile version 11-jul-2023
Global Tsys min/max = 157.443 387.956
Tsys[0] = 228.58 +/- 10.2429
Tsys[1] = 186.372 +/- 7.53352
Tsys[2] = 237.356 +/- 11.3537
Tsys[3] = 176.799 +/- 7.17621
Tsys[4] = 196.994 +/- 10.7623
Tsys[5] = 189.769 +/- 8.25308
Tsys[6] = 246.279 +/- 13.8708
Tsys[7] = 191.478 +/- 5.46294
Tsys[8] = 180.613 +/- 9.61018
Tsys[9] = 190.855 +/- 9.23923
Tsys[10] = 187.981 +/- 9.43315
Tsys[11] = 189.746 +/- 7.27313
Tsys[12] = 169.172 +/- 6.84111
Tsys[13] = 205.224 +/- 8.25166
Tsys[14] = 195.98 +/- 8.06468
Tsys[15] = 198.31 +/- 10.424
First/Last beam: [mapcoord=0] 287.776 1.65713 -91.4461 334.485
X-range: -359.687 331.804 Y-range: -356.455 340.486 arcsec
X-ramp: -276.287 248.404 Y-ramp: -273.055 257.086 arcsec
MapSize: 691.491 x 696.941 arcsec
79448
IRC+10216
DATE-OBS 2018-11-16T06:48:52
DeltaFreq=390625 DeltaTime=0.1 : sqrt(Df.Dt) = 197.642
jinc filter
CF.n_cells= 7 delta=0.129412 cell=5.5 resol=11 otf_select=1
C-AXIS 2 starts at 880.446
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
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=91816 rms_cutoff=-4
Cube Completed, 84857/91816 Spectra accepted = 0.924
Weighting Completed, fuzzy_edge=0
Weight at 0.000000 0.000000 is 16.974379
write to IRC+10216_79448.fits
write weights to IRC+10216_79448.wt.fits
write expected RMS to radiometer.rms.fits
STDERR ****************************
Exit Code: 0
LMTOY>> grid_data smooth
@todo this needs a new option to bin in channel space
{'--cell': '5.5',
'--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.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': '22',
'--rmax': '3',
'--rms_cut': '-4',
'--sample': '-1',
'--weight': 'IRC+10216_79448.wt.fits.fits',
'--x_extent': '400',
'--y_extent': '400'}
STDOUT ****************************
spec_driver_fits LMTSLR 4-dec-2022
CMD: spec_driver_fits -i IRC+10216_79448.nc -o IRC+10216_79448.fits.fits -w IRC+10216_79448.wt.fits.fits -l 22 -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_79448.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_79448.nc"
file: IRC+10216_79448.nc nspec= 91816 nchan= 1317 npix=16 ncal=1
SpecFile version 11-jul-2023
Global Tsys min/max = 157.443 387.956
Tsys[0] = 228.58 +/- 10.2429
Tsys[1] = 186.372 +/- 7.53352
Tsys[2] = 237.356 +/- 11.3537
Tsys[3] = 176.799 +/- 7.17621
Tsys[4] = 196.994 +/- 10.7623
Tsys[5] = 189.769 +/- 8.25308
Tsys[6] = 246.279 +/- 13.8708
Tsys[7] = 191.478 +/- 5.46294
Tsys[8] = 180.613 +/- 9.61018
Tsys[9] = 190.855 +/- 9.23923
Tsys[10] = 187.981 +/- 9.43315
Tsys[11] = 189.746 +/- 7.27313
Tsys[12] = 169.172 +/- 6.84111
Tsys[13] = 205.224 +/- 8.25166
Tsys[14] = 195.98 +/- 8.06468
Tsys[15] = 198.31 +/- 10.424
First/Last beam: [mapcoord=0] 287.776 1.65713 -91.4461 334.485
X-range: -359.687 331.804 Y-range: -356.455 340.486 arcsec
X-ramp: -276.287 248.404 Y-ramp: -273.055 257.086 arcsec
MapSize: 691.491 x 696.941 arcsec
79448
IRC+10216
DATE-OBS 2018-11-16T06:48:52
DeltaFreq=390625 DeltaTime=0.1 : sqrt(Df.Dt) = 197.642
jinc filter
CF.n_cells= 13 delta=0.258824 cell=5.5 resol=22 otf_select=1
C-AXIS 2 starts at 880.446
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
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=91816 rms_cutoff=-4
Cube Completed, 84857/91816 Spectra accepted = 0.924
Weighting Completed, fuzzy_edge=0
Weight at 0.000000 0.000000 is 69.648529
write to IRC+10216_79448.fits.fits
write weights to IRC+10216_79448.wt.fits.fits
write expected RMS to radiometer.rms.fits
STDERR ****************************
Exit Code: 0
LMTOY>> Some NEMO post-processing
@todo rid the code of fits header stealing
### nemo Debug Info: FITS file: Image size 147 147 1317
### nemo Debug Info: CUNIT1 deg
### nemo Debug Info: CUNIT2 deg
### nemo Debug Info: CUNIT3 m/s
### Warning [fitsccd]: Setting map minimum from 0 to -675.63
### Warning [fitsccd]: Setting map maximum from 0 to 1488.1
### nemo Debug Info: There were 16906329 blank values in the image (59.405803 %)
### nemo Debug Info: FITS file: Image size 147 147
### nemo Debug Info: CUNIT1 deg
### nemo Debug Info: CUNIT2 deg
### Warning [fitsccd]: Setting map minimum from 0 to -3.35573
### Warning [fitsccd]: Setting map maximum from 0 to 53.311
### 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
### nemo Debug Info: WCS Corner: 146.989 13.2788 880446
QAC_STATS: on 0.00110446 0.220483 -16.5098 20.1344 92342 0.0641871 6912134
LMTOY>> STATS IRC+10216_79448.wt.ccd centerbox robust
### Warning [ccdsub]: Sampling will be done
### nemo Debug Info: WCS Corner: 146.989 13.2788 0
QAC_STATS: wt 17.351 2.46149 -1.34482 53.311 91823 0.999576 4582
### Warning [ccdsub]: Sampling will be done
### nemo Debug Info: WCS Corner: 146.989 13.2788 880446
### nemo Debug Info: Keeping 147*147*1317 cube
SLABS: [[-920,-120],[80,880]] == -920,-120,80,880
### Warning [ccdslice]: axis=1 not fully supported yet
### nemo Debug Info: Axis: 1317 0 880446 -1367.42
### nemo Debug Info: Axis min/max: 880446 -919078
### nemo Debug Info: slab -920 -120 -> 1316 731
### nemo Debug Info: slab 80 880 -> 585 0
### nemo Debug Info: Total planes: 1172
### nemo Debug Info: Z_SLICE: 732 1
### nemo Debug Info: FITS file: Image size 147 147
### nemo Debug Info: CUNIT1 deg
### nemo Debug Info: CUNIT2 deg
### Warning [fitsccd]: Setting map minimum from 0 to -872.872
### Warning [fitsccd]: Setting map maximum from 0 to 549433
### nemo Debug Info: There were 10038 blank values in the image (46.452867 %)
### nemo Debug Info: FITS file: Image size 147 147 1317
### nemo Debug Info: CUNIT1 deg
### nemo Debug Info: CUNIT2 deg
### nemo Debug Info: CUNIT3 m/s
### Warning [fitsccd]: Setting map minimum from 0 to -16.1588
### Warning [fitsccd]: Setting map maximum from 0 to 17.3618
### nemo Debug Info: There were 16870770 blank values in the image (59.280855 %)
### nemo Debug Info: Averaging map 1 * 1 * 4 pixels; mapsize 147 * 147 * 1317
### Warning [ccdsub]: Z binning not even
### Warning [ccdslice]: axis=1 not fully supported yet
### Warning [ccdslice]: Stepping assumed to be regular as 4 for WCS correction
### nemo Debug Info: Z_SLICE: 1 4
### Warning [ccdsub]: Sampling will be done
### nemo Debug Info: WCS Corner: 146.989 13.2788 880446
### Warning [ccdsub]: Sampling will be done
### nemo Debug Info: WCS Corner: 146.989 13.2788 0
### Warning [ccdsub]: Sampling will be done
### nemo Debug Info: WCS Corner: 146.989 13.2788 880446
[0;31mQAC_STATS: IRC+10216_79448-full 0.00100673 0.244634 -675.63 1488.1 91148.1 0.0291053 11053099[0m
[0;31mQAC_STATS: IRC+10216_79448-cent 0.00110446 0.220483 -16.5098 20.1344 92342 0.0641871 6912134[0m
[0;31mQAC_STATS: RMS/radiometer 1.98187 0.196105 0.888413 28.5834 11207.2 1 4947[0m
### nemo Debug Info: FITS file: Image size 147 147 330
### nemo Debug Info: CUNIT1 deg
### nemo Debug Info: CUNIT2 deg
### nemo Debug Info: CUNIT3 m/s
### nemo Debug Info: There were no blank values set in the image
cubespec : # 73 73 1317 min/mean/sig/max -1.21648 0.283584 1.96331 15.2313 -0.00012564 0.432642 0.155821
cubespecs: # 73 73 330 min/mean/sig/max -0.137512 0.152597 1.03677 7.92212 0.000359187 0.451607 0.308009
PJT TABPLOT IRC+10216_79448 -1 ,
Traceback (most recent call last):
File "/home/teuben/LMT/lmtoy/bin/tab_plot.py", line 65, in
xmin = float(cw[0])
ValueError: could not convert string to float: ''
mv: cannot stat 'tab_plot.png': No such file or directory
### nemo Debug Info: axis=1: new xmin,ymin=147.101 13.1672
### nemo Debug Info: Image minmax: -67.7187 288.176
### nemo Debug Info: axis=1: new xmin,ymin=147.101 13.1672
### nemo Debug Info: Image minmax: 0 23639.7
### nemo Debug Info: axis=1: new xmin,ymin=147.101 13.1672
### nemo Debug Info: Image minmax: -2.30473e+06 7.70142e+06
### nemo Debug Info: axis=1: new xmin,ymin=147.101 13.1672
### nemo Debug Info: Image minmax: 0 157198
### nemo Debug Info: axis=1: new xmin,ymin=147.101 13.1672
### nemo Debug Info: Image minmax: -3.35573 53.311
### nemo Debug Info: axis=1: new xmin,ymin=147.101 13.1672
### nemo Debug Info: Image minmax: 0 147.901
### nemo Debug Info: axis=1: new xmin,ymin=147.101 13.1672
### nemo Debug Info: Image minmax: 0 148.805
### nemo Debug Info: axis=1: new xmin,ymin=147.101 13.1672
### nemo Debug Info: Image minmax: 0 1
### nemo Debug Info: axis=1: new xmin,ymin=147.101 13.1672
### nemo Debug Info: Image minmax: 0 1.08237
### nemo Debug Info: axis=1: new xmin,ymin=147.101 13.1672
### nemo Debug Info: Image minmax: 0 0.00118117
### nemo Debug Info: axis=1: new xmin,ymin=147.101 13.1672
### nemo Debug Info: Image minmax: 0 785.235
WARNING: FITSFixedWarning: RADECSYS= 'FK5 ' /
the RADECSYS keyword is deprecated, use RADESYSa. [astropy.wcs.wcs]
WARNING: FITSFixedWarning: 'datfix' made the change 'Set MJD-OBS to 60137.505440 from DATE-OBS'. [astropy.wcs.wcs]
Data min/max/mean/sig: -67.7188 288.176 8.58246 19.5163
Data min/max: -49.9664 67.1313
BINS: [-49.9664278 -46.18908048 -42.41173317 -38.63438585 -34.85703853
-31.07969121 -27.30234389 -23.52499657 -19.74764925 -15.97030194
-12.19295462 -8.4156073 -4.63825998 -0.86091266 2.91643466
6.69378198 10.47112929 14.24847661 18.02582393 21.80317125
25.58051857 29.35786589 33.13521321 36.91256052 40.68990784
44.46725516 48.24460248 52.0219498 55.79929712 59.57664444
63.35399175 67.13133907]
INFO: Auto-setting vmin to -1.599e+01 [aplpy.core]
INFO: Auto-setting vmax to 1.130e+02 [aplpy.core]
INFO: Auto-setting vmin to -1.681e+01 [aplpy.core]
INFO: Auto-setting vmax to 1.229e+02 [aplpy.core]
Writing IRC+10216_79448.mom0.png
WARNING: FITSFixedWarning: RADECSYS= 'FK5 ' /
the RADECSYS keyword is deprecated, use RADESYSa. [astropy.wcs.wcs]
WARNING: FITSFixedWarning: 'datfix' made the change 'Set MJD-OBS to 60137.505440 from DATE-OBS'. [astropy.wcs.wcs]
Data min/max/mean/sig: 248.069 23639.7 1145.19 972.546
Data min/max: -1772.45 4062.83
BINS: [-1772.44665527 -1584.21187075 -1395.97708622 -1207.74230169
-1019.50751717 -831.27273264 -643.03794812 -454.80316359
-266.56837906 -78.33359454 109.90118999 298.13597451
486.37075904 674.60554357 862.84032809 1051.07511262
1239.30989715 1427.54468167 1615.7794662 1804.01425072
1992.24903525 2180.48381978 2368.7186043 2556.95338883
2745.18817336 2933.42295788 3121.65774241 3309.89252693
3498.12731146 3686.36209599 3874.59688051 4062.83166504]
INFO: Auto-setting vmin to -5.776e+02 [aplpy.core]
INFO: Auto-setting vmax to 6.411e+03 [aplpy.core]
INFO: Auto-setting vmin to -5.826e+02 [aplpy.core]
INFO: Auto-setting vmax to 6.467e+03 [aplpy.core]
Writing IRC+10216_79448.peak.png
WARNING: FITSFixedWarning: RADECSYS= 'FK5 ' /
the RADECSYS keyword is deprecated, use RADESYSa. [astropy.wcs.wcs]
WARNING: FITSFixedWarning: 'datfix' made the change 'Set MJD-OBS to 60137.505440 from DATE-OBS'. [astropy.wcs.wcs]
Data min/max/mean/sig: -2.30473e+06 7.70142e+06 1246.08 108677
Data min/max: -324784 327276
BINS: [-324784.30847168 -303750.08972168 -282715.87097168 -261681.65222168
-240647.43347168 -219613.21472168 -198578.99597168 -177544.77722168
-156510.55847168 -135476.33972168 -114442.12097168 -93407.90222168
-72373.68347168 -51339.46472168 -30305.24597168 -9271.02722168
11763.19152832 32797.41027832 53831.62902832 74865.84777832
95900.06652832 116934.28527832 137968.50402832 159002.72277832
180036.94152832 201071.16027832 222105.37902832 243139.59777832
264173.81652832 285208.03527832 306242.25402832 327276.47277832]
INFO: Auto-setting vmin to -5.625e+04 [aplpy.core]
INFO: Auto-setting vmax to 7.682e+04 [aplpy.core]
INFO: Auto-setting vmin to -8.467e+04 [aplpy.core]
INFO: Auto-setting vmax to 6.105e+04 [aplpy.core]
Writing IRC+10216_79448.mom1.png
WARNING: FITSFixedWarning: RADECSYS= 'FK5 ' /
the RADECSYS keyword is deprecated, use RADESYSa. [astropy.wcs.wcs]
WARNING: FITSFixedWarning: 'datfix' made the change 'Set MJD-OBS to 60137.505440 from DATE-OBS'. [astropy.wcs.wcs]
Data min/max/mean/sig: 118.944 157198 353.945 2025.76
Data min/max: -5723.33 6431.22
BINS: [-5723.32803345 -5331.24589982 -4939.16376618 -4547.08163255
-4154.99949892 -3762.91736529 -3370.83523166 -2978.75309803
-2586.67096439 -2194.58883076 -1802.50669713 -1410.4245635
-1018.34242987 -626.26029624 -234.17816261 157.90397103
549.98610466 942.06823829 1334.15037192 1726.23250555
2118.31463918 2510.39677282 2902.47890645 3294.56104008
3686.64317371 4078.72530734 4470.80744097 4862.8895746
5254.97170824 5647.05384187 6039.1359755 6431.21810913]
INFO: Auto-setting vmin to -1.291e+02 [aplpy.core]
INFO: Auto-setting vmax to 1.433e+03 [aplpy.core]
INFO: Auto-setting vmin to -1.641e+02 [aplpy.core]
INFO: Auto-setting vmax to 1.822e+03 [aplpy.core]
Writing IRC+10216_79448.rms.png
Data min/max/mean/sig: -3.35573 53.311 11.597 8.5637
Data min/max: -14.0941 37.2881
BINS: [-14.09413528 -12.43664538 -10.77915549 -9.12166559 -7.46417569
-5.80668579 -4.14919589 -2.49170599 -0.83421609 0.82327381
2.48076371 4.13825361 5.79574351 7.45323341 9.11072331
10.76821321 12.42570311 14.08319301 15.74068291 17.39817281
19.05566271 20.71315261 22.37064251 24.02813241 25.68562231
27.34311221 29.00060211 30.65809201 32.31558191 33.97307181
35.63056171 37.28805161]
INFO: Auto-setting vmin to -5.220e+00 [aplpy.core]
INFO: Auto-setting vmax to 4.471e+01 [aplpy.core]
INFO: Auto-setting vmin to -4.866e+00 [aplpy.core]
INFO: Auto-setting vmax to 4.050e+01 [aplpy.core]
Writing IRC+10216_79448.wt.png
WARNING: FITSFixedWarning: RADECSYS= 'FK5 ' /
the RADECSYS keyword is deprecated, use RADESYSa. [astropy.wcs.wcs]
WARNING: FITSFixedWarning: 'datfix' made the change 'Set MJD-OBS to 60137.505301 from DATE-OBS'. [astropy.wcs.wcs]
Data min/max/mean/sig: 0.120473 147.901 0.342059 1.93623
Data min/max: -5.46663 6.15075
BINS: [-5.46663493 -5.09188047 -4.71712601 -4.34237155 -3.96761709 -3.59286263
-3.21810817 -2.84335371 -2.46859925 -2.09384479 -1.71909033 -1.34433587
-0.96958141 -0.59482695 -0.22007249 0.15468197 0.52943642 0.90419088
1.27894534 1.6536998 2.02845426 2.40320872 2.77796318 3.15271764
3.5274721 3.90222656 4.27698102 4.65173548 5.02648994 5.4012444
5.77599886 6.15075332]
INFO: Setting slices=[0] [aplpy.core]
INFO: Auto-setting vmin to -1.469e-01 [aplpy.core]
INFO: Auto-setting vmax to 1.631e+00 [aplpy.core]
INFO: Auto-setting vmin to -1.603e-01 [aplpy.core]
INFO: Auto-setting vmax to 1.780e+00 [aplpy.core]
Writing IRC+10216_79448.wt2.png
WARNING: FITSFixedWarning: RADECSYS= 'FK5 ' /
the RADECSYS keyword is deprecated, use RADESYSa. [astropy.wcs.wcs]
WARNING: FITSFixedWarning: 'datfix' made the change 'Set MJD-OBS to 60137.505312 from DATE-OBS'. [astropy.wcs.wcs]
Data min/max/mean/sig: 0.122217 148.805 0.343418 1.95035
Data min/max: -5.50764 6.19448
BINS: [-5.50764081 -5.13015315 -4.75266549 -4.37517783 -3.99769017 -3.62020251
-3.24271485 -2.86522719 -2.48773953 -2.11025187 -1.73276421 -1.35527655
-0.9777889 -0.60030124 -0.22281358 0.15467408 0.53216174 0.9096494
1.28713706 1.66462472 2.04211238 2.41960004 2.7970877 3.17457536
3.55206302 3.92955068 4.30703834 4.684526 5.06201366 5.43950132
5.81698897 6.19447663]
INFO: Setting slices=[0] [aplpy.core]
INFO: Auto-setting vmin to -1.640e-01 [aplpy.core]
INFO: Auto-setting vmax to 1.821e+00 [aplpy.core]
INFO: Auto-setting vmin to -1.568e-01 [aplpy.core]
INFO: Auto-setting vmax to 1.741e+00 [aplpy.core]
Writing IRC+10216_79448.wt3.png
WARNING: FITSFixedWarning: RADECSYS= 'FK5 ' /
the RADECSYS keyword is deprecated, use RADESYSa. [astropy.wcs.wcs]
WARNING: FITSFixedWarning: 'datfix' made the change 'Set MJD-OBS to 60137.505440 from DATE-OBS'. [astropy.wcs.wcs]
Data min/max/mean/sig: 2.83583e-07 1 0.448799 0.198321
Data min/max: -0.146163 1.04376
BINS: [-0.14616305 -0.10777843 -0.06939381 -0.03100919 0.00737543 0.04576005
0.08414467 0.1225293 0.16091392 0.19929854 0.23768316 0.27606778
0.3144524 0.35283702 0.39122164 0.42960626 0.46799088 0.5063755
0.54476012 0.58314474 0.62152936 0.65991398 0.6982986 0.73668322
0.77506784 0.81345246 0.85183708 0.8902217 0.92860632 0.96699094
1.00537556 1.04376018]
INFO: Auto-setting vmin to -8.332e-02 [aplpy.core]
INFO: Auto-setting vmax to 9.249e-01 [aplpy.core]
INFO: Auto-setting vmin to -8.411e-02 [aplpy.core]
INFO: Auto-setting vmax to 9.336e-01 [aplpy.core]
Writing IRC+10216_79448.wtn.png
WARNING: FITSFixedWarning: RADECSYS= 'FK5 ' /
the RADECSYS keyword is deprecated, use RADESYSa. [astropy.wcs.wcs]
WARNING: FITSFixedWarning: 'datfix' made the change 'Set MJD-OBS to 60137.505312 from DATE-OBS'. [astropy.wcs.wcs]
Data min/max/mean/sig: 0.900381 1.08237 1.00286 0.0205756
Data min/max: 0.941134 1.06459
BINS: [0.9411337 0.94511607 0.94909845 0.95308083 0.95706321 0.96104558
0.96502796 0.96901034 0.97299272 0.97697509 0.98095747 0.98493985
0.98892223 0.9929046 0.99688698 1.00086936 1.00485173 1.00883411
1.01281649 1.01679887 1.02078124 1.02476362 1.028746 1.03272838
1.03671075 1.04069313 1.04467551 1.04865789 1.05264026 1.05662264
1.06060502 1.0645874 ]
INFO: Setting slices=[0] [aplpy.core]
INFO: Auto-setting vmin to -1.048e-01 [aplpy.core]
INFO: Auto-setting vmax to 1.164e+00 [aplpy.core]
INFO: Auto-setting vmin to -1.053e-01 [aplpy.core]
INFO: Auto-setting vmax to 1.169e+00 [aplpy.core]
Writing IRC+10216_79448.wtr.png
WARNING: FITSFixedWarning: RADECSYS= 'FK5 ' /
the RADECSYS keyword is deprecated, use RADESYSa. [astropy.wcs.wcs]
WARNING: FITSFixedWarning: 'datfix' made the change 'Set MJD-OBS to 60137.505312 from DATE-OBS'. [astropy.wcs.wcs]
Data min/max/mean/sig: 0.000770751 0.00118117 0.000970385 4.03201e-05
Data min/max: 0.000849424 0.00109134
BINS: [0.00084942 0.00085723 0.00086503 0.00087284 0.00088064 0.00088844
0.00089625 0.00090405 0.00091186 0.00091966 0.00092746 0.00093527
0.00094307 0.00095087 0.00095868 0.00096648 0.00097429 0.00098209
0.00098989 0.0009977 0.0010055 0.00101331 0.00102111 0.00102891
0.00103672 0.00104452 0.00105233 0.00106013 0.00106793 0.00107574
0.00108354 0.00109134]
INFO: Setting slices=[0] [aplpy.core]
INFO: Auto-setting vmin to -1.045e-04 [aplpy.core]
INFO: Auto-setting vmax to 1.160e-03 [aplpy.core]
INFO: Auto-setting vmin to -1.044e-04 [aplpy.core]
INFO: Auto-setting vmax to 1.158e-03 [aplpy.core]
Writing IRC+10216_79448.wtr3.png
WARNING: FITSFixedWarning: RADECSYS= 'FK5 ' /
the RADECSYS keyword is deprecated, use RADESYSa. [astropy.wcs.wcs]
WARNING: FITSFixedWarning: 'datfix' made the change 'Set MJD-OBS to 60137.505312 from DATE-OBS'. [astropy.wcs.wcs]
Data min/max/mean/sig: 0.380349 785.235 2.39865 10.4532
Data min/max: -28.961 33.7584
BINS: [-28.96104956 -26.93784285 -24.91463614 -22.89142942 -20.86822271
-18.845016 -16.82180929 -14.79860258 -12.77539587 -10.75218916
-8.72898245 -6.70577574 -4.68256903 -2.65936232 -0.63615561
1.38705111 3.41025782 5.43346453 7.45667124 9.47987795
11.50308466 13.52629137 15.54949808 17.57270479 19.5959115
21.61911821 23.64232492 25.66553164 27.68873835 29.71194506
31.73515177 33.75835848]
INFO: Setting slices=[0] [aplpy.core]
INFO: Auto-setting vmin to -7.623e-01 [aplpy.core]
INFO: Auto-setting vmax to 8.462e+00 [aplpy.core]
INFO: Auto-setting vmin to -7.147e-01 [aplpy.core]
INFO: Auto-setting vmax to 7.933e+00 [aplpy.core]
Writing IRC+10216_79448.wtr4.png
LMTOY>> Created IRC+10216_79448.nf.fits and IRC+10216_79448.nfs.fits
LMTOY>> skipping ADMIT post-processing
LMTOY>> Running maskmoment vlsr=-20
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
YES IRC+10216_79448.fits
NO IRC+10216_79448__0.fits
NO IRC+10216_79448__1.fits
IRC+10216_79448
IRC+10216_79448_specviews
IRC+10216_79448_specpoint
Writing index_admit.html
BASE: IRC+10216_79448 IRC+10216_79448_specviews IRC+10216_79448_specpoint
Writing index_mm.html
Writing index_pars.html
Writing index_log.html
RC: ./lmtoy_79448.rc
Writing index_pipeline.html 115.271204
BASE: IRC+10216_79448 IRC+10216_79448_specviews IRC+10216_79448_specpoint
Found IRC+10216_79448.fits
Writing index.html
BASE: IRC+10216_79448
BASE: IRC+10216_79448__0
BASE: IRC+10216_79448__1
LMTOY>> Processed 1 bands