bmarechal / gnuplot_scripts

Download as
Browse Code »

Commit 7047b0001eaada1e95b448d30b8b6aa0ef020200

Authored by bmarechal
0 parents
Exists in master

initial commit

Showing 10 changed files with 516 additions and 0 deletions Side-by-side Diff

OFC-FRR/temp_mon01.gnu
Diff comments   View file @ 7047b00
... ... @@ -0,0 +1,88 @@
  1 +set grid
  2 +# init graph
  3 +#set log y
  4 +#set format y "%.1e"
  5 +set timefmt "%s"
  6 +set format x "%y/%m/%d - %H:%M:%S"
  7 +set xdata time
  8 +set xtics rotate by 90 right
  9 +
  10 +
  11 +
  12 +
  13 +file1 = '/home/sicav/server/data/2017/2017-03/20170308-150806-AG34972A.dat'
  14 +file2 = '/home/sicav/server/data/2017/2017-03/20170309-150807-AG34972A.dat'
  15 +file3 = '/home/sicav/server/data/2017/2017-03/20170310-150807-AG34972A.dat'
  16 +file4 = '/home/sicav/server/data/2017/2017-03/20170311-150807-AG34972A.dat'
  17 +file5 = '/home/sicav/server/data/2017/2017-03/20170312-150808-AG34972A.dat'
  18 +file6 = '/home/sicav/server/data/2017/2017-03/20170313-150808-AG34972A.dat'
  19 +file7 = '/home/sicav/server/data/2017/2017-03/20170314-150809-AG34972A.dat'
  20 +
  21 +N_dec =1
  22 +Tamb = 273.15
  23 +
  24 +
  25 +# formula for NTC 10k: T = 1 / ( A + B * ln(R) + C * ln(R)^3)
  26 +A = 2.108508173e-3
  27 +B = 0.7979204727e-4
  28 +C = 6.535076315e-7
  29 +T(R) = 1 / (A + B*log(R) + C*log(R)**(3) ) - Tamb
  30 +
  31 +#formula for PT 100 Ohms in °C:
  32 +kpt = 0.003916
  33 +Tpt(R) = ((R-100)/(100*kpt))-0.9264
  34 +
  35 +#plot file1 u 1:(Tpt($3)) every N_dec title 'Labo (ch 101)' lc 'red' w l, file2 u 1:(Tpt($3)) every N_dec title '' lc 'red' w l, file3 u 1:(Tpt($3)) every N_dec title '' lc 'red' w l
  36 +
  37 +#replot file1 u 1:(Tpt($4)) every N_dec title 'Temp side - monitor (ch 104)' lc 'blue' w l, file2 u 1:(Tpt($4)) every N_dec title '' lc 'blue' w l, file3 u 1:(Tpt($4)) every N_dec title '' lc 'blue' w l
  38 +
  39 +#replot file1 u 1:(T($5)) every N_dec title 'Control side - monitor (ch 105)' lc 'purple' w l, file2 u 1:(T($5)) every N_dec title '' lc 'purple' w l, file3 u 1:(T($5)) every N_dec title '' lc 'purple' w l
  40 +
  41 +
  42 +#plot file1 u 1:(Tpt($3)) every N_dec title 'Labo (ch 101)' lc 'red' w l, file2 u 1:(Tpt($3)) every N_dec title '' lc 'red' w l, file3 u 1:(Tpt($3)) every N_dec title '' lc 'red' w l, file4 u 1:(Tpt($3)) every 10 title '' lc 'red' w l
  43 +
  44 +#replot file1 u 1:($4+22.21) every N_dec title 'PD DC Out+22.21 - monitor (ch 103)' lc 'black' w l, file2 u 1:($4+22.21) every N_dec title '' lc 'black' w l, file3 u 1:($4+22.21) every N_dec title '' lc 'black' w l, file4 u 1:($4+22.21) every 10 title '' lc 'black' w l
  45 +
  46 +#replot file1 u 1:(Tpt($5)) every N_dec title 'Temp side - monitor (ch 104)' lc 'blue' w l, file2 u 1:(Tpt($5)) every N_dec title '' lc 'blue' w l, file3 u 1:(Tpt($5)) every N_dec title '' lc 'blue' w l, file4 u 1:(Tpt($5)) every 10 title '' lc 'blue' w l
  47 +
  48 +#replot file1 u 1:(T($6)) every N_dec title 'Control side - monitor (ch 105)' lc 'purple' w l, file2 u 1:(T($6)) every N_dec title '' lc 'purple' w l, file3 u 1:(T($6)) every N_dec title '' lc 'purple' w l, file4 u 1:(T($6)) every 10 title '' lc 'purple' w l
  49 +
  50 +
  51 +#plot file16 u 1:(Tpt($3)) every 120 title 'Labo (ch 101)' lc 'red' w l, file17 u 1:(Tpt($3)) every 120 title '' lc 'red' w l, file18 u 1:(Tpt($3)) every 120 title '' lc 'red' w l, file19 u 1:(Tpt($3)) every 120 title '' lc 'red' w l, file20 u 1:(Tpt($3)) title '' lc 'red' w l
  52 +
  53 +#replot file16 u 1:($4+22.21) every 120 title 'PD DC Out+22.21 - monitor (ch 103)' lc 'black' w l, file17 u 1:($4+22.21) every 120 title '' lc 'black' w l, file18 u 1:($4+22.21) every 120 title '' lc 'black' w l, file19 u 1:($4+22.21) every 120 title '' lc 'black' w l, file20 u 1:($4+22.21) title '' lc 'black' w l
  54 +
  55 +#replot file16 u 1:(Tpt($5)) every 120 title 'Temp side - monitor (ch 104)' lc 'blue' w l, file17 u 1:(Tpt($5)) every 120 title '' lc 'blue' w l, file18 u 1:(Tpt($5)) every 120 title '' lc 'blue' w l, file19 u 1:(Tpt($5)) every 120 title '' lc 'blue' w l, file20 u 1:(Tpt($5)) title '' lc 'blue' w l
  56 +
  57 +#replot file16 u 1:(T($6)) every 120 title 'Control side - monitor (ch 105)' lc 'purple' w l, file17 u 1:(T($6)) every 120 title '' lc 'purple' w l, file18 u 1:(T($6)) every 120 title '' lc 'purple' w l, file19 u 1:(T($6)) every 120 title '' lc 'purple' w l, file20 u 1:(T($6)) title '' lc 'purple' w l
  58 +
  59 +#plot file1 u 1:(Tpt($3)) title 'Labo (ch 101)' lc 'red' w l, file2 u 1:(Tpt($3)) title '' lc 'red' w l, file3 u 1:(Tpt($3)) title '' lc 'red' w l
  60 +
  61 +#replot file1 u 1:($4+21.41) title 'PD DC Out+21.41 - monitor (ch 103)' lc 'black' w l, file2 u 1:($4+21.41) title '' lc 'black' w l, file3 u 1:($4+21.41) title '' lc 'black' w l
  62 +
  63 +#replot file1 u 1:(Tpt($5)) title 'Temp side - monitor (ch 104)' lc 'blue' w l, file2 u 1:(Tpt($5)) title '' lc 'blue' w l, file3 u 1:(Tpt($5)) title '' lc 'blue' w l
  64 +
  65 +#replot file1 u 1:(T($6)) title 'Control side - monitor (ch 105)' lc 'purple' w l, file2 u 1:(T($6)) title '' lc 'purple' w l, file3 u 1:(T($6)) title '' lc 'purple' w l
  66 +
  67 +
  68 +#plot file1 u 1:(Tpt($3)) every N_dec title 'Labo (ch 101)' lc 'red' w l, file2 u 1:(Tpt($3)) every N_dec title '' lc 'red' w l, file3 u 1:(Tpt($3)) every N_dec title '' lc 'red' w l, file4 u 1:(Tpt($3)) every N_dec title '' lc 'red' w l, file5 u 1:(Tpt($3)) every N_dec title '' lc 'red' w l, file6 u 1:(Tpt($3)) every N_dec title '' lc 'red' w l, file7 u 1:(Tpt($3)) every N_dec title '' lc 'red' w l, file8 u 1:(Tpt($3)) every N_dec title '' lc 'red' w l, file9 u 1:(Tpt($3)) every N_dec title '' lc 'red' w l, file10 u 1:(Tpt($3)) every N_dec title '' lc 'red' w l, file11 u 1:(Tpt($3)) every N_dec title '' lc 'red' w l, file12 u 1:(Tpt($3)) every N_dec title '' lc 'red' w l, file13 u 1:(Tpt($3)) every N_dec title '' lc 'red' w l
  69 +
  70 +#replot file1 u 1:(Tpt($4)-0.3579+0.7) every N_dec title 'Temp side 1 - monitor (ch 102)' lc 'green' w l, file2 u 1:(Tpt($4)-0.3579+0.7) every N_dec title '' lc 'green' w l, file3 u 1:(Tpt($4)-0.3579+0.7) every N_dec title '' lc 'green' w l, file4 u 1:(Tpt($4)-0.3579+0.7) every N_dec title '' lc 'green' w l, file5 u 1:(Tpt($4)-0.3579+0.7) every N_dec title '' lc 'green' w l, file6 u 1:(Tpt($4)-0.3579+0.7) every N_dec title '' lc 'green' w l, file7 u 1:(Tpt($4)-0.3579+0.7) every N_dec title '' lc 'green' w l, file8 u 1:(Tpt($4)-0.3579+0.7) every N_dec title '' lc 'green' w l, file9 u 1:(Tpt($4)-0.3579+0.7) every N_dec title '' lc 'green' w l, file10 u 1:(Tpt($4)-0.3579+0.7) every N_dec title '' lc 'green' w l, file11 u 1:(Tpt($4)-0.3579+0.7) every N_dec title '' lc 'green' w l, file12 u 1:(Tpt($4)-0.3579+0.7) every N_dec title '' lc 'green' w l, file13 u 1:(Tpt($4)-0.3579+0.7) every N_dec title '' lc 'green' w l
  71 +
  72 +#replot file1 u 1:($5+21.75) every N_dec title 'PD DC Out+21.75 - monitor (ch 103)' lc 'black' w l, file2 u 1:($5+21.75) every N_dec title '' lc 'black' w l, file3 u 1:($5+21.75) every N_dec title '' lc 'black' w l, file4 u 1:($5+21.75) every N_dec title '' lc 'black' w l, file5 u 1:($5+21.75) every N_dec title '' lc 'black' w l, file6 u 1:($5+21.75) every N_dec title '' lc 'black' w l, file7 u 1:($5+21.75) every N_dec title '' lc 'black' w l, file8 u 1:($5+21.75) every N_dec title '' lc 'black' w l, file9 u 1:($5+21.75) every N_dec title '' lc 'black' w l, file10 u 1:($5+21.75) every N_dec title '' lc 'black' w l, file11 u 1:($5+21.75) every N_dec title '' lc 'black' w l, file12 u 1:($5+21.75) every N_dec title '' lc 'black' w l, file13 u 1:($5+21.75) every N_dec title '' lc 'black' w l
  73 +
  74 +#replot file1 u 1:(Tpt($6)-0.2128+0.7) every N_dec title 'Temp side 2 - monitor (ch 104)' lc 'blue' w l, file2 u 1:(Tpt($6)-0.2128+0.7) every N_dec title '' lc 'blue' w l, file3 u 1:(Tpt($6)-0.2128+0.7) every N_dec title '' lc 'blue' w l, file4 u 1:(Tpt($6)-0.2128+0.7) every N_dec title '' lc 'blue' w l, file5 u 1:(Tpt($6)-0.2128+0.7) every N_dec title '' lc 'blue' w l, file6 u 1:(Tpt($6)-0.2128+0.7) every N_dec title '' lc 'blue' w l, file7 u 1:(Tpt($6)-0.2128+0.7) every N_dec title '' lc 'blue' w l, file8 u 1:(Tpt($6)-0.2128+0.7) every N_dec title '' lc 'blue' w l, file9 u 1:(Tpt($6)-0.2128+0.7) every N_dec title '' lc 'blue' w l, file10 u 1:(Tpt($6)-0.2128+0.7) every N_dec title '' lc 'blue' w l, file11 u 1:(Tpt($6)-0.2128+0.7) every N_dec title '' lc 'blue' w l, file12 u 1:(Tpt($6)-0.2128+0.7) every N_dec title '' lc 'blue' w l, file13 u 1:(Tpt($6)-0.2128+0.7) every N_dec title '' lc 'blue' w l
  75 +
  76 +#replot file1 u 1:(T($7)+0.7) every N_dec title 'Control side - monitor (ch 105)' lc 'purple' w l, file2 u 1:(T($7)+0.7) every N_dec title '' lc 'purple' w l, file3 u 1:(T($7)+0.7) every N_dec title '' lc 'purple' w l, file4 u 1:(T($7)+0.7) every N_dec title '' lc 'purple' w l, file5 u 1:(T($7)+0.7) every N_dec title '' lc 'purple' w l, file6 u 1:(T($7)+0.7) every N_dec title '' lc 'purple' w l, file7 u 1:(T($7)+0.7) every N_dec title '' lc 'purple' w l, file8 u 1:(T($7)+0.7) every N_dec title '' lc 'purple' w l, file9 u 1:(T($7)+0.7) every N_dec title '' lc 'purple' w l, file10 u 1:(T($7)+0.7) every N_dec title '' lc 'purple' w l, file11 u 1:(T($7)+0.7) every N_dec title '' lc 'purple' w l, file12 u 1:(T($7)+0.7) every N_dec title '' lc 'purple' w l, file13 u 1:(T($7)+0.7) every N_dec title '' lc 'purple' w l
  77 +
  78 +
  79 +#plot file16 u 1:(Tpt($3)) every N_dec title 'Labo (ch 101)' lc 'red' w l
  80 +
  81 +#plot file16 u 1:(Tpt($4)-0.3579+0.7) every N_dec title 'Temp side 1 - monitor (ch 102)' lc 'green' w l
  82 +
  83 +#replot file16 u 1:(Tpt($5)-0.2128+0.7) every N_dec title 'Temp side 2 - monitor (ch 104)' lc 'blue' w l
  84 +
  85 +#replot file16 u 1:(T($6)+0.7) every N_dec title 'Control side - monitor (ch 105)' lc 'purple' w l
  86 +
  87 +
  88 +plot file1 u 1:(T($6)+0.7) every N_dec title 'Control side - monitor (ch 105)' lc 'purple' w l, file2 u 1:(T($6)+0.7) every N_dec title '' lc 'purple' w l, file3 u 1:(T($6)+0.7) every N_dec title '' lc 'purple' w l, file4 u 1:(T($6)+0.7) every N_dec title '' lc 'purple' w l, file5 u 1:(T($6)+0.7) every N_dec title '' lc 'purple' w l, file6 u 1:(T($6)+0.7) every N_dec title '' lc 'purple' w l, file7 u 1:(T($6)+0.7) every N_dec title '' lc 'purple' w l
