iccsa-21-wind

Wind Simulation Using High-Frequency Velocity Component Measurements
git clone https://git.igankevich.com/iccsa-21-wind.git
Log | Files | Refs

daily-stats.gnuplot (1157B)


      1 load 'gnuplot/paired.pal'
      2 set terminal svg size 500,400 font 'Free Serif, 10' enhanced dynamic round
      3 set xtics nomirror out rotate 90 offset 0,0.25
      4 set ytics nomirror out offset 0.5,0
      5 set border 1+2
      6 set key top center outside maxrows 1
      7 
      8 stats 'build/daily-stats' using 6
      9 
     10 # time axis
     11 set xdata time
     12 set timefmt "%s"
     13 set format x "%Y-%m-%d"
     14 set xtics 0,60*60*24 
     15 set mxtics 1
     16 localtime(t) = t+3*60*60
     17 
     18 set output 'build/daily-stats.svg'
     19 
     20 do for [timestamp in "1614546000 1614978000 1615496400 1616187600 1616274000 1616360400 1616533200"] {
     21     t = int(timestamp)
     22     set obj t rectangle behind from first localtime(t), first 0 to first localtime(t)+24*60*60, graph 1 back fillstyle solid 1.0 fillcolor '#f0f0c0' lw 0
     23 }
     24 
     25 plot \
     26          'build/daily-stats' using (localtime($2)):3 with lines lc '#c04040' title 'υ_x',\
     27          'build/daily-stats' using (localtime($2)):4 with lines lc '#40c040' title 'υ_y',\
     28          'build/daily-stats' using (localtime($2)):5 with lines lc '#4040c0' title 'υ_z',\
     29          'build/daily-stats' using (localtime($2)):6 with lines lc '#404040' title 'υ',\
     30          STATS_mean with lines lc '#808080' dashtype 2 notitle
     31