09 - b - Uranus - Feasibility Charts - Drag

[1]:

from AMAT.planet import Planet
from AMAT.vehicle import Vehicle

import numpy as np
from scipy import interpolate

import matplotlib.pyplot as plt
from matplotlib import rcParams
from matplotlib.patches import Polygon
import os

[2]:

# Create a planet object for Titan
planet=Planet("URANUS")
planet.h_skip = 1000e3
planet.h_trap = 50e3

# Load an nominal atmospheric profile with height, temp, pressure, density data
planet.loadAtmosphereModel('../atmdata/Uranus/uranus-gram-avg.dat', 0 , 1 , 2, 3, heightInKmFlag=True)

[3]:

beta1 = 20.0

runID = 'uranus-drag-'

[9]:

vinf_kms_array  = np.linspace( 0.0,   30.0,   11)
betaRatio_array = np.linspace( 1.0,   41.0 ,  11)

[10]:

v0_kms_array    = np.zeros(len(vinf_kms_array))
v0_kms_array[:] = np.sqrt(1.0*(vinf_kms_array[:]*1E3)**2.0 + 2*np.ones(len(vinf_kms_array))*\
                          planet.GM/(planet.RP+1000.0*1.0E3))/1.0E3

overShootLimit_array  = np.zeros((len(v0_kms_array),len(betaRatio_array)))
underShootLimit_array = np.zeros((len(v0_kms_array),len(betaRatio_array)))
exitflag_os_array     = np.zeros((len(v0_kms_array),len(betaRatio_array)))
exitflag_us_array     = np.zeros((len(v0_kms_array),len(betaRatio_array)))
TCW_array             = np.zeros((len(v0_kms_array),len(betaRatio_array)))

[11]:

for i in range(0,len(v0_kms_array)):
    for j in range(0,len(betaRatio_array)):
        vehicle=Vehicle('DMVehicle', 300.0, beta1, 0.0, 3.1416, 0.0, 0.10, planet)
        vehicle.setInitialState(1000.0,0.0,0.0,v0_kms_array[i],0.0,-4.5,0.0,0.0)
        vehicle.setSolverParams(1E-6)
        vehicle.setDragModulationVehicleParams(beta1,betaRatio_array[j])

        underShootLimit_array[i,j], exitflag_us_array[i,j] = vehicle.findUnderShootLimitD(2400.0, 1.0, -20.0,-4.0,1E-10,4000)
        overShootLimit_array[i,j] , exitflag_os_array[i,j] = vehicle.findOverShootLimitD (2400.0, 1.0, -20.0,-4.0,1E-10,4000)

        TCW_array[i,j]    = overShootLimit_array[i,j] - underShootLimit_array[i,j]

        print('VINF: '+str(vinf_kms_array[i])+' km/s, BETA RATIO: '+str(betaRatio_array[j])+' TCW: '+str(TCW_array[i,j])+' deg.')

np.savetxt('../data/jsr-paper/uranus/'+runID+'vinf_kms_array.txt',vinf_kms_array)
np.savetxt('../data/jsr-paper/uranus/'+runID+'v0_kms_array.txt',v0_kms_array)
np.savetxt('../data/jsr-paper/uranus/'+runID+'betaRatio_array.txt',betaRatio_array)
np.savetxt('../data/jsr-paper/uranus/'+runID+'overShootLimit_array.txt',overShootLimit_array)
np.savetxt('../data/jsr-paper/uranus/'+runID+'exitflag_os_array.txt',exitflag_os_array)
np.savetxt('../data/jsr-paper/uranus/'+runID+'underShootLimit_array.txt',underShootLimit_array)
np.savetxt('../data/jsr-paper/uranus/'+runID+'exitflag_us_array.txt',exitflag_us_array)
np.savetxt('../data/jsr-paper/uranus/'+runID+'TCW_array.txt',TCW_array)

