diff --git a/src/cooling/grackle/cooling.h b/src/cooling/grackle/cooling.h index 07741c97997caddc24c63cd4d68244a41fb4d200..43d20d918a4679314964735a6b51e0e5ffe36380 100644 --- a/src/cooling/grackle/cooling.h +++ b/src/cooling/grackle/cooling.h @@ -491,10 +491,12 @@ __attribute__((always_inline)) INLINE static gr_float cooling_rate( const struct cosmology* restrict cosmo, const struct cooling_function_data* restrict cooling, const struct part* restrict p, struct xpart* restrict xp, double dt) { - - error( - "Check cosmology factors (physical vs. co-moving and drifted vs. " - "un-drifted)!"); + + if (cosmo->Omega_m != 0. || cosmo->Omega_r != 0. || cosmo->Omega_k != 0. || + cosmo->Omega_lambda != 0. || cosmo->Omega_b != 0.) + error( + "Check cosmology factors (physical vs. co-moving and drifted vs. " + "un-drifted)!"); /* set current time */ code_units units = cooling->units; @@ -646,9 +648,11 @@ __attribute__((always_inline)) INLINE static void cooling_cool_part( struct part* restrict p, struct xpart* restrict xp, double dt, double dt_therm) { - error( - "Check cosmology factors (physical vs. co-moving and drifted vs. " - "un-drifted)!"); + if (cosmo->Omega_m != 0. || cosmo->Omega_r != 0. || cosmo->Omega_k != 0. || + cosmo->Omega_lambda != 0. || cosmo->Omega_b != 0.) + error( + "Check cosmology factors (physical vs. co-moving and drifted vs. " + "un-drifted)!"); if (dt == 0.) return;