diff --git a/DESCRIPTION b/DESCRIPTION
index 90b5e63dfa08b122ea259ab99f03e4b8c469effa..3295a7a25087c56b5854db43588eb7818005a2e5 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -1,7 +1,7 @@
 Package: SWATdoctR
 Type: Package
 Title: Finding the right diagnoses and treatments for SWAT+ models
-Version: 0.1.20
+Version: 0.1.21
 Author: c(person("Christoph", "Schürz",
              email = "christoph.schuerz@ufz.de",
              role = c("aut", "cre")),
diff --git a/R/update_landuse_labels.R b/R/update_landuse_labels.R
index 011de7dfc7fb93c73f1f1df0ec3d06cf44f7aa31..26e32427d4898e7f1536097a1404484adaae5891 100644
--- a/R/update_landuse_labels.R
+++ b/R/update_landuse_labels.R
@@ -153,8 +153,13 @@ remove_consonants <- function(tbl) {
     n <- tbl$n_chr[1] - 12
     txt_split <- str_split(tbl$lu_mgt_upd[1], '', simplify = TRUE)
     alph_pos <- which(str_detect(txt_split, '[:alpha:]'))
-    pos_rmv <- sample(alph_pos[alph_pos > 1], n)
-    txt <- paste(txt_split[-pos_rmv], collapse = '')
+    pos_rmv <- seq(2, length(alph_pos), 2)
+    pos_rmv <- pos_rmv[1:min(length(pos_rmv), n)]
+    if (length(pos_rmv) < n) {
+      n_diff <- n - length(pos_rmv)
+      pos_rmv <- c(pos_rmv, seq(3, length(alph_pos), length.out = n_diff))
+    }
+    txt <- paste(txt_split[-alph_pos[pos_rmv]], collapse = '')
     tbl$lu_mgt_upd <- txt
   }
   return(tbl)