plasmonic_fdtd/plot.p at master · danieljg/plasmonic_fdtd · 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
set term png enhanced size 800,800

red = 1.0/20.0
yellow = 1.0/6.0
green = 1.0/3.0
blue = 2.0/3.0

#set size ratio -1
set size ratio 0
#set xrange [0:1.620]
#set yrange [0:0.645]
set xlabel 'x ({/Symbol m}m)'
set ylabel 'z ({/Symbol m}m)'
set pm3d
unset surface
set view map


 #print ' plotting index '.idx.' of '.snapshots
 set palette model HSV defined ( 0 blue 0.7 0.3, 1 blue 0.7 1, 2 green 0.5 0.6, 3 green 0.7 0.3 )
outfile = sprintf('tmp/H%04.0f.png',idx)
 set output outfile
 set title 'Magnetic field'
 set cbrange [-1.:1.]
 unset logscale zcb
 splot 'tmp/H.dat' u ($1*1e4):($2*1e4):($3) noti
 set output

 set palette model HSV defined ( 0 blue 0.7 0.3, 1.5 green 0.5 0.6, 2 yellow 0.5 0.7, 2.5 red 0.7 1, 3.5 red 0.8 0.9, 5 red 0.7 0.7 )
 outfile = sprintf('tmp/E%04.0f.png',idx)
 set title 'Electric field intensity'
 set output outfile
 set cbrange [0:1]
 splot 'tmp/E.dat' u ($1*1e4):($2*1e4):($3) noti
 set output

 set cbrange [0:0.001]
 #set cbrange [0:0.018]
 #set autoscale cb
 set size ratio 0
 outfile = sprintf('tmp/Psi%04.0f.png',idx)
 set title 'Material response'
 set output outfile
 set xlabel 'x (nm)'
 splot 'tmp/Psi.dat' u ($1*1e7):($2*1e4):($3) noti
 set output

#system('rm animation/*.eps')
#system('rm animation/*.pdf')
#system('mv *.png animation/')