diff --git a/extern/SoilMoistureProfiles/SoilMoistureProfiles b/extern/SoilMoistureProfiles/SoilMoistureProfiles
index a02a463fee..832d5a70d6 160000
--- a/extern/SoilMoistureProfiles/SoilMoistureProfiles
+++ b/extern/SoilMoistureProfiles/SoilMoistureProfiles
@@ -1 +1 @@
-Subproject commit a02a463feec7a3feea4ad528bb354041b09fbeec
+Subproject commit 832d5a70d632d9978dec40a480e970a6b59720c5
diff --git a/include/core/DomainLayer.hpp b/include/core/DomainLayer.hpp
index 29d0b098e8..c2aa6d5c8f 100644
--- a/include/core/DomainLayer.hpp
+++ b/include/core/DomainLayer.hpp
@@ -46,6 +46,20 @@ namespace ngen
             formulation->write_output("Time Step,""Time,"+formulation->get_output_header_line(",")+"\n");
         }
 
+        /**
+         * @brief Get a list of output variables names for this layer
+         * 
+         * @return vector of output variable names
+        */
+
+        virtual std::vector<std::string>&& get_output_variable_names() 
+        { 
+            std::vector<std::string> output_variables;
+
+
+            return std::move(output_variables);
+        }
+
         /***
          * @brief Run one simulation timestep for this model associated with the domain
          * 
diff --git a/include/core/Layer.hpp b/include/core/Layer.hpp
index 0610bc79f7..19e068d46e 100644
--- a/include/core/Layer.hpp
+++ b/include/core/Layer.hpp
@@ -87,6 +87,17 @@ namespace ngen
         */
         virtual int class_id() { return LayerClass::kLayer; }
 
+        /**
+         * @brief Get a list of output variables names for this layer
+         * 
+         * @return vector of output variable names
+        */
+
+        virtual std::vector<std::string>&& get_output_variable_names() 
+        { 
+            return std::move(std::vector<std::string>());
+        }
+
         /***
          * @brief Return the next timestep that will be processed by this layer in epoch time units
         */