%% modper.m calculates the polar n(H+) profile, based on Persoon 
%% et al., JGR (1983). The fit used being taken as 
%% n(H+)=490*1e6*(dismat/re)^(-3.85).
%%

function eden=modper(dismat)

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

re=6371;                                %Earth radi
eden=490*1e6*(dismat/re).^(-3.85);      %Electron density