import planetaryconstants as constants from sim import Params import numpy as np import matplotlib.pyplot as plt import matplotlib as mpl from scipy import stats plt.close('all') mpl.rcParams["figure.autolayout"] = True mpl.rcParams.update({'font.size': 15}) # ============================================================================= # Load Mars constants # ============================================================================= params = Params() params.p = constants.MARS # ============================================================================= # Load data # ============================================================================= filename = '../results/AeroDrop_dispersed_500_0302193435.npz' data = np.load(filename, allow_pickle=True) xxvecArr_O = data['xxvecArr_O'] tvecArr_O = data['tvecArr_O'] sigdvecArr_O = data['sigdvecArr_O'] tsvecArr_O = data['tsvecArr_O'] BCList_O = data['BCList_O'] L_DList_O = data['L_DList_O'] gam0List = data['gam0List']
import numpy as np import time import matplotlib.pyplot as plt import sys import scipy.interpolate as interp import planetaryconstants as constants import ODEs from sim import Params, Outs, mainAD from atm import getRho_from_table tic = time.time() plt.close('all') ### CREATE params INPUT CLASS params = Params() params.p = constants.MARS params.returnTimeVectors = True params.atmMod = 'nom' ### INPUT ATM TABLE - GET ATM TABLE FROM GRAM DATA FILE params.dMode = 'table' filename = '../data/dens_Mars_nom.txt' # filename = '../data/dat_raw_Earth_nom.txt' atmdata = np.genfromtxt(filename, names=True) # atmdata.sort(order='Var_X') # put in ascending altitude order params.atmdat = np.array([atmdata['Var_X'], atmdata['DENSAV']]) # atmdata.sort(order='Hgtkm') # put in ascending altitude order # params.atmdat = np.array([atmdata['Hgtkm'], atmdata['DensMean']]) # alter density if requested
import numpy as np import time import matplotlib.pyplot as plt from datetime import datetime import sys import planetaryconstants as constants import ODEs from sim import Params, Outs, mainAD tic = time.time() plt.close('all') ### CREATE params INPUT CLASS params = Params() params.p = constants.TITAN params.returnTimeVectors = False # ============================================================================= # comment/uncomment the below blocks of code for desired scenario at Titan # ============================================================================= #### Make all BC and efpa ranges the same #### BCList = np.arange(10, 200, 2.5) efpaList = np.arange(-24.9, -35.4, -0.1) ### Lift-up, nominal atmosphere params.atmMod = 'nom' params.LD = 0.25 params.bank = 0 # deg modestr = 'lift-up, nominal atmosphere'
Created on Wed Feb 10 18:24:45 2021 @author: Samuel Albert """ from sim import Params import planetaryconstants as constants from atm import getRho_from_table import numpy as np import scipy.interpolate as interp import time ### CREATE params INPUT CLASS FOR NOMINAL VALUES params = Params() params.p = constants.MARS params.returnTimeVectors = True ### INPUT ATM TABLE - GET ATM TABLE FROM GRAM DATA FILE params.dMode = 'table' filename = '../data/dens_Mars_nom.txt' atmdata = np.genfromtxt(filename, names=True) params.atmdat = np.array([atmdata['Var_X'], atmdata['DENSAV']]) scipyfun = interp.interp1d(params.atmdat[0,:], params.atmdat[1,:]) def myinterp(h, params): ''' h is altitude in km. assumes atm data sorted in ascending order! also, will not work if h is outside of data range!
import numpy as np import time import matplotlib.pyplot as plt from datetime import datetime import sys import planetaryconstants as constants import ODEs from sim import Params, Outs, mainAD tic = time.time() plt.close('all') ### CREATE params INPUT CLASS params = Params() params.p = constants.EARTH params.returnTimeVectors = False # ============================================================================= # comment/uncomment the below blocks of code for desired scenario at Mars # ============================================================================= #### Make all BC and efpa ranges the same #### BCList = np.arange(10, 200, 2.5) efpaList = np.arange(-3.7, -7.4, -0.02) # efpaList = np.arange(-4.4, -8.2, -0.02) # for 12 km/s entry # ### Lift-up, nominal atmosphere # params.atmMod = 'nom' # params.LD = 0.25 # params.bank = 0 # deg
def getQoIParams(params): paramsQ = Params() paramsQ.CD = params.CD paramsQ.BC = params.BC paramsQ.m = params.m paramsQ.efpaWR = params.efpaWR paramsQ.vmagWR = params.vmagWR paramsQ.m_Y = params.m_Y paramsQ.CD_Y = params.CD_Y paramsQ.efpa_Y = params.efpa_Y paramsQ.vmag_Y = params.vmag_Y paramsQ.atm_Ys = params.atm_Ys return paramsQ
return paramsQ tic = time.time() datestring = datetime.now().strftime('%m%d%H%M%S') plt.close('all') # ============================================================================= # Number of Monte Carlo trials # ============================================================================= Nmc = 10000 # ============================================================================= # Create Params input class for Mars # ============================================================================= params = Params() params.p = constants.MARS params.returnTimeVectors = True # ============================================================================= # Generate eigenvalues, eigenvectors to use in KLE function # ============================================================================= params.dMode = 'fun' filename = '../data/Mars_0.1_50000.npz' alpha = 0.95 pfact = 1 # use GRAM dispersions directly evals, evecs, densSampMean, d, h = getEproblem(filename, alpha, pfact) params.evals = evals params.evecs = evecs
@author: Samuel Albert """ import numpy as np import matplotlib.pyplot as plt import planetaryconstants as constants from sim import Params from UQ import getEproblem, getKLEdensfun plt.close('all') # ============================================================================= # Create Params input class for Mars # ============================================================================= params = Params() params.p = constants.MARS params.returnTimeVectors = True # ============================================================================= # Load MarsGRAM density data # ============================================================================= filename = '../data/Mars_0.1_50000.npz' densdata = np.load(filename) densTot = densdata['densTot'] h = densdata['h'] densMean = densdata['densMean'] densCentered = densTot - densMean[:, None] # ============================================================================= # Generate eigenvalues, eigenvectors to use in KLE function
import ODEs from sim import Params, Outs, mainAD from atm import getMarsGRAMDensTable from conversions import RV2LLAEHV, VN2Vinf tic = time.time() plt.close('all') mpl.rcParams["figure.autolayout"] = True mpl.rcParams.update({'font.size': 16}) outsList = [] # ============================================================================= # Create Params input class for Mars # ============================================================================= params = Params() params.p = constants.MARS params.returnTimeVectors = True # ============================================================================= # Generate atm table for Nmc profiles # ============================================================================= filename = './../data/Mars_0.1_5000.txt' Nmc = 1 densAll, densMean, h = getMarsGRAMDensTable(filename, Nmc) # densAll, densMean, h = getMarsGRAMDensTableAll(filename, Nmc) params.dMode = 'table' # ============================================================================= # Load atm table for this run # =============================================================================