obsname = names[3] + '_' + names[4] + '_' + names[5]
return obsname, srcname
params = {
'legend.fontsize': 'small',
'axes.titlesize': 'small',
'axes.labelsize': 'x-small',
'xtick.labelsize': 'xx-small',
'ytick.labelsize': 'xx-small',
'image.interpolation': 'none'
}
rcParams.update(params)
savefile = '/data/jrv/burst_paper/all_burst_dynspec.npy'
ds_list = load_dict(savefile)
ds_dir = '/data/jrv/burst_paper/ds/' # where to save ds plots
if not os.path.exists(ds_dir):
os.system('mkdir ' + ds_dir)
obs = ds_list.keys()[0]
ds = ds_list[obs]
close('all')
### PLOT INDIVIDUAL OBSERVATIONS ###
srclist = [('ADLeo', '15A-416'), ('UVCet', '13A-423'), ('UVCet', '15A-416'),
('EQPeg', '15A-416'), ('YZCMi', '15A-416')]
#srclist = [('UVCet','15A-416')]
from dynspec import load_dict
from dynspec.pipeline_utils import load_burst_filelist
import numpy as np
import pickle
def get_band(filename):
# get band from dynspec filename
return filename.split('/')[6]
savefile = '/data/jrv/burst_paper/all_burst_epoch_dynspec_LSband.npy'
savefileP = '/data/jrv/burst_paper/all_burst_epoch_dynspec_Pband.npy'
filelist = load_burst_filelist()
try:
ds_list = load_dict(savefile)
except:
ds_list = {}
try:
dsP_list = load_dict(savefileP)
except:
dsP_list = {}
reload_obslist = filelist.keys()
#reload_obslist = ['/data/jrv/15A-416/EQPeg/2'] # change this to only reload one observation
for obs in reload_obslist:
ds_files = filelist[obs]
ds_obs = None
band_list = [get_band(f) for f in ds_files]
for f in ds_files:
band = get_band(f)
'ds.phase does not update if you bin or truncate the time list, so make sure to do',\
'this last before plotting.'
params = {
'legend.fontsize': 'small',
'axes.titlesize': 'small',
'axes.labelsize': 'small',
'xtick.labelsize': 'x-small',
'ytick.labelsize': 'x-small',
'image.interpolation': 'nearest'
}
rcParams.update(params)
loadfile = '/data/jrv/burst_paper/all_burst_epoch_dynspec_LSband.npy'
ds_list = load_dict(loadfile)
loadfileP = '/data/jrv/burst_paper/all_burst_epoch_dynspec_Pband.npy'
dsP_list = load_dict(loadfileP)
close('all')
# params for ds plots
fig_params = {
'tint_P': 300,
'tint_LS': 60,
'df_MHz_P': 16,
'df_MHz_LS': 16,
'smax_P': None,
'smax_LS': 0.1,
'pixflag_sigfacP': 7.,
'pixflag_sigfacLS': 10.,
import numpy as np
import pickle
def get_band(filename):
# get band from dynspec filename
return filename.split('/')[6]
savefile = '/data/jrv/burst_paper/all_burst_epoch_dynspec_LSband.npy'
savefileP = '/data/jrv/burst_paper/all_burst_epoch_dynspec_Pband.npy'
filelist = load_burst_filelist()
try:
ds_list = load_dict(savefile)
except:
ds_list = {}
try:
dsP_list = load_dict(savefileP)
except:
dsP_list = {}
reload_obslist = filelist.keys()
#reload_obslist = ['/data/jrv/15A-416/EQPeg/2'] # change this to only reload one observation
for obs in reload_obslist:
ds_files = filelist[obs]
ds_obs = None
band_list = [get_band(f) for f in ds_files]
for f in ds_files:
band = get_band(f)
# convert to obs name such as '15A-416_YZCMi_1' and srcname 'YZCMi'
names = obsfile.split('/')
srcname = names[4]
obsname = names[3]+'_'+names[4]+'_'+names[5]
return obsname,srcname
params = {'legend.fontsize': 'small',
'axes.titlesize': 'small',
'axes.labelsize': 'x-small',
'xtick.labelsize': 'xx-small',
'ytick.labelsize': 'xx-small',
'image.interpolation': 'hanning'}
rcParams.update(params)
savefile = '/data/jrv/burst_paper/ds/all_burst_Pdynspec.npy'
ds_list = load_dict(savefile)
ds_dir = '/data/jrv/burst_paper/ds/' # where to save ds plots
if not os.path.exists(ds_dir):
os.system('mkdir '+ds_dir)
obs = ds_list.keys()[0]
ds = ds_list[obs]
close('all')
### PLOT INDIVIDUAL OBSERVATIONS ###
srclist = ['ADLeo','UVCet','EQPeg']
#srclist = ['UVCet']
# and file directory ('/data/jrv/15A-416/YZCMi/1')
names = obsname.split('_')
srcname = names[1]
obsfile = '/data/jrv/'+names[0]+'/'+names[1]+'/'+names[2]
return obsfile, srcname
params = {'legend.fontsize': 'small',
'axes.titlesize': 'small',
'axes.labelsize': 'small',
'xtick.labelsize': 'x-small',
'ytick.labelsize': 'x-small',
'image.interpolation': 'nearest'}
rcParams.update(params)
loadfile = '/data/jrv/burst_paper/all_burst_epoch_dynspec_LSband.npy'
ds_list = load_dict(loadfile)
loadfileP = '/data/jrv/burst_paper/all_burst_epoch_dynspec_Pband.npy'
dsP_list = load_dict(loadfileP)
ds_dir = '/data/jrv/burst_paper/ds/all_burst_dynspec/' # where to save ds plots
if not os.path.exists(ds_dir):
os.system('mkdir '+ds_dir)
close('all')
# note: throughout, "LS" can also include C band, I initially wrote this code for 2015 data (which only has LS band)
# but it works for the 2013 data with LSC band
# params that can be changed are listed in default_fig_params
default_fig_params = {
'tint_P': 300,
'ds.phase does not update if you bin or truncate the time list, so make sure to do',\
'this last before plotting.'
params = {
'legend.fontsize': 'small',
'axes.titlesize': 'small',
'axes.labelsize': 'small',
'xtick.labelsize': 'x-small',
'ytick.labelsize': 'x-small',
'image.interpolation': 'nearest'
}
rcParams.update(params)
loadfile = '/users/jvillads/data_1/UVCet_paper/all_burst_epoch_dynspec_LSband.npy'
ds_list = load_dict(loadfile)
#loadfileP = '/data/jrv/burst_paper/all_burst_epoch_dynspec_Pband.npy'
#dsP_list = load_dict(loadfileP)
close('all')
# note: throughout, "LS" can also include C band, I initially wrote this code for 2015 data (which only has LS band)
# but it works for the 2013 data with LSC band
# params that can be changed are listed in default_fig_params
default_fig_params = {
'tint_P': 300,
'tint_LS': 60,
'df_MHz_P': 16,
'df_MHz_LS': 16,
'smax_P': None,