Platon Technologies
not logged in Login Registration
EnglishSlovak
open source software development celebrating 10 years of open source development! Thursday, March 28, 2024

File: [Platon] / doc / diplomova-praca-rajo / images / vdrift_T-300K-x2-model.gpi (download)

Revision 1.2, Sun Jan 16 15:01:31 2005 UTC (19 years, 2 months ago) by rajo


Changes since 1.1: +2 -2 lines

Fix: gnuplot 4.0 has different syntax (incompatible with gnuplot 3.7.x)
Change lines from:
    set linestyle 1 linetype 7 linewidth 1
to:
    set line style 1 linetype 7 linewidth 1

#!/usr/bin/gnuplot

reset

#set terminal postscript enhanced
#set output "vdrift_T-300K-x2-model.ps"
set terminal postscript eps enhanced 16
set output "vdrift_T-300K-x2-model.eps"

set line style 1 linetype 7 linewidth 1
set line style 2 linetype 1 linewidth 3

set multiplot
set size 1.0,1.0
set lmargin 8
set rmargin 1


set mytics 2
set mxtics 2

set format y "%2.2f"

set grid

# ------------------------------------------------------------------------
#     +-----+  +-----+
#     |  X  |  |     |
#     +-----+  +-----+
#     +-----+  +-----+
#     |     |  |     |
#     +-----+  +-----+
set origin 0.0,0.5
set size 0.5,0.5
set key bottom
set ytics 0.2

set xlabel "{/Symbol t}\n a)"
set ylabel "u_d({/Symbol t})"

set label 1     "{/Symbol a} = 0.1" at 5,0.10
set label 2     "{/Symbol a} = 0.2" at 5,0.35
set label 5     "{/Symbol a} = 0.5" at 5,0.50
set label 10    "{/Symbol a} = 1.0"  at 5,0.60
set label 10000 "{/Symbol a} = 10^3" at 5,0.85

plot [:8][0:0.95] \
    "x2-model_N-5e5_T-300_alpha-1000_speed-0.dat"    using 2:6 notitle with lines 1, \
    "x2-model_N-5e5_T-300_alpha-1_speed-0.dat"        using 2:6 notitle with lines 1, \
    "x2-model_N-5e5_T-300_alpha-0.5_speed-0.dat"    using 2:6 notitle with lines 1, \
    "x2-model_N-5e5_T-300_alpha-0.2_speed-0.dat"    using 2:6 notitle with lines 1, \
    "x2-model_N-5e5_T-300_alpha-0.1_speed-0.dat"    using 2:6 notitle with lines 1



# ------------------------------------------------------------------------
#     +-----+  +-----+
#     |     |  |     |
#     +-----+  +-----+
#     +-----+  +-----+
#     |  X  |  |     |
#     +-----+  +-----+
set origin 0.0,0.0
set size 0.5,0.5
set key bottom
set ytics 0.01
set xlabel "{/Symbol t}\n c)"
set ylabel "S({/Symbol t})"

set nolabel
set label 1 "-{/Symbol S}_i {/Symbol f}_i ln{/Symbol f}_i({/Symbol t})" at 2.2, 1.207

plot [][1.185:] "x2-model_N-5e5_T-300_alpha-1_beta-0_speed-0.ENTROPY.dat" using 2:3 \
        notitle with l
    


# ------------------------------------------------------------------------
#     +-----+  +-----+
#     |     |  |  X  |
#     +-----+  +-----+
#     +-----+  +-----+
#     |     |  |     |
#     +-----+  +-----+

set origin 0.5,0.5
set size 0.5,0.5
set lmargin 9
set rmargin 0
set key top
#set ytics 0.005
set xlabel "{/Symbol a}\n b)"
set ylabel ""
set nolabel
set label 1 "v_{drift} [m.s^{-1}]" at graph -0.2,graph 0.5 center rotate

set logscale xy
set key bottom
set xtics ("10^{-4}" 1e-4, "10^{-3}" 1e-3, "10^{-2}" 1e-2, "10^{-1}" 1e-1, "1" 1, "10^1" 10, "10^2" 1e2, "10^3" 1e3)
set ytics (    "10^{-3}" 1e-3,    "10^{-2}" 1e-2,    "10^{-1}" 1e-1,    \
            "10^{0 }" 1,    "10^{1 }" 10,    "10^{2 }" 1e2,    \
            "10^{3 }" 1e3,    "10^{4 }" 1e4,    "10^{5 }" 1e5    )

k = 891.053540967816
#a               = 1.194            +/- 0.007405
a = 1.194

plot [1e-3:1e3] [1e-1:1e5] \
    "zavislost-vdrift-E_T-300-x2-model.dat" title "MC" with points 1, \
    k * sqrt(x) title "~E^{0.50}" with lines 2, \
    k * x        title "~E^{1.00}" with lines 2, \
    k * x**a title "~E^{1.19}" with lines 1


set nologscale xy
set format y "%2.2f"

# ------------------------------------------------------------------------
#     +-----+  +-----+
#     |     |  |     |
#     +-----+  +-----+
#     +-----+  +-----+
#     |     |  |  X  |
#     +-----+  +-----+
# 
set origin 0.5,0.0
set size 0.5,0.5
set lmargin 9
set rmargin 0
set key top
#set xlabel "N\n d)"
#set ylabel "\n\n{/Symbol e}"
set xlabel "{/Symbol t}\n d)"
set ylabel "u_d({/Symbol t})"

set nolabel

set nomultiplot

reset

set terminal postscript eps enhanced 16
set output "x2-model-distrib-3D_projection.eps"

set view 50,50;
#set view 90,0;
set ticslevel 0

set ylabel "{/Symbol t}"
set xlabel "u"
set zlabel "{/Symbol f}(u, {/Symbol t})"
set format z "%2.2f"

v_T = 1116.77
Nx  = 500
f_max_v          = 5.000000
fi_norm   =   2* Nx / (f_max_v * 4 * pi / (sqrt(pi))**3 )

set border 1+2+4+8+16+32+64+128+256+512+1024+2048

splot [-4:4][:4][] \
    "x2-model-distrib-N-5e5_T-300_alpha-1_beta-0_speed-0.dat" using ($6 / v_T):2:(fi_norm * $7) notitle with l 1, \
    "" using (-$6/ v_T):2:(fi_norm * $8) notitle with lines 1, 0 notitle 



Platon Group <platon@platon.org> http://platon.org/
Copyright © 2002-2006 Platon Group
Site powered by Metafox CMS
Go to Top