diff --git a/src/homogenization_mech_RGC.f90 b/src/homogenization_mech_RGC.f90
index 57b573506..5003f7c24 100644
--- a/src/homogenization_mech_RGC.f90
+++ b/src/homogenization_mech_RGC.f90
@@ -1125,7 +1125,16 @@ subroutine mech_RGC_results(instance,group)
 
   outputsLoop: do o = 1,size(prm%outputID)
     select case(prm%outputID(o))
-
+    
+      case (constitutivework_ID)
+        call results_writeDataset(group,stt%work,'W',&
+                                  'work density','J/m³')
+      case (magnitudemismatch_ID)
+        call results_writeDataset(group,dst%mismatch,'M',&
+                                  'average mismatch tensor','tbd')
+      case (penaltyenergy_ID)
+        call results_writeDataset(group,stt%penaltyEnergy,'R',&
+                                  'mismatch penalty density','tbd')
 
     end select
   enddo outputsLoop