#!/usr/bin/gnuplot # jeans_unik.plt # Jeansuv unik castic z atmosfery. # Miroslav Broz (miroslav.broz@email.cz), Mar 1st 2010 k = 1.38e-23 # J/K; Boltzmannova konstanta m_u = 1.66e-27 # kg; atomova hmotnostni jednotka G = 6.67e-11 # gravitacni konstanta M_Z = 5.e24 # kg; hmotnost Zeme R_Z = 6378.e3 # m; polomer Zeme z = 100.e3 # m n_exo = 1.e10 # m^-3; koncentrace castic na exobazi T = 1000. # K; teplota exosfery v_esc = sqrt(2.*G*M_Z/(R_Z+z)) v_T(m) = sqrt(2.*k*T/m) lambda_esc(m) = (v_esc/v_T(m))**2 Phi_J(m) = 1./sqrt(pi) * n_exo * v_T(m) * (1+lambda_esc(m)) * exp(-lambda_esc(m)) print "v_esc = ", v_esc, " m/s = ", v_esc/1.e3, " km/s" ######################################################################## set xl "m / m_u" set yl "Phi_J / castic/m^2" set xr [1:20] set logscale y p Phi_J(x*m_u) pa -1 set term png small set out "jeans_unik.png" rep