saisons.gnuplot

set terminal pngcairo size 800,600
set output "/tmp/temps.png"
stats "time_era5_t2m_France_metropolitan_mon12_ave3_dump0.txt" using ($3+$4/12-2000):($2) name "st0"
w0(x)=st0_intercept+st0_slope*x
title0=sprintf("hiver (régr. %.2f°+%.4f×(Y−2000))", st0_intercept, st0_slope)
stats "time_era5_t2m_France_metropolitan_mon3_ave3_dump0.txt" using ($3+$4/12-2000):($2) name "st1"
w1(x)=st1_intercept+st1_slope*x
title1=sprintf("printemps (régr. %.2f°+%.4f×(Y−2000))", st1_intercept, st1_slope)
stats "time_era5_t2m_France_metropolitan_mon6_ave3_dump0.txt" using ($3+$4/12-2000):($2) name "st2"
w2(x)=st2_intercept+st2_slope*x
title2=sprintf("été (régr. %.2f°+%.4f×(Y−2000))", st2_intercept, st2_slope)
stats "time_era5_t2m_France_metropolitan_mon9_ave3_dump0.txt" using ($3+$4/12-2000):($2) name "st3"
w3(x)=st3_intercept+st3_slope*x
title3=sprintf("automne (régr. %.2f°+%.4f×(Y−2000))", st3_intercept, st3_slope)
plot [1950:2022] w0(x-2000) lt rgb "#8080ff" dashtype 2 notitle, [1950:2022] w1(x-2000) lt rgb "#80ff80" dashtype 2 notitle, [1950:2022] w2(x-2000) lt rgb "#ff8080" dashtype 2 notitle, [1950:2022] w3(x-2000) lt rgb "#ffc080" dashtype 2 notitle, "time_era5_t2m_France_metropolitan_mon12_ave3_dump0.txt" using ($3+$4/12):($2) with lines lt rgb "#8080ff" title title0, "time_era5_t2m_France_metropolitan_mon3_ave3_dump0.txt" using ($3+$4/12):($2) with lines lt rgb "#80ff80" title title1, "time_era5_t2m_France_metropolitan_mon6_ave3_dump0.txt" using ($3+$4/12):($2) with lines lt rgb "#ff8080" title title2, "time_era5_t2m_France_metropolitan_mon9_ave3_dump0.txt" using ($3+$4/12):($2) with lines lt rgb "#ffc080" title title3

Voir https://twitter.com/gro_tsen/status/1558640335617609728 pour explications et commentaires