3 path1 =
'../Execute/3D_Without_0.1hr/';
5 [time, Kp, Bf, Lpp, Bw2] = load_plt([path1,
'Output/Out1d.dat']);
8 f1 = figure(
'Units',
'inches',
'Position', [0.5 0.5 8.5 2], ...
9 'DefaultAxesFontName',
'times new roman',
'DefaultAxesFontWeight',
'bold',
'DefaultAxesFontSize', 7,
'DefaultTextFontWeight',
'bold',
'DefaultTextFontSize', 7,...
10 'PaperOrientation',
'portrait',
'PaperPositionMode',
'auto',
'DefaultLineLineWidth', 2);%,
'PaperPosition', [1 1 3.8 9]);%
12 xmin = 0.05; ymin = 0.2;
13 xmax = 1; ymax = 0.95;
14 subplot1(1, 3,
'Min', [xmin ymin],
'Max', [xmax ymax],
'Gap', [0.06 0.01],
'YTickL',
'All')
17 plot(time.arr, Kp.arr, 'k');
21 text(min(xlim),
max(ylim), [' (A) '], 'VerticalAlignment', 'top', 'HorizontalAlignment', 'left', 'FontWeight', 'bold', 'Color', 'black');
24 plot(time.arr, Bf.arr, 'k');
27 ylabel('Boundary flux');
28 text(min(xlim),
max(ylim), [' (B) '], 'VerticalAlignment', 'top', 'HorizontalAlignment', 'left', 'FontWeight', 'bold', 'Color', 'black');
31 plot(time.arr, Lpp.arr, 'k');
34 ylabel('Plasmapause location');
35 text(min(xlim),
max(ylim), [' (C) '], 'VerticalAlignment', 'top', 'HorizontalAlignment', 'left', 'FontWeight', 'bold', 'Color', 'black');
37 save_picture(
gcf, [path1, 'Figure_4_bw'], 'jpg', 'pdf', 'eps');
double max(double v1, double v2)
Return maximum.
void gcf(double *gammcf, double a, double x, double *gln)
Returns the incomplete gamma function Q(a, x) evaluated by its continued fraction representation as g...