'''
# ============= Define the Mission ====================== #
This section is rather convoluted and should be cleaned up.
'''
# If you do a propagation, then you need to worry about dtval!
propagationParamFile = [] # Points to thrust profile for doing propagations
precomputedParamFile = 'AntaresNominalTrajectory.txt' # Points to file with precomputed profile for nominal trajectory
pathToMissionFiles = curFilePath # Kind of a holdover from a previous file structure
# Planet info
omegaE = 7.2921158494529352e-05 # rad/s
planetModel = 0 # 0 means spherical, 1 means elliptical
curMission = dict(propagationParamFile = propagationParamFile, precomputedParamFile = precomputedParamFile,
pathToMissionFiles = pathToMissionFiles, omegaE = omegaE, planetModel = planetModel)
curMission = TJC.InitializeMission(curMission)
TJC.SetupOutputFolders(curMission, tempDir, outputDir, vehicleName, launchLocation)
# If we're propagating trajectories, this needs to be set. Should do it in the input file, but this is fine for now.
curMission['useLoverD'] = False
curMission['loverd'] = 0.
# These hold files that need to be read in
# curMission['debrisCatPath'] = curMission['pathToMissionFiles'] + 'DebrisCatalog/'
# curMission['debrisCatFile'] = 'testAntaresCatalog.txt'
# curMission['debrisCatFile'] = 'Halcon9_1stNEW.txt'
curMission['debrisCatPath'] = debrisPath + 'Falcon9/'
curMission['debrisCatFile'] = 'Falcon9_FTS.txt'
curMission['atmospherePickle'] = rootDir + "data/AtmoProfiles/WallopsOverOcean.pkl"
'''
# If you do a propagation, then you need to worry about dtval!
propagationParamFile = [] # Points to thrust profile for doing propagations
precomputedParamFile = 'HTHL_Abridged.txt' # Points to file with precomputed profile for nominal trajectory
pathToMissionFiles = curFilePath # Kind of a holdover from a previous file structure
# Planet info
omegaE = 7.2921158494529352e-05 # rad/s
planetModel = 0 # 0 means spherical, 1 means elliptical
curMission = dict(propagationParamFile=propagationParamFile,
precomputedParamFile=precomputedParamFile,
pathToMissionFiles=pathToMissionFiles,
omegaE=omegaE,
planetModel=planetModel)
curMission = TJC.InitializeMission(curMission)
TJC.SetupOutputFolders(curMission, tempDir, outputDir, vehicleName,
launchLocation)
# If we're propagating trajectories, this needs to be set. Should do it in the input file, but this is fine for now.
curMission['useLoverD'] = False
curMission['loverd'] = 0.
# # These hold files that need to be read in
curMission[
'debrisCatPath'] = curMission['pathToMissionFiles'] + 'DebrisCatalog/'
curMission['debrisCatFile'] = 'LynxDebrisCatalog.txt'
curMission[
'atmospherePickle'] = rootDir + "data/AtmoProfiles/SpaceportAmerica.pkl"
'''
PROPAGATION and PROBABILITY PARAMETERS
# If you do a propagation, then you need to worry about dtval!
propagationParamFile = 'nominalParam_new.txt' # Points to thrust profile for doing propagations
precomputedParamFile = [
] # Points to file with precomputed profile for nominal trajectory
pathToMissionFiles = curFilePath # Kind of a holdover from a previous file structure
# Planet info
omegaE = 7.2921158494529352e-05 # rad/s
planetModel = 0 # 0 means spherical, 1 means elliptical
curMission = dict(propagationParamFile=propagationParamFile,
precomputedParamFile=precomputedParamFile,
pathToMissionFiles=pathToMissionFiles,
omegaE=omegaE,
planetModel=planetModel)
curMission = TJC.InitializeMission(curMission)
TJC.SetupOutputFolders(curMission, tempDir, outputDir, vehicleName,
launchLocation)
# If we're propagating trajectories, this needs to be set. Should do it in the input file, but this is fine for now.
curMission['useLoverD'] = False
curMission['loverd'] = 0.
# These hold files that need to be read in
curMission[
'debrisCatPath'] = curMission['pathToMissionFiles'] + 'DebrisCatalog/'
curMission['debrisCatFile'] = 'Halcon9_1stNEW.txt'
curMission['atmospherePickle'] = rootDir + "data/AtmoProfiles/Cape.pkl"
'''
PROPAGATION and PROBABILITY PARAMETERS
Set parameters related to:
'''
# ============= Define the Mission ====================== #
This section is rather convoluted and should be cleaned up.
'''
# If you do a propagation, then you need to worry about dtval!
propagationParamFile = [] # Points to thrust profile for doing propagations
precomputedParamFile = 'HTHLwCarrier_Abridged.txt' # Points to file with precomputed profile for nominal trajectory
pathToMissionFiles = curFilePath # Kind of a holdover from a previous file structure
# Planet info
omegaE = 7.2921158494529352e-05 # rad/s
planetModel = 0 # 0 means spherical, 1 means elliptical
curMission = dict(propagationParamFile = propagationParamFile, precomputedParamFile = precomputedParamFile,
pathToMissionFiles = pathToMissionFiles, omegaE = omegaE, planetModel = planetModel)
curMission = TJC.InitializeMission(curMission)
TJC.SetupOutputFolders(curMission, tempDir, outputDir, vehicleName, launchLocation)
# If we're propagating trajectories, this needs to be set. Should do it in the input file, but this is fine for now.
curMission['useLoverD'] = False
curMission['loverd'] = 0.
# # These hold files that need to be read in
curMission['debrisCatPath'] = curMission['pathToMissionFiles'] + 'DebrisCatalog/'
curMission['debrisCatFile'] = 'LynxDebrisCatalog.txt'
curMission['atmospherePickle'] = rootDir + "data/AtmoProfiles/SpaceportAmerica.pkl"
rootDir = os.path.abspath(curFilePath + "../../../") + "/"
outputDir = rootDir + "outputs/" # Where to store results, gitignored
# tempDir = rootDir + ct.tempFolderName # temp files here, gitignored
debrisPath = rootDir + "src/python/packages/DebrisCatalogs/"
from Simulation import TJC
curMission = dict()
curMission['atmospherePickle'] = rootDir + "data/AtmoProfiles/Cape.pkl"
curMission['numTrajSamples'] = 0
curMission[
'numWindSamples'] = 100 # Best results if this is a multiple of the number of nodes you're running on.
atmStorage, stateVecStorage, thetagStorage, tfailStorage = \
TJC.GenerateWindTrajProfiles(curMission, curMission['numTrajSamples'], curMission['numWindSamples'])
import numpy as np
from bokeh.plotting import figure, output_notebook, output_file, show
from bokeh.palettes import Spectral11
mypalette = Spectral11 * 10
# output to static HTML file
output_file("lines.html", title="line plot example")
p = figure(plot_width=1000, plot_height=1000, y_range=(0, 60), x_range=(0, 40))
p.title = "Winds"
p.title_text_font_size = "40pt"
for ix, curProfile in enumerate(atmStorage):
# Initialize the mission
propagationParamFile = [] # Points to thrust profile for doing propagations
precomputedParamFile = 'PM2_WSands_40kmBasedOnSonda3_COARSENED.txt' # Points to file with precomputed profile for nominal trajectory
pathToMissionFiles = './' # Kind of a holdover from a previous file structure
# Planet info
omegaE = 7.2921158494529352e-05 # rad/s
planetModel = 0 # 0 means spherical, 1 means elliptical
curMission = dict(propagationParamFile=propagationParamFile,
precomputedParamFile=precomputedParamFile,
pathToMissionFiles=pathToMissionFiles,
omegaE=omegaE,
planetModel=planetModel)
curMission = TJC.InitializeMission(curMission)
# If we're propagating trajectories, this needs to be set. Should do it in the input file, but this is fine for now.
curMission['useLoverD'] = False
curMission['loverd'] = 0.
'''
Defines the file structure for folders / files that are specific to this vehicle.
In theory, all vehicles will have the same file structure and this section would
be the same across all main files. Also, should bring pathToMissionFiles
into this section
'''
# These hold output files, create them in a moment if they don't already exist
curMission['GeneratedFilesFolder'] = curMission[
'pathToMissionFiles'] + 'GeneratedFiles/'
curMission['debrisPickleFolder'] = curMission[
'GeneratedFilesFolder'] + 'debrisPickleFolder'
# Initialize the mission
propagationParamFile = [] # Points to thrust profile for doing propagations
precomputedParamFile = 'ColumbiaMainPieceTrajectory.txt' # Points to file with precomputed profile for nominal trajectory
pathToMissionFiles = curFilePath # Kind of a holdover from a previous file structure
# Planet info
omegaE = 7.2921158494529352e-05 # rad/s
planetModel = initVec.planetModel # 0 means spherical, 1 means elliptical
curMission = dict(propagationParamFile=propagationParamFile,
precomputedParamFile=precomputedParamFile,
pathToMissionFiles=pathToMissionFiles,
omegaE=omegaE,
planetModel=planetModel)
curMission = TJC.InitializeMission(curMission)
TJC.SetupOutputFolders(curMission, tempDir, outputDir, vehicleName,
launchLocation)
# Use the same parameters that were used in the initVec that generated the main trajectories
if initVec.cloption == 0:
curMission['useLoverD'] = True
else:
curMission['useLoverD'] = False
curMission['loverd'] = initVec.loverd
# These hold files that need to be read in
curMission['debrisCatPath'] = debrisPath + 'Columbia/'
# curMission['debrisCatFile'] = 'testFileDistributed.txt'
# curMission['debrisCatFile'] = 'debugDistributed.txt'
curMission['debrisCatFile'] = 'debugColumbiaMarch16.txt'
# If you do a propagation, then you need to worry about dtval!
propagationParamFile = 'nominalParam_new.txt' # Points to thrust profile for doing propagations
precomputedParamFile = [
] # Points to file with precomputed profile for nominal trajectory
pathToMissionFiles = curFilePath # Kind of a holdover from a previous file structure
# Planet info
omegaE = 7.2921158494529352e-05 # rad/s
planetModel = 0 # 0 means spherical, 1 means elliptical
curMission = dict(propagationParamFile=propagationParamFile,
precomputedParamFile=precomputedParamFile,
pathToMissionFiles=pathToMissionFiles,
omegaE=omegaE,
planetModel=planetModel)
curMission = TJC.InitializeMission(curMission)
TJC.SetupOutputFolders(curMission, tempDir, outputDir, vehicleName,
launchLocation)
# If we're propagating trajectories, this needs to be set. Should do it in the input file, but this is fine for now.
curMission['useLoverD'] = False
curMission['loverd'] = 0.
# These hold files that need to be read in
curMission['debrisCatPath'] = debrisPath + 'Falcon9/'
curMission['debrisCatFile'] = 'Halcon9_1stNEW.txt'
curMission['atmospherePickle'] = rootDir + "data/AtmoProfiles/Cape.pkl"
'''
PROPAGATION and PROBABILITY PARAMETERS
Set parameters related to:
* ASH / density estimation
if len(sys.argv) != 2:
print "ERROR: Only one argument allowed. Must pass in .txt file containing atmo data"
sys.exit()
atmoFile = sys.argv[1]
baseName = atmoFile.split('.')[0]
pickleName = baseName + ".pkl"
from Simulation import TJC
# Note: You have to import all the frisco stuff because TJC demands to know about it.
# Create the atmosphere pickle. Should only need to be done once ever.
# atmoFile = 'Cape.txt'
# pickleName = 'Cape.pkl'
# atmoFile = 'WSands.txt'
# pickleName = 'WSands.pkl'
# atmoFile = 'FrontRange.txt'
# pickleName = 'FrontRange.pkl'
# atmoFile = 'SpaceportAmerica.txt'
# pickleName = 'SpaceportAmerica.pkl'
# atmoFile = 'WestTexas.txt'
# pickleName = 'WestTexas.pkl'
atmoFolder = './'
TJC.createAtmoPickle(atmoFolder, atmoFile, pickleName)
# Initialize the mission
propagationParamFile = [] # Points to thrust profile for doing propagations
precomputedParamFile = 'ReentryToOKwithLDpt3.txt' # Points to file with precomputed profile for nominal trajectory
pathToMissionFiles = curFilePath # Kind of a holdover from a previous file structure
# Planet info
omegaE = 7.2921158494529352e-05 # rad/s
planetModel = initVec.planetModel # 0 means spherical, 1 means elliptical
curMission = dict(propagationParamFile=propagationParamFile,
precomputedParamFile=precomputedParamFile,
pathToMissionFiles=pathToMissionFiles,
omegaE=omegaE,
planetModel=planetModel)
curMission = TJC.InitializeMission(curMission)
TJC.SetupOutputFolders(curMission, tempDir, outputDir, vehicleName,
launchLocation)
# Use the same parameters that were used in the initVec that generated the main trajectories
if initVec.cloption == 0:
curMission['useLoverD'] = True
else:
curMission['useLoverD'] = False
curMission['loverd'] = initVec.loverd
# These hold files that need to be read in
curMission[
'debrisCatPath'] = curMission['pathToMissionFiles'] + 'DebrisCatalog/'
curMission['debrisCatFile'] = 'testFileDistributed.txt'
curMission['atmospherePickle'] = rootDir + "data/AtmoProfiles/WestTexas.pkl"