VINF: 0.0 km/s, BETA RATIO: 1.0 TCW: 0.0 deg.
VINF: 0.0 km/s, BETA RATIO: 5.0 TCW: 0.4300069945747964 deg.
VINF: 0.0 km/s, BETA RATIO: 9.0 TCW: 0.548236197559163 deg.
VINF: 0.0 km/s, BETA RATIO: 13.0 TCW: 0.6287807369953953 deg.
VINF: 0.0 km/s, BETA RATIO: 17.0 TCW: 0.6904995632939972 deg.
VINF: 0.0 km/s, BETA RATIO: 21.0 TCW: 0.7407522390130907 deg.
VINF: 0.0 km/s, BETA RATIO: 25.0 TCW: 0.7830369498115033 deg.
VINF: 0.0 km/s, BETA RATIO: 29.0 TCW: 0.8199181346572004 deg.
VINF: 0.0 km/s, BETA RATIO: 33.0 TCW: 0.8523810709593818 deg.
VINF: 0.0 km/s, BETA RATIO: 37.0 TCW: 0.8816194402170368 deg.
VINF: 0.0 km/s, BETA RATIO: 41.0 TCW: 0.9082424457883462 deg.
VINF: 3.0 km/s, BETA RATIO: 1.0 TCW: 0.0 deg.
VINF: 3.0 km/s, BETA RATIO: 5.0 TCW: 0.4306469949078746 deg.
VINF: 3.0 km/s, BETA RATIO: 9.0 TCW: 0.5533987928647548 deg.
VINF: 3.0 km/s, BETA RATIO: 13.0 TCW: 0.6368987987516448 deg.
VINF: 3.0 km/s, BETA RATIO: 17.0 TCW: 0.7001530813286081 deg.
VINF: 3.0 km/s, BETA RATIO: 21.0 TCW: 0.7520165814785287 deg.
VINF: 3.0 km/s, BETA RATIO: 25.0 TCW: 0.7956012652721256 deg.
VINF: 3.0 km/s, BETA RATIO: 29.0 TCW: 0.8335674138506874 deg.
VINF: 3.0 km/s, BETA RATIO: 33.0 TCW: 0.8669485827558674 deg.
VINF: 3.0 km/s, BETA RATIO: 37.0 TCW: 0.8972552987979725 deg.
VINF: 3.0 km/s, BETA RATIO: 41.0 TCW: 0.9245952178025618 deg.
VINF: 6.0 km/s, BETA RATIO: 1.0 TCW: 0.0 deg.
VINF: 6.0 km/s, BETA RATIO: 5.0 TCW: 0.43009640561649576 deg.
VINF: 6.0 km/s, BETA RATIO: 9.0 TCW: 0.5665283860289492 deg.
VINF: 6.0 km/s, BETA RATIO: 13.0 TCW: 0.6570447093108669 deg.
VINF: 6.0 km/s, BETA RATIO: 17.0 TCW: 0.7259380482137203 deg.
VINF: 6.0 km/s, BETA RATIO: 21.0 TCW: 0.7817554083303548 deg.
VINF: 6.0 km/s, BETA RATIO: 25.0 TCW: 0.829021891055163 deg.
VINF: 6.0 km/s, BETA RATIO: 29.0 TCW: 0.8698718391242437 deg.
VINF: 6.0 km/s, BETA RATIO: 33.0 TCW: 0.9059941348386928 deg.
VINF: 6.0 km/s, BETA RATIO: 37.0 TCW: 0.9384503098553978 deg.
VINF: 6.0 km/s, BETA RATIO: 41.0 TCW: 0.9679333599633537 deg.
VINF: 9.0 km/s, BETA RATIO: 1.0 TCW: 0.0 deg.
VINF: 9.0 km/s, BETA RATIO: 5.0 TCW: 0.4336301910225302 deg.
VINF: 9.0 km/s, BETA RATIO: 9.0 TCW: 0.5847309437813237 deg.
VINF: 9.0 km/s, BETA RATIO: 13.0 TCW: 0.6850503110326827 deg.
VINF: 9.0 km/s, BETA RATIO: 17.0 TCW: 0.7613813919597305 deg.
VINF: 9.0 km/s, BETA RATIO: 21.0 TCW: 0.8230569479637779 deg.
VINF: 9.0 km/s, BETA RATIO: 25.0 TCW: 0.8752881792606786 deg.
VINF: 9.0 km/s, BETA RATIO: 29.0 TCW: 0.9202921183314174 deg.
VINF: 9.0 km/s, BETA RATIO: 33.0 TCW: 0.9601127228233963 deg.
VINF: 9.0 km/s, BETA RATIO: 37.0 TCW: 0.9956420034868643 deg.
VINF: 9.0 km/s, BETA RATIO: 41.0 TCW: 1.0279079706524499 deg.
VINF: 12.0 km/s, BETA RATIO: 1.0 TCW: 0.0 deg.
VINF: 12.0 km/s, BETA RATIO: 5.0 TCW: 0.43706142978044227 deg.
VINF: 12.0 km/s, BETA RATIO: 9.0 TCW: 0.6043902909150347 deg.
VINF: 12.0 km/s, BETA RATIO: 13.0 TCW: 0.7161780579481274 deg.
VINF: 12.0 km/s, BETA RATIO: 17.0 TCW: 0.8005300595541485 deg.
VINF: 12.0 km/s, BETA RATIO: 21.0 TCW: 0.8689311987836845 deg.
VINF: 12.0 km/s, BETA RATIO: 25.0 TCW: 0.9263564213179052 deg.
VINF: 12.0 km/s, BETA RATIO: 29.0 TCW: 0.9758633595192805 deg.
VINF: 12.0 km/s, BETA RATIO: 33.0 TCW: 1.0195248335949145 deg.
VINF: 12.0 km/s, BETA RATIO: 37.0 TCW: 1.0581424981355667 deg.
VINF: 12.0 km/s, BETA RATIO: 41.0 TCW: 1.0930620817816816 deg.
VINF: 15.0 km/s, BETA RATIO: 1.0 TCW: 0.0 deg.
VINF: 15.0 km/s, BETA RATIO: 5.0 TCW: 0.4454347515129484 deg.
VINF: 15.0 km/s, BETA RATIO: 9.0 TCW: 0.6298921164707281 deg.
VINF: 15.0 km/s, BETA RATIO: 13.0 TCW: 0.7525014239363372 deg.
VINF: 15.0 km/s, BETA RATIO: 17.0 TCW: 0.8450111636775546 deg.
VINF: 15.0 km/s, BETA RATIO: 21.0 TCW: 0.9195918379700743 deg.
VINF: 15.0 km/s, BETA RATIO: 25.0 TCW: 0.982135801226832 deg.
VINF: 15.0 km/s, BETA RATIO: 29.0 TCW: 1.035611552069895 deg.
VINF: 15.0 km/s, BETA RATIO: 33.0 TCW: 1.0823744336958043 deg.
VINF: 15.0 km/s, BETA RATIO: 37.0 TCW: 1.1237224540673196 deg.
VINF: 15.0 km/s, BETA RATIO: 41.0 TCW: 1.1603765798499808 deg.
VINF: 18.0 km/s, BETA RATIO: 1.0 TCW: 0.0 deg.
VINF: 18.0 km/s, BETA RATIO: 5.0 TCW: 0.47190874099032953 deg.
VINF: 18.0 km/s, BETA RATIO: 9.0 TCW: 0.6728750808397308 deg.
VINF: 18.0 km/s, BETA RATIO: 13.0 TCW: 0.8054167323280126 deg.
VINF: 18.0 km/s, BETA RATIO: 17.0 TCW: 0.9050693272729404 deg.
VINF: 18.0 km/s, BETA RATIO: 21.0 TCW: 0.9851783454068936 deg.
VINF: 18.0 km/s, BETA RATIO: 25.0 TCW: 1.0515629358706065 deg.
VINF: 18.0 km/s, BETA RATIO: 29.0 TCW: 1.1083436924382113 deg.
VINF: 18.0 km/s, BETA RATIO: 33.0 TCW: 1.1569818182615563 deg.
VINF: 18.0 km/s, BETA RATIO: 37.0 TCW: 1.1988971406244673 deg.
VINF: 18.0 km/s, BETA RATIO: 41.0 TCW: 1.2348786691436544 deg.
VINF: 21.0 km/s, BETA RATIO: 1.0 TCW: 0.0 deg.
VINF: 21.0 km/s, BETA RATIO: 5.0 TCW: 0.5037200312362984 deg.
VINF: 21.0 km/s, BETA RATIO: 9.0 TCW: 0.7183635247056372 deg.
VINF: 21.0 km/s, BETA RATIO: 13.0 TCW: 0.8598193866200745 deg.
VINF: 21.0 km/s, BETA RATIO: 17.0 TCW: 0.9657282422995195 deg.
VINF: 21.0 km/s, BETA RATIO: 21.0 TCW: 1.0499617995228618 deg.
VINF: 21.0 km/s, BETA RATIO: 25.0 TCW: 1.1193449369748123 deg.
VINF: 21.0 km/s, BETA RATIO: 29.0 TCW: 1.177060860558413 deg.
VINF: 21.0 km/s, BETA RATIO: 33.0 TCW: 1.2252492255647667 deg.
VINF: 21.0 km/s, BETA RATIO: 37.0 TCW: 1.2646268653916195 deg.
VINF: 21.0 km/s, BETA RATIO: 41.0 TCW: 1.2969450384844095 deg.
VINF: 24.0 km/s, BETA RATIO: 1.0 TCW: 0.0 deg.
VINF: 24.0 km/s, BETA RATIO: 5.0 TCW: 0.5357726637739688 deg.
VINF: 24.0 km/s, BETA RATIO: 9.0 TCW: 0.7629505083314143 deg.
VINF: 24.0 km/s, BETA RATIO: 13.0 TCW: 0.9118634675396606 deg.
VINF: 24.0 km/s, BETA RATIO: 17.0 TCW: 1.0225540153332986 deg.
VINF: 24.0 km/s, BETA RATIO: 21.0 TCW: 1.1099410562310368 deg.
VINF: 24.0 km/s, BETA RATIO: 25.0 TCW: 1.1800493606715463 deg.
VINF: 24.0 km/s, BETA RATIO: 29.0 TCW: 1.236652796738781 deg.
VINF: 24.0 km/s, BETA RATIO: 33.0 TCW: 1.2806346639990807 deg.
VINF: 24.0 km/s, BETA RATIO: 37.0 TCW: 1.3166232906514779 deg.
VINF: 24.0 km/s, BETA RATIO: 41.0 TCW: 1.34690492216032 deg.
VINF: 27.0 km/s, BETA RATIO: 1.0 TCW: 0.0 deg.
VINF: 27.0 km/s, BETA RATIO: 5.0 TCW: 0.5671427182387561 deg.
VINF: 27.0 km/s, BETA RATIO: 9.0 TCW: 0.8048952496028505 deg.
VINF: 27.0 km/s, BETA RATIO: 13.0 TCW: 0.9600149451871403 deg.
VINF: 27.0 km/s, BETA RATIO: 17.0 TCW: 1.0744621621561237 deg.
VINF: 27.0 km/s, BETA RATIO: 21.0 TCW: 1.1627207996789366 deg.
VINF: 27.0 km/s, BETA RATIO: 25.0 TCW: 1.2315627416828647 deg.
VINF: 27.0 km/s, BETA RATIO: 29.0 TCW: 1.283228982181754 deg.
VINF: 27.0 km/s, BETA RATIO: 33.0 TCW: 1.3243737865122966 deg.
VINF: 27.0 km/s, BETA RATIO: 37.0 TCW: 1.358186729310546 deg.
VINF: 27.0 km/s, BETA RATIO: 41.0 TCW: 1.3864998817443848 deg.
VINF: 30.0 km/s, BETA RATIO: 1.0 TCW: 0.0 deg.
VINF: 30.0 km/s, BETA RATIO: 5.0 TCW: 0.5939707551151514 deg.
VINF: 30.0 km/s, BETA RATIO: 9.0 TCW: 0.8407557697501034 deg.
VINF: 30.0 km/s, BETA RATIO: 13.0 TCW: 1.0008656184654683 deg.
VINF: 30.0 km/s, BETA RATIO: 17.0 TCW: 1.1174983243690804 deg.
VINF: 30.0 km/s, BETA RATIO: 21.0 TCW: 1.2049253375153057 deg.
VINF: 30.0 km/s, BETA RATIO: 25.0 TCW: 1.2690666004200466 deg.
VINF: 30.0 km/s, BETA RATIO: 29.0 TCW: 1.3177535697468556 deg.
VINF: 30.0 km/s, BETA RATIO: 33.0 TCW: 1.3566493302350864 deg.
VINF: 30.0 km/s, BETA RATIO: 37.0 TCW: 1.388413718319498 deg.
VINF: 30.0 km/s, BETA RATIO: 41.0 TCW: 1.4149620372336358 deg.

