bmarechal / gnuplot_scripts

Blame view

disp_local_3days_Regul_LS350.gnu 3.43 KB
edit raw normal view history
7047b0001   bmarechal   initial commit 
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
51
52
53
54
55
56
57
58
59
 
  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 ""