In [ ]:

    




    

In [ ]:

    

#%run '/Users/admin/s3_handling_poc.py'
import sys
sys.path.append('../../python')
import pySimE as sim
import pySimE.space as space
import numpy as np
import matplotlib.pyplot as plt

'''
Ts = np.linspace( 300,3000 )
Ps = sim.thermalRadiativePower( temperature = Ts )
plt.plot( Ps )

plt.show()
'''

#print dir(space)


plt.figure(figsize=(16,16))

Rs = np.linspace( 100.0,10000.0,1000 )

plt.subplot(2,2,1)
plt.title( 'Escape velocity from asteroides (density=2.5kg/l)'); 
vs = space.escapeVelocity( radius=Rs, density=2500.0 )
plt.plot( Rs, vs )
plt.xlabel("Radius[m]"); plt.ylabel("Escape velocity [m/s]"); 
plt.grid()

plt.subplot(2,2,3)
plt.title( 'Surface gravity from asteroides (density=2.5kg/l)'); 
Gs = space.gravitationalAcceleration( radius=Rs, density=2500.0 )
plt.plot( Rs, Gs )
plt.xlabel("Radius[m]"); plt.ylabel("Surface gravity [m/s^2]"); 
plt.grid()

Rs = np.linspace( 10.0,1.0e+9, 1000 )
plt.subplot(2,2,2)
plt.title("Escape velocity from bodies log-scale (density=2.5kg/l)"); 
vs = space.escapeVelocity( radius=Rs, density=2500.0 )
plt.plot( Rs, vs )
plt.xlabel("Radius[m]"); plt.ylabel("Escape velocity [m/s]"); 
plt.xscale('log'); plt.yscale('log'); plt.grid()

plt.subplot(2,2,4)
plt.title("Surface gravity for bodies log-scale (density=2.5kg/l)"); 
Gs = space.gravitationalAcceleration( radius=Rs, density=2500.0 )
plt.plot( Rs, Gs )
plt.xlabel("Radius[m]"); plt.ylabel("Surface gravity [m/s^2]"); 
plt.xscale('log'); plt.yscale('log'); plt.grid()

plt.show()