[12]:

acc_net_g_max_array       = np.zeros((len(v0_kms_array),len(betaRatio_array)))
stag_pres_atm_max_array   = np.zeros((len(v0_kms_array),len(betaRatio_array)))
q_stag_total_max_array    = np.zeros((len(v0_kms_array),len(betaRatio_array)))
heatload_max_array        = np.zeros((len(v0_kms_array),len(betaRatio_array)))


for i in range(0,len(v0_kms_array)):
    for j in range(0,len(betaRatio_array)):
        vehicle=Vehicle('DMVehicle', 300, beta1, 0.0, 3.1416, 0.0, 0.10, planet)
        vehicle.setInitialState(1000.0,0.0,0.0,v0_kms_array[i],0.0,overShootLimit_array[i,j],0.0,0.0)
        vehicle.setSolverParams(1E-6)

        vehicle.propogateEntry (2400.0, 1.0, 0.0)

        # Extract and save variables to plot
        t_min_os         = vehicle.t_minc
        h_km_os          = vehicle.h_kmc
        acc_net_g_os     = vehicle.acc_net_g
        q_stag_con_os    = vehicle.q_stag_con
        q_stag_rad_os    = vehicle.q_stag_rad
        rc_os            = vehicle.rc
        vc_os            = vehicle.vc
        stag_pres_atm_os = vehicle.computeStagPres(rc_os,vc_os)/(1.01325E5)
        heatload_os      = vehicle.heatload


        vehicle=Vehicle('DMVehicle', 300.0, beta1, 0.0, 3.1416, 0.0, 0.10, planet)
        vehicle.setInitialState(1000.0,0.0,0.0,v0_kms_array[i],0.0,underShootLimit_array[i,j],0.0,0.0)
        vehicle.setSolverParams( 1E-6)

        vehicle.propogateEntry (2400.0, 1.0, 0.0)

        # Extract and save variable to plot
        t_min_us         = vehicle.t_minc
        h_km_us          = vehicle.h_kmc
        acc_net_g_us     = vehicle.acc_net_g
        q_stag_con_us    = vehicle.q_stag_con
        q_stag_rad_us    = vehicle.q_stag_rad
        rc_us            = vehicle.rc
        vc_us            = vehicle.vc
        stag_pres_atm_us = vehicle.computeStagPres(rc_us,vc_us)/(1.01325E5)
        heatload_us      = vehicle.heatload

        q_stag_total_os  = q_stag_con_os + q_stag_rad_os
        q_stag_total_us  = q_stag_con_us + q_stag_rad_us

        acc_net_g_max_array[i,j]      = max(max(acc_net_g_os),max(acc_net_g_os))
        stag_pres_atm_max_array[i,j]  = max(max(stag_pres_atm_os),max(stag_pres_atm_os))
        q_stag_total_max_array[i,j]   = max(max(q_stag_total_os),max(q_stag_total_us))
        heatload_max_array[i,j]       = max(max(heatload_os),max(heatload_os))

        print("V_infty: "+str(vinf_kms_array[i])+" km/s"+", BR: "+str(betaRatio_array[j])+" G_MAX: "+str(acc_net_g_max_array[i,j])+" QDOT_MAX: "+str(q_stag_total_max_array[i,j])+" J_MAX: "+str(heatload_max_array[i,j])+" STAG. PRES: "+str(stag_pres_atm_max_array[i,j]))

np.savetxt('../data/jsr-paper/uranus/'+runID+'acc_net_g_max_array.txt',acc_net_g_max_array)
np.savetxt('../data/jsr-paper/uranus/'+runID+'stag_pres_atm_max_array.txt',stag_pres_atm_max_array)
np.savetxt('../data/jsr-paper/uranus/'+runID+'q_stag_total_max_array.txt',q_stag_total_max_array)
np.savetxt('../data/jsr-paper/uranus/'+runID+'heatload_max_array.txt',heatload_max_array)