0 89 \ No newline at end of file
... ... @@ -0,0 +1,43 @@
  1 +import Gnuplot, time
  2 +
  3 +g= Gnuplot.Gnuplot(persist=1)
  4 +
  5 +# init graph
  6 +g('set grid')
  7 +#g('set log y')
  8 +#g('set format y "%.1e"')
  9 +g('set timefmt "%s"')
  10 +g('set format x "%y/%m/%d - %H:%M:%S"')
  11 +g('set xdata time')
  12 +g('set xtics rotate by 90 right')
  13 +
  14 +# set res2temp functions
  15 +g('res2temp627(R) = 0.399341181655472610 + 10.8420092277810909 * (1000./R) + -26.4597939187660813 * (1000./R)**2 + 245.9828566655493379 * (1000./R)**3 + -668.069876596331596 * (1000./R)**4 + 1001.69882618263364 * (1000./R) **5 + -267.272089680656791 * (1000./R)**6')
  16 +g('res2temp625(R) = 0.333548856582638109 + 11.7361551595386118 * (1000./R) + -31.32988932320903987 * (1000./R)**2 + 262.878643524833024 * (1000./R)**3 + -704.163538021035492 * (1000./R)**4 + 1056.6040485650301 * (1000./R)**5 + -307.057196729816496 * (1000./R)**6')
  17 +g('res2temp628(R) = 0.463200932294057566 + 13.5049710820894688 * (1000./R) + -30.5191222755238414 * (1000./R)**2 + 231.098593852017075 * (1000./R)**3 + -550.122691885568202 * (1000./R)**4 + 806.038547554984689 * (1000./R)**5 + -198.510489917360246 * (1000./R)**6')
  18 +g('res2temp16945(R) = 10**(3.2497 * (log10(1000./R))**2 + 5.1777 * log10(1000./R) + 2.499)')
  19 +g('res2temp16943(R) = 10**(3.4738 * (log10(1000./R))**2 + 5.1198 * log10(1000./R) + 2.3681)')
  20 +g('res2temp16944(R) = 10**(3.3674 * (log10(1000./R))**2 + 5.2874 * log10(1000./R) + 2.5165)')
  21 +g('res2temp16941(R) = 10**(2.9486 * (log10(1000./R))**2 + 4.5862 * log10(1000./R) + 2.266)')
  22 +g('res2temp16947(R) = 10**(3.4597 * (log10(1000./R))**2 + 5.2422 * log10(1000./R) + 2.4169)')
  23 +g('pt100(R) = (R/100.-1)/0.003850+273.15')
  24 +
  25 +# plot
  26 +g('plot for [file in system("find ../../../data/2016/ -name 20161*T7*.dat | sort")] file u 1:(pt100($3)) w l lt 1 title "pt100 lab"')
  27 +g('replot for [file in system("find ../../../data/2016/ -name 20161*T7*.dat | sort")] file u 1:(res2temp628($4)) w l lt 2 title "table exp 628"')
  28 +g('replot for [file in system("find ../../../data/2016/ -name 20161*T7*.dat | sort")] file u 1:(res2temp16947($5)) w l lt 3 title "4K 16947"')
  29 +g('replot for [file in system("find ../../../data/2016/ -name 20161*T7*.dat | sort")] file u 1:(res2temp16941($6)) w l lt 4 title "Ecran PT1 16941"')
  30 +g('replot for [file in system("find ../../../data/2016/ -name 20161*LS*.dat | sort")] file u 1:(res2temp16944($3)) w l lt 5 title "Regul 16944"')
  31 +g('replot for [file in system("find ../../../data/2016/ -name 20161*LS*.dat | sort")] file u 1:(res2temp16945($4)) w l lt 6 title "Regul 16945"')
  32 +g('replot for [file in system("find ../../../data/2016/ -name 20161*LS*.dat | sort")] file u 1:(res2temp16943($5)) w l lt 7 title "Regul 16943"')
  33 +g('replot for [file in system("find ../../../data/2016/ -name 20161*LS*.dat | sort")] file u 1:(res2temp627($6)) w l lt 8 title "Regul 627"')
  34 +
  35 +# update
  36 +try:
  37 + while True:
  38 + time.sleep(60)
  39 + print('updating plot\n')
  40 + g('set autoscale')
  41 + g('rep')
  42 +except KeyboardInterrupt:
  43 + g.close()
