#!/usr/bin/gnuplot # dipol.plt # Intenzita magnetickeho pole dipolu. # Miroslav Broz (miroslav.broz@email.cz), May 31st 2011 deg = pi/180. R_Z = 6378.e3 # m M_B = 7.9e25 * 1e-4 * 1e-6 # G cm^3 -> T m^3 theta = 0.*deg r = R_Z B_r = -2.*M_B/r**3 * sin(theta) B_theta = -M_B/r**3 * cos(theta) print "M_B = ", M_B, " T m^3" print "B_r = ", B_r, " T" print "B_theta = ", B_theta, " T"