V_infty: 0.0 km/s, BR: 1.0 G_MAX: 1.8278605124542677 QDOT_MAX: 271.63896100022583 J_MAX: 251148.25410870728 STAG. PRES: 0.003553624041011706
V_infty: 0.0 km/s, BR: 5.0 G_MAX: 1.8278605124542677 QDOT_MAX: 435.71346317409103 J_MAX: 251148.25410870728 STAG. PRES: 0.003553624041011706
V_infty: 0.0 km/s, BR: 9.0 G_MAX: 1.8278605124542677 QDOT_MAX: 479.64582633787853 J_MAX: 251148.25410870728 STAG. PRES: 0.003553624041011706
V_infty: 0.0 km/s, BR: 13.0 G_MAX: 1.8278605124542677 QDOT_MAX: 504.58188910669907 J_MAX: 251148.25410870728 STAG. PRES: 0.003553624041011706
V_infty: 0.0 km/s, BR: 17.0 G_MAX: 1.8278605124542677 QDOT_MAX: 520.1699448478237 J_MAX: 251148.25410870728 STAG. PRES: 0.003553624041011706
V_infty: 0.0 km/s, BR: 21.0 G_MAX: 1.8278605124542677 QDOT_MAX: 531.6403839661771 J_MAX: 251148.25410870728 STAG. PRES: 0.003553624041011706
V_infty: 0.0 km/s, BR: 25.0 G_MAX: 1.8278605124542677 QDOT_MAX: 540.1466279922666 J_MAX: 251148.25410870728 STAG. PRES: 0.003553624041011706
V_infty: 0.0 km/s, BR: 29.0 G_MAX: 1.8278605124542677 QDOT_MAX: 546.7376408682949 J_MAX: 251148.25410870728 STAG. PRES: 0.003553624041011706
V_infty: 0.0 km/s, BR: 33.0 G_MAX: 1.8278605124542677 QDOT_MAX: 552.6315994156959 J_MAX: 251148.25410870728 STAG. PRES: 0.003553624041011706
V_infty: 0.0 km/s, BR: 37.0 G_MAX: 1.8278605124542677 QDOT_MAX: 557.110194331954 J_MAX: 251148.25410870728 STAG. PRES: 0.003553624041011706
V_infty: 0.0 km/s, BR: 41.0 G_MAX: 1.8278605124542677 QDOT_MAX: 561.3700996211214 J_MAX: 251148.25410870728 STAG. PRES: 0.003553624041011706
V_infty: 3.0 km/s, BR: 1.0 G_MAX: 1.813509125837342 QDOT_MAX: 285.27881234617496 J_MAX: 254310.6625966577 STAG. PRES: 0.0035257034148754485
V_infty: 3.0 km/s, BR: 5.0 G_MAX: 1.813509125837342 QDOT_MAX: 453.72407122074287 J_MAX: 254310.6625966577 STAG. PRES: 0.0035257034148754485
V_infty: 3.0 km/s, BR: 9.0 G_MAX: 1.813509125837342 QDOT_MAX: 500.0958186744359 J_MAX: 254310.6625966577 STAG. PRES: 0.0035257034148754485
V_infty: 3.0 km/s, BR: 13.0 G_MAX: 1.813509125837342 QDOT_MAX: 525.7400282935042 J_MAX: 254310.6625966577 STAG. PRES: 0.0035257034148754485
V_infty: 3.0 km/s, BR: 17.0 G_MAX: 1.813509125837342 QDOT_MAX: 541.6476177013972 J_MAX: 254310.6625966577 STAG. PRES: 0.0035257034148754485
V_infty: 3.0 km/s, BR: 21.0 G_MAX: 1.813509125837342 QDOT_MAX: 553.0037740228105 J_MAX: 254310.6625966577 STAG. PRES: 0.0035257034148754485
V_infty: 3.0 km/s, BR: 25.0 G_MAX: 1.813509125837342 QDOT_MAX: 561.5656671054957 J_MAX: 254310.6625966577 STAG. PRES: 0.0035257034148754485
V_infty: 3.0 km/s, BR: 29.0 G_MAX: 1.813509125837342 QDOT_MAX: 568.412971091998 J_MAX: 254310.6625966577 STAG. PRES: 0.0035257034148754485
V_infty: 3.0 km/s, BR: 33.0 G_MAX: 1.813509125837342 QDOT_MAX: 574.2790466667215 J_MAX: 254310.6625966577 STAG. PRES: 0.0035257034148754485
V_infty: 3.0 km/s, BR: 37.0 G_MAX: 1.813509125837342 QDOT_MAX: 578.989663885992 J_MAX: 254310.6625966577 STAG. PRES: 0.0035257034148754485
V_infty: 3.0 km/s, BR: 41.0 G_MAX: 1.813509125837342 QDOT_MAX: 583.235087172835 J_MAX: 254310.6625966577 STAG. PRES: 0.0035257034148754485
V_infty: 6.0 km/s, BR: 1.0 G_MAX: 1.9239973748466965 QDOT_MAX: 327.2763452091848 J_MAX: 264002.34462654754 STAG. PRES: 0.003740432746181789
V_infty: 6.0 km/s, BR: 5.0 G_MAX: 1.9239973748466965 QDOT_MAX: 509.3556477521291 J_MAX: 264002.34462654754 STAG. PRES: 0.003740432746181789
V_infty: 6.0 km/s, BR: 9.0 G_MAX: 1.9239973748466965 QDOT_MAX: 562.4487767763521 J_MAX: 264002.34462654754 STAG. PRES: 0.003740432746181789
V_infty: 6.0 km/s, BR: 13.0 G_MAX: 1.9239973748466965 QDOT_MAX: 589.3801338713673 J_MAX: 264002.34462654754 STAG. PRES: 0.003740432746181789
V_infty: 6.0 km/s, BR: 17.0 G_MAX: 1.9239973748466965 QDOT_MAX: 605.8638333199464 J_MAX: 264002.34462654754 STAG. PRES: 0.003740432746181789
V_infty: 6.0 km/s, BR: 21.0 G_MAX: 1.9239973748466965 QDOT_MAX: 617.6453709454064 J_MAX: 264002.34462654754 STAG. PRES: 0.003740432746181789
V_infty: 6.0 km/s, BR: 25.0 G_MAX: 1.9239973748466965 QDOT_MAX: 626.4498205844337 J_MAX: 264002.34462654754 STAG. PRES: 0.003740432746181789
V_infty: 6.0 km/s, BR: 29.0 G_MAX: 1.9239973748466965 QDOT_MAX: 633.2704677685279 J_MAX: 264002.34462654754 STAG. PRES: 0.003740432746181789
V_infty: 6.0 km/s, BR: 33.0 G_MAX: 1.9239973748466965 QDOT_MAX: 639.2913984541158 J_MAX: 264002.34462654754 STAG. PRES: 0.003740432746181789
V_infty: 6.0 km/s, BR: 37.0 G_MAX: 1.9239973748466965 QDOT_MAX: 644.1799308185472 J_MAX: 264002.34462654754 STAG. PRES: 0.003740432746181789
V_infty: 6.0 km/s, BR: 41.0 G_MAX: 1.9239973748466965 QDOT_MAX: 648.0977508314228 J_MAX: 264002.34462654754 STAG. PRES: 0.003740432746181789
V_infty: 9.0 km/s, BR: 1.0 G_MAX: 2.0258130250579387 QDOT_MAX: 400.25377250048905 J_MAX: 278489.1582496572 STAG. PRES: 0.003938281787666398
V_infty: 9.0 km/s, BR: 5.0 G_MAX: 2.0258130250579387 QDOT_MAX: 609.3306408965296 J_MAX: 278489.1582496572 STAG. PRES: 0.003938281787666398
V_infty: 9.0 km/s, BR: 9.0 G_MAX: 2.0258130250579387 QDOT_MAX: 668.5018256552937 J_MAX: 278489.1582496572 STAG. PRES: 0.003938281787666398
V_infty: 9.0 km/s, BR: 13.0 G_MAX: 2.0258130250579387 QDOT_MAX: 696.3639900479758 J_MAX: 278489.1582496572 STAG. PRES: 0.003938281787666398
V_infty: 9.0 km/s, BR: 17.0 G_MAX: 2.0258130250579387 QDOT_MAX: 713.3785806026733 J_MAX: 278489.1582496572 STAG. PRES: 0.003938281787666398
V_infty: 9.0 km/s, BR: 21.0 G_MAX: 2.0258130250579387 QDOT_MAX: 724.9975795128324 J_MAX: 278489.1582496572 STAG. PRES: 0.003938281787666398
V_infty: 9.0 km/s, BR: 25.0 G_MAX: 2.0258130250579387 QDOT_MAX: 734.2923884518979 J_MAX: 278489.1582496572 STAG. PRES: 0.003938281787666398
V_infty: 9.0 km/s, BR: 29.0 G_MAX: 2.0258130250579387 QDOT_MAX: 740.974411219486 J_MAX: 278489.1582496572 STAG. PRES: 0.003938281787666398
V_infty: 9.0 km/s, BR: 33.0 G_MAX: 2.0258130250579387 QDOT_MAX: 746.8216917769777 J_MAX: 278489.1582496572 STAG. PRES: 0.003938281787666398
V_infty: 9.0 km/s, BR: 37.0 G_MAX: 2.0258130250579387 QDOT_MAX: 752.1778747450783 J_MAX: 278489.1582496572 STAG. PRES: 0.003938281787666398
V_infty: 9.0 km/s, BR: 41.0 G_MAX: 2.0258130250579387 QDOT_MAX: 756.2828367440743 J_MAX: 278489.1582496572 STAG. PRES: 0.003938281787666398
V_infty: 12.0 km/s, BR: 1.0 G_MAX: 2.5189076325782014 QDOT_MAX: 506.5215997712668 J_MAX: 296835.90644155553 STAG. PRES: 0.004886858063109952
V_infty: 12.0 km/s, BR: 5.0 G_MAX: 2.5189076325782014 QDOT_MAX: 756.9382577972025 J_MAX: 296835.90644155553 STAG. PRES: 0.004886858063109952
V_infty: 12.0 km/s, BR: 9.0 G_MAX: 2.5189076325782014 QDOT_MAX: 819.0144172542351 J_MAX: 296835.90644155553 STAG. PRES: 0.004886858063109952
V_infty: 12.0 km/s, BR: 13.0 G_MAX: 2.5189076325782014 QDOT_MAX: 847.0458882792574 J_MAX: 296835.90644155553 STAG. PRES: 0.004886858063109952
V_infty: 12.0 km/s, BR: 17.0 G_MAX: 2.5189076325782014 QDOT_MAX: 864.3742549780326 J_MAX: 296835.90644155553 STAG. PRES: 0.004886858063109952
V_infty: 12.0 km/s, BR: 21.0 G_MAX: 2.5189076325782014 QDOT_MAX: 876.6277585268331 J_MAX: 296835.90644155553 STAG. PRES: 0.004886858063109952
V_infty: 12.0 km/s, BR: 25.0 G_MAX: 2.5189076325782014 QDOT_MAX: 886.3051663733083 J_MAX: 296835.90644155553 STAG. PRES: 0.004886858063109952
V_infty: 12.0 km/s, BR: 29.0 G_MAX: 2.5189076325782014 QDOT_MAX: 893.9000127633334 J_MAX: 296835.90644155553 STAG. PRES: 0.004886858063109952
V_infty: 12.0 km/s, BR: 33.0 G_MAX: 2.5189076325782014 QDOT_MAX: 900.1261414641633 J_MAX: 296835.90644155553 STAG. PRES: 0.004886858063109952
V_infty: 12.0 km/s, BR: 37.0 G_MAX: 2.5189076325782014 QDOT_MAX: 905.897535259384 J_MAX: 296835.90644155553 STAG. PRES: 0.004886858063109952
V_infty: 12.0 km/s, BR: 41.0 G_MAX: 2.5189076325782014 QDOT_MAX: 912.1255446132999 J_MAX: 296835.90644155553 STAG. PRES: 0.004886858063109952
V_infty: 15.0 km/s, BR: 1.0 G_MAX: 3.6836563647030673 QDOT_MAX: 645.2381706851423 J_MAX: 317919.7825151407 STAG. PRES: 0.007145096565839003
V_infty: 15.0 km/s, BR: 5.0 G_MAX: 3.6836563647030673 QDOT_MAX: 953.6745584519241 J_MAX: 317919.7825151407 STAG. PRES: 0.007145096565839003
V_infty: 15.0 km/s, BR: 9.0 G_MAX: 3.6836563647030673 QDOT_MAX: 1017.1843569444363 J_MAX: 317919.7825151407 STAG. PRES: 0.007145096565839003
V_infty: 15.0 km/s, BR: 13.0 G_MAX: 3.6836563647030673 QDOT_MAX: 1046.1243324464817 J_MAX: 317919.7825151407 STAG. PRES: 0.007145096565839003
V_infty: 15.0 km/s, BR: 17.0 G_MAX: 3.6836563647030673 QDOT_MAX: 1064.9488358779624 J_MAX: 317919.7825151407 STAG. PRES: 0.007145096565839003
V_infty: 15.0 km/s, BR: 21.0 G_MAX: 3.6836563647030673 QDOT_MAX: 1077.9082122746656 J_MAX: 317919.7825151407 STAG. PRES: 0.007145096565839003
V_infty: 15.0 km/s, BR: 25.0 G_MAX: 3.6836563647030673 QDOT_MAX: 1089.8843744095432 J_MAX: 317919.7825151407 STAG. PRES: 0.007145096565839003
V_infty: 15.0 km/s, BR: 29.0 G_MAX: 3.6836563647030673 QDOT_MAX: 1101.1681341272536 J_MAX: 317919.7825151407 STAG. PRES: 0.007145096565839003
V_infty: 15.0 km/s, BR: 33.0 G_MAX: 3.6836563647030673 QDOT_MAX: 1110.858375608081 J_MAX: 317919.7825151407 STAG. PRES: 0.007145096565839003
V_infty: 15.0 km/s, BR: 37.0 G_MAX: 3.6836563647030673 QDOT_MAX: 1119.3099028026904 J_MAX: 317919.7825151407 STAG. PRES: 0.007145096565839003
V_infty: 15.0 km/s, BR: 41.0 G_MAX: 3.6836563647030673 QDOT_MAX: 1126.680541361289 J_MAX: 317919.7825151407 STAG. PRES: 0.007145096565839003
V_infty: 18.0 km/s, BR: 1.0 G_MAX: 5.426536661833298 QDOT_MAX: 846.7541409867274 J_MAX: 341452.1706464557 STAG. PRES: 0.010522977092040508
V_infty: 18.0 km/s, BR: 5.0 G_MAX: 5.426536661833298 QDOT_MAX: 1202.7713958257084 J_MAX: 341452.1706464557 STAG. PRES: 0.010522977092040508
V_infty: 18.0 km/s, BR: 9.0 G_MAX: 5.426536661833298 QDOT_MAX: 1269.3041890068848 J_MAX: 341452.1706464557 STAG. PRES: 0.010522977092040508
V_infty: 18.0 km/s, BR: 13.0 G_MAX: 5.426536661833298 QDOT_MAX: 1300.5053341718078 J_MAX: 341452.1706464557 STAG. PRES: 0.010522977092040508
V_infty: 18.0 km/s, BR: 17.0 G_MAX: 5.426536661833298 QDOT_MAX: 1322.5482784091985 J_MAX: 341452.1706464557 STAG. PRES: 0.010522977092040508
V_infty: 18.0 km/s, BR: 21.0 G_MAX: 5.426536661833298 QDOT_MAX: 1343.0843053141464 J_MAX: 341452.1706464557 STAG. PRES: 0.010522977092040508
V_infty: 18.0 km/s, BR: 25.0 G_MAX: 5.426536661833298 QDOT_MAX: 1359.6247359709578 J_MAX: 341452.1706464557 STAG. PRES: 0.010522977092040508
V_infty: 18.0 km/s, BR: 29.0 G_MAX: 5.426536661833298 QDOT_MAX: 1373.619805964247 J_MAX: 341452.1706464557 STAG. PRES: 0.010522977092040508
V_infty: 18.0 km/s, BR: 33.0 G_MAX: 5.426536661833298 QDOT_MAX: 1385.2337491766125 J_MAX: 341452.1706464557 STAG. PRES: 0.010522977092040508
V_infty: 18.0 km/s, BR: 37.0 G_MAX: 5.426536661833298 QDOT_MAX: 1395.200072332318 J_MAX: 341452.1706464557 STAG. PRES: 0.010522977092040508
V_infty: 18.0 km/s, BR: 41.0 G_MAX: 5.426536661833298 QDOT_MAX: 1403.3071515720953 J_MAX: 341452.1706464557 STAG. PRES: 0.010522977092040508
V_infty: 21.0 km/s, BR: 1.0 G_MAX: 7.573289638134486 QDOT_MAX: 1109.2358111938918 J_MAX: 368407.0119218065 STAG. PRES: 0.014682577162996762
V_infty: 21.0 km/s, BR: 5.0 G_MAX: 7.573289638134486 QDOT_MAX: 1516.6279209877707 J_MAX: 368407.0119218065 STAG. PRES: 0.014682577162996762
V_infty: 21.0 km/s, BR: 9.0 G_MAX: 7.573289638134486 QDOT_MAX: 1589.5153466287386 J_MAX: 368407.0119218065 STAG. PRES: 0.014682577162996762
V_infty: 21.0 km/s, BR: 13.0 G_MAX: 7.573289638134486 QDOT_MAX: 1627.0660644783588 J_MAX: 368407.0119218065 STAG. PRES: 0.014682577162996762
V_infty: 21.0 km/s, BR: 17.0 G_MAX: 7.573289638134486 QDOT_MAX: 1658.4778851537048 J_MAX: 368407.0119218065 STAG. PRES: 0.014682577162996762
V_infty: 21.0 km/s, BR: 21.0 G_MAX: 7.573289638134486 QDOT_MAX: 1683.5658170619477 J_MAX: 368407.0119218065 STAG. PRES: 0.014682577162996762
V_infty: 21.0 km/s, BR: 25.0 G_MAX: 7.573289638134486 QDOT_MAX: 1703.7330065643607 J_MAX: 368407.0119218065 STAG. PRES: 0.014682577162996762
V_infty: 21.0 km/s, BR: 29.0 G_MAX: 7.573289638134486 QDOT_MAX: 1720.0423825436449 J_MAX: 368407.0119218065 STAG. PRES: 0.014682577162996762
V_infty: 21.0 km/s, BR: 33.0 G_MAX: 7.573289638134486 QDOT_MAX: 1733.48294170645 J_MAX: 368407.0119218065 STAG. PRES: 0.014682577162996762
V_infty: 21.0 km/s, BR: 37.0 G_MAX: 7.573289638134486 QDOT_MAX: 1743.9092978600181 J_MAX: 368407.0119218065 STAG. PRES: 0.014682577162996762
V_infty: 21.0 km/s, BR: 41.0 G_MAX: 7.573289638134486 QDOT_MAX: 1752.617310778698 J_MAX: 368407.0119218065 STAG. PRES: 0.014682577162996762
V_infty: 24.0 km/s, BR: 1.0 G_MAX: 10.164674157032687 QDOT_MAX: 1445.2283791763298 J_MAX: 401761.6388641296 STAG. PRES: 0.019702790245935957
V_infty: 24.0 km/s, BR: 5.0 G_MAX: 10.164674157032687 QDOT_MAX: 1923.6465719753949 J_MAX: 401761.6388641296 STAG. PRES: 0.019702790245935957
V_infty: 24.0 km/s, BR: 9.0 G_MAX: 10.164674157032687 QDOT_MAX: 2012.320255331466 J_MAX: 401761.6388641296 STAG. PRES: 0.019702790245935957
V_infty: 24.0 km/s, BR: 13.0 G_MAX: 10.164674157032687 QDOT_MAX: 2057.3052051223267 J_MAX: 401761.6388641296 STAG. PRES: 0.019702790245935957
V_infty: 24.0 km/s, BR: 17.0 G_MAX: 10.164674157032687 QDOT_MAX: 2093.4961750458774 J_MAX: 401761.6388641296 STAG. PRES: 0.019702790245935957
V_infty: 24.0 km/s, BR: 21.0 G_MAX: 10.164674157032687 QDOT_MAX: 2123.3810490858928 J_MAX: 401761.6388641296 STAG. PRES: 0.019702790245935957
V_infty: 24.0 km/s, BR: 25.0 G_MAX: 10.164674157032687 QDOT_MAX: 2147.1521790826096 J_MAX: 401761.6388641296 STAG. PRES: 0.019702790245935957
V_infty: 24.0 km/s, BR: 29.0 G_MAX: 10.164674157032687 QDOT_MAX: 2165.6502942315856 J_MAX: 401761.6388641296 STAG. PRES: 0.019702790245935957
V_infty: 24.0 km/s, BR: 33.0 G_MAX: 10.164674157032687 QDOT_MAX: 2179.7459047020434 J_MAX: 401761.6388641296 STAG. PRES: 0.019702790245935957
V_infty: 24.0 km/s, BR: 37.0 G_MAX: 10.164674157032687 QDOT_MAX: 2190.964736799414 J_MAX: 401761.6388641296 STAG. PRES: 0.019702790245935957
V_infty: 24.0 km/s, BR: 41.0 G_MAX: 10.164674157032687 QDOT_MAX: 2200.379990292354 J_MAX: 401761.6388641296 STAG. PRES: 0.019702790245935957
V_infty: 27.0 km/s, BR: 1.0 G_MAX: 13.231949270350407 QDOT_MAX: 1905.7084344396844 J_MAX: 447065.25571461295 STAG. PRES: 0.025643897044314443
V_infty: 27.0 km/s, BR: 5.0 G_MAX: 13.231949270350407 QDOT_MAX: 2494.8637498622566 J_MAX: 447065.25571461295 STAG. PRES: 0.025643897044314443
V_infty: 27.0 km/s, BR: 9.0 G_MAX: 13.231949270350407 QDOT_MAX: 2622.9345131712353 J_MAX: 447065.25571461295 STAG. PRES: 0.025643897044314443
V_infty: 27.0 km/s, BR: 13.0 G_MAX: 13.231949270350407 QDOT_MAX: 2685.7417146694183 J_MAX: 447065.25571461295 STAG. PRES: 0.025643897044314443
V_infty: 27.0 km/s, BR: 17.0 G_MAX: 13.231949270350407 QDOT_MAX: 2724.0702343817607 J_MAX: 447065.25571461295 STAG. PRES: 0.025643897044314443
V_infty: 27.0 km/s, BR: 21.0 G_MAX: 13.231949270350407 QDOT_MAX: 2756.6469392507915 J_MAX: 447065.25571461295 STAG. PRES: 0.025643897044314443
V_infty: 27.0 km/s, BR: 25.0 G_MAX: 13.231949270350407 QDOT_MAX: 2775.92697999241 J_MAX: 447065.25571461295 STAG. PRES: 0.025643897044314443
V_infty: 27.0 km/s, BR: 29.0 G_MAX: 13.231949270350407 QDOT_MAX: 2792.0420524520086 J_MAX: 447065.25571461295 STAG. PRES: 0.025643897044314443
V_infty: 27.0 km/s, BR: 33.0 G_MAX: 13.231949270350407 QDOT_MAX: 2803.6854909597773 J_MAX: 447065.25571461295 STAG. PRES: 0.025643897044314443
V_infty: 27.0 km/s, BR: 37.0 G_MAX: 13.231949270350407 QDOT_MAX: 2815.738320993501 J_MAX: 447065.25571461295 STAG. PRES: 0.025643897044314443
V_infty: 27.0 km/s, BR: 41.0 G_MAX: 13.231949270350407 QDOT_MAX: 2825.663474147248 J_MAX: 447065.25571461295 STAG. PRES: 0.025643897044314443
V_infty: 30.0 km/s, BR: 1.0 G_MAX: 16.90302063923796 QDOT_MAX: 2655.1890245179093 J_MAX: 519761.8784462273 STAG. PRES: 0.03275324229726046
V_infty: 30.0 km/s, BR: 5.0 G_MAX: 16.90302063923796 QDOT_MAX: 3436.725170579263 J_MAX: 519761.8784462273 STAG. PRES: 0.03275324229726046
V_infty: 30.0 km/s, BR: 9.0 G_MAX: 16.90302063923796 QDOT_MAX: 3636.966863267076 J_MAX: 519761.8784462273 STAG. PRES: 0.03275324229726046
V_infty: 30.0 km/s, BR: 13.0 G_MAX: 16.90302063923796 QDOT_MAX: 3756.883255846687 J_MAX: 519761.8784462273 STAG. PRES: 0.03275324229726046
V_infty: 30.0 km/s, BR: 17.0 G_MAX: 16.90302063923796 QDOT_MAX: 3828.3357624603673 J_MAX: 519761.8784462273 STAG. PRES: 0.03275324229726046
V_infty: 30.0 km/s, BR: 21.0 G_MAX: 16.90302063923796 QDOT_MAX: 3875.2080696198727 J_MAX: 519761.8784462273 STAG. PRES: 0.03275324229726046
V_infty: 30.0 km/s, BR: 25.0 G_MAX: 16.90302063923796 QDOT_MAX: 3891.8761523320354 J_MAX: 519761.8784462273 STAG. PRES: 0.03275324229726046
V_infty: 30.0 km/s, BR: 29.0 G_MAX: 16.90302063923796 QDOT_MAX: 3919.5713578059385 J_MAX: 519761.8784462273 STAG. PRES: 0.03275324229726046
V_infty: 30.0 km/s, BR: 33.0 G_MAX: 16.90302063923796 QDOT_MAX: 3939.793314640587 J_MAX: 519761.8784462273 STAG. PRES: 0.03275324229726046
V_infty: 30.0 km/s, BR: 37.0 G_MAX: 16.90302063923796 QDOT_MAX: 3967.916036764199 J_MAX: 519761.8784462273 STAG. PRES: 0.03275324229726046
V_infty: 30.0 km/s, BR: 41.0 G_MAX: 16.90302063923796 QDOT_MAX: 3967.6619222539325 J_MAX: 519761.8784462273 STAG. PRES: 0.03275324229726046