... ... @@ -0,0 +1,21 @@
  1 +import Gnuplot, time
  2 +
  3 +g= Gnuplot.Gnuplot(persist=1)
  4 +g('set grid')
  5 +g('set log y')
  6 +g('set format y "%.1e"')
  7 +g('set timefmt "%s"')
  8 +g('set format x "%y/%m/%d - %H:%M:%S"')
  9 +g('set xdata time')
  10 +g('set xtics rotate by 90 right')
  11 +
  12 +g('plot for [file in system("find ~/server/data/2016/ -name 20161*TPG*.dat | sort")] file u 1:3 w l lt 1 title file')
  13 +
  14 +try:
  15 + while True:
  16 + time.sleep(60)
  17 + print('updating plot\n')
  18 + g('set autoscale')
  19 + g('rep')
  20 +except KeyboardInterrupt:
  21 + g.close()
disp_local_3days_Regul_LS350.gnu
Diff comments   View file @ 7047b00
... ... @@ -0,0 +1,60 @@
  1 +set grid
  2 +# init graph
  3 +#set log y
  4 +#set format y "%.1e"
  5 +set timefmt "%s"
  6 +set format x "%y/%m/%d - %H:%M:%S"
  7 +set xdata time
  8 +set xtics rotate by 90 right
  9 +
  10 +n = 10*10
  11 +
  12 +filename1 = "20161202"
  13 +filename2 = "20161203"
  14 +filename3 = "20161204"
  15 +
  16 +
  17 +# set res2temp functions
  18 +res2temp627(R) = 0.399341181655472610 + 10.8420092277810909 * (1000./R) + -26.4597939187660813 * (1000./R)**2 + 245.9828566655493379 * (1000./R)**3 + -668.069876596331596 * (1000./R)**4 + 1001.69882618263364 * (1000./R) **5 + -267.272089680656791 * (1000./R)**6
  19 +res2temp625(R) = 0.333548856582638109 + 11.7361551595386118 * (1000./R) + -31.32988932320903987 * (1000./R)**2 + 262.878643524833024 * (1000./R)**3 + -704.163538021035492 * (1000./R)**4 + 1056.6040485650301 * (1000./R)**5 + -307.057196729816496 * (1000./R)**6
  20 +res2temp628(R) = 0.463200932294057566 + 13.5049710820894688 * (1000./R) + -30.5191222755238414 * (1000./R)**2 + 231.098593852017075 * (1000./R)**3 + -550.122691885568202 * (1000./R)**4 + 806.038547554984689 * (1000./R)**5 + -198.510489917360246 * (1000./R)**6
  21 +res2temp16945(R) = 10**(3.2497 * (log10(1000./R))**2 + 5.1777 * log10(1000./R) + 2.499)
  22 +res2temp16943(R) = 10**(3.4738 * (log10(1000./R))**2 + 5.1198 * log10(1000./R) + 2.3681)
  23 +res2temp16944(R) = 10**(3.3674 * (log10(1000./R))**2 + 5.2874 * log10(1000./R) + 2.5165)
  24 +res2temp16941(R) = 10**(2.9486 * (log10(1000./R))**2 + 4.5862 * log10(1000./R) + 2.266)
  25 +res2temp16947(R) = 10**(3.4597 * (log10(1000./R))**2 + 5.2422 * log10(1000./R) + 2.4169)
  26 +pt100(R) = (R/100.-1)/0.003850+273.15
  27 +
  28 +
  29 +############################ regul probe 1
  30 +
  31 +plot for [file in system("find ../../../data/2016/ -name ".filename1."*LS*.dat | sort")] file every n u 1:(res2temp16944($3)) w l lt 5 title "Regul 16944"
  32 +
  33 +replot for [file in system("find ../../../data/2016/ -name ".filename2."*LS*.dat | sort")] file every n u 1:(res2temp16944($3)) w l lt 5 title ""
  34 +
  35 +replot for [file in system("find ../../../data/2016/ -name ".filename3."*LS*.dat | sort")] file every n u 1:(res2temp16944($3)) w l lt 5 title ""
  36 +
  37 +############################ regul probe 2
  38 +
  39 +replot for [file in system("find ../../../data/2016/ -name ".filename1."*LS*.dat | sort")] file every n u 1:(res2temp16945($4)) w l lt 6 title "Regul 16945"
  40 +
  41 +replot for [file in system("find ../../../data/2016/ -name ".filename2."*LS*.dat | sort")] file every n u 1:(res2temp16945($4)) w l lt 6 title ""
  42 +
  43 +replot for [file in system("find ../../../data/2016/ -name ".filename3."*LS*.dat | sort")] file every n u 1:(res2temp16945($4)) w l lt 6 title ""
  44 +
  45 +############################ regul probe 3
  46 +
  47 +replot for [file in system("find ../../../data/2016/ -name ".filename1."*LS*.dat | sort")] file every n u 1:(res2temp16943($5)) w l lt 7 title "Regul 16943"
  48 +
  49 +replot for [file in system("find ../../../data/2016/ -name ".filename2."*LS*.dat | sort")] file every n u 1:(res2temp16943($5)) w l lt 7 title ""
  50 +
  51 +replot for [file in system("find ../../../data/2016/ -name ".filename3."*LS*.dat | sort")] file every n u 1:(res2temp16943($5)) w l lt 7 title ""
  52 +
  53 +############################ regul probe 4 (error signal)
  54 +
  55 +replot for [file in system("find ../../../data/2016/ -name ".filename1."*LS*.dat | sort")] file every n u 1:(res2temp627($6)) w l lt 8 title "Regul 627"
  56 +
  57 +replot for [file in system("find ../../../data/2016/ -name ".filename2."*LS*.dat | sort")] file every n u 1:(res2temp627($6)) w l lt 8 title ""
  58 +
  59 +replot for [file in system("find ../../../data/2016/ -name ".filename3."*LS*.dat | sort")] file every n u 1:(res2temp627($6)) w l lt 8 title ""
  60 +
