move prep table makkink to dev/bln_tables

Author: gerardhros

dev/bln_tables.R

@@ -49,10 +49,12 @@
                                    crop_name = B_LU_NAME,
                                    crop_cat1 = B_LU_CULTCAT4,
                                    bln_country ="NL",
+                                   B_LU = B_LU,
                                    B_LU_EOM,
                                    B_LU_EOM_RESIDUE,
                                    B_LU_HC,
-                                   B_LU_WATERSTRESS_OBIC
+                                   B_LU_WATERSTRESS_OBIC,
+                                   B_LU_MAKKINK
                                    )]
 
   # replace löss with loess
@@ -67,6 +69,23 @@
   # save updated crop table
   usethis::use_data(bln_crops,overwrite = TRUE)
 
+# make makkink table
+
+  # load csv file with makkink per crop per month
+  dt <- fread('dev/bln_makkink.csv',header = TRUE)
+
+  # colnames for the 12 months
+  cols <- colnames(dt)[-1]
+
+  # replace all missing values with the value 0.36 for braak soil
+  dt[,c(cols) := lapply(.SD,function(x) fifelse(is.na(x),0.36,x)),.SDcols = cols]
+
+  # write file to data
+  bln_makkink <- copy(dt)
+
+  # write file to data
+  usethis::use_data(bln_makkink, overwrite = TRUE)
+
 # make LSW table with the averaged properties across whole NL
 
   # set LSW to averaged values across whole NL