VERB_code_2.3
Figure_4.m
1 clear all
2 
3 path1 = '../Execute/3D_Without_0.1hr/';
4 
5 [time, Kp, Bf, Lpp, Bw2] = load_plt([path1, 'Output/Out1d.dat']);
6 
7 
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]);%
11 
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')
15 
16 subplot1(1)
17 plot(time.arr, Kp.arr, 'k');
18 ylim([0 7])
19 xlabel('Time, days');
20 ylabel('Kp');
21 text(min(xlim), max(ylim), [' (A) '], 'VerticalAlignment', 'top', 'HorizontalAlignment', 'left', 'FontWeight', 'bold', 'Color', 'black');
22 
23 subplot1(2)
24 plot(time.arr, Bf.arr, 'k');
25 ylim([0 1.2])
26 xlabel('Time, days');
27 ylabel('Boundary flux');
28 text(min(xlim), max(ylim), [' (B) '], 'VerticalAlignment', 'top', 'HorizontalAlignment', 'left', 'FontWeight', 'bold', 'Color', 'black');
29 
30 subplot1(3)
31 plot(time.arr, Lpp.arr, 'k');
32 ylim([0 7])
33 xlabel('Time, days');
34 ylabel('Plasmapause location');
35 text(min(xlim), max(ylim), [' (C) '], 'VerticalAlignment', 'top', 'HorizontalAlignment', 'left', 'FontWeight', 'bold', 'Color', 'black');
36 
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...
Definition: erf.cpp:103