disp_local_3days_allprobe_LS350-T7Pro.gnu
Diff comments   View file @ 7047b00
... ... @@ -0,0 +1,93 @@
  1 +# init graph
  2 +set grid
  3 +#set log y
  4 +#set format y "%.1e"
  5 +set timefmt "%s"
  6 +set format x "%y/%m/%d - %H:%M:%S"
  7 +set xdata time
  8 +set xtics rotate by 90 right
  9 +
  10 +n = 10*60
  11 +
  12 +filename1 = "20170202"
  13 +filename2 = "20170203"
  14 +filename3 = "20170204"
  15 +
  16 +
  17 +# set res2temp functions
  18 +res2temp627(R) = 0.399341181655472610 + 10.8420092277810909 * (1000./R) + -26.4597939187660813 * (1000./R)**2 + 245.9828566655493379 * (1000./R)**3 + -668.069876596331596 * (1000./R)**4 + 1001.69882618263364 * (1000./R) **5 + -267.272089680656791 * (1000./R)**6
  19 +res2temp625(R) = 0.333548856582638109 + 11.7361551595386118 * (1000./R) + -31.32988932320903987 * (1000./R)**2 + 262.878643524833024 * (1000./R)**3 + -704.163538021035492 * (1000./R)**4 + 1056.6040485650301 * (1000./R)**5 + -307.057196729816496 * (1000./R)**6
  20 +res2temp628(R) = 0.463200932294057566 + 13.5049710820894688 * (1000./R) + -30.5191222755238414 * (1000./R)**2 + 231.098593852017075 * (1000./R)**3 + -550.122691885568202 * (1000./R)**4 + 806.038547554984689 * (1000./R)**5 + -198.510489917360246 * (1000./R)**6
  21 +res2temp16945(R) = 10**(3.2497 * (log10(1000./R))**2 + 5.1777 * log10(1000./R) + 2.499)
  22 +res2temp16943(R) = 10**(3.4738 * (log10(1000./R))**2 + 5.1198 * log10(1000./R) + 2.3681)
  23 +res2temp16944(R) = 10**(3.3674 * (log10(1000./R))**2 + 5.2874 * log10(1000./R) + 2.5165)
  24 +res2temp16941(R) = 10**(2.9486 * (log10(1000./R))**2 + 4.5862 * log10(1000./R) + 2.266)
  25 +res2temp16947(R) = 10**(3.4597 * (log10(1000./R))**2 + 5.2422 * log10(1000./R) + 2.4169)
  26 +pt100(R) = (R/100.-1)/0.003850+273.15
  27 +
  28 +
  29 +############################ Labo
  30 +
  31 +#plot for [file in system("find ../../../data/2017/ -name ".filename1."*T7*.dat | sort")] file every n u 1:(pt100($3)-270) w l lt 1 title "pt100 lab-270K"
  32 +
  33 +#replot for [file in system("find ../../../data/2017/ -name ".filename2."*T7*.dat | sort")] file every n u 1:(pt100($3)-270) w l lt 1 title ""
  34 +
  35 +#replot for [file in system("find ../../../data/2017/ -name ".filename3."*T7*.dat | sort")] file every n u 1:(pt100($3)-270) w l lt 1 title ""
  36 +
  37 +############################ Table exp
  38 +
  39 +plot for [file in system("find ../../../data/2017/ -name ".filename1."*T7*.dat | sort")] file every n u 1:(res2temp628($4)) w l lt 2 title "table exp 628"
  40 +
  41 +replot for [file in system("find ../../../data/2017/ -name ".filename2."*T7*.dat | sort")] file every n u 1:(res2temp628($4)) w l lt 2 title ""
  42 +
  43 +replot for [file in system("find ../../../data/2017/ -name ".filename3."*T7*.dat | sort")] file every n u 1:(res2temp628($4)) w l lt 2 title ""
  44 +
  45 +############################ Tranfer plate 4K
  46 +
  47 +replot for [file in system("find ../../../data/2017/ -name ".filename1."*T7*.dat | sort")] file every n u 1:(res2temp16947($5)) w l lt 3 title "4K 16947"
  48 +
  49 +replot for [file in system("find ../../../data/2017/ -name ".filename2."*T7*.dat | sort")] file every n u 1:(res2temp16947($5)) w l lt 3 title ""
  50 +
  51 +replot for [file in system("find ../../../data/2017/ -name ".filename3."*T7*.dat | sort")] file every n u 1:(res2temp16947($5)) w l lt 3 title ""
  52 +
  53 +############################ Ecran PT1
  54 +
  55 +replot for [file in system("find ../../../data/2017/ -name ".filename1."*T7*.dat | sort")] file every n u 1:(res2temp16941($6)) w l lt 4 title "Ecran PT1 16941"
  56 +
  57 +replot for [file in system("find ../../../data/2017/ -name ".filename2."*T7*.dat | sort")] file every n u 1:(res2temp16941($6)) w l lt 4 title ""
  58 +
  59 +replot for [file in system("find ../../../data/2017/ -name ".filename3."*T7*.dat | sort")] file every n u 1:(res2temp16941($6)) w l lt 4 title ""
  60 +
  61 +
  62 +############################ regul probe 1
  63 +
  64 +#replot for [file in system("find ../../../data/2017/ -name ".filename1."*LS*.dat | sort")] file every n u 1:(res2temp16944($3)) w l lt 5 title "Regul 16944"
  65 +
  66 +#replot for [file in system("find ../../../data/2017/ -name ".filename2."*LS*.dat | sort")] file every n u 1:(res2temp16944($3)) w l lt 5 title ""
  67 +
  68 +#replot for [file in system("find ../../../data/2017/ -name ".filename3."*LS*.dat | sort")] file every n u 1:(res2temp16944($3)) w l lt 5 title ""
  69 +
  70 +############################ regul probe 2
  71 +
  72 +#replot for [file in system("find ../../../data/2017/ -name ".filename1."*LS*.dat | sort")] file every n u 1:(res2temp16945($4)) w l lt 6 title "Regul 16945"
  73 +
  74 +#replot for [file in system("find ../../../data/2017/ -name ".filename2."*LS*.dat | sort")] file every n u 1:(res2temp16945($4)) w l lt 6 title ""
  75 +
  76 +#replot for [file in system("find ../../../data/2017/ -name ".filename3."*LS*.dat | sort")] file every n u 1:(res2temp16945($4)) w l lt 6 title ""
  77 +
  78 +############################ regul probe 3
  79 +
  80 +replot for [file in system("find ../../../data/2017/ -name ".filename1."*LS*.dat | sort")] file every n u 1:(res2temp16943($5)) w l lt 7 title "Regul 16943"
  81 +
  82 +replot for [file in system("find ../../../data/2017/ -name ".filename2."*LS*.dat | sort")] file every n u 1:(res2temp16943($5)) w l lt 7 title ""
  83 +
  84 +replot for [file in system("find ../../../data/2017/ -name ".filename3."*LS*.dat | sort")] file every n u 1:(res2temp16943($5)) w l lt 7 title ""
  85 +
  86 +############################ regul probe 4 (error signal)
  87 +
  88 +replot for [file in system("find ../../../data/2017/ -name ".filename1."*LS*.dat | sort")] file every n u 1:(res2temp627($6)) w l lt 8 title "Regul 627 (error sig)"
  89 +
  90 +replot for [file in system("find ../../../data/2017/ -name ".filename2."*LS*.dat | sort")] file every n u 1:(res2temp627($6)) w l lt 8 title ""
  91 +
  92 +replot for [file in system("find ../../../data/2017/ -name ".filename3."*LS*.dat | sort")] file every n u 1:(res2temp627($6)) w l lt 8 title ""
  93 +
