My Project
(2004)_README Description

The folder contains all of the functions needed to reproduce the plots observed in the 2004 paper (same name as this directory).

There are only two main functions here:

All other functions are helper functions used in these codes. The steady_state.m file plots the time-independent PSD; the 1D Radial Diffusion folder holds the code needed to reproduce pcolor plots of flux versus L versus time (please see README in that folder).

Reproduction Tests

Please run the following code segments to ensure that the output is in line with expectations.

1 L = linspace(2,8);
2 Kp = 2;
3 plot(L, 1./RadialDiffCoeff(Kp,L))
4 set(gca,'ylim',[10e-1 10e4],'xlim',[2 8],'yscale','log');

The code above will reproduce our diffusion coefficients vs. L-star plots.

1 L = linspace(2,8);
2 Kp = 2;
3 DLL = RadialDiffCoeff(Kp, L);
4 tau = 10;
5 plot(L, steady_state(L, DLL, tau))
6 set(gca,'ylim',[10e-1 10e4],'xlim',[2 8],'yscale','log');

The code above will reproduce our phase space density vs. L-star plots.

1 L = linspace(2,8);
2 Kp = 2;
3 tau = 10;
4 plot(L, ElectronFlux(L, Kp, tau))
5 set(gca,'ylim',[10e-1 10e4],'xlim',[2 8],'yscale','log');

The code above will reproduce our electron flux vs. L-star plots.