diff --git a/docs/source/_static/parameter-files/parameter_LEVEL2.prm b/docs/source/_static/parameter-files/parameter_LEVEL2.prm
index 7e671f7a..03c37688 100644
--- a/docs/source/_static/parameter-files/parameter_LEVEL2.prm
+++ b/docs/source/_static/parameter-files/parameter_LEVEL2.prm
@@ -156,6 +156,10 @@ MULTI_SCATTERING = TRUE
 # daily tables do not exist or if there is no valid daily value.
 # Type: full directory path
 DIR_WVPLUT = NULL
+# This indicates if the water vapor climatology tables are permitted (FALSE),
+# or if only daily tables are allowed (TRUE)
+# Type: Logical. Valid values: {TRUE,FALSE}
+STRICT_WATER_VAPOR = FALSE
 # This specifies a global value for atmospheric water vapor content in
 # g cm-2. This parameter can be a dummy value to quickly process an image
 # without needing to generate a water vapor database. Note that especially
diff --git a/docs/source/_static/parameter-files/parameter_LEVEL3.prm b/docs/source/_static/parameter-files/parameter_LEVEL3.prm
index 979abd8b..e20b7c7c 100644
--- a/docs/source/_static/parameter-files/parameter_LEVEL3.prm
+++ b/docs/source/_static/parameter-files/parameter_LEVEL3.prm
@@ -222,6 +222,13 @@ Y_FACTOR = 0.75
 # Should selection or weighting be performed?
 # Type: Logical. Valid values: {TRUE,FALSE}
 SELECT = TRUE
+# How to combine individual scores?
+# In the default behaviour, the individual scores are summed up.
+# You can also choose to multiply the scores. This applies a 
+# strong penalty on an observation that performs poor in a single score.
+# The weigting of scores is applied in both cases.
+# Type: Character. Valid values: {ADDITIVE,MULTIPLICATIVE}
+COMBINE_SCORES = ADDITIVE
 # These parameters specify the function values used for fitting the DOY
 # scoring functions. The function type is automatically chosen from the 
 # given values, i.e.
@@ -251,6 +258,16 @@ DOY_STATIC = 120 180 240
 # parts of the year are selected.
 # Type: Logical. Valid values: {TRUE,FALSE}
 OFF_SEASON = FALSE
+# This parameter specifies whether observations with a cloud score of less than 1%
+# should be candidates. On one hand, this reduces nodata gaps in the composite. On the other hand,
+# it will include clouds if there is no better observation.
+# Type: Logical. Valid values: {TRUE,FALSE}
+USE_CLOUDY = FALSE
+# This parameter specifies whether observations with a haze score of less than 1%
+# should be candidates. On one hand, this reduces nodata gaps in the composite. On the other hand,
+# it will include haze if there is no better observation.
+# Type: Logical. Valid values: {TRUE,FALSE}
+USE_HAZY = FALSE
 # This parameter controls the strength of the DOY score.
 # 0 disables the use of this score.
 # Type: Float. Valid values: [0,1]
diff --git a/docs/source/_static/parameter-files/parameter_TSA.prm b/docs/source/_static/parameter-files/parameter_TSA.prm
index 46f00bdd..8428660e 100644
--- a/docs/source/_static/parameter-files/parameter_TSA.prm
+++ b/docs/source/_static/parameter-files/parameter_TSA.prm
@@ -256,7 +256,7 @@ OUTPUT_RMS = FALSE
 # Radial Basis Function or harmonic Interpolation.
 # Harmonic interpolation can be used as a simple near-real time monitoring component.
 # Type: Character. Valid values: {NONE,LINEAR,MOVING,RBF,HARMONIC}
-INTERPOLATE = RBF
+INTERPOLATE = NONE
 # Max temporal distance for the moving average filter in days. For each
 # interpolation date, MOVING_MAX days before and after are considered.
 # Type: Integer. Valid range: [1,365]
diff --git a/misc/force-version.txt b/misc/force-version.txt
index a436e30c..5878221f 100644
--- a/misc/force-version.txt
+++ b/misc/force-version.txt
@@ -1 +1 @@
-3.7.12-dev:::2024-05-23_13:00:29
+3.7.12-dev:::2024-07-11_06:19:54
diff --git a/src/aux-level/param-aux.c b/src/aux-level/param-aux.c
index d610333e..2a1306d3 100755
--- a/src/aux-level/param-aux.c
+++ b/src/aux-level/param-aux.c
@@ -1499,7 +1499,7 @@ void write_par_hl_tsi(FILE *fp, bool verbose){
     fprintf(fp, "# Harmonic interpolation can be used as a simple near-real time monitoring component.\n");
     fprintf(fp, "# Type: Character. Valid values: {NONE,LINEAR,MOVING,RBF,HARMONIC}\n");
   }
-  fprintf(fp, "INTERPOLATE = RBF\n");
+  fprintf(fp, "INTERPOLATE = NONE\n");
 
   if (verbose){
     fprintf(fp, "# Max temporal distance for the moving average filter in days. For each\n");