disp_local_5days_labo_T7Pro.gnu
Diff comments   View file @ 7047b00
... ... @@ -0,0 +1,43 @@
  1 +# init graph
  2 +set grid
  3 +#set log y
  4 +#set format y "%.1e"
  5 +set timefmt "%s"
  6 +set format x "%y/%m/%d - %H:%M:%S"
  7 +set xdata time
  8 +set xtics rotate by 90 right
  9 +
  10 +n = 10*10
  11 +
  12 +filename1 = "20161130"
  13 +filename2 = "20161201"
  14 +filename3 = "20161202"
  15 +filename4 = "20161203"
  16 +filename5 = "20161204"
  17 +
  18 +
  19 +# set res2temp functions
  20 +res2temp627(R) = 0.399341181655472610 + 10.8420092277810909 * (1000./R) + -26.4597939187660813 * (1000./R)**2 + 245.9828566655493379 * (1000./R)**3 + -668.069876596331596 * (1000./R)**4 + 1001.69882618263364 * (1000./R) **5 + -267.272089680656791 * (1000./R)**6
  21 +res2temp625(R) = 0.333548856582638109 + 11.7361551595386118 * (1000./R) + -31.32988932320903987 * (1000./R)**2 + 262.878643524833024 * (1000./R)**3 + -704.163538021035492 * (1000./R)**4 + 1056.6040485650301 * (1000./R)**5 + -307.057196729816496 * (1000./R)**6
  22 +res2temp628(R) = 0.463200932294057566 + 13.5049710820894688 * (1000./R) + -30.5191222755238414 * (1000./R)**2 + 231.098593852017075 * (1000./R)**3 + -550.122691885568202 * (1000./R)**4 + 806.038547554984689 * (1000./R)**5 + -198.510489917360246 * (1000./R)**6
  23 +res2temp16945(R) = 10**(3.2497 * (log10(1000./R))**2 + 5.1777 * log10(1000./R) + 2.499)
  24 +res2temp16943(R) = 10**(3.4738 * (log10(1000./R))**2 + 5.1198 * log10(1000./R) + 2.3681)
  25 +res2temp16944(R) = 10**(3.3674 * (log10(1000./R))**2 + 5.2874 * log10(1000./R) + 2.5165)
  26 +res2temp16941(R) = 10**(2.9486 * (log10(1000./R))**2 + 4.5862 * log10(1000./R) + 2.266)
  27 +res2temp16947(R) = 10**(3.4597 * (log10(1000./R))**2 + 5.2422 * log10(1000./R) + 2.4169)
  28 +pt100(R) = (R/100.-1)/0.003850+273.15
  29 +
  30 +
  31 +############################ Labo
  32 +
  33 +plot for [file in system("find ../../../data/2016/ -name ".filename1."*T7*.dat | sort")] file every n u 1:(pt100($3)-273.15) w l lt 1 title "pt100 lab in °C"
  34 +
  35 +replot for [file in system("find ../../../data/2016/ -name ".filename2."*T7*.dat | sort")] file every n u 1:(pt100($3)-273.15) w l lt 1 title ""
  36 +
  37 +replot for [file in system("find ../../../data/2016/ -name ".filename3."*T7*.dat | sort")] file every n u 1:(pt100($3)-273.15) w l lt 1 title ""
  38 +
  39 +replot for [file in system("find ../../../data/2016/ -name ".filename4."*T7*.dat | sort")] file every n u 1:(pt100($3)-273.15) w l lt 1 title ""
  40 +
  41 +replot for [file in system("find ../../../data/2016/ -name ".filename5."*T7*.dat | sort")] file every n u 1:(pt100($3)-273.15) w l lt 1 title ""
  42 +
  43 +
