diff --git a/examples/parameter_example.yml b/examples/parameter_example.yml
index 845cfb74fe6d9f412eddb28c6f398702ac8881cd..e9ec53dc044d1267acf88790180406e68dd58b4e 100644
--- a/examples/parameter_example.yml
+++ b/examples/parameter_example.yml
@@ -292,6 +292,10 @@ SchayeSF:
   A:              2.5e-4  # (optional) The normalization of the Kennicutt-Schmidt law
   nks:            1.4     # (optional) The power law of the Kennicutt-Schmidt law
   gamma:          1.6667  # (optional) The heat capacity ratio (gamma)
+  Schaye2004:     1       # (optional) whether to use the metallicity dependent critical star formation of Schaye (2004) (1) or not (0).
+  norm_ncrit:     .1      # (optional) Critical sf normalization to use (is not a normalization when Schaye2004=0, than it is the value.
+  Z0:             0.002   # (optional) Scale metallicity to use for the equation (not used for Schaye2004=0)
+  n_Z0:           -0.64   # (optional) Scaling of the critical density with the metallicity (not used for Schaye2004=0)
 
 # Structure finding options (requires velociraptor)
 StructureFinding:
diff --git a/src/starformation/schaye08/starformation.h b/src/starformation/schaye08/starformation.h
index 4c22c3b51b9e48aa7d427f581f40b33d9d353be7..d5c3a2fb7a27eb8cc87856809f4adbe8169fe309 100644
--- a/src/starformation/schaye08/starformation.h
+++ b/src/starformation/schaye08/starformation.h
@@ -268,11 +268,15 @@ static void starformation_init_backend(
     staform->n_Z0 = parser_get_opt_param_double(
     parameter_file, "SchayeSF:n_Z0", powerlawZ_default);
   }
+  /* Conversion of number density from cgs */
+  static const float dimension_numb_den = {0, -3, 0, 0, 0};
+  static const double conversion_numb_density = 1/
+  units_general_cgs_conversion_factor(us, dimension_numb_den);
 
   /* Calculate the prefactor that is always common */
   /* !!!DONT FORGET TO DO THE CORRECT UNIT CONVERSION!!!*/
   starform->den_crit_star = starform->den_crit / pow(starform->Z0,
-  starform->n_Z0);
+  starform->n_Z0) * conversion_numb_density;