streamer_1d/emissionTestImproved.gnuplot at master · MD-CWI/streamer_1d · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
set terminal pngcairo size 1280,1024 enhanced linewidth 2 font "Helvetica, 16"

set xrange [1.95 : 2]
do for [i=1:500] {
        set output sprintf("output/emissionLEATest/emissionManyOp_%06d.png", i)
        set multiplot layout 3,2 columnsfirst title sprintf("t = %.2f ns", (i-1)*0.1)
        # Electron and Ion density
        set ylabel "density m-3"
        set yrange [-1 : 5e5]
        #set log y
        set key bottom left
        plot sprintf("output/emissionLEATest/emissionManyOutputs_fluid_000%03.0f.txt", i)\
        using ($1*1e+3):($3/1e15) with lines lt 1 lw 3 title "LFA ne",\
        "" using ($1*1e+3):($4/1e15) with points ps 0.5 pt 4 title "LFA ni",\
        sprintf("output/emissionLEATest/emissionManyOutputsLEA_fluid_000%03.0f.txt", i) \
        using ($1*1e3):($3/1e15) with lines lt 1 lw 2 lc 6 title "LEA ne",\
        "" using ($1*1e3):($4/1e15) with points ps 0.5 pt 5 title "LEA ni",\
        sprintf("output/emissionLEATest/emission_particle_000%03.0f.txt", i) \
        using ($1*1e3):($3/1e15) with lines lt 1 lw 2 lc 2 title "part ne",\
        "" using ($1*1e3):($4/1e15) with points ps 0.5 pt 3 title "part ni"
        #unset yrange
        #unset log y
        # Electric  field
        set key top center
        set yrange [-1e6:4e7]
        set ylabel "E(V/m)"
        plot sprintf("output/emissionLEATest/emissionManyOutputs_fluid_000%03.0f.txt", i)\
        using ($1*1e+3):2 with lines lt 1 lw 3 title "LFA",\
        sprintf("output/emissionLEATest/emissionManyOutputsLEA_fluid_000%03.0f.txt", i) \
        using ($1*1e3):2 with lines lt 1 lw 2 lc 6 title "LEA",\
        sprintf("output/emissionLEATest/emission_particle_000%03.0f.txt", i) \
        using ($1*1e3):2 with lines lt 1 lw 2 lc 2 title "part"
        unset xlabel
        unset yrange
        # Net charge
        set ylabel "density m-3"
        set yrange [-1 : 1e5]
        #set log y
        set key top center
        plot sprintf("output/emissionLEATest/emissionManyOutputs_fluid_000%03.0f.txt", i)\
        using ($1*1e+3):($4-$3)/1e15 with lines lt 1 lw 3 title "LFA ne-ni",\
        sprintf("output/emissionLEATest/emissionManyOutputsLEA_fluid_000%03.0f.txt", i) \
        using ($1*1e3):($4-$3)/1e15 with lines lt 1 lw 2 lc 6 title "LEA ne-ni",\
        sprintf("output/emissionLEATest/emission_particle_000%03.0f.txt", i) \
        using ($1*1e3):($4-$3)/1e15 with lines lt 1 lw 2 lc 2 title "part ne-ni"
        unset yrange
        # Fluxes
        set ylabel "Fluxes"
        set yrange [-1e9 : 1e9]
        #set log y
        set key bottom left
        plot sprintf("output/emissionLEATest/emissionManyOutputs_fluid_000%03.0f.txt", i)\
        using ($1*1e+3):($9/1e15) with lines lt 1 lw 3  lc 1title "LFA advective",\
        "" using ($1*1e+3):($10/1e15) with lines lt 1 lw 3 lc 2 title "LFA diff",\
        sprintf("output/emissionLEATest/emissionManyOutputsLEA_fluid_000%03.0f.txt", i) \
        using ($1*1e3):($9/1e15) with lines lt 1 lw 2 lc 3 title "LEA advective",\
        "" using ($1*1e3):($10/1e15) with lines lt 1 lw 2 lc 4 title "LEA diff"
        unset yrange
        # Source term
        set key top center
        set yrange [-1:1e6]
        set ylabel "Source term"
        plot sprintf("output/emissionLEATest/emission_fluid_000%03.0f.txt", i)\
        using ($1*1e+3):($6/1e15) with lines lt 1 lw 3 title "LFA",\
        sprintf("output/emissionLEATest/emissionLEA_fluid_000%03.0f.txt", i) \
        using ($1*1e3):($6/1e15) with lines lt 1 lw 2 lc 6 title "LEA",\
        sprintf("output/emissionLEATest/emission_particle_000%03.0f.txt", i) \
        using ($1*1e3):($6/1e15) with lines lt 1 lw 2 lc 2 title "part"
        unset xlabel
        unset yrange
        unset multiplot
}