disp_local_LS350-T7Pro.gnu
Diff comments   View file @ 7047b00
... ... @@ -0,0 +1,31 @@
  1 +# init graph
  2 +set grid
  3 +#set log y
  4 +#set format y "%.1e"
  5 +set timefmt "%s"
  6 +set format x "%y/%m/%d - %H:%M:%S"
  7 +set xdata time
  8 +set xtics rotate by 90 right
  9 +
  10 +n = 10*60
  11 +
  12 +# set res2temp functions
  13 +res2temp627(R) = 0.399341181655472610 + 10.8420092277810909 * (1000./R) + -26.4597939187660813 * (1000./R)**2 + 245.9828566655493379 * (1000./R)**3 + -668.069876596331596 * (1000./R)**4 + 1001.69882618263364 * (1000./R) **5 + -267.272089680656791 * (1000./R)**6
  14 +res2temp625(R) = 0.333548856582638109 + 11.7361551595386118 * (1000./R) + -31.32988932320903987 * (1000./R)**2 + 262.878643524833024 * (1000./R)**3 + -704.163538021035492 * (1000./R)**4 + 1056.6040485650301 * (1000./R)**5 + -307.057196729816496 * (1000./R)**6
  15 +res2temp628(R) = 0.463200932294057566 + 13.5049710820894688 * (1000./R) + -30.5191222755238414 * (1000./R)**2 + 231.098593852017075 * (1000./R)**3 + -550.122691885568202 * (1000./R)**4 + 806.038547554984689 * (1000./R)**5 + -198.510489917360246 * (1000./R)**6
  16 +res2temp16945(R) = 10**(3.2497 * (log10(1000./R))**2 + 5.1777 * log10(1000./R) + 2.499)
  17 +res2temp16943(R) = 10**(3.4738 * (log10(1000./R))**2 + 5.1198 * log10(1000./R) + 2.3681)
  18 +res2temp16944(R) = 10**(3.3674 * (log10(1000./R))**2 + 5.2874 * log10(1000./R) + 2.5165)
  19 +res2temp16941(R) = 10**(2.9486 * (log10(1000./R))**2 + 4.5862 * log10(1000./R) + 2.266)
  20 +res2temp16947(R) = 10**(3.4597 * (log10(1000./R))**2 + 5.2422 * log10(1000./R) + 2.4169)
  21 +pt100(R) = (R/100.-1)/0.003850+273.15
  22 +
  23 +# plot
  24 +plot for [file in system("find ../../../data/2016/ -name 20161*T7*.dat | sort")] file every n u 1:(pt100($3)) w l lt 1 title "pt100 lab"
  25 +replot for [file in system("find ../../../data/2016/ -name 20161*T7*.dat | sort")] file every n u 1:(res2temp628($4)) w l lt 2 title "table exp 628"
  26 +replot for [file in system("find ../../../data/2016/ -name 20161*T7*.dat | sort")] file every n u 1:(res2temp16947($5)) w l lt 3 title "4K 16947"
  27 +replot for [file in system("find ../../../data/2016/ -name 20161*T7*.dat | sort")] file every n u 1:(res2temp16941($6)) w l lt 4 title "Ecran PT1 16941"
  28 +replot for [file in system("find ../../../data/2016/ -name 20161*LS*.dat | sort")] file every n u 1:(res2temp16944($3)) w l lt 5 title "Regul 16944"
  29 +replot for [file in system("find ../../../data/2016/ -name 20161*LS*.dat | sort")] file every n u 1:(res2temp16945($4)) w l lt 6 title "Regul 16945"
  30 +replot for [file in system("find ../../../data/2016/ -name 20161*LS*.dat | sort")] file every n u 1:(res2temp16943($5)) w l lt 7 title "Regul 16943"
  31 +replot for [file in system("find ../../../data/2016/ -name 20161*LS*.dat | sort")] file every n u 1:(res2temp627($6)) w l lt 8 title "Regul 627"
