#!/usr/bin/gnuplot # destnik_io.plt # Vulkanicky destnik na Io, pocitany v HOMOGENNIM gravitacnim poli. # Miroslav Broz (miroslav.broz@email.cz), Feb 11th 2009 km = 1.e3 # m G = 6.67e-11 # N m^2/kg^2; gravitacni konstanta k = 1.38e-23 # J/K; Boltzmannova konstanta m_u = 1.67e-27 # kg; atomova hmotnostni konstanta M = 8.9319e22 # kg; hmotnost Io R = 1821.3*km # m; polomer g_Z = 9.81 # m/s^2; tihove zrychleni na Zemi v_0 = 800. # m/s; rychlost plynu tryskajiciho z gejziru T = 1600. # K; teplota plynu mu = 32.+2*16. # ; relativni molekulova hmotnost SO2 g_Io = sqrt(G*M/R**2) # tihove zrychleni na Io v_esc = sqrt(2*G*M/R) # unikova rychlost z Io v_th = sqrt(3.*k*T/(mu*m_u)) # termalni rychlost plynu print "v_0 = ", v_0, " m/s" print "v_th = ", v_th, " m/s" print "v_esc = ", v_esc, " m/s" print "g_Io = ", g_Io, " m/s^2 = ", g_Io/g_Z, " g_Z" ######################################################################## set xl "x [km]" offset 0,+0.25 set yl "y [km]" offset +1 set tr [0:1] tmp=500 set xr [-tmp:tmp] set yr [-50:300] set xtics 100 set ytics 100 set parametric set size ratio -1 set ang deg set nokey set samples 1000 lw=3 set style line 1 lt 9 lw 3 set style line 2 lt 1 lw 3 set style line 3 lt 7 lw 6 set style line 4 lt 7 lw 0.5 a = g_Io h = 3600. x(phi,t) = v_0*cos(phi)*(t*h)/km y_tmp(phi,t) = (v_0*sin(phi)*(t*h) - 1./2.*a*(t*h)**2)/km y(phi,t) = y_tmp(phi,t) > 0 ? y_tmp(phi,t) : NaN p \ x( 0,t), y( 0,t) ls 1,\ x( 10,t), y( 10,t) ls 1,\ x( 20,t), y( 20,t) ls 1,\ x( 30,t), y( 30,t) ls 1,\ x( 40,t), y( 40,t) ls 1,\ x( 50,t), y( 50,t) ls 1,\ x( 60,t), y( 60,t) ls 1,\ x( 70,t), y( 70,t) ls 1,\ x( 80,t), y( 80,t) ls 1,\ x(100,t), y(100,t) ls 1,\ x(110,t), y(110,t) ls 1,\ x(120,t), y(120,t) ls 1,\ x(130,t), y(130,t) ls 1,\ x(140,t), y(140,t) ls 1,\ x(150,t), y(150,t) ls 1,\ x(160,t), y(160,t) ls 1,\ x(170,t), y(170,t) ls 1,\ x(180,t), y(180,t) ls 1,\ x( 45,t), y( 45,t) ls 2,\ x( 90,t), y( 90,t) ls 2,\ R/km*cos(t*360.), R/km*(sin(t*360.)-1) ls 3,\ tmp*2.*(t-0.5), 0. ls 4 #pa -1 set term post eps enh color solid "Helvetica" 18 set out "destnik_io.eps" set size 0.85,0.5 rep