LMTOY>> seq_pipeline: 2-feb-2022
LMTOY>> READING lmtoy_functions 24-dec-2021 via /lma1/teuben/lmtoy/examples/seq_pipeline.sh
Working directory ./2018-S1-MU-46/85776
LMTOY>> Hang on, creating a bootstrap lmtoy_85776.rc from path=/n/lma1/lmt/data_lmt/
#
# version=2-feb-2022
# ifproc="/n/lma1/lmt/data_lmt//ifproc/ifproc_2019-10-31_085776_00_0001.nc"
date_obs="2019-10-31T03:01:57"
# skytime=1471.28 sec
# inttime=1297.48 sec
obsnum="85776"
calobsnum="85775"
obspgm="Map"
obsgoal="Science"
ProjectId="2018-S1-MU-46"
# SkyOff=0 0
# HPBW=12 arcsec
# XYLength=300 300 arcsec
# XYRamp=0 0 arcsec
# XYoff=0 0 arcsec
# XYstep=1.25 0.4
numbands=1
vlsr=-296 # km/s
skyfreq=115.38505778665652 # GHz
restfreq=115.271204 # Ghz
src="Region_J-K"
resolution=11 # arcsec
cell=5.5 # arcsec
x_extent=300 # arcsec
y_extent=300 # arcsec
instrument="SEQ"
#
LMTOY>> this is your startup lmtoy_85776.rc file:
# seq_pipeline: 2-feb-2022
# DATE: 2022-02-03T21:38:32.109911125
# obsnum=85776
# Using ifproc=/n/lma1/lmt/data_lmt//ifproc/ifproc_2019-10-31_085776_00_0001.nc
# path=/n/lma1/lmt/data_lmt/
#
# version=2-feb-2022
# ifproc="/n/lma1/lmt/data_lmt//ifproc/ifproc_2019-10-31_085776_00_0001.nc"
date_obs="2019-10-31T03:01:57"
# skytime=1471.28 sec
# inttime=1297.48 sec
obsnum="85776"
calobsnum="85775"
obspgm="Map"
obsgoal="Science"
ProjectId="2018-S1-MU-46"
# SkyOff=0 0
# HPBW=12 arcsec
# XYLength=300 300 arcsec
# XYRamp=0 0 arcsec
# XYoff=0 0 arcsec
# XYstep=1.25 0.4
numbands=1
vlsr=-296 # km/s
skyfreq=115.38505778665652 # GHz
restfreq=115.271204 # Ghz
src="Region_J-K"
resolution=11 # arcsec
cell=5.5 # arcsec
x_extent=300 # arcsec
y_extent=300 # arcsec
instrument="SEQ"
#
# based on vlsr=-296, dv=100, dw=250
b_order=0
b_regions=[[-646,-396],[-196,54]]
l_regions=[[-396,-196]]
slice=[-646,54]
v_range=-396,-196
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_85776.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 536, 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 667, 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: divide by zero encountered in true_divide
self.tsys_spectrum = 280 * sky_spectrum / (hot_spectrum -
/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: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': '1',
'--b_regions': '[[-646,-396],[-196,54]]',
'--bank': '0',
'--eliminate_list': '0',
'--help': False,
'--l_regions': '[[-396,-196]]',
'--obsnum': '85776',
'--output': 'Region_J-K_85776.nc',
'--path': '/n/lma1/lmt/data_lmt/',
'--pix_list': '0,1,2,4,5,6,7,8,9,10,11,12,13,14,15',
'--restfreq': None,
'--sample': None,
'--save_tsys': False,
'--slice': '[-646,54]',
'--stype': '2',
'--tsys': '250.0',
'--use_cal': True,
'--use_otf_cal': False,
'--x_axis': 'VLSR'}
program process_otf_map options
obsnum = 85776
path = /n/lma1/lmt/data_lmt/
output = Region_J-K_85776.nc
use_cal = True
use_otf_cal = False
bank = 0
x_axis = VLSR
b_regions = [[-646, -396], [-196, 54]]
eliminate_list = [0]
l_regions = [[-396, -196]]
data_path = /n/lma1/lmt/data_lmt/
b_order = 1
stype = 2
pix_list = [0, 1, 2, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15]
slice = [-646, 54]
tsys = 250.0
output_file_name = Region_J-K_85776.nc
2019-10-31T03:01:57 begin /n/lma1/lmt/data_lmt/ifproc/ifproc_2019-10-31_085776_00_0001.nc
2019-10-31T03:26:28 end /n/lma1/lmt/data_lmt/ifproc/ifproc_2019-10-31_085776_00_0001.nc
BufPos 3 4.2 sec 85776
BufPos 2 4.0 sec 85776
BufPos 1 10.6 sec 85776
BufPos 0 718.0 sec 85776
BufPos 1 10.6 sec 85776
BufPos 0 718.3 sec 85776
BufPos 1 5.6 sec 85776
before read npix
from pixels npix = 16
from xlen npix = 16
TRACKING Sequoia PIXEL 10
Map Parameters: Ra Continuous
HPBW= 12.0 XLength= 300.0 YLength= 300.0 XStep= 1.25 YStep= 0.40
/n/lma1/lmt/data_lmt/ifproc/ifproc_2019-10-31_085776_00_0001.nc does not have bs parameters
85776 is a Map observation
get chop
'Data.Msip1mm.BeamChopperActPos'
no chop
2019-10-31T03:01:32 begin /n/lma1/lmt/data_lmt/ifproc/ifproc_2019-10-31_085775_00_0001.nc
2019-10-31T03:01:38 end /n/lma1/lmt/data_lmt/ifproc/ifproc_2019-10-31_085775_00_0001.nc
BufPos 3 4.2 sec 85775
BufPos 2 2.6 sec 85775
before read npix
from pixels npix = 16
from xlen npix = 16
TRACKING Sequoia PIXEL 10
/n/lma1/lmt/data_lmt/ifproc/ifproc_2019-10-31_085775_00_0001.nc does not have map parameters
/n/lma1/lmt/data_lmt/ifproc/ifproc_2019-10-31_085775_00_0001.nc does not have bs parameters
85775 is a Cal observation
get chop cal
'Data.Msip1mm.BeamChopperActPos'
no chop cal
read_roach /n/lma1/lmt/data_lmt/spectrometer/roach0/roach0_85776_0_1_Region_J-K_2019-10-31_070149.nc nspec,nchan=57824,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
read_roach /n/lma1/lmt/data_lmt/spectrometer/roach1/roach1_85776_0_1_Region_J-K_2019-10-31_070149.nc nspec,nchan=57832,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_85776_0_1_Region_J-K_2019-10-31_070149.nc nspec,nchan=57852,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_85776_0_1_Region_J-K_2019-10-31_070149.nc nspec,nchan=57860,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_85775_0_1_Region_J-K_2019-10-31_070131.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
read_roach /n/lma1/lmt/data_lmt/spectrometer/roach1/roach1_85775_0_1_Region_J-K_2019-10-31_070131.nc nspec,nchan=232,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_85775_0_1_Region_J-K_2019-10-31_070131.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_85775_0_1_Region_J-K_2019-10-31_070131.nc nspec,nchan=240,2048
r:3 inp:0 pix:12 time_offset:-0.030000
r:3 inp:1 pix:13 time_offset:-0.030000
r:3 inp:2 pix:14 time_offset:-0.030000
r:3 inp:3 pix:15 time_offset:-0.030000
TSYS[0] = 273.223 +/- 16.1213 (2048 channels)
TSYS[1] = 228.846 +/- 16.7252 (2048 channels)
TSYS[2] = 286.859 +/- 21.7043 (2048 channels)
TSYS[4] = 235.844 +/- 15.7642 (2048 channels)
TSYS[5] = 228.369 +/- 13.6459 (2048 channels)
TSYS[6] = 279.71 +/- 22.8392 (2048 channels)
TSYS[7] = 229.166 +/- 15.9674 (2048 channels)
TSYS[8] = 216.69 +/- 13.6095 (2048 channels)
TSYS[9] = 228.28 +/- 12.6938 (2048 channels)
TSYS[10] = 224.335 +/- 12.6476 (2048 channels)
TSYS[11] = 227.001 +/- 14.1093 (2048 channels)
TSYS[12] = 211.449 +/- 12.028 (2048 channels)
TSYS[13] = 251.519 +/- 22.6493 (2048 channels)
TSYS[14] = 238.724 +/- 12.2556 (2048 channels)
TSYS[15] = 239.218 +/- 14.4866 (2048 channels)
Spectral Band velocity range: -1335.8 743.795 km/s
Total Number of OTF Spectra = 213432
NSPEC: 213432
NCHAN: 689
FAST_NC: True
Looping over 15 pixel list [0, 1, 2, 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: 272.454 (6.45201) minmax: 260.734 285.758
0 14228 0.030 1.457 1.426 -9.211 11.926 1.429 2.738 *M 3.0
TSYS[1] slice: 230.15 (6.81687) minmax: 216.916 241.434
1 14228 -0.016 1.293 1.213 -13.063 10.682 1.231 3.469 *M 7.1
TSYS[2] slice: 289.298 (8.88673) minmax: 271.389 304.781
2 14228 0.007 1.552 1.516 -10.198 14.028 1.529 3.604 *M 4.9
TSYS[4] slice: 235.854 (5.61154) minmax: 224.217 246.073
4 14229 0.019 1.287 1.238 -9.342 14.002 1.241 2.613 *M 3.6
TSYS[5] slice: 228.077 (4.93491) minmax: 218.338 236.471
5 14229 -0.007 1.219 1.195 -9.478 8.651 1.196 2.569 *M 3.6
TSYS[6] slice: 279.547 (8.62068) minmax: 262.185 295.214
6 14229 0.018 1.513 1.463 -9.237 16.300 1.469 3.071 *M 3.7
TSYS[7] slice: 229.672 (5.46014) minmax: 217.48 240.052
7 14229 -0.003 1.203 1.196 -8.235 7.995 1.197 2.236 *M 2.2
TSYS[8] slice: 216.237 (4.8358) minmax: 207.009 224.71
8 14229 0.014 1.165 1.135 -8.000 12.013 1.138 2.569 *M 3.6
TSYS[9] slice: 227.167 (5.16018) minmax: 218.023 237.007
9 14229 -0.019 1.239 1.194 -11.164 10.268 1.197 2.946 *M 4.2
TSYS[10] slice: 223.816 (5.41305) minmax: 214.017 233.764
10 14229 0.005 1.180 1.168 -7.755 9.570 1.170 2.144 *M 2.6
TSYS[11] slice: 226.906 (5.23418) minmax: 216.192 236.9
11 14229 -0.004 1.192 1.182 -8.561 8.061 1.183 2.291 *M 2.7
TSYS[12] slice: 211.688 (3.91914) minmax: 203.808 219.665
12 14229 0.009 1.129 1.108 -6.700 10.835 1.114 2.418 *M 3.6
TSYS[13] slice: 251.151 (4.45858) minmax: 243.572 260.409
13 14229 -0.012 1.364 1.322 -11.027 10.461 1.331 3.406 *M 4.6
TSYS[14] slice: 238.938 (3.64413) minmax: 231.789 246.136
14 14229 0.002 1.280 1.252 -9.046 10.835 1.261 3.058 *M 4.6
TSYS[15] slice: 238.164 (5.3363) minmax: 228.165 249.234
15 14229 0.003 1.260 1.247 -8.113 8.703 1.253 2.425 *M 3.1
CPU TIME: 116.651 sec
Traceback (most recent call last):
File "/lma1/teuben/lmtoy/anaconda3/bin/process_otf_map2.py", line 7, in
exec(compile(f.read(), __file__, 'exec'))
File "/lma1/teuben/lmtoy/SpectralLineReduction/bin/process_otf_map2.py", line 110, in
main(sys.argv[1:])
File "/lma1/teuben/lmtoy/SpectralLineReduction/bin/process_otf_map2.py", line 105, in main
specfile.write_ncdata()
File "/lma1/teuben/lmtoy/SpectralLineReduction/lmtslr/spec/specfile.py", line 325, in write_ncdata
self._create_nc_data()
File "/lma1/teuben/lmtoy/SpectralLineReduction/lmtslr/spec/specfile.py", line 286, in _create_nc_data
nc_tsys[:,:,:] = tmp_tsys[:,:,:]
File "src/netCDF4/_netCDF4.pyx", line 4918, in netCDF4._netCDF4.Variable.__setitem__
File "/lma1/teuben/lmtoy/anaconda3/lib/python3.8/site-packages/netCDF4/utils.py", line 365, in _StartCountStride
datashape = broadcasted_shape(shape, datashape)
File "/lma1/teuben/lmtoy/anaconda3/lib/python3.8/site-packages/netCDF4/utils.py", line 973, in broadcasted_shape
return np.broadcast(a, b).shape
ValueError: shape mismatch: objects cannot be broadcast to a single shape
LMTOY>> view_spec_file
{'--help': False,
'--input': 'Region_J-K_85776.nc',
'--pix_list': '0,1,2,4,5,6,7,8,9,10,11,12,13,14,15',
'--plot_range': '-1,3',
'--plots': 'Region_J-K_85776_specviews',
'--rms_cut': '-4',
'--show_pixel': None,
'--skip_tsys': False}
SLR Plots initialized with 1 method=Region_J-K_85776_specviews
Max sequence=14229
Plots saving in Region_J-K_85776_specviews.1.png
Plots saving in Region_J-K_85776_specviews.2.png
Plots saving in Region_J-K_85776_specviews.3.png
Plots saving in Region_J-K_85776_specviews.4.png
Plots saving in Region_J-K_85776_specviews.5.png
NEW TSYS PIC
SHAPE (1, 15, 689)
PIXEL: 0 [ 0 1 2 ... 14225 14226 14227]
PIXEL: 1 [14228 14229 14230 ... 28453 28454 28455]
PIXEL: 2 [28456 28457 28458 ... 42681 42682 42683]
PIXEL: 4 [42684 42685 42686 ... 56910 56911 56912]
PIXEL: 5 [56913 56914 56915 ... 71139 71140 71141]
PIXEL: 6 [71142 71143 71144 ... 85368 85369 85370]
PIXEL: 7 [85371 85372 85373 ... 99597 99598 99599]
PIXEL: 8 [ 99600 99601 99602 ... 113826 113827 113828]
PIXEL: 9 [113829 113830 113831 ... 128055 128056 128057]
PIXEL: 10 [128058 128059 128060 ... 142284 142285 142286]
PIXEL: 11 [142287 142288 142289 ... 156513 156514 156515]
PIXEL: 12 [156516 156517 156518 ... 170742 170743 170744]
PIXEL: 13 [170745 170746 170747 ... 184971 184972 184973]
PIXEL: 14 [184974 184975 184976 ... 199200 199201 199202]
PIXEL: 15 [199203 199204 199205 ... 213429 213430 213431]
Plots saving in Region_J-K_85776_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=Region_J-K_85776_specpoint,png,1
MAD: 0 1.410631 0.04698525716977158 1.5985719700700776
MAD: 1 1.1916258 0.039566597468212915 1.3498922233842041
MAD: 2 1.4956046 0.05008978288716834 1.6959637658336466
MAD: 4 1.2221187 0.04042820475934093 1.3838315543507793
MAD: 5 1.1790614 0.03881828048962091 1.3343345347697628
MAD: 6 1.4467978 0.047192749879470926 1.6355688472656864
MAD: 7 1.1874303 0.041880653727236324 1.354952877474558
MAD: 8 1.1216302 0.03748318685327608 1.2715629392160828
MAD: 9 1.1771525 0.03962280077458804 1.3356437175560547
MAD: 10 1.1585048 0.03797540763395127 1.310406474247658
MAD: 11 1.1717888 0.03865161471002425 1.3263952705237518
MAD: 12 1.097261 0.0361083267265715 1.2416942589021356
MAD: 13 1.3046552 0.044874275121518654 1.4841522947686063
MAD: 14 1.2376969 0.03982057278664388 1.3969791772091975
MAD: 15 1.235385 0.04121964628747863 1.4002635262509888
Plots saving in Region_J-K_85776_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=Region_J-K_85776_specpoint,png,2
MAD: 0 1.410631 0.04698525716977158 1.5985719700700776
MAD: 1 1.1916258 0.039566597468212915 1.3498922233842041
MAD: 2 1.4956046 0.05008978288716834 1.6959637658336466
MAD: 4 1.2221187 0.04042820475934093 1.3838315543507793
MAD: 5 1.1790614 0.03881828048962091 1.3343345347697628
MAD: 6 1.4467978 0.047192749879470926 1.6355688472656864
MAD: 7 1.1874303 0.041880653727236324 1.354952877474558
MAD: 8 1.1216302 0.03748318685327608 1.2715629392160828
MAD: 9 1.1771525 0.03962280077458804 1.3356437175560547
MAD: 10 1.1585048 0.03797540763395127 1.310406474247658
MAD: 11 1.1717888 0.03865161471002425 1.3263952705237518
MAD: 12 1.097261 0.0361083267265715 1.2416942589021356
MAD: 13 1.3046552 0.044874275121518654 1.4841522947686063
MAD: 14 1.2376969 0.03982057278664388 1.3969791772091975
MAD: 15 1.235385 0.04121964628747863 1.4002635262509888
Plots saving in Region_J-K_85776_specpoint.2.png
LMTOY>> make_spec_fits (waterfall)
pix [0, 1, 2, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15]
Min nsamp 14228
Binning [1, 1]
FITS file will be 14228 x 689 x 15
Written waterfall cube to Region_J-K_85776.wf.fits
pix [0, 1, 2, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15]
Min nsamp 14228
Binning [10, 1]
FITS file will be 1422 x 689 x 15
New shape: (15, 689, 1422)
Written waterfall cube to Region_J-K_85776.wf10.fits
LMTOY>> grid_data
{'--cell': '5.5',
'--edge': '0',
'--help': False,
'--input': 'Region_J-K_85776.nc',
'--n_samples': '256',
'--noise_sigma': '1',
'--otf_a': '1.1',
'--otf_b': '1',
'--otf_c': '2',
'--otf_select': '2',
'--output': 'Region_J-K_85776.fits',
'--pix_list': '0,1,2,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': 'Region_J-K_85776.wt.fits',
'--x_extent': '250',
'--y_extent': '250'}
STDOUT ****************************
spec_driver_fits LMTSLR 26-jan-2022
CMD: spec_driver_fits -i Region_J-K_85776.nc -o Region_J-K_85776.fits -w Region_J-K_85776.wt.fits -l 11 -c 5.5 -u [0,1,2,4,5,6,7,8,9,10,11,12,13,14,15] -z -4 -s 1 -x 250 -y 250 -f 2 -r 3 -n 256 -0 1.1 -1 1 -2 2 -b -1 -e 0
file list=Region_J-K_85776.nc
DEBUG:5.5
sample_list=-1
Warning: sample list is not a multiple of 3
nsamples=1
Processing 1 SpecFiles:
Opening SpecFile file Region_J-K_85776.nc
file: Region_J-K_85776.nc nspec= 213432 nchan= 689
SpecFile version 6-mar-2021
X-range: -216.22 192.529 Y-range: -185.827 224.666 arcsec
X-ramp: -132.82 109.129 Y-ramp: -102.427 141.266 arcsec
MapSize: 408.748 x 410.493 arcsec
85776
Region_J-K
DATE-OBS 2019-10-31T03:01:57
CF.n_cells= 7 delta=0.129412 cell=5.5 resol=11 oft_select=2
C-AXIS 2 starts at 53.9848
C-AXIS 1 starts at -247.5
C-AXIS 0 starts at -247.5
P-AXIS 1 starts at -247.5
P-AXIS 0 starts at -247.5
P-AXIS 1 starts at -247.5
P-AXIS 0 starts at -247.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 0 1 1 1 1 1 1 1 1 1 1 1 1
rms_stats: n=213432 rms_cutoff=-4
Cube Completed, 208941/213432 Spectra accepted = 0.979
Weighting Completed, fuzzy_edge=0
Weight at 0.000000 0.000000 is 189.436752
write to Region_J-K_85776.fits
write weights to Region_J-K_85776.wt.fits
STDERR ****************************
Exit Code: 0
LMTOY>> Some NEMO post-processing
### nemo Debug Info: FITS file: Image size 91 91 689
### Warning [fitsccd]: Setting map minimum from 0 to -1.11521
### Warning [fitsccd]: Setting map maximum from 0 to 1.33357
### nemo Debug Info: There were 1962961 blank values in the image (34.404057 %)
### nemo Debug Info: FITS file: Image size 91 91
### Warning [fitsccd]: Setting map maximum from 0 to 216.069
### nemo Debug Info: There were no blank values set in the image
LMTOY>> STATS Region_J-K_85776.ccd centerbox robust
### Warning [ccdsub]: Sampling will be done
### Warning [ccdsub]: Attempting to fix the WCS
### nemo Debug Info: WCS Corner: 10.6151 41.5473 53984.8
Number of points : 1457924
Min and Max : -0.332266 1.333569
Mean and dispersion : 0.002518 0.071444
Skewness and kurtosis : 1.504429 14.158853
Sum and Sum*Dx*Dy*Dz* : 3670.647218 8.704084
Points per beam (map) : 1
N Robust : 1439336
Mean Robust : 0.000439
Sigma Robust : 0.062382
Median Robust : 0.000203
0/1457924 out-of-range points discarded
LMTOY>> STATS Region_J-K_85776.wt.ccd centerbox robust
### Warning [ccdsub]: Sampling will be done
### Warning [ccdsub]: Attempting to fix the WCS
### nemo Debug Info: WCS Corner: 10.6151 41.5473 0
Number of points : 2116
Min and Max : 98.285919 216.068710
Mean and dispersion : 181.754615 17.794531
Skewness and kurtosis : -1.678943 2.682912
Sum and Sum*Dx*Dy* : 384592.764648 -0.897680
Points per beam (map) : 1
N Robust : 1683
Mean Robust : 187.853586
Sigma Robust : 5.674239
Median Robust : 188.463165
0/2116 out-of-range points discarded
### nemo Debug Info: Keeping 91*91*689 cube
### Warning [ccdfits]: Replacing header from Region_J-K_85776.wt.fits
SLABS: [[-646,-396],[-196,54]] == -646,-396,-196,54
### Warning [ccdslice]: axis=1 not fully supported yet
### nemo Debug Info: Axis: 689 0 53984.8 -1015.92
### nemo Debug Info: Axis min/max: 53984.8 -644969
### nemo Debug Info: slab -646 -396 -> 689 442
### nemo Debug Info: slab -196 54 -> 246 0
### nemo Debug Info: Total planes: 495
### Warning [ccdslice]: Code not converted to fix reference pixel value
### Warning [ccdfits]: Replacing header from Region_J-K_85776.wt.fits
### Warning [ccdfits]: Replacing header from Region_J-K_85776.wt.fits
### Warning [ccdfits]: Replacing header from Region_J-K_85776.fits
[0;31mQAC_STATS: Region_J-K_85776-full 0.00146799 0.11 -1.11521 1.33357 5494.16 0.0187528[0m
### Warning [ccdsub]: Sampling will be done
### Warning [ccdsub]: Attempting to fix the WCS
### nemo Debug Info: WCS Corner: 10.6151 41.5473 53984.8
[0;31mQAC_STATS: Region_J-K_85776-cent 0.00251772 0.0714442 -0.332266 1.33357 3670.65 0.0474508[0m
### nemo Debug Info: FITS file: Image size 91 91 689
### Warning [fitsccd]: Setting map minimum from 0 to -0.0902669
### Warning [fitsccd]: Setting map maximum from 0 to 0.746021
### nemo Debug Info: There were no blank values set in the image
spectab : # 45 45 689 min/mean/sig/max -0.177768 0.00193127 0.0608632 0.162234 -7.97393e-05 0.0842558 0.978881
specstab: # 45 45 689 min/mean/sig/max -0.0471451 0.00209688 0.0194742 0.0730862 -2.17013e-05 0.00644006 0.233838
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: Region_J-K_85776.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: Region_J-K_85776.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: Region_J-K_85776.mom2.fits
INFO: Auto-setting vmin to -2.109e+01 [aplpy.core]
INFO: Auto-setting vmax to 2.341e+02 [aplpy.core]
INFO: Auto-setting vmin to -2.109e+01 [aplpy.core]
INFO: Auto-setting vmax to 2.341e+02 [aplpy.core]
Writing Region_J-K_85776.wt.png
INFO: Auto-setting vmin to -3.396e+01 [aplpy.core]
INFO: Auto-setting vmax to 3.770e+02 [aplpy.core]
INFO: Auto-setting vmin to -3.396e+01 [aplpy.core]
INFO: Auto-setting vmax to 3.770e+02 [aplpy.core]
Writing Region_J-K_85776.wt2.png
INFO: Auto-setting vmin to -4.141e+03 [aplpy.core]
INFO: Auto-setting vmax to 4.596e+04 [aplpy.core]
INFO: Auto-setting vmin to -4.141e+03 [aplpy.core]
INFO: Auto-setting vmax to 4.596e+04 [aplpy.core]
Writing Region_J-K_85776.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: Region_J-K_85776.wtn.fits
INFO: Auto-setting vmin to -1.083e-01 [aplpy.core]
INFO: Auto-setting vmax to 1.202e+00 [aplpy.core]
INFO: Auto-setting vmin to -1.083e-01 [aplpy.core]
INFO: Auto-setting vmax to 1.202e+00 [aplpy.core]
Writing Region_J-K_85776.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 -5.451e+03 [aplpy.core]
INFO: Auto-setting vmax to 1.538e+04 [aplpy.core]
INFO: Auto-setting vmin to -5.451e+03 [aplpy.core]
INFO: Auto-setting vmax to 1.538e+04 [aplpy.core]
Writing Region_J-K_85776.mom0.png
LMTOY>> Created Region_J-K_85776.nf.fits and Region_J-K_85776.nfs.fits
LMTOY>> skipping ADMIT post-processing
LMTOY>> Created Region_J-K_85776.fits and Region_J-K_85776.wt.fits
LMTOY>> Parameter file used: lmtoy_85776.rc
LMTOY>> Making summary index.html:
Making index.html for obsnum=85776
Region_J-K_85776_?
Region_J-K_85776_?_specviews
Region_J-K_85776_?_specpoint
Writing index.html
Writing index_pipeline.html
cp: cannot stat 'Region_J-K_85776_?_specviews.1.png': No such file or directory