disp_local_LS350-T7Pro.py
Diff comments   View file @ 7047b00
... ... @@ -0,0 +1,45 @@
  1 +import Gnuplot, time
  2 +
  3 +g= Gnuplot.Gnuplot(persist=1)
  4 +
  5 +# init graph
  6 +g('set grid')
  7 +#g('set log y')
  8 +#g('set format y "%.1e"')
  9 +g('set timefmt "%s"')
  10 +g('set format x "%y/%m/%d - %H:%M:%S"')
  11 +g('set xdata time')
  12 +g('set xtics rotate by 90 right')
  13 +
  14 +g('n = 540')
  15 +
  16 +# set res2temp functions
  17 +g('res2temp627(R) = 0.399341181655472610 + 10.8420092277810909 * (1000./R) + -26.4597939187660813 * (1000./R)**2 + 245.9828566655493379 * (1000./R)**3 + -668.069876596331596 * (1000./R)**4 + 1001.69882618263364 * (1000./R) **5 + -267.272089680656791 * (1000./R)**6')
  18 +g('res2temp625(R) = 0.333548856582638109 + 11.7361551595386118 * (1000./R) + -31.32988932320903987 * (1000./R)**2 + 262.878643524833024 * (1000./R)**3 + -704.163538021035492 * (1000./R)**4 + 1056.6040485650301 * (1000./R)**5 + -307.057196729816496 * (1000./R)**6')
  19 +g('res2temp628(R) = 0.463200932294057566 + 13.5049710820894688 * (1000./R) + -30.5191222755238414 * (1000./R)**2 + 231.098593852017075 * (1000./R)**3 + -550.122691885568202 * (1000./R)**4 + 806.038547554984689 * (1000./R)**5 + -198.510489917360246 * (1000./R)**6')
  20 +g('res2temp16945(R) = 10**(3.2497 * (log10(1000./R))**2 + 5.1777 * log10(1000./R) + 2.499)')
  21 +g('res2temp16943(R) = 10**(3.4738 * (log10(1000./R))**2 + 5.1198 * log10(1000./R) + 2.3681)')
  22 +g('res2temp16944(R) = 10**(3.3674 * (log10(1000./R))**2 + 5.2874 * log10(1000./R) + 2.5165)')
  23 +g('res2temp16941(R) = 10**(2.9486 * (log10(1000./R))**2 + 4.5862 * log10(1000./R) + 2.266)')
  24 +g('res2temp16947(R) = 10**(3.4597 * (log10(1000./R))**2 + 5.2422 * log10(1000./R) + 2.4169)')
  25 +g('pt100(R) = (R/100.-1)/0.003850+273.15')
  26 +
  27 +# plot
  28 +g('plot for [file in system("find ../../../data/2016/ -name 20161*T7*.dat | sort")] file every n u 1:(pt100($3)) w l lt 1 title "pt100 lab"')
  29 +g('replot for [file in system("find ../../../data/2016/ -name 20161*T7*.dat | sort")] file every n u 1:(res2temp628($4)) w l lt 2 title "table exp 628"')
  30 +g('replot for [file in system("find ../../../data/2016/ -name 20161*T7*.dat | sort")] file every n u 1:(res2temp16947($5)) w l lt 3 title "4K 16947"')
  31 +g('replot for [file in system("find ../../../data/2016/ -name 20161*T7*.dat | sort")] file every n u 1:(res2temp16941($6)) w l lt 4 title "Ecran PT1 16941"')
  32 +g('replot for [file in system("find ../../../data/2016/ -name 20161*LS*.dat | sort")] file every n u 1:(res2temp16944($3)) w l lt 5 title "Regul 16944"')
  33 +g('replot for [file in system("find ../../../data/2016/ -name 20161*LS*.dat | sort")] file every n u 1:(res2temp16945($4)) w l lt 6 title "Regul 16945"')
  34 +g('replot for [file in system("find ../../../data/2016/ -name 20161*LS*.dat | sort")] file every n u 1:(res2temp16943($5)) w l lt 7 title "Regul 16943"')
  35 +g('replot for [file in system("find ../../../data/2016/ -name 20161*LS*.dat | sort")] file every n u 1:(res2temp627($6)) w l lt 8 title "Regul 627"')
  36 +
  37 +# update
  38 +try:
  39 + while True:
  40 + time.sleep(60)
  41 +# print('updating plot\n')
  42 +# g('set autoscale')
  43 +# g('rep')
  44 +except KeyboardInterrupt:
  45 + g.close()
