Use as_apsim_met() to import SILO data and add comments back

R/get_data_drill_apsim.R

@@ -61,6 +61,8 @@
 #'
 #' @section Saving objects:
 #' To save \dQuote{met} objects, please use [apsimx::write_apsim_met()].
+#'   Note that when saving, comments from SILO will be included, but these will
+#'   not be printed as a part of the resulting `met` object in your \R session.
 #'
 #' @return An \CRANpkg{apsimx} object of class \sQuote{met} with attributes.
 #'

R/get_patched_point_apsim.R

@@ -72,6 +72,8 @@
 #'
 #' @section Saving objects:
 #' To save \dQuote{met} objects, please use [apsimx::write_apsim_met()].
+#'   Note that when saving, comments from SILO will be included, but these will
+#'   not be printed as a part of the resulting `met` object in your \R session.
 #'
 #' @return An \CRANpkg{apsimx} object of class \sQuote{met} with attributes.
 #'

R/query_silo_api.R

@@ -123,7 +123,35 @@
     met_file_path <- file.path(tempdir(), "apsim.met")
     on.exit(unlink(met_file_path))
     writeLines(text = response$parse("UTF8"), con = met_file_path)
-    return(apsimx::read_apsim_met(file = "apsim.met", src.dir = tempdir()))
+    apsim <- readLines(met_file_path)
+    data <- data.table::fread(file = met_file_path,
+                              skip = grep(pattern = "^year",
+                                          x = apsim) + 2)
+    apsim_met <- apsimx::as_apsim_met(
+      x = data,
+      filename = sprintf("station.%s.met",
+                         trimws(sub("^.*=", "",
+                                    apsim[grep(pattern = "^!station number =",
+                                               x = apsim)]))),
+      site = trimws(sub("^.*=", "",
+                        apsim[grep(pattern = "^!station name =", x = apsim)])),
+      latitude = trimws(sub("^.*=", "",
+                            apsim[grep(pattern = "^latitude =", x = apsim)])),
+      longitude = trimws(sub("^.*=", "",
+                             apsim[grep(pattern = "^longitude =", x = apsim)])),
+      colnames = scan(
+        text = apsim[grep(pattern = "year", x = apsim)],
+        what = '',
+        quiet = TRUE
+      ),
+      units = scan(
+        text = apsim[grep(pattern = "\\(\\)", x = apsim)],
+        what = '',
+        quiet = TRUE
+      ),
+      comments = apsim[grep(pattern = "^!", x = apsim)]
+    )
+    return(apsim_met)
   }
 
   response_data <- data.table::fread(response$parse("UTF8"))