[19]:

x = np.loadtxt('../data/jsr-paper/uranus/'+runID+'vinf_kms_array.txt')
y = np.loadtxt('../data/jsr-paper/uranus/'+runID+'betaRatio_array.txt')

Z1 = np.loadtxt('../data/jsr-paper/uranus/'+runID+'TCW_array.txt')
G1 = np.loadtxt('../data/jsr-paper/uranus/'+runID+'acc_net_g_max_array.txt')
Q1 = np.loadtxt('../data/jsr-paper/uranus/'+runID+'q_stag_total_max_array.txt')
H1 = np.loadtxt('../data/jsr-paper/uranus/'+runID+'heatload_max_array.txt')
S1 = np.loadtxt('../data/jsr-paper/uranus/'+runID+'stag_pres_atm_max_array.txt')


f1 = interpolate.interp2d(x, y, np.transpose(Z1), kind='cubic')
g1 = interpolate.interp2d(x, y, np.transpose(G1), kind='cubic')
q1 = interpolate.interp2d(x, y, np.transpose(Q1), kind='cubic')
h1 = interpolate.interp2d(x, y, np.transpose(H1), kind='cubic')
s1 = interpolate.interp2d(x, y, np.transpose(S1), kind='cubic')


x_new =  np.linspace( 0.0,   30,  110)
y_new =  np.linspace( 0.0,   41 ,110)

