From cb08390a07bb8200d6e4f7bf5a06bcc2b7db0bd3 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Roland=20Holl=C3=B3s?= Date: Tue, 29 Sep 2020 15:37:59 +0200 Subject: [PATCH] debug soilGrid function --- RBBGCMuso/R/soilQuery.R | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/RBBGCMuso/R/soilQuery.R b/RBBGCMuso/R/soilQuery.R index 2da0615..a81a0a1 100644 --- a/RBBGCMuso/R/soilQuery.R +++ b/RBBGCMuso/R/soilQuery.R @@ -50,13 +50,13 @@ createSoilFile <- function(lat,lon, } soilDepth <- tryCatch(getMeanSoil(rest,"bdod")/100,error=function(e){stop("There is no data for the given coordinates")}) - outFile[47] <- sprintf("%s (m) soil depth",paste(soilDepth,collapse="\t")) outFile[55] <- sprintf("%s (%%) percentage of sand by volume in rock free soil", paste(createMusoLayers(getMeanSoil(rest,"sand")/10), collapse="\t")) outFile[56] <- sprintf("%s (%%) percentage of silt by volume in rock free soil", paste(createMusoLayers(getMeanSoil(rest,"silt")/10), collapse="\t")) outFile[57] <- sprintf("%s (dimless) soil PH", paste(createMusoLayers(getMeanSoil(rest,"phh2o")/10), collapse="\t")) + outFile[58] <- sprintf("%s (%%) bulk density",paste(createMusoLayers(soilDepth),collapse="\t")) writeLines(outFile,outputFile) } # createSoilFile(60,50)