Wellcome to the new version of BEST methods

The BEST-version-003 is functionning similarly to the BEST-version-002.
When you excute the main.sce script, the choice of "BEST-WR" (Di Prima et al., 2021)
and "BEST-WR-3T" (Yilmaz et al., 2022) methods are available for Beerkan infiltration data treatment.
"Use this version if your data is influenced by water repellency (WR)" 

Note: This version of BEST WR methods is adapted for manually recorded data. However for the automated recorded
beerkan data and specifically for the choice of BEST-WR methods, the user have to select the steady-state points.
Because for the optimization procedure of S (sorptivity) and Ks, BEST-WR method consider  the points before 
the assumed steady state.
For BEST-WR-3T method all the data (included steady state points) set can be use for the optimization procedure.

The data set need to organized with two column as follow:

1.27	0.00       // dry bulk density kg/dm^3 (M/L^3) 
2.65	0.00	   // specific density kg/dm^3 (M/L^3) 
0.17	0.00	   // volumetric initial water content
52.00	0.00	   // Beerkan ring radius (L) in the example [mm]
0.002	0.2        // particle fraction below 0.002 mm and % of fraction<0.002 mm/100
0.05	0.35	   // particle fraction below 0.050 mm and % of fraction<0.050 mm/100
2.00	1.00	   // the complete fraction [0-2 mm] is equal to 100 %/100 =1 
0.00	0.00       // set 0 and 0 to setup the beginning of the cumulated infiltration data
563.00	5.89	   //     t1 (T)  and Icum1 (L)    in the example t(s) and Icum (mm)
1017.00	11.78      //     t2      and Icum2
1498.00	17.67		   .		.
1919.00	23.56		   .		.
2147.00	29.44		   .		.
2374.00	35.33		   .		.
2607.00	41.22		   .		.
2856.00	47.11		   .		.
3107.00	53.00		   .		.
3357.00	58.89	   //     tfinal  and Icum_final
0.00	0.00       //  set 0 and 0 to endup the experimental data