cngf-pf-exp1

commit 50440e41d31728921faf47fd9cbc5cb53352b054
parent af37119c92536999d088c37d99fe6764322051d4
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Wed, 18 Dec 2019 20:59:01 +0100

Plot rate flux as m^2/year

Diffstat:
M
fig-hysteresis_rate/fig.gp
 | 
18
+++++++++++-------

1 file changed, 11 insertions(+), 7 deletions(-)
diff --git a/fig-hysteresis_rate/fig.gp b/fig-hysteresis_rate/fig.gp
@@ -40,15 +40,19 @@ plot "" u ($2/1000):(8.0 - $6) w l lt 4 t ""
 
 set xtics format "%g"
 set xlabel "{/:Normal Effective normal stress [kPa]}"
-set ytics 3e-5
-set yrange [0:1e-4]
+to_year=60*60*24*365
+#set ytics 3e-5*to_year
+#set yrange [0:1e-4*to_year]
+set ytics 1000
+set yrange [0:3000]
 unset arrow
-set arrow from 135,1.5e-5 to 100,1e-5 head
-set arrow from 110,4.0e-5 to 140,5.5e-5 head
-set arrow from 165.3,4e-5 to 165,1.5e-5 head
+set arrow from 135,1.5e-5*to_year to 100,1e-5*to_year head
+set arrow from 110,4.0e-5*to_year to 140,5.5e-5*to_year head
+set arrow from 165.3,4e-5*to_year to 165,1.5e-5*to_year head
 unset y2tics
 unset y2label
 unset label
 set label "{/:Bold f}" at graph 0.01,0.93
-set ylabel "{/:Normal Specific sediment flux [m^2/s]}"
-plot "" u ($2/1000):8 w l lt 6 t ""
+#set ylabel "{/:Normal Specific sediment flux [m^2/s]}"
+set ylabel "{/:Normal Specific sediment flux [m^2/a]}"
+plot "" u ($2/1000):($8*to_year) w l lt 6 t ""