From af92d4c3e3a120ff8ca7081c3fa937fc61c7a9bb Mon Sep 17 00:00:00 2001
From: lhausamm <loic_hausammann@hotmail.com>
Date: Fri, 3 Nov 2017 17:03:21 +0100
Subject: [PATCH] Code works and seems to give good results

---
 src/cooling/grackle/grackle_wrapper.c | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/cooling/grackle/grackle_wrapper.c b/src/cooling/grackle/grackle_wrapper.c
index c88efc918a..34647ea932 100644
--- a/src/cooling/grackle/grackle_wrapper.c
+++ b/src/cooling/grackle/grackle_wrapper.c
@@ -134,6 +134,8 @@ int wrap_get_cooling_time(double rho, double u, double Z, double a_now,
 
 int wrap_do_cooling(double rho, double *u, double dt, double Z, double a_now) {
 
+  GRACKLE_ASSERT(FIELD_SIZE == 1);
+
   gr_float den_factor = 1.0;
   gr_float u_factor = 1.0;
 
@@ -149,8 +151,6 @@ int wrap_do_cooling(double rho, double *u, double dt, double Z, double a_now) {
   int grid_start[3] = {0, 0, 0};
   int grid_end[3] = {0, 0, 0};
 
-  GRACKLE_ASSERT(FIELD_SIZE == 1);
-
   if (solve_chemistry_table(&my_units, a_now, dt, grid_rank, grid_dimension,
                             grid_start, grid_end, density, energy, x_velocity,
                             y_velocity, z_velocity, metal_density) == 0) {
-- 
GitLab