z1_new =  np.zeros((len(x_new),len(y_new)))
g1_new =  np.zeros((len(x_new),len(y_new)))
q1_new =  np.zeros((len(x_new),len(y_new)))
h1_new =  np.zeros((len(x_new),len(y_new)))
s1_new =  np.zeros((len(x_new),len(y_new)))

for i in range(0,len(x_new)):
    for j in range(0,len(y_new)):

        z1_new[i,j] = f1(x_new[i],y_new[j])
        g1_new[i,j] = g1(x_new[i],y_new[j])
        q1_new[i,j] = q1(x_new[i],y_new[j])
        h1_new[i,j] = h1(x_new[i],y_new[j])
        s1_new[i,j] = s1(x_new[i],y_new[j])


Z1 = z1_new
G1 = g1_new
Q1 = q1_new
S1 = s1_new
H1 = h1_new/1000.0

X, Y = np.meshgrid(x_new, y_new)

Zlevels = np.array([0.5, 1.0])

Glevels = np.array([8])
Qlevels = np.array([1200, 1600])
Hlevels = np.array([400])
#Slevels = np.array([0.8])


fig = plt.figure()
fig.set_size_inches([3.25,3.25])
rcParams['font.family'] = 'sans-serif'
rcParams['font.sans-serif'] = ['DejaVu Sans']

