%% modvance.m calculates the polar n(H+) profile, based on Vance
%% Henize IMAGE RPI derived polar profile. The fit used being taken as 
%% n(H+)=3.5e4*1e6*(dismat/re)^(-6.4).
%%

function eden=modvance(dismat)

% dismat = radial distance of point of interest (km)
% eden = electron density (e-/m^3)

re=6371;                                %Earth radi
eden=3.5e10*((dismat/re).^(-6.4));      %electron density
