import hbt ########## # Select which of the sequences we want to read in ########## # (NB: The spellling / capitalization is inconsistent in SAPNAME vs. VISITNAM. I have standardized it here.) sequence = 'O_RING_OC3' #sequence = 'O_RING_OC2' fs = 15 # Font size dir_images = '/Users/throop/Data/NH_Alice_Ring/' + sequence + '/data/pluto/level2/ali/all' file_tm = '/Users/throop/gv/dev/gv_kernels_new_horizons.txt' cspice.furnsh(file_tm) #============================================================================== # Read the Alice data #============================================================================== file_list = glob.glob(dir_images + '/*fit') met_all = [] # A long array with a list of all of the timestamps, one per 4 ms (i.e., at 250 hz) count_rate_fits_all = [] # The count rate as read from the COUNT_RATE extension directly count_rate_all = [] # Count rate computed from the PIXEL_LIST_TABLE. Should match that in COUNT_RATE extension count_rate_target_all = [] # Count rate for the target only, extracted by spatially filtering the PIXEL_LIST_TABLE d_target_summed = np.zeros((5,541)) d_summed = np.zeros((32,1024))
# HBT 15-Jun-2015
import cspice
print "Spice imported"
cspice.furnsh("/Users/throop/gv/dev/gv_kernels_new_horizons.txt")
UT = "2015 1 Jun 00:00:00"
et = cspice.utc2et(UT) # Works OK
print "ET = " + repr(et)
utc = cspice.et2utc(et, "C", 1) # Crashes the python kernel
# Works now after recompiling the kernel!
print "UTC = " + utc
#
# Syntax:
# python lorri_production_jupiter.py path/to/NHJULO_2001_supplemental_index.lbl
################################################################################
import oops
import oops.inst.nh.lorri as lorri
import numpy as np
import os, sys, traceback
import metadata_jupiter as meta
PLANET = 'JUPITER'
import cspice
cspice.furnsh('/Resources/SPICE/General/PCK/pck00010_edit_v01.tpc')
# This reload last overrides some bug that produces this error message:
# SPICE(INSUFFICIENTANGLES) -- TISBOD: Insufficient number of nutation/precession angles for body 599 at time 2.8472754946623E+07.
############################################
# Key parameters of run
############################################
execfile('COLUMNS_JUPITER.py') # define the columns, tiles, etc.
SAMPLING = 8 # pixel sampling density
SELECTION = 'S' # summary files only
############################################
# Construct the meshgrid for each LORRI FOV
############################################
plt.title('Shifted, dx=' + repr(dx) + ', dy = ' + repr(dy))
plt.show()
return t0
hbt.set_plot_defaults()
d2r = np.pi /180.
r2d = 1. / d2r
dir_data = '/Users/throop/data/NH_Jring/data/jupiter/level2/lor/all'
# Start up SPICE
file_tm = "/Users/throop/gv/dev/gv_kernels_new_horizons.txt"
cspice.furnsh(file_tm)
# Get the full list of files
file_list = glob.glob(dir_data + '/*fit')
files = np.array(file_list)
indices = np.argsort(file_list)
files = files[indices]
DO_FAST = True
if (DO_FAST):
files = files[0:100]
# Read the JD from each file. Then sort the files based on JD.
jd = []
# HBT 15-Jun-2015
import cspice
print "Spice imported"
cspice.furnsh("/Users/throop/gv/dev/gv_kernels_new_horizons.txt")
UT = "2015 1 Jun 00:00:00"
et = cspice.utc2et(UT) # Works OK
print "ET = " + repr(et)
utc = cspice.et2utc(et, "C", 1) # Crashes the python kernel
# Works now after recompiling the kernel!
print "UTC = " + utc