plt.xlim([0.0,30.0])
plt.ylim([1.0,41.0])


ZCS1 = plt.contour(X, Y, np.transpose(Z1), levels=Zlevels, colors='black')




plt.clabel(ZCS1, inline=1, fontsize=10, colors='black',fmt='%.2f',inline_spacing=1)
ZCS1.collections[0].set_linewidths(1.5)
ZCS1.collections[1].set_linewidths(1.5)
ZCS1.collections[0].set_label(r'$TCW, deg$')


GCS1 = plt.contour(X, Y, np.transpose(G1), levels=Glevels, colors='blue',linestyles='dashed')

Glabels=plt.clabel(GCS1, inline=1, fontsize=10, colors='blue',fmt='%d',inline_spacing=0)
GCS1.collections[0].set_linewidths(1.5)
GCS1.collections[0].set_linewidths(1.5)

GCS1.collections[0].set_label(r'$g$'+r'-load')





QCS1 = plt.contour(X, Y, np.transpose(Q1), levels=Qlevels, colors='red',linestyles='dotted')

plt.clabel(QCS1, inline=1, fontsize=10, colors='red',fmt='%d',inline_spacing=0)
QCS1.collections[0].set_linewidths(1.5)
QCS1.collections[1].set_linewidths(1.5)
QCS1.collections[0].set_label(r'$\dot{q}$'+', '+r'$W/cm^2$')


HCS1 = plt.contour(X, Y, np.transpose(H1), levels=Hlevels, colors='xkcd:brown',linestyles='dashdot')

Hlabels=plt.clabel(HCS1, inline=1, fontsize=10, colors='xkcd:brown',fmt='%d',inline_spacing=0)
HCS1.collections[0].set_linewidths(1.5)

HCS1.collections[0].set_label(r'$Q$'+', '+r'$kJ/cm^2$')


#GCS1.collections[0].set_label(r'$Peak$'+r' '+r'$g-load$')
#plt.grid(True,linestyle='dotted', linewidth=0.3)
params = {'mathtext.default': 'regular' }
plt.rcParams.update(params)
plt.ylabel(r'$\beta_2$'+' / '+r'$ \beta_1 $' ,fontsize=10)
plt.xlabel("Arrival "+r'$V_\infty$'+r', km/s' ,fontsize=10)
plt.xticks(fontsize=10)
plt.yticks(fontsize=10)
ax = plt.gca()
ax.tick_params(direction='in')
ax.yaxis.set_ticks_position('both')
ax.xaxis.set_ticks_position('both')
plt.legend(loc='lower right', fontsize=8)

for l in Hlabels:
    l.set_rotation(-90)
for l in Glabels:
    l.set_rotation(-90)

dat0 = ZCS1.allsegs[1][0]

x1,y1=dat0[:,0],dat0[:,1]
F1 = interpolate.interp1d(x1, y1, kind='linear',fill_value='extrapolate', bounds_error=False)


dat2 = QCS1.allsegs[0][0]
x3,y3= dat2[:,0],dat2[:,1]
F3 = interpolate.interp1d(x3, y3, kind='linear',fill_value='extrapolate', bounds_error=False)

dat0a = ZCS1.allsegs[0][0]
x1a,y1a=dat0a[:,0],dat0a[:,1]
F1a = interpolate.interp1d(x1a, y1a, kind='linear',fill_value='extrapolate', bounds_error=False)


x4 = np.linspace(0,30,301)
y4 = F1(x4)
y4a =F1a(x4)


y6 = F3(x4)

y7 = y6
y8 = np.minimum(y4,y6)



