saisons.gnuplot

set terminal pngcairo size 1024,768
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
res0_stddev=sqrt(1-st0_correlation**2)*st0_stddev_y
title0=sprintf("hiver (régr. %.2f° + %.2f×(Y−2000)/100 ; σ=%.2f)", st0_intercept, st0_slope*100, res0_stddev)
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
res1_stddev=sqrt(1-st1_correlation**2)*st1_stddev_y
title1=sprintf("printemps (régr. %.2f° + %.2f×(Y−2000)/100 ; σ=%.2f)", st1_intercept, st1_slope*100, res1_stddev)
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
res2_stddev=sqrt(1-st2_correlation**2)*st2_stddev_y
title2=sprintf("été (régr. %.2f° + %.2f×(Y−2000)/100 ; σ=%.2f)", st2_intercept, st2_slope*100, res2_stddev)
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
res3_stddev=sqrt(1-st3_correlation**2)*st3_stddev_y
title3=sprintf("automne (régr. %.2f° + %.2f×(Y−2000)/100 ; σ=%.2f)", st3_intercept, st3_slope*100, res3_stddev)
plot [1950:2022] "+" using 1:(w0($1-2000)-res0_stddev):(w0($1-2000)+res0_stddev) with filledcurves lc rgb "#8080ff" fs transparent solid 0.08 notitle, [1950:2022] "+" using 1:(w1($1-2000)-res1_stddev):(w1($1-2000)+res1_stddev) with filledcurves lc rgb "#80ff80" fs transparent solid 0.08 notitle, [1950:2022] "+" using 1:(w2($1-2000)-res2_stddev):(w2($1-2000)+res2_stddev) with filledcurves lc rgb "#ff8080" fs transparent solid 0.08 notitle, [1950:2022] "+" using 1:(w3($1-2000)-res3_stddev):(w3($1-2000)+res3_stddev) with filledcurves lc rgb "#ffc080" fs transparent solid 0.08 notitle, [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