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/85824
LMTOY>> Hang on, creating a bootstrap lmtoy_85824.rc from path=/n/lma1/lmt/data_lmt/
#
# version=2-feb-2022
# ifproc="/n/lma1/lmt/data_lmt//ifproc/ifproc_2019-11-01_085824_00_0001.nc"
date_obs="2019-11-01T03:09:48"
# skytime=1488.4 sec
# inttime=1297.48 sec
obsnum="85824"
calobsnum="85823"
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.38488685307446 # 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_85824.rc file:
# seq_pipeline: 2-feb-2022
# DATE: 2022-02-03T21:48:34.544887453
# obsnum=85824
# Using ifproc=/n/lma1/lmt/data_lmt//ifproc/ifproc_2019-11-01_085824_00_0001.nc
# path=/n/lma1/lmt/data_lmt/
#
# version=2-feb-2022
# ifproc="/n/lma1/lmt/data_lmt//ifproc/ifproc_2019-11-01_085824_00_0001.nc"
date_obs="2019-11-01T03:09:48"
# skytime=1488.4 sec
# inttime=1297.48 sec
obsnum="85824"
calobsnum="85823"
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.38488685307446 # 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_85824.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': '85824',
'--output': 'Region_J-K_85824.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
output = Region_J-K_85824.nc
use_cal = True
b_regions = [[-646, -396], [-196, 54]]
stype = 2
path = /n/lma1/lmt/data_lmt/
pix_list = [0, 1, 2, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15]
b_order = 1
output_file_name = Region_J-K_85824.nc
eliminate_list = [0]
l_regions = [[-396, -196]]
bank = 0
tsys = 250.0
slice = [-646, 54]
x_axis = VLSR
obsnum = 85824
use_otf_cal = False
data_path = /n/lma1/lmt/data_lmt/
2019-11-01T03:09:48 begin /n/lma1/lmt/data_lmt/ifproc/ifproc_2019-11-01_085824_00_0001.nc
2019-11-01T03:34:37 end /n/lma1/lmt/data_lmt/ifproc/ifproc_2019-11-01_085824_00_0001.nc
BufPos 3 5.2 sec 85824
BufPos 2 4.0 sec 85824
BufPos 1 11.3 sec 85824
BufPos 0 725.7 sec 85824
BufPos 1 10.7 sec 85824
BufPos 0 726.0 sec 85824
BufPos 1 5.5 sec 85824
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-11-01_085824_00_0001.nc does not have bs parameters
85824 is a Map observation
get chop
'Data.Msip1mm.BeamChopperActPos'
no chop
2019-11-01T03:09:20 begin /n/lma1/lmt/data_lmt/ifproc/ifproc_2019-11-01_085823_00_0001.nc
2019-11-01T03:09:27 end /n/lma1/lmt/data_lmt/ifproc/ifproc_2019-11-01_085823_00_0001.nc
BufPos 3 4.3 sec 85823
BufPos 2 2.6 sec 85823
before read npix
from pixels npix = 16
from xlen npix = 16
TRACKING Sequoia PIXEL 10
/n/lma1/lmt/data_lmt/ifproc/ifproc_2019-11-01_085823_00_0001.nc does not have map parameters
/n/lma1/lmt/data_lmt/ifproc/ifproc_2019-11-01_085823_00_0001.nc does not have bs parameters
85823 is a Cal observation
get chop cal
'Data.Msip1mm.BeamChopperActPos'
no chop cal
read_roach /n/lma1/lmt/data_lmt/spectrometer/roach0/roach0_85824_0_1_Region_J-K_2019-11-01_070938.nc nspec,nchan=58396,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_85824_0_1_Region_J-K_2019-11-01_070938.nc nspec,nchan=58412,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_85824_0_1_Region_J-K_2019-11-01_070938.nc nspec,nchan=58420,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_85824_0_1_Region_J-K_2019-11-01_070938.nc nspec,nchan=58440,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_85823_0_1_Region_J-K_2019-11-01_070920.nc nspec,nchan=228,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_85823_0_1_Region_J-K_2019-11-01_070920.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_85823_0_1_Region_J-K_2019-11-01_070920.nc nspec,nchan=236,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_85823_0_1_Region_J-K_2019-11-01_070920.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] = 310.917 +/- 17.9693 (2048 channels)
TSYS[1] = 259.91 +/- 19.6778 (2048 channels)
TSYS[2] = 323.08 +/- 23.3755 (2048 channels)
TSYS[4] = 268.243 +/- 17.3265 (2048 channels)
TSYS[5] = 261.767 +/- 15.1227 (2048 channels)
TSYS[6] = 315.264 +/- 25.2557 (2048 channels)
TSYS[7] = 260.886 +/- 18.333 (2048 channels)
TSYS[8] = 250.407 +/- 15.649 (2048 channels)
TSYS[9] = 261.689 +/- 14.2288 (2048 channels)
TSYS[10] = 256.424 +/- 14.2224 (2048 channels)
TSYS[11] = 259.034 +/- 15.6052 (2048 channels)
TSYS[12] = 242.229 +/- 12.9186 (2048 channels)
TSYS[13] = 284.374 +/- 23.9932 (2048 channels)
TSYS[14] = 272.825 +/- 13.5293 (2048 channels)
TSYS[15] = 274.394 +/- 16.2792 (2048 channels)
Spectral Band velocity range: -1335.8 743.795 km/s
Total Number of OTF Spectra = 215587
NSPEC: 215587
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: 310.428 (7.37492) minmax: 297.063 325.249
0 14373 0.022 1.654 1.622 -9.974 13.449 1.627 3.187 *M 3.0
TSYS[1] slice: 261.541 (6.93753) minmax: 247.878 272.734
1 14373 -0.018 1.460 1.370 -15.389 12.272 1.388 3.963 *M 6.5
TSYS[2] slice: 325.599 (9.05173) minmax: 307.652 340.31
2 14373 0.013 1.739 1.702 -11.890 15.523 1.715 3.980 *M 4.7
TSYS[4] slice: 268.201 (6.33614) minmax: 255.435 279.531
4 14372 0.031 1.460 1.404 -10.228 16.254 1.409 3.097 *M 3.7
TSYS[5] slice: 261.476 (5.44509) minmax: 250.411 270.84
5 14372 -0.012 1.399 1.369 -10.711 9.910 1.372 3.021 *M 3.7
TSYS[6] slice: 315.148 (9.53852) minmax: 296.342 332.186
6 14372 0.030 1.705 1.647 -10.325 19.713 1.651 3.446 *M 3.5
TSYS[7] slice: 261.473 (6.16351) minmax: 249.098 273.49
7 14372 -0.002 1.370 1.362 -7.919 8.609 1.363 2.562 *M 2.4
TSYS[8] slice: 249.793 (5.32861) minmax: 239.423 260.29
8 14373 0.021 1.342 1.306 -10.057 14.815 1.310 3.045 *M 3.7
TSYS[9] slice: 260.392 (5.73754) minmax: 248.926 270.907
9 14373 -0.022 1.421 1.367 -12.948 11.432 1.374 3.519 *M 4.6
TSYS[10] slice: 255.94 (6.16668) minmax: 244.02 268.202
10 14373 0.007 1.350 1.334 -9.948 11.079 1.336 2.586 *M 2.7
TSYS[11] slice: 258.762 (5.65284) minmax: 247.835 268.786
11 14373 -0.007 1.362 1.349 -9.471 9.406 1.351 2.635 *M 2.8
TSYS[12] slice: 242.256 (4.35996) minmax: 233.213 251.137
12 14372 0.012 1.288 1.264 -8.360 13.131 1.270 2.831 *M 3.8
TSYS[13] slice: 283.862 (5.21147) minmax: 273.504 293.977
13 14372 -0.011 1.541 1.490 -13.171 12.096 1.501 3.813 *M 4.9
TSYS[14] slice: 272.797 (4.33609) minmax: 263.161 281.438
14 14372 0.006 1.462 1.428 -9.929 13.590 1.438 3.540 *M 4.6
TSYS[15] slice: 273.337 (5.97683) minmax: 262.045 285.107
15 14372 -0.002 1.446 1.428 -9.877 10.846 1.435 2.852 *M 3.5
CPU TIME: 119.247 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_85824.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_85824_specviews',
'--rms_cut': '-4',
'--show_pixel': None,
'--skip_tsys': False}
SLR Plots initialized with 1 method=Region_J-K_85824_specviews
Max sequence=14373
Plots saving in Region_J-K_85824_specviews.1.png
Plots saving in Region_J-K_85824_specviews.2.png
Plots saving in Region_J-K_85824_specviews.3.png
Plots saving in Region_J-K_85824_specviews.4.png
Plots saving in Region_J-K_85824_specviews.5.png
NEW TSYS PIC
SHAPE (1, 15, 689)
PIXEL: 0 [ 0 1 2 ... 14370 14371 14372]
PIXEL: 1 [14373 14374 14375 ... 28743 28744 28745]
PIXEL: 2 [28746 28747 28748 ... 43116 43117 43118]
PIXEL: 4 [43119 43120 43121 ... 57488 57489 57490]
PIXEL: 5 [57491 57492 57493 ... 71860 71861 71862]
PIXEL: 6 [71863 71864 71865 ... 86232 86233 86234]
PIXEL: 7 [ 86235 86236 86237 ... 100604 100605 100606]
PIXEL: 8 [100607 100608 100609 ... 114977 114978 114979]
PIXEL: 9 [114980 114981 114982 ... 129350 129351 129352]
PIXEL: 10 [129353 129354 129355 ... 143723 143724 143725]
PIXEL: 11 [143726 143727 143728 ... 158096 158097 158098]
PIXEL: 12 [158099 158100 158101 ... 172468 172469 172470]
PIXEL: 13 [172471 172472 172473 ... 186840 186841 186842]
PIXEL: 14 [186843 186844 186845 ... 201212 201213 201214]
PIXEL: 15 [201215 201216 201217 ... 215584 215585 215586]
Plots saving in Region_J-K_85824_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_85824_specpoint,png,1
MAD: 0 1.609211 0.05241753649287122 1.818881113989063
MAD: 1 1.3486049 0.04468534010732053 1.5273462779555271
MAD: 2 1.6812463 0.05537245183653672 1.9027360880163129
MAD: 4 1.3870989 0.04626371629468848 1.5721537736031315
MAD: 5 1.3507328 0.04517384934889547 1.5314282007403084
MAD: 6 1.6274322 0.053527905273694144 1.8415438482294813
MAD: 7 1.351484 0.04486031266490345 1.5309251917377997
MAD: 8 1.2915685 0.04320372904649465 1.4643834338709152
MAD: 9 1.3474652 0.044929594728108 1.527183536421282
MAD: 10 1.3231797 0.04363886282100894 1.497735173116311
MAD: 11 1.3371965 0.044240132155248434 1.5141569979279395
MAD: 12 1.2507411 0.04129670490879798 1.4159279437883292
MAD: 13 1.4695039 0.04942992425713214 1.6672235765756478
MAD: 14 1.4105041 0.04649542237851798 1.5964857922209812
MAD: 15 1.4142089 0.04683246547681155 1.6015387509148145
Plots saving in Region_J-K_85824_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_85824_specpoint,png,2
MAD: 0 1.609211 0.05241753649287122 1.818881113989063
MAD: 1 1.3486049 0.04468534010732053 1.5273462779555271
MAD: 2 1.6812463 0.05537245183653672 1.9027360880163129
MAD: 4 1.3870989 0.04626371629468848 1.5721537736031315
MAD: 5 1.3507328 0.04517384934889547 1.5314282007403084
MAD: 6 1.6274322 0.053527905273694144 1.8415438482294813
MAD: 7 1.351484 0.04486031266490345 1.5309251917377997
MAD: 8 1.2915685 0.04320372904649465 1.4643834338709152
MAD: 9 1.3474652 0.044929594728108 1.527183536421282
MAD: 10 1.3231797 0.04363886282100894 1.497735173116311
MAD: 11 1.3371965 0.044240132155248434 1.5141569979279395
MAD: 12 1.2507411 0.04129670490879798 1.4159279437883292
MAD: 13 1.4695039 0.04942992425713214 1.6672235765756478
MAD: 14 1.4105041 0.04649542237851798 1.5964857922209812
MAD: 15 1.4142089 0.04683246547681155 1.6015387509148145
Plots saving in Region_J-K_85824_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 14372
Binning [1, 1]
FITS file will be 14372 x 689 x 15
Written waterfall cube to Region_J-K_85824.wf.fits
pix [0, 1, 2, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15]
Min nsamp 14372
Binning [10, 1]
FITS file will be 1437 x 689 x 15
New shape: (15, 689, 1437)
Written waterfall cube to Region_J-K_85824.wf10.fits
LMTOY>> grid_data
{'--cell': '5.5',
'--edge': '0',
'--help': False,
'--input': 'Region_J-K_85824.nc',
'--n_samples': '256',
'--noise_sigma': '1',
'--otf_a': '1.1',
'--otf_b': '4.75',
'--otf_c': '2',
'--otf_select': '1',
'--output': 'Region_J-K_85824.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_85824.wt.fits',
'--x_extent': '250',
'--y_extent': '250'}
STDOUT ****************************
spec_driver_fits LMTSLR 26-jan-2022
CMD: spec_driver_fits -i Region_J-K_85824.nc -o Region_J-K_85824.fits -w Region_J-K_85824.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 1 -r 3 -n 256 -0 1.1 -1 4.75 -2 2 -b -1 -e 0
file list=Region_J-K_85824.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_85824.nc
file: Region_J-K_85824.nc nspec= 215587 nchan= 689
SpecFile version 6-mar-2021
X-range: -216.685 194.43 Y-range: -185.993 223.181 arcsec
X-ramp: -133.285 111.03 Y-ramp: -102.593 139.781 arcsec
MapSize: 411.114 x 409.173 arcsec
85824
Region_J-K
DATE-OBS 2019-11-01T03:09:48
CF.n_cells= 7 delta=0.129412 cell=5.5 resol=11 oft_select=1
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=215587 rms_cutoff=-4
Cube Completed, 211256/215587 Spectra accepted = 0.980
Weighting Completed, fuzzy_edge=0
Weight at 0.000000 0.000000 is 53.315315
write to Region_J-K_85824.fits
write weights to Region_J-K_85824.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 -24.579
### Warning [fitsccd]: Setting map maximum from 0 to 16.9879
### nemo Debug Info: There were 1980875 blank values in the image (34.718029 %)
### nemo Debug Info: FITS file: Image size 91 91
### Warning [fitsccd]: Setting map minimum from 0 to -1.73376
### Warning [fitsccd]: Setting map maximum from 0 to 66.5245
### nemo Debug Info: There were no blank values set in the image
LMTOY>> STATS Region_J-K_85824.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.777011 1.653386
Mean and dispersion : 0.004103 0.141541
Skewness and kurtosis : 0.253195 1.443946
Sum and Sum*Dx*Dy*Dz* : 5981.783294 14.184405
Points per beam (map) : 1
N Robust : 1443157
Mean Robust : 0.002612
Sigma Robust : 0.133632
Median Robust : 0.002528
0/1457924 out-of-range points discarded
LMTOY>> STATS Region_J-K_85824.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 : 23.789110 66.524490
Mean and dispersion : 48.748379 5.335928
Skewness and kurtosis : -0.994461 2.113953
Sum and Sum*Dx*Dy* : 103151.569948 -0.240767
Points per beam (map) : 1
N Robust : 1848
Mean Robust : 49.651785
Sigma Robust : 3.073833
Median Robust : 49.895718
0/2116 out-of-range points discarded
### nemo Debug Info: Keeping 91*91*689 cube
### Warning [ccdfits]: Replacing header from Region_J-K_85824.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_85824.wt.fits
### Warning [ccdfits]: Replacing header from Region_J-K_85824.wt.fits
### Warning [ccdfits]: Replacing header from Region_J-K_85824.fits
[0;31mQAC_STATS: Region_J-K_85824-full 0.00368699 0.315649 -24.579 16.9879 13733.1 0.0204856[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_85824-cent 0.00410295 0.141541 -0.777011 1.65339 5981.78 0.0369215[0m
### nemo Debug Info: FITS file: Image size 91 91 689
### Warning [fitsccd]: Setting map minimum from 0 to -0.0908337
### Warning [fitsccd]: Setting map maximum from 0 to 0.77721
### nemo Debug Info: There were no blank values set in the image
spectab : # 45 45 689 min/mean/sig/max -0.375437 0.00218038 0.128436 0.369485 -1.90449e-05 0.176907 0.973961
specstab: # 45 45 689 min/mean/sig/max -0.0730288 0.0013163 0.0246846 0.0653572 2.65203e-06 0.00761745 0.218207
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_85824.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_85824.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_85824.mom2.fits
INFO: Auto-setting vmin to -7.058e+00 [aplpy.core]
INFO: Auto-setting vmax to 6.739e+01 [aplpy.core]
INFO: Auto-setting vmin to -7.058e+00 [aplpy.core]
INFO: Auto-setting vmax to 6.739e+01 [aplpy.core]
Writing Region_J-K_85824.wt.png
INFO: Auto-setting vmin to -7.255e+00 [aplpy.core]
INFO: Auto-setting vmax to 8.053e+01 [aplpy.core]
INFO: Auto-setting vmin to -7.255e+00 [aplpy.core]
INFO: Auto-setting vmax to 8.053e+01 [aplpy.core]
Writing Region_J-K_85824.wt2.png
INFO: Auto-setting vmin to -2.717e+03 [aplpy.core]
INFO: Auto-setting vmax to 3.016e+04 [aplpy.core]
INFO: Auto-setting vmin to -2.717e+03 [aplpy.core]
INFO: Auto-setting vmax to 3.016e+04 [aplpy.core]
Writing Region_J-K_85824.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_85824.wtn.fits
INFO: Auto-setting vmin to -1.107e-01 [aplpy.core]
INFO: Auto-setting vmax to 1.229e+00 [aplpy.core]
INFO: Auto-setting vmin to -1.107e-01 [aplpy.core]
INFO: Auto-setting vmax to 1.229e+00 [aplpy.core]
Writing Region_J-K_85824.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 -6.959e+03 [aplpy.core]
INFO: Auto-setting vmax to 1.791e+04 [aplpy.core]
INFO: Auto-setting vmin to -6.959e+03 [aplpy.core]
INFO: Auto-setting vmax to 1.791e+04 [aplpy.core]
Writing Region_J-K_85824.mom0.png
LMTOY>> Created Region_J-K_85824.nf.fits and Region_J-K_85824.nfs.fits
LMTOY>> ADMIT post-processing
LMTOY>> lmtoy_admit: 31-oct-2021
/lma1/teuben/lmtoy/examples/lmtoy_admit.sh: line 33: export: `Region_J-K_85824.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
Found 33 Xvfb
1256822 Sl Xvfb
1256146 Sl Xvfb
1254824 Sl Xvfb
1247243 Sl Xvfb
### (1/1) : processing Region_J-K_85824.nfs.fits
Region_J-K_85824.nfs.fits [91, 91, 689] Region_J-K 10.64875 41.5137 5.398476e-05 0.021
/lma1/teuben/lmtoy/admit/admit/test/admit1.py --basename x Region_J-K_85824.nfs.fits
117.474u 3.035s 2:16.29 88.4% 0+0k 416+202952io 0pf+0w
Logfile: Region_J-K_85824.nfs.fits.log
LMTOY>> lmtoy_admit: 31-oct-2021
/lma1/teuben/lmtoy/examples/lmtoy_admit.sh: line 33: export: `Region_J-K_85824.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 Region_J-K_85824.nf.fits
Region_J-K_85824.nf.fits [91, 91, 689] Region_J-K 10.64875 41.5137 5.398476e-05 0.021
/lma1/teuben/lmtoy/admit/admit/test/admit1.py --basename x Region_J-K_85824.nf.fits
66.847u 2.204s 1:21.24 84.9% 0+0k 312+140864io 0pf+0w
Logfile: Region_J-K_85824.nf.fits.log
LMTOY>> Created Region_J-K_85824.fits and Region_J-K_85824.wt.fits
LMTOY>> Parameter file used: lmtoy_85824.rc
LMTOY>> Making summary index.html:
Making index.html for obsnum=85824
Region_J-K_85824_?
Region_J-K_85824_?_specviews
Region_J-K_85824_?_specpoint
Writing index.html
Writing index_pipeline.html
cp: cannot stat 'Region_J-K_85824_?_specviews.1.png': No such file or directory