plt.fill_between(x4, y4, y7, where=y4<=y7,color='xkcd:neon green')

plt.fill_between(x4, y4a, y8, where=y4a<=y8,color='xkcd:bright yellow')







plt.savefig('../data/jsr-paper/uranus/uranus-drag-small.png', dpi= 300,bbox_inches='tight')
plt.savefig('../data/jsr-paper/uranus/uranus-drag-small.pdf', dpi=300,bbox_inches='tight')
plt.savefig('../data/jsr-paper/uranus/uranus-drag-small.eps', dpi=300,bbox_inches='tight')

plt.show()

C:\Users\AthulGirija\anaconda3\envs\env1\lib\site-packages\scipy\interpolate\interpolate.py:630: RuntimeWarning: divide by zero encountered in true_divide
  slope = (y_hi - y_lo) / (x_hi - x_lo)[:, None]
The PostScript backend does not support transparency; partially transparent artists will be rendered opaque.
The PostScript backend does not support transparency; partially transparent artists will be rendered opaque.
../_images/jsr-notebooks_09-b-uranus-feasibility-drag_8_1.png 
[25]:

x = np.loadtxt('../data/jsr-paper/uranus/'+runID+'vinf_kms_array.txt')
y = np.loadtxt('../data/jsr-paper/uranus/'+runID+'betaRatio_array.txt')

Z1 = np.loadtxt('../data/jsr-paper/uranus/'+runID+'TCW_array.txt')
G1 = np.loadtxt('../data/jsr-paper/uranus/'+runID+'acc_net_g_max_array.txt')
Q1 = np.loadtxt('../data/jsr-paper/uranus/'+runID+'q_stag_total_max_array.txt')
H1 = np.loadtxt('../data/jsr-paper/uranus/'+runID+'heatload_max_array.txt')
S1 = np.loadtxt('../data/jsr-paper/uranus/'+runID+'stag_pres_atm_max_array.txt')

f1 = interpolate.interp2d(x, y, np.transpose(Z1), kind='cubic')
g1 = interpolate.interp2d(x, y, np.transpose(G1), kind='cubic')
q1 = interpolate.interp2d(x, y, np.transpose(Q1), kind='cubic')
h1 = interpolate.interp2d(x, y, np.transpose(H1), kind='cubic')

x_new =  np.linspace( 0.0,   30,  310)
y_new =  np.linspace( 1.0,   41 , 110)

z1_new =  np.zeros((len(x_new),len(y_new)))
g1_new =  np.zeros((len(x_new),len(y_new)))
q1_new =  np.zeros((len(x_new),len(y_new)))
h1_new =  np.zeros((len(x_new),len(y_new)))
#s1_new =  np.zeros((len(x_new),len(y_new)))

for i in range(0,len(x_new)):
    for j in range(0,len(y_new)):

        z1_new[i,j] = f1(x_new[i],y_new[j])
        g1_new[i,j] = g1(x_new[i],y_new[j])
        q1_new[i,j] = q1(x_new[i],y_new[j])
        h1_new[i,j] = h1(x_new[i],y_new[j])

Z1 = z1_new
G1 = g1_new
Q1 = q1_new
H1 = h1_new/1000.0

X, Y = np.meshgrid(x_new, y_new)

Zlevels = np.array([0.4, 0.6 ,0.7, 0.8,0.9, 1.0, 1.1])

Glevels = np.array([2.0,  4.0, 5.0 , 8.0])
Qlevels = np.array([700, 900.0, 1000.0 , 1200.0, 1500.0 ])
Hlevels = np.array([300, 350, 400, 450])
#Slevels = np.array([0.8])

fig = plt.figure()
fig.set_size_inches([6.5,6.5])
rcParams['font.family'] = 'sans-serif'
rcParams['font.sans-serif'] = ['DejaVu Sans']

ZCS1 = plt.contour(X, Y, np.transpose(Z1), levels=Zlevels, colors='black')




plt.clabel(ZCS1, inline=1, fontsize=10, colors='black',fmt='%.1f',inline_spacing=1)
ZCS1.collections[0].set_linewidths(1.5)
ZCS1.collections[1].set_linewidths(1.5)
ZCS1.collections[2].set_linewidths(1.5)
ZCS1.collections[3].set_linewidths(1.5)
ZCS1.collections[4].set_linewidths(1.5)
ZCS1.collections[5].set_linewidths(1.5)
ZCS1.collections[6].set_linewidths(1.5)


ZCS1.collections[0].set_label(r'$TCW, deg$')


GCS1 = plt.contour(X, Y, np.transpose(G1), levels=Glevels, colors='blue',linestyles='dashed')

Glabels=plt.clabel(GCS1, inline=1, fontsize=10, colors='blue',fmt='%d',inline_spacing=0)
GCS1.collections[0].set_linewidths(1.5)
GCS1.collections[1].set_linewidths(1.5)
GCS1.collections[2].set_linewidths(1.5)

GCS1.collections[0].set_label(r'$g$'+r'-load')


for l in Glabels:
    l.set_rotation(-90)

QCS1 = plt.contour(X, Y, np.transpose(Q1), levels=Qlevels, colors='red',linestyles='dotted')

plt.clabel(QCS1, inline=1, fontsize=10, colors='red',fmt='%d',inline_spacing=0)
QCS1.collections[0].set_linewidths(1.5)
QCS1.collections[1].set_linewidths(1.5)
QCS1.collections[2].set_linewidths(1.5)
QCS1.collections[3].set_linewidths(1.5)
QCS1.collections[4].set_linewidths(1.5)

QCS1.collections[0].set_label(r'$\dot{q}$'+', '+r'$W/cm^2$')


HCS1 = plt.contour(X, Y, np.transpose(H1), levels=Hlevels, colors='xkcd:brown',linestyles='dashdot')

Hlabels=plt.clabel(HCS1, inline=1, fontsize=10, colors='xkcd:brown',fmt='%d',inline_spacing=0)
HCS1.collections[0].set_linewidths(1.75)
HCS1.collections[1].set_linewidths(1.75)
HCS1.collections[2].set_linewidths(1.75)
HCS1.collections[3].set_linewidths(1.75)

HCS1.collections[0].set_label(r'$Q$'+', '+r'$kJ/cm^2$')

for l in Hlabels:
    l.set_rotation(-90)

plt.ylim(1,11)
#plt.grid(True,linestyle='dotted', linewidth=0.3)
params = {'mathtext.default': 'regular' }
plt.rcParams.update(params)
plt.ylabel(r'$\beta_2$'+' / '+r'$ \beta_1 $' ,fontsize=12)
plt.xlabel("Arrival "+r'$V_\infty$'+r', km/s' ,fontsize=12)
plt.xticks(np.array([ 0.0, 5, 10, 15, 20, 25, 30]),fontsize=12)
plt.yticks(np.array([ 1, 10, 20, 30, 40]),fontsize=12)
ax = plt.gca()
ax.tick_params(direction='in')
ax.yaxis.set_ticks_position('both')
ax.xaxis.set_ticks_position('both')
plt.legend(loc='upper right', fontsize=12)

plt.savefig('../data/jsr-paper/uranus/uranus-drag-large.png', dpi= 300,bbox_inches='tight')
plt.savefig('../data/jsr-paper/uranus/uranus-drag-large.pdf', dpi=300,bbox_inches='tight')
plt.savefig('../data/jsr-paper/uranus/uranus-drag-large.eps', dpi=300,bbox_inches='tight')


plt.show()


The PostScript backend does not support transparency; partially transparent artists will be rendered opaque.
The PostScript backend does not support transparency; partially transparent artists will be rendered opaque.
../_images/jsr-notebooks_09-b-uranus-feasibility-drag_9_1.png