diff --git a/cara/mc-output-publication.py b/cara/mc-output-publication.py
index a12bc6aa..f65fd2cd 100644
--- a/cara/mc-output-publication.py
+++ b/cara/mc-output-publication.py
@@ -5,5 +5,26 @@ Code version: <code version>
 Availability: <where it's located> """
 
 from cara.montecarlo import *
-from cara.model_scenarios import *
+from cara.model_scenarios_publication import *
+
+# qR table:
+#TODO table with qR values
+
+# qR values for Breathing, light activity:
+print('\n<<<<<<<<<<< qR for Breathing, light activity >>>>>>>>>>>')
+present_qR_values(qR_models[0].concentration_model)
+# qR values for Speaking, light activity:
+print('\n<<<<<<<<<<< qR for Speaking, light activity >>>>>>>>>>>')
+present_qR_values(qR_models[1].concentration_model)
+# qR values for Shouting, light activity:
+print('\n<<<<<<<<<<< qR for Shouting / Singing, light activity >>>>>>>>>>>')
+present_qR_values(qR_models[2].concentration_model)
+
+# qR values for shouting and light activity - qID=100, 500 and 1000
+print('\n<<<<<<<<<<< shouting and light activity - qID=100 >>>>>>>>>>>')
+present_qR_values(qR_models_shout_light[0].concentration_model)
+print('\n<<<<<<<<<<< shouting and light activity - qID=500 >>>>>>>>>>>')
+present_qR_values(qR_models_shout_light[1].concentration_model)
+print('\n<<<<<<<<<<< shouting and light activity - qID=1000 >>>>>>>>>>>')
+present_qR_values(qR_models_shout_light[2].concentration_model)