Ejemplo n.º 1
0
#     for inc in incs:
#         rpar, rperp = find_Rpar_Rperp(beta, inc)
#         rpars.append(R0(beta)*rpar)
#         rperps.append(R0(beta)*rperp)

#     plt.plot(rpars, rperps, "--", label="beta = {:.4f}".format(beta))

# Second, plot solutions as function of beta for constant inc
incs_deg = np.linspace(0.0, 90.0, 18, endpoint=False)
betas = np.logspace(-6.0, -0.2, 300)
for inc_deg in incs_deg:
    inc = np.radians(inc_deg)
    print "inc = ", inc_deg
    rpars, rperps = list(), list()
    for beta in betas:
        thlim = theta_lim(beta)
        inc_lim = (1.0 - tol)*(thlim - 0.5*np.pi)
        if inc > inc_lim:
            break
        rpar, rperp = find_Rpar_Rperp(beta, inc)
        rpars.append(R0(beta)*rpar)
        rperps.append(R0(beta)*rperp)

    plt.plot(rpars, rperps, "-", label="inc = {:.0f}".format(inc_deg))
    

plt.xlim(0.0, 0.35)
plt.ylim(0.0, 1.0)
plt.xlabel("R_par")
plt.ylabel("R_perp")
plt.legend(loc="upper left", ncol=2, prop=dict(size="xx-small"))
Ejemplo n.º 2
0
import sys
sys.path.append("projected")
import bowfuncs as bow
import numpy as np
import matplotlib.pyplot as plt

#print bowfuncs.xt(1.5, 0.1, 0.0)

#steps:
# 0 Array for theta, beta and i
# 1 Create bowshock for different i and beta
# 2 select R0 and R90 for each i
# 3 Make the plot

N = 100
beta = np.array([0.15, 0.05, 0.02, 0.005,0.2,0.3])
i = np.linspace(0,30*np.pi/180.,5)

for b in beta:
    theta = np.linspace(0.01,bow.theta_lim(b),100)
    for j in i:
        xth,yth = np.array([bow.xt(t,b,j) for t in theta ]), np.array([bow.yt(t,b,j) for t in theta ])
        plt.plot(xth,yth,label='i={}'.format(j))
    plt.axis([-1.0,2.0,-0.05,2.0])
    plt.legend()
    plt.xlabel('z')
    plt.ylabel('r')
    plt.title('Bowshock shapes for isotropic wind and beta = {}'.format(b))
    plt.savefig('isotropic-will-wind-plane-sky-beta-{}.png'.format(b))
    plt.clf()