diff --git a/src/starformation/schaye08/starformation.h b/src/starformation/schaye08/starformation.h
index 02bf188a85edcaa4affc9e99400a968a75ec4827..13f2b71cfac4a5e675b353f555a529978430db45 100644
--- a/src/starformation/schaye08/starformation.h
+++ b/src/starformation/schaye08/starformation.h
@@ -23,53 +23,95 @@
 /* Starformation struct */
 struct star_formation {
   
-  /* Normalization of the KS star formation law */
+  /*! Normalization of the KS star formation law */
   double A;
 
-  /* Slope of the KS law */
+  /*! Slope of the KS law */
   double nks;
 
-  /* Critical density */
-  double rho_crit;
+  /*! Critical density */
+  double Delta_crit;
 
-  /* Critical temperature */
+  /*! Critical temperature */
   double T_crit;
 
-  /* Ratio of the specific heats */
+  /*! Ratio of the specific heats */
   double gamma;
 
-  /* gas fraction */
+  /*! gas fraction */
   double fg;
 
-  /* Star formation law slope */
+  /*! Star formation law slope */
   double nstar;
 
-  /* star formation normalization of schaye+08 */
+  /*! star formation normalization of schaye+08 */
   double Astar;
   
 };
 
 
-/* Brief initialization of the star formation law 
- * 
+/* 
+ * Brief initialization of the star formation law 
+ *
+ * @param parameter_file The parsed parameter file
+ * @param phys_const Physical constants in internal units
+ * @param us The current internal system of units
+ * @param starform the star formation law properties to initialize
+ *
  * */
 static void starformation_init_backend(
-    
-    ) {
+  struct swift_params* parameter_file, const struct phys_const* phys_const,
+  const struct unit_system* us, const struct space* s,
+  const struct star_formation* starform) {
+  
   /* Default values for the normalization and the power law */
   static const double normalization_default = 2.5e-4;
-  static const double KS_power_law = 1.4;
+  static const double KS_power_law_default = 1.4;
+
+  /* Default value for the heat capacity ratio gamma */
+  static const double gamma_default = 5.f/3.f;
+
+  /* Read the critical density contrast from the parameter file*/
+  starform->Delta_crit = parser_get_param_double(parameter_file, 
+  "SchayeSF:Delta_crit");
+
+  /* Read the critical temperature from the parameter file */
+  starform->T_crit = parser_get_param_double(parameter_file,
+  "SchayeSF:T_crit");
+
+  /* Read the gas fraction from the file */
+  starform->fg = parser_get_param_double(parameter_file,
+  "SchayeSF:fg");
+
+  /* Read the normalization */
+  const double normalization = parser_get_opt_param_double(
+  parameter_file, "SchayeSF:A", normalization_default);
+
+  /* Read the Kennicutt-Schmidt power law exponent */
+  starform->nks = parser_get_opt_param_double(
+  parameter_file, "SchayeSF:nks", KS_power_law_default);
+
+  /* Read the heat capacity ratio gamma */
+  starform->gamma = parser_get_opt_param_double(
+  parameter_file, "SchayeSF:gamma", gamma_default); 
+
+  /* Calculate the power law of the star formation */
+  starform->nstar = (starform->nks - 1.f)/2.f;
 
+  /* Get the appropriate constant to calculate the 
+   * star formation constant */ 
+  const double G_newton = phys_const->const_newton_G;
 
   
 }
 
 /* @brief Prints the used parameters of the star formation law 
  *
+ * @param starform the star formation law properties.
  * */
 static void starformation_print_backend(
+    const struct star_formation* starform){ 
 
-    ){ 
   message("Star formation law is Schaye and Dalla Vecchia (2008)"
   " with properties, normalization = %e, slope of the Kennicutt"
   "-Schmidt law = %e, gamma = %e, gas fraction = %e, critical "