% Init script for foreshock.m to adjust independent physical variables

global Te Ti Matot nhnc ecc
% global theta
global kappa
global vsw
global Ne
global Rp Rrange xrange

Te = 1.5e5;     % electron temperature [K]
Ti = 5.0e4;     % ion temperature [K]
Ne = 8e6;     % electron number density [m^-3]
vsw = 4.0e5;    % Solar wind velocity (upstream wrt shock) [m/s]
theta_deg = 300;    % angle between -vsw & tangential B-field (set between 330 = 30 and 270 = 90)
theta = theta_deg;
valfven = 4.0e5/9;
Matot = vsw/valfven;         % Alfven Mach number
nhnc = 0.05;    
kappa = 3;
Rp = 6.3780e6;  % Planet's radius
Rrange = 50*Rp;   % spatial extent of source region
xrange = 1.25*Rp;

ecc = 1.00;     % eccentricity of conic bow shock  

vr = 300001;              % beams array size

RR = 250;       % dimensions wrt B-field
xR = 250;
