diff --git a/tardis/model/parse_input.py b/tardis/model/parse_input.py
index 745def52e1b..b8c3ef1eb81 100644
--- a/tardis/model/parse_input.py
+++ b/tardis/model/parse_input.py
@@ -678,6 +678,13 @@ def parse_csvy_radiation_field_state(
             geometry, packet_source
         )
 
+    if np.any(t_radiative < 1000 * u.K):
+        logging.critical(
+            "Radiative temperature is too low in some of the shells, temperatures below 1000K "
+            f"(e.g., T_rad = {t_radiative[np.argmin(t_radiative)]} in shell {np.argmin(t_radiative)} in your model) "
+            "are not accurately handled by TARDIS.",
+        )
+
     if hasattr(csvy_model_data, "columns") and (
         "dilution_factor" in csvy_model_data.columns
     ):