#######################       LOOPER
reset

set output outputfile(i)

set multiplot
set pm3d
unset surf
set style data lines
set palette defined ( 0 'blue', 250 'cyan', 500 'green',750 'orange', 1000 'red')
set view map
set key noauto

set xtics nomirror
set ytics nomirror
unset clabel

unset ztics
set ticslevel 0

t_marge_tmp=t_marge
####time_mark_tmp=`sed -e "$i!d" ../SD2_ne9/Toutput.dat | awk '{print $1}'`
######set label @title center at screen l_marge,t_marge font @police
#str=eval(TIME(i))

time=system(VAR_FILE(i,1))
dt=system(VAR_FILE(i,2))
ua=system(VAR_FILE(i,4))

#; tmp2=system(TIME(tmp)) 
 ###; time_label=sprintf("t=%d s",tmp)
#tmp=system(VAR_FILE(i,2)); dt_label=sprintf("t=%d s",tmp)
#tmp=system(VAR_FILE(i,4)); ua_label=sprintf("t=%d V",tmp)

#set label system(TIME(i,1)) at screen 0.70, screen 0.96 font "arial,20"
#set label str  at screen 0.70, screen 0.96 font "arial,20"

##########
set label "time=".time." s"  at screen l_marge, t_marge+0.05
set label "dt=".dt." s" at screen l_marge + 0.15, t_marge+0.05 
set label "ua=".ua." V" at screen l_marge + 0.3, t_marge+0.05 
#########

set lmargin screen lm2D
set rmargin screen rm2D
set tmargin screen tm2D
set bmargin screen bm2D

set colorbox horizontal user origin lm2D,tm2D+decv_cb size L_cb, H_cb 

set yrange [0.6:1.4]
set xrange [-0.2:0.2]
set cblabel "Electric field V/cm" offset screen 0,3*H_cb @police
set xlabel " r (cm)" offset screen 0, 0.01 @police
set xtics offset screen 0, 0.007 0.1 @police 
set ylabel " x (cm)" offset screen (rm2D-lm2D)+0.11,0 @police 
set format y "%.2f"
set ytics offset  screen (rm2D-lm2D)+0.05,0 @police 
set ytics mirror scale 2 @police
set cbrange [0:60]
set logscale z

set arrow from -0.25,0.8,1 to -0.2,0.8,1 lw 4 lt -1
set arrow from -0.25,0.9,1 to -0.2,0.9,1 lw 4 lt 1
set arrow from -0.25,1.0,1 to -0.2,1.0,1 lw 4 lt 2
set arrow from -0.25,1.1,1 to -0.2,1.1,1 lw 4 lt 3
set arrow from -0.25,1.2,1 to -0.2,1.2,1 lw 4 lt 4

splot EF(i) binary u 2:1:(0.001*$3)  lt 1,\
          "" binary u (-$2):1:(0.001*$3) lt 1 

set style data lines

set ylabel offset screen 0,0
set ytics offset screen 0,0 @police
set xlabel offset screen 0,0
set xtics offset screen 0,0 @police

set format y "%.0f"

unset label
unset arrow
set mxtics 10
###########     FIRST IMAGE FROM TOP

set lmargin screen l_marge
set rmargin screen l_marge+L_image
set tmargin screen t_marge_tmp
set bmargin screen t_marge_tmp-H_image
t_marge_tmp=t_marge_tmp-H_image-int_v

unset colorbox
set autoscale y
set autoscale x
set grid
set ylabel "Electric Field (kV/cm)" @police
set xlabel ""
set format x ""
set xrange[-0.2:0.2]
set yrange[:60]
set key inside left top vertical Right noreverse enhanced autotitles box spacing 1.25

plot   EF(i) binary u 2:(0.001*$3) every ::430::430 lt -1 lw 4 title "0,8cm" axis x1y1,\
       EF(i) binary u 2:(0.001*$3) every ::630::630 lt 1 lw 4 title "0,9cm" axis x1y1,\
       EF(i) binary u 2:(0.001*$3) every ::830::830 lt 2 lw 4 title "1,0cm" axis x1y1,\
       EF(i) binary u 2:(0.001*$3) every ::1030::1030 lt 3 lw 4 title "1,1cm" axis x1y1,\
       EF(i) binary u 2:(0.001*$3) every ::1230::1230 lt 4 lw 4 title "1,2cm" axis x1y1,\
       EF(i) binary u (-$2):(0.001*$3) every ::430::430 lt -1 lw 4 title "" axis x1y1,\
       EF(i) binary u (-$2):(0.001*$3) every ::630::630 lt 1 lw 4 title "" axis x1y1,\
       EF(i) binary u (-$2):(0.001*$3) every ::830::830 lt 2 lw 4 title "" axis x1y1,\
       EF(i) binary u (-$2):(0.001*$3) every ::1030::1030 lt 3 lw 4 title "" axis x1y1,\
       EF(i) binary u (-$2):(0.001*$3) every ::1230::1230 lt 4 lw 4 title "" axis x1y1

       
unset logscale y


######################### THIRD IMAGE FROM TOP        

set lmargin screen l_marge
set rmargin screen l_marge+L_image
set tmargin screen t_marge_tmp
set bmargin screen t_marge_tmp-H_image
t_marge_tmp=t_marge_tmp-H_image-int_v

set ylabel "Electron density (cm^{-3})" @police 
unset yrange
set yrange [1e8:1e15]
set logscale y
unset key
set xlabel "r (cm)" @police
set format y "10^{%L}"
set format x "%.1f" 
plot   Ne(i) binary u 2:3 every ::430::430 lt -1 lw 4 title "" axis x1y1,\
       Ne(i) binary u 2:3 every ::630::630 lt 1 lw 4 title "" axis x1y1,\
       Ne(i) binary u 2:3 every ::830::830 lt 2 lw 4 title "" axis x1y1,\
       Ne(i) binary u 2:3 every ::1030::1030 lt 3 lw 4 title "" axis x1y1,\
       Ne(i) binary u 2:3 every ::1230::1230 lt 4 lw 4 title "" axis x1y1,\
       Ne(i) binary u (-$2):3 every ::430::430 lt -1 lw 4 title "" axis x1y1,\
       Ne(i) binary u (-$2):3 every ::630::630 lt 1 lw 4 title "" axis x1y1,\
       Ne(i) binary u (-$2):3 every ::830::830 lt 2 lw 4 title "" axis x1y1,\
       Ne(i) binary u (-$2):3 every ::1030::1030 lt 3 lw 4 title "" axis x1y1,\
       Ne(i) binary u (-$2):3 every ::1230::1230 lt 4 lw 4 title "" axis x1y1

unset multiplot
######
show output
i=i+10
if(i < 262) reread

########  END LOOPER


