diff --git a/src/cooling/grackle/cooling_io.h b/src/cooling/grackle/cooling_io.h
index 83f0bcbac711fe186f670ed91b961c957bce0a75..3905cafd05fb8e15ddf33f4ea688d6144698df73 100644
--- a/src/cooling/grackle/cooling_io.h
+++ b/src/cooling/grackle/cooling_io.h
@@ -61,8 +61,6 @@ __attribute__((always_inline)) INLINE static int cooling_write_particles(
 
   int num = 0;
 
-  if (cooling->output_mode == 0) return num;
-
 #if COOLING_GRACKLE_MODE >= 1
   /* List what we want to write */
   list[0] = io_make_output_field("HI", FLOAT, 1, UNIT_CONV_NO_UNITS, xparts,
@@ -86,8 +84,6 @@ __attribute__((always_inline)) INLINE static int cooling_write_particles(
   num += 6;
 #endif
 
-  if (cooling->output_mode == 1) return num;
-
 #if COOLING_GRACKLE_MODE >= 2
   list += num;
 
@@ -103,8 +99,6 @@ __attribute__((always_inline)) INLINE static int cooling_write_particles(
   num += 3;
 #endif
 
-  if (cooling->output_mode == 2) return num;
-
 #if COOLING_GRACKLE_MODE >= 3
   list += num;
 
@@ -153,9 +147,6 @@ __attribute__((always_inline)) INLINE static void cooling_read_parameters(
   cooling->self_shielding_method = parser_get_opt_param_int(
       parameter_file, "GrackleCooling:SelfShieldingMethod", 0);
 
-  cooling->output_mode =
-      parser_get_opt_param_int(parameter_file, "GrackleCooling:OutputMode", 0);
-
   cooling->max_step = parser_get_opt_param_int(
       parameter_file, "GrackleCooling:MaxSteps", 10000);
 
diff --git a/src/cooling/grackle/cooling_struct.h b/src/cooling/grackle/cooling_struct.h
index 2f1b187761a49a4e394d2b3d2a0f759fc388dbbb..66c385234cccf6532100e86f2c16a508b1112baa 100644
--- a/src/cooling/grackle/cooling_struct.h
+++ b/src/cooling/grackle/cooling_struct.h
@@ -59,9 +59,6 @@ struct cooling_function_data {
   /* Self shielding method (<= 3) means grackle modes */
   int self_shielding_method;
 
-  /* Output mode (correspond to primordial chemistry mode */
-  int output_mode;
-
   /* convergence limit for first init */
   float convergence_limit;