From 7047b0001eaada1e95b448d30b8b6aa0ef020200 Mon Sep 17 00:00:00 2001 From: bma Date: Thu, 12 Oct 2017 08:55:19 +0200 Subject: [PATCH] initial commit --- OFC-FRR/temp_mon01.gnu | 88 +++++++++++++++++++++++++++++ disp_LS350-T7Pro.py | 43 ++++++++++++++ disp_TPG.py | 21 +++++++ disp_local_3days_Regul_LS350.gnu | 60 ++++++++++++++++++++ disp_local_3days_allprobe_LS350-T7Pro.gnu | 93 +++++++++++++++++++++++++++++++ disp_local_5days_labo_T7Pro.gnu | 43 ++++++++++++++ disp_local_LS350-T7Pro.gnu | 31 +++++++++++ disp_local_LS350-T7Pro.py | 45 +++++++++++++++ disp_local_singleday_LS350-T7Pro.gnu | 58 +++++++++++++++++++ plotAll.gnu | 34 +++++++++++ 10 files changed, 516 insertions(+) create mode 100644 OFC-FRR/temp_mon01.gnu create mode 100644 disp_LS350-T7Pro.py create mode 100644 disp_TPG.py create mode 100644 disp_local_3days_Regul_LS350.gnu create mode 100644 disp_local_3days_allprobe_LS350-T7Pro.gnu create mode 100644 disp_local_5days_labo_T7Pro.gnu create mode 100644 disp_local_LS350-T7Pro.gnu create mode 100644 disp_local_LS350-T7Pro.py create mode 100644 disp_local_singleday_LS350-T7Pro.gnu create mode 100755 plotAll.gnu diff --git a/OFC-FRR/temp_mon01.gnu b/OFC-FRR/temp_mon01.gnu new file mode 100644 index 0000000..ed847ed --- /dev/null +++ b/OFC-FRR/temp_mon01.gnu @@ -0,0 +1,88 @@ +set grid +# init graph +#set log y +#set format y "%.1e" +set timefmt "%s" +set format x "%y/%m/%d - %H:%M:%S" +set xdata time +set xtics rotate by 90 right + + + + +file1 = '/home/sicav/server/data/2017/2017-03/20170308-150806-AG34972A.dat' +file2 = '/home/sicav/server/data/2017/2017-03/20170309-150807-AG34972A.dat' +file3 = '/home/sicav/server/data/2017/2017-03/20170310-150807-AG34972A.dat' +file4 = '/home/sicav/server/data/2017/2017-03/20170311-150807-AG34972A.dat' +file5 = '/home/sicav/server/data/2017/2017-03/20170312-150808-AG34972A.dat' +file6 = '/home/sicav/server/data/2017/2017-03/20170313-150808-AG34972A.dat' +file7 = '/home/sicav/server/data/2017/2017-03/20170314-150809-AG34972A.dat' + +N_dec =1 +Tamb = 273.15 + + +# formula for NTC 10k: T = 1 / ( A + B * ln(R) + C * ln(R)^3) +A = 2.108508173e-3 +B = 0.7979204727e-4 +C = 6.535076315e-7 +T(R) = 1 / (A + B*log(R) + C*log(R)**(3) ) - Tamb + +#formula for PT 100 Ohms in °C: +kpt = 0.003916 +Tpt(R) = ((R-100)/(100*kpt))-0.9264 + +#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 + +#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 + +#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 + + +#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 + +#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 + +#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 + +#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 + + +#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 + +#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 + +#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 + +#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 + +#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 + +#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 + +#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 + +#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 + + +#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 + +#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 + +#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 + +#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 + +#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 + + +#plot file16 u 1:(Tpt($3)) every N_dec title 'Labo (ch 101)' lc 'red' w l + +#plot file16 u 1:(Tpt($4)-0.3579+0.7) every N_dec title 'Temp side 1 - monitor (ch 102)' lc 'green' w l + +#replot file16 u 1:(Tpt($5)-0.2128+0.7) every N_dec title 'Temp side 2 - monitor (ch 104)' lc 'blue' w l + +#replot file16 u 1:(T($6)+0.7) every N_dec title 'Control side - monitor (ch 105)' lc 'purple' w l + + +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 \ No newline at end of file diff --git a/disp_LS350-T7Pro.py b/disp_LS350-T7Pro.py new file mode 100644 index 0000000..6fe5a9d --- /dev/null +++ b/disp_LS350-T7Pro.py @@ -0,0 +1,43 @@ +import Gnuplot, time + +g= Gnuplot.Gnuplot(persist=1) + +# init graph +g('set grid') +#g('set log y') +#g('set format y "%.1e"') +g('set timefmt "%s"') +g('set format x "%y/%m/%d - %H:%M:%S"') +g('set xdata time') +g('set xtics rotate by 90 right') + +# set res2temp functions +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') +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') +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') +g('res2temp16945(R) = 10**(3.2497 * (log10(1000./R))**2 + 5.1777 * log10(1000./R) + 2.499)') +g('res2temp16943(R) = 10**(3.4738 * (log10(1000./R))**2 + 5.1198 * log10(1000./R) + 2.3681)') +g('res2temp16944(R) = 10**(3.3674 * (log10(1000./R))**2 + 5.2874 * log10(1000./R) + 2.5165)') +g('res2temp16941(R) = 10**(2.9486 * (log10(1000./R))**2 + 4.5862 * log10(1000./R) + 2.266)') +g('res2temp16947(R) = 10**(3.4597 * (log10(1000./R))**2 + 5.2422 * log10(1000./R) + 2.4169)') +g('pt100(R) = (R/100.-1)/0.003850+273.15') + +# plot +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"') +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"') +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"') +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"') +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"') +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"') +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"') +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"') + +# update +try: + while True: + time.sleep(60) + print('updating plot\n') + g('set autoscale') + g('rep') +except KeyboardInterrupt: + g.close() diff --git a/disp_TPG.py b/disp_TPG.py new file mode 100644 index 0000000..5fdcd93 --- /dev/null +++ b/disp_TPG.py @@ -0,0 +1,21 @@ +import Gnuplot, time + +g= Gnuplot.Gnuplot(persist=1) +g('set grid') +g('set log y') +g('set format y "%.1e"') +g('set timefmt "%s"') +g('set format x "%y/%m/%d - %H:%M:%S"') +g('set xdata time') +g('set xtics rotate by 90 right') + +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') + +try: + while True: + time.sleep(60) + print('updating plot\n') + g('set autoscale') + g('rep') +except KeyboardInterrupt: + g.close() diff --git a/disp_local_3days_Regul_LS350.gnu b/disp_local_3days_Regul_LS350.gnu new file mode 100644 index 0000000..cbd3c55 --- /dev/null +++ b/disp_local_3days_Regul_LS350.gnu @@ -0,0 +1,60 @@ +set grid +# init graph +#set log y +#set format y "%.1e" +set timefmt "%s" +set format x "%y/%m/%d - %H:%M:%S" +set xdata time +set xtics rotate by 90 right + +n = 10*10 + +filename1 = "20161202" +filename2 = "20161203" +filename3 = "20161204" + + +# set res2temp functions +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 +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 +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 +res2temp16945(R) = 10**(3.2497 * (log10(1000./R))**2 + 5.1777 * log10(1000./R) + 2.499) +res2temp16943(R) = 10**(3.4738 * (log10(1000./R))**2 + 5.1198 * log10(1000./R) + 2.3681) +res2temp16944(R) = 10**(3.3674 * (log10(1000./R))**2 + 5.2874 * log10(1000./R) + 2.5165) +res2temp16941(R) = 10**(2.9486 * (log10(1000./R))**2 + 4.5862 * log10(1000./R) + 2.266) +res2temp16947(R) = 10**(3.4597 * (log10(1000./R))**2 + 5.2422 * log10(1000./R) + 2.4169) +pt100(R) = (R/100.-1)/0.003850+273.15 + + +############################ regul probe 1 + +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" + +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 "" + +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 "" + +############################ regul probe 2 + +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" + +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 "" + +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 "" + +############################ regul probe 3 + +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" + +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 "" + +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 "" + +############################ regul probe 4 (error signal) + +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" + +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 "" + +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 "" + diff --git a/disp_local_3days_allprobe_LS350-T7Pro.gnu b/disp_local_3days_allprobe_LS350-T7Pro.gnu new file mode 100644 index 0000000..8968312 --- /dev/null +++ b/disp_local_3days_allprobe_LS350-T7Pro.gnu @@ -0,0 +1,93 @@ +# init graph +set grid +#set log y +#set format y "%.1e" +set timefmt "%s" +set format x "%y/%m/%d - %H:%M:%S" +set xdata time +set xtics rotate by 90 right + +n = 10*60 + +filename1 = "20170202" +filename2 = "20170203" +filename3 = "20170204" + + +# set res2temp functions +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 +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 +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 +res2temp16945(R) = 10**(3.2497 * (log10(1000./R))**2 + 5.1777 * log10(1000./R) + 2.499) +res2temp16943(R) = 10**(3.4738 * (log10(1000./R))**2 + 5.1198 * log10(1000./R) + 2.3681) +res2temp16944(R) = 10**(3.3674 * (log10(1000./R))**2 + 5.2874 * log10(1000./R) + 2.5165) +res2temp16941(R) = 10**(2.9486 * (log10(1000./R))**2 + 4.5862 * log10(1000./R) + 2.266) +res2temp16947(R) = 10**(3.4597 * (log10(1000./R))**2 + 5.2422 * log10(1000./R) + 2.4169) +pt100(R) = (R/100.-1)/0.003850+273.15 + + +############################ Labo + +#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" + +#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 "" + +#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 "" + +############################ Table exp + +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" + +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 "" + +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 "" + +############################ Tranfer plate 4K + +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" + +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 "" + +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 "" + +############################ Ecran PT1 + +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" + +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 "" + +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 "" + + +############################ regul probe 1 + +#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" + +#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 "" + +#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 "" + +############################ regul probe 2 + +#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" + +#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 "" + +#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 "" + +############################ regul probe 3 + +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" + +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 "" + +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 "" + +############################ regul probe 4 (error signal) + +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)" + +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 "" + +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 "" + diff --git a/disp_local_5days_labo_T7Pro.gnu b/disp_local_5days_labo_T7Pro.gnu new file mode 100644 index 0000000..d4c16a7 --- /dev/null +++ b/disp_local_5days_labo_T7Pro.gnu @@ -0,0 +1,43 @@ +# init graph +set grid +#set log y +#set format y "%.1e" +set timefmt "%s" +set format x "%y/%m/%d - %H:%M:%S" +set xdata time +set xtics rotate by 90 right + +n = 10*10 + +filename1 = "20161130" +filename2 = "20161201" +filename3 = "20161202" +filename4 = "20161203" +filename5 = "20161204" + + +# set res2temp functions +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 +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 +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 +res2temp16945(R) = 10**(3.2497 * (log10(1000./R))**2 + 5.1777 * log10(1000./R) + 2.499) +res2temp16943(R) = 10**(3.4738 * (log10(1000./R))**2 + 5.1198 * log10(1000./R) + 2.3681) +res2temp16944(R) = 10**(3.3674 * (log10(1000./R))**2 + 5.2874 * log10(1000./R) + 2.5165) +res2temp16941(R) = 10**(2.9486 * (log10(1000./R))**2 + 4.5862 * log10(1000./R) + 2.266) +res2temp16947(R) = 10**(3.4597 * (log10(1000./R))**2 + 5.2422 * log10(1000./R) + 2.4169) +pt100(R) = (R/100.-1)/0.003850+273.15 + + +############################ Labo + +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" + +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 "" + +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 "" + +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 "" + +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 "" + + diff --git a/disp_local_LS350-T7Pro.gnu b/disp_local_LS350-T7Pro.gnu new file mode 100644 index 0000000..b9d766f --- /dev/null +++ b/disp_local_LS350-T7Pro.gnu @@ -0,0 +1,31 @@ +# init graph +set grid +#set log y +#set format y "%.1e" +set timefmt "%s" +set format x "%y/%m/%d - %H:%M:%S" +set xdata time +set xtics rotate by 90 right + +n = 10*60 + +# set res2temp functions +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 +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 +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 +res2temp16945(R) = 10**(3.2497 * (log10(1000./R))**2 + 5.1777 * log10(1000./R) + 2.499) +res2temp16943(R) = 10**(3.4738 * (log10(1000./R))**2 + 5.1198 * log10(1000./R) + 2.3681) +res2temp16944(R) = 10**(3.3674 * (log10(1000./R))**2 + 5.2874 * log10(1000./R) + 2.5165) +res2temp16941(R) = 10**(2.9486 * (log10(1000./R))**2 + 4.5862 * log10(1000./R) + 2.266) +res2temp16947(R) = 10**(3.4597 * (log10(1000./R))**2 + 5.2422 * log10(1000./R) + 2.4169) +pt100(R) = (R/100.-1)/0.003850+273.15 + +# plot +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" +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" +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" +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" +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" +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" +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" +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" diff --git a/disp_local_LS350-T7Pro.py b/disp_local_LS350-T7Pro.py new file mode 100644 index 0000000..6121568 --- /dev/null +++ b/disp_local_LS350-T7Pro.py @@ -0,0 +1,45 @@ +import Gnuplot, time + +g= Gnuplot.Gnuplot(persist=1) + +# init graph +g('set grid') +#g('set log y') +#g('set format y "%.1e"') +g('set timefmt "%s"') +g('set format x "%y/%m/%d - %H:%M:%S"') +g('set xdata time') +g('set xtics rotate by 90 right') + +g('n = 540') + +# set res2temp functions +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') +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') +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') +g('res2temp16945(R) = 10**(3.2497 * (log10(1000./R))**2 + 5.1777 * log10(1000./R) + 2.499)') +g('res2temp16943(R) = 10**(3.4738 * (log10(1000./R))**2 + 5.1198 * log10(1000./R) + 2.3681)') +g('res2temp16944(R) = 10**(3.3674 * (log10(1000./R))**2 + 5.2874 * log10(1000./R) + 2.5165)') +g('res2temp16941(R) = 10**(2.9486 * (log10(1000./R))**2 + 4.5862 * log10(1000./R) + 2.266)') +g('res2temp16947(R) = 10**(3.4597 * (log10(1000./R))**2 + 5.2422 * log10(1000./R) + 2.4169)') +g('pt100(R) = (R/100.-1)/0.003850+273.15') + +# plot +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"') +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"') +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"') +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"') +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"') +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"') +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"') +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"') + +# update +try: + while True: + time.sleep(60) +# print('updating plot\n') +# g('set autoscale') +# g('rep') +except KeyboardInterrupt: + g.close() diff --git a/disp_local_singleday_LS350-T7Pro.gnu b/disp_local_singleday_LS350-T7Pro.gnu new file mode 100644 index 0000000..5766f0c --- /dev/null +++ b/disp_local_singleday_LS350-T7Pro.gnu @@ -0,0 +1,58 @@ +# init graph +set grid +#set log y +#set format y "%.1e" +set timefmt "%s" +set format x "%y/%m/%d - %H:%M:%S" +set xdata time +set xtics rotate by 90 right + +n = 10*1 + +filename1 = "20170129" + + +# set res2temp functions +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 +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 +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 +res2temp16945(R) = 10**(3.2497 * (log10(1000./R))**2 + 5.1777 * log10(1000./R) + 2.499) +res2temp16943(R) = 10**(3.4738 * (log10(1000./R))**2 + 5.1198 * log10(1000./R) + 2.3681) +res2temp16944(R) = 10**(3.3674 * (log10(1000./R))**2 + 5.2874 * log10(1000./R) + 2.5165) +res2temp16941(R) = 10**(2.9486 * (log10(1000./R))**2 + 4.5862 * log10(1000./R) + 2.266) +res2temp16947(R) = 10**(3.4597 * (log10(1000./R))**2 + 5.2422 * log10(1000./R) + 2.4169) +pt100(R) = (R/100.-1)/0.003850+273.15 + + +############################ Labo (shifted) + +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" + +############################ Table exp + +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" + +############################ Tranfer plate 4K + +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" + +############################ Ecran PT1 + +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" + +############################ regul probe 1 + +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" + +############################ regul probe 2 + +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" + +############################ regul probe 3 + +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" + +############################ regul probe 4 (error signal) + +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" + diff --git a/plotAll.gnu b/plotAll.gnu new file mode 100755 index 0000000..fdee954 --- /dev/null +++ b/plotAll.gnu @@ -0,0 +1,34 @@ +clear + +if (!exists("name")) print('No files to plot') exit +if (!exists("unit")) unit='Y' +if (!exists("path")) path='.' +if (!exists("ncol")) ncol=1 +if (!exists("sample")) sample=1 +if (exists("logy")) set log y +if (exists("logx")) set log x +if (exists("logxy")) set log x y + +set grid x y +unset key +#set key outside + +set timefmt "%s" +set format x "%y/%m/%d - %H:%M:%S" +set xdata time +set xtics rotate by 90 right +set xlabel "GMT time" +set ylabel unit +set ytics nomirror +set tics out +set autoscale y + +list = system('find '.path.' -name '.name.' | sort -d' ) +print(list) +plot for [file in list] for [col=3:(ncol+2)] file u 1:col every sample w l ls col-2 + +pause -1 + + +# set format y "%1.1e" + -- 2.16.4