disp_local_singleday_LS350-T7Pro.gnu
Diff comments   View file @ 7047b00
... ... @@ -0,0 +1,58 @@
  1 +# init graph
  2 +set grid
  3 +#set log y
  4 +#set format y "%.1e"
  5 +set timefmt "%s"
  6 +set format x "%y/%m/%d - %H:%M:%S"
  7 +set xdata time
  8 +set xtics rotate by 90 right
  9 +
  10 +n = 10*1
  11 +
  12 +filename1 = "20170129"
  13 +
  14 +
  15 +# set res2temp functions
  16 +res2temp627(R) = 0.399341181655472610 + 10.8420092277810909 * (1000./R) + -26.4597939187660813 * (1000./R)**2 + 245.9828566655493379 * (1000./R)**3 + -668.069876596331596 * (1000./R)**4 + 1001.69882618263364 * (1000./R) **5 + -267.272089680656791 * (1000./R)**6
  17 +res2temp625(R) = 0.333548856582638109 + 11.7361551595386118 * (1000./R) + -31.32988932320903987 * (1000./R)**2 + 262.878643524833024 * (1000./R)**3 + -704.163538021035492 * (1000./R)**4 + 1056.6040485650301 * (1000./R)**5 + -307.057196729816496 * (1000./R)**6
  18 +res2temp628(R) = 0.463200932294057566 + 13.5049710820894688 * (1000./R) + -30.5191222755238414 * (1000./R)**2 + 231.098593852017075 * (1000./R)**3 + -550.122691885568202 * (1000./R)**4 + 806.038547554984689 * (1000./R)**5 + -198.510489917360246 * (1000./R)**6
  19 +res2temp16945(R) = 10**(3.2497 * (log10(1000./R))**2 + 5.1777 * log10(1000./R) + 2.499)
  20 +res2temp16943(R) = 10**(3.4738 * (log10(1000./R))**2 + 5.1198 * log10(1000./R) + 2.3681)
  21 +res2temp16944(R) = 10**(3.3674 * (log10(1000./R))**2 + 5.2874 * log10(1000./R) + 2.5165)
  22 +res2temp16941(R) = 10**(2.9486 * (log10(1000./R))**2 + 4.5862 * log10(1000./R) + 2.266)
  23 +res2temp16947(R) = 10**(3.4597 * (log10(1000./R))**2 + 5.2422 * log10(1000./R) + 2.4169)
  24 +pt100(R) = (R/100.-1)/0.003850+273.15
  25 +
  26 +
  27 +############################ Labo (shifted)
  28 +
  29 +plot for [file in system("find ../../../data/ -name ".filename1."*T7*.dat | sort")] file every n u 1:(pt100($3)) w l lt 1 title "pt100 lab"
  30 +
  31 +############################ Table exp
  32 +
  33 +replot for [file in system("find ../../../data/ -name ".filename1."*T7*.dat | sort")] file every n u 1:(res2temp628($4)) w l lt 2 title "table exp 628"
  34 +
  35 +############################ Tranfer plate 4K
  36 +
  37 +replot for [file in system("find ../../../data/ -name ".filename1."*T7*.dat | sort")] file every n u 1:(res2temp16947($5)) w l lt 3 title "4K 16947"
  38 +
  39 +############################ Ecran PT1
  40 +
  41 +replot for [file in system("find ../../../data/ -name ".filename1."*T7*.dat | sort")] file every n u 1:(res2temp16941($6)) w l lt 4 title "Ecran PT1 16941"
  42 +
  43 +############################ regul probe 1
  44 +
  45 +replot for [file in system("find ../../../data/ -name ".filename1."*LS*.dat | sort")] file every n u 1:(res2temp16944($3)) w l lt 5 title "Regul 16944"
  46 +
  47 +############################ regul probe 2
  48 +
  49 +replot for [file in system("find ../../../data/ -name ".filename1."*LS*.dat | sort")] file every n u 1:(res2temp16945($4)) w l lt 6 title "Regul 16945"
  50 +
  51 +############################ regul probe 3
  52 +
  53 +replot for [file in system("find ../../../data/ -name ".filename1."*LS*.dat | sort")] file every n u 1:(res2temp16943($5)) w l lt 7 title "Regul 16943"
  54 +
  55 +############################ regul probe 4 (error signal)
  56 +
  57 +replot for [file in system("find ../../../data/ -name ".filename1."*LS*.dat | sort")] file every n u 1:(res2temp627($6)) w l lt 8 title "Regul 627"
  58 +
... ... @@ -0,0 +1,34 @@
  1 +clear
  2 +
  3 +if (!exists("name")) print('No files to plot') exit
  4 +if (!exists("unit")) unit='Y'
  5 +if (!exists("path")) path='.'
  6 +if (!exists("ncol")) ncol=1
  7 +if (!exists("sample")) sample=1
  8 +if (exists("logy")) set log y
  9 +if (exists("logx")) set log x
  10 +if (exists("logxy")) set log x y
  11 +
  12 +set grid x y
  13 +unset key
  14 +#set key outside
  15 +
  16 +set timefmt "%s"
  17 +set format x "%y/%m/%d - %H:%M:%S"
  18 +set xdata time
  19 +set xtics rotate by 90 right
  20 +set xlabel "GMT time"
  21 +set ylabel unit
  22 +set ytics nomirror
  23 +set tics out
  24 +set autoscale y
  25 +
  26 +list = system('find '.path.' -name '.name.' | sort -d' )
  27 +print(list)
  28 +plot for [file in list] for [col=3:(ncol+2)] file u 1:col every sample w l ls col-2
  29 +
  30 +pause -1
  31 +
  32 +
  33 +# set format y "%1.1e"
  34 +