#!/usr/bin/gnuplot # vitr.plt # Ztrata hmoty pusobenim korpuskularniho slunecniho vetru. # Miroslav Broz (miroslav.broz@email.cz), May 25th 2011 m_u = 1.66e-27 # kg AU = 1.5e11 # m M_S = 2.e30 # kg yr = 365.25*86400. # s n = 1.e7 # castic/m^3 v = 400.e3 # m/s r = 1.*AU # m rho = m_u * n dM_dt = 4*pi*r**2 * v * rho print "rho = ", rho, " kg/m^3" print "dM/dt = ", dM_dt, " kg/s = ", dM_dt/M_S*yr, " M_S/yr" dt = 4.5e9*yr dM = dM_dt * dt print "dt = ", dt, " s = ", dt/(1.e9*yr), " 10^9 yr" print "dM = ", dM, " kg = ", dM/M_S, " M_S"