Python cspice 예제들

예제 #1

파일: extract_scs.py 프로젝트: vmerkin/LFM-helio

Omega = 2*pi/27.27/24./3600.   #  note, the sun rotation rate used in
                               #  the runs was 25.38 days. But here we
                               #  have to use 27.27 because this is
                               #  the rate at which Earth moves
                               #  relative to the frame rotating with
                               #  the sun
####################################################################################################
home = os.path.expanduser("~");
if os.path.join(home,'GitHub/Python/LFM-helio/lib/') not in sys.path: sys.path.append(os.path.join(home,'GitHub/Python/LFM-helio/lib/'))
if os.path.join(home,'GitHub/Python/LFM-helio/analysis/') not in sys.path: sys.path.append(os.path.join(home,'GitHub/Python/LFM-helio/analysis/'))
import lfmhlib,ace,cspice


# get ACE position
ACE_dati,ACE_x,ACE_y,ACE_z,ACE_r,ACE_theta,ACE_phi=cspice.cspice(ace_loc_file)

# GET THIS FOR OTHER SCS FROM MAP_SCS.PY
# get STEREO A position
STA_dati,STA_x,STA_y,STA_z,STA_r,STA_theta,STA_phi=cspice.cspice(sta_loc_file)

STB_dati,STB_x,STB_y,STB_z,STB_r,STB_theta,STB_phi=cspice.cspice(stb_loc_file)

# get MSGR position
MSGR_dati,MSGR_x,MSGR_y,MSGR_z,MSGR_r,MSGR_theta,MSGR_phi=cspice.cspice(msgr_loc_file)
####################################################################################################
lfmFiles = sorted(glob.glob(os.path.join(lfmFileDir,'*.hdf')))
start_ind = [start_step==int(f[-11:-4]) for f in lfmFiles].index(True)
dates = []
vr = []
br = []

예제 #2

파일: map_ace_2066.py 프로젝트: vmerkin/LFM-helio

if os.path.join(home,'GitHub/Python/LFM-helio/analysis/') not in sys.path: sys.path.append(os.path.join(home,'GitHub/Python/LFM-helio/analysis/'))
import lfmhlib,ace,cspice
import summary

####################################################################################################
########### CR 2066  ###############################################################################
####################################################################################################
lfmFileDir = '/glade/p/ujhb0003/LFM-helio/ADAPT/2066rf'
lfmFiles = sorted(glob.glob(os.path.join(lfmFileDir,'*.hdf')))
start_step = 25000
start_dati = datetime.datetime(2008,1,1,23,0,0)
nk=192
Omega = 2*pi/27.27/24./3600.

# get ACE position
ACE_dati,ACE_x,ACE_y,ACE_z,ACE_r,ACE_theta,ACE_phi=cspice.cspice('/glade/u/home/vgm/LFM-helio_2.0/SCDATA/ACE/Earth_cspice_HCI_position_2008_001_060')

# get MSGR position
MSGR_dati,MSGR_x,MSGR_y,MSGR_z,MSGR_r,MSGR_theta,MSGR_phi=cspice.cspice('/glade/u/home/vgm/LFM-helio_2.0/SCDATA/MSGR/MSGR_cspice_HCI_position_2008_001_060')

# get STEREO A position
STA_dati,STA_x,STA_y,STA_z,STA_r,STA_theta,STA_phi=cspice.cspice('/glade/u/home/vgm/LFM-helio_2.0/SCDATA/STEREO/STA_scpice_HCI_position_2008_001_060')

# get STEREO B position
STB_dati,STB_x,STB_y,STB_z,STB_r,STB_theta,STB_phi=cspice.cspice('/glade/u/home/vgm/LFM-helio_2.0/SCDATA/STEREO/STB_scpice_HCI_position_2008_001_060')

# get ACE data
ACE_mf_time,ACE_bx,ACE_by,ACE_bz=ace.cdaweb_mfi('/glade/u/home/vgm/LFM-helio_2.0/SCDATA/ACE/2066/AC_H2_MFI_20155.txt')
ACE_plasma_time,ACE_v,ACE_n=ace.cdaweb_swe('/glade/u/home/vgm/LFM-helio_2.0/SCDATA/ACE/2066/AC_H2_SWE_20155.txt')
####################################################################################################

예제 #3

파일: map_scs_islice.py 프로젝트: vmerkin/LFM-helio

Omega = 2*pi/27.27/24./3600.   #  note, the sun rotation rate used in
                               #  the runs was 25.38 days. But here we
                               #  have to use 27.27 because this is
                               #  the rate at which Earth moves
                               #  relative to the frame rotating with
                               #  the sun
####################################################################################################
home = os.path.expanduser("~");
if os.path.join(home,'GitHub/Python/LFM-helio/lib/') not in sys.path: sys.path.append(os.path.join(home,'GitHub/Python/LFM-helio/lib/'))
if os.path.join(home,'GitHub/Python/LFM-helio/analysis/') not in sys.path: sys.path.append(os.path.join(home,'GitHub/Python/LFM-helio/analysis/'))
import lfmhlib,ace,cspice

# get ACE position
ACE_dati,ACE_x,ACE_y,ACE_z,ACE_r,ACE_theta,ACE_phi=cspice.cspice(ace_loc_file)

# GET THIS FOR OTHER SCS FROM MAP_SCS.PY
# get STEREO A position
STA_dati,STA_x,STA_y,STA_z,STA_r,STA_theta,STA_phi=cspice.cspice(sta_loc_file)

STB_dati,STB_x,STB_y,STB_z,STB_r,STB_theta,STB_phi=cspice.cspice(stb_loc_file)

# get MSGR position
MSGR_dati,MSGR_x,MSGR_y,MSGR_z,MSGR_r,MSGR_theta,MSGR_phi=cspice.cspice(msgr_loc_file)
####################################################################################################
lfmFiles = sorted(glob.glob(os.path.join(lfmFileDir,'*.hdf')))
start_ind = [start_step==int(f[-11:-4]) for f in lfmFiles].index(True)

plt.style.use('publication')
plt.close('all')

예제 #4

import lfmhlib,ace,cspice

####################################################################################################
########### CR 2094 ################################################################################
####################################################################################################

lfmFileDir = '/glade/scratch/vgm/LFM-helio/ADAPT/2094rf'
lfmFiles = sorted(glob.glob(os.path.join(lfmFileDir,'*v1*.hdf')))
start_step = 25000
#start_file = 'cr2094_192x96x192_v1_mhd_0025000.hdf'
start_dati = datetime.datetime(2010,1,1,23,0,0)
nk=192
Omega = 2*pi/27.27/24./3600.

# get ACE position
ACE_dati,ACE_x,ACE_y,ACE_z,ACE_r,ACE_theta,ACE_phi=cspice.cspice('/glade/u/home/vgm/LFM-helio_2.0/SCDATA/ACE/Earth_cspice_HCI_position_2010_001_060')

# get MSGR position
MSGR_dati,MSGR_x,MSGR_y,MSGR_z,MSGR_r,MSGR_theta,MSGR_phi=cspice.cspice('/glade/u/home/vgm/LFM-helio_2.0/SCDATA/MSGR/MSGR_cspice_HCI_position_2010_001_060')

# get ACE data
ACE_mf_time,ACE_bx,ACE_by,ACE_bz=ace.cdaweb_mfi('/glade/u/home/vgm/LFM-helio_2.0/SCDATA/ACE/2094/AC_H2_MFI_14914.txt')
ACE_plasma_time,ACE_v,ACE_n=ace.cdaweb_swe('/glade/u/home/vgm/LFM-helio_2.0/SCDATA/ACE/2094/AC_H2_SWE_14914.txt')

####################################################################################################

####################################################################################################
########### CR 2066  ###############################################################################
####################################################################################################
# lfmFileDir = '/glade/scratch/vgm/LFM-helio/ADAPT/2066rf'
# lfmFiles = sorted(glob.glob(os.path.join(lfmFileDir,'*.hdf')))

예제 #5

파일: map_ace_2066.py 프로젝트: vmerkin/LFM-helio

import lfmhlib, ace, cspice
import summary

####################################################################################################
########### CR 2066  ###############################################################################
####################################################################################################
lfmFileDir = '/glade/p/ujhb0003/LFM-helio/ADAPT/2066rf'
lfmFiles = sorted(glob.glob(os.path.join(lfmFileDir, '*.hdf')))
start_step = 25000
start_dati = datetime.datetime(2008, 1, 1, 23, 0, 0)
nk = 192
Omega = 2 * pi / 27.27 / 24. / 3600.

# get ACE position
ACE_dati, ACE_x, ACE_y, ACE_z, ACE_r, ACE_theta, ACE_phi = cspice.cspice(
    '/glade/u/home/vgm/LFM-helio_2.0/SCDATA/ACE/Earth_cspice_HCI_position_2008_001_060'
)

# get MSGR position
MSGR_dati, MSGR_x, MSGR_y, MSGR_z, MSGR_r, MSGR_theta, MSGR_phi = cspice.cspice(
    '/glade/u/home/vgm/LFM-helio_2.0/SCDATA/MSGR/MSGR_cspice_HCI_position_2008_001_060'
)

# get STEREO A position
STA_dati, STA_x, STA_y, STA_z, STA_r, STA_theta, STA_phi = cspice.cspice(
    '/glade/u/home/vgm/LFM-helio_2.0/SCDATA/STEREO/STA_scpice_HCI_position_2008_001_060'
)

# get STEREO B position
STB_dati, STB_x, STB_y, STB_z, STB_r, STB_theta, STB_phi = cspice.cspice(
    '/glade/u/home/vgm/LFM-helio_2.0/SCDATA/STEREO/STB_scpice_HCI_position_2008_001_060'