Update initialize_solver.inl

Split up InitializeVariables() functions
NCAR · Aug 4, 2024 · ec26a2c · ec26a2c
1 parent 0f82e21
commit ec26a2c
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/include/tuvx/radiative_transfer/solvers/initialize_solver.inl b/include/tuvx/radiative_transfer/solvers/initialize_solver.inl
@@ -30,6 +30,12 @@ namespace tuvx
     // Check for consistency between the grids and profiles.
     assert(vertical_grid.NumberOfColumns() == number_of_columns);
     assert(wavelength_grid.NumberOfColumns() == 1);
+
+    // Scale variables for the Delta-Eddington approximation
+    ScaleVariables(solver_variables);
+
+    // Generate functions that define the source variables (C functions from the paper)
+    BuildSourceFunctions(solver_variables, source_functions); 
   }
 
   template<typename T>