suggested changes from MR !530

examples/config/homogenization/mechanical/RGC_8grains.yaml

@@ -1,4 +1,4 @@
-# For Relaxed Grain Cluster homogenization, requires N_constituents = 8
+# For Relaxed Grain Cluster homogenization, requires N_constituents = prod(cluster_size)
 type: RGC
 D_alpha: [4.0e-06, 4.0e-06, 2.0e-06]
 a_g: [0.0, 0.0, 0.0]

src/homogenization_mechanical.f90

@@ -204,7 +204,7 @@ module subroutine mechanical_results(group_base,ho)
                                   'deformation gradient','1')
       case('P')
         call results_writeDataset(reshape(homogenization_P,[3,3,discretization_nCells]),group,'P', &
-                                  'deformation gradient','1')
+                                  'first Piola-Kirchhoff stress','Pa')
     end select
   end do
 

src/homogenization_thermal_pass.f90

@@ -11,7 +11,7 @@ module subroutine pass_init()
   print'(/,1x,a)', '<<<+-  homogenization:thermal:pass init  -+>>>'
 
   if (homogenization_Nconstituents(1) /= 1) &
-    call IO_error(211,ext_msg='N_constituents (pass)')
+    call IO_error(211,ext_msg='(pass) with N_constituents !=1')
 
 end subroutine pass_init
 

src/phase_thermal.f90

@@ -388,7 +388,7 @@ end function thermal_active
 
 
 !----------------------------------------------------------------------------------------------
-!< @brief writes damage sources results to HDF5 output file
+!< @brief writes thermal sources results to HDF5 output file
 !----------------------------------------------------------------------------------------------
 module subroutine thermal_results(group,ph)
 
@@ -398,18 +398,16 @@ module subroutine thermal_results(group,ph)
 
   integer :: ou
 
-  if (allocated(param(ph)%output)) then
+  if (.not. allocated(param(ph)%output)) return
+
   call results_closeGroup(results_addGroup(group//'thermal'))
-  else
-    return
-  endif
 
   do ou = 1, size(param(ph)%output)
 
     select case(trim(param(ph)%output(ou)))
 
       case ('T')
-        call results_writeDataset(current(ph)%T,group//'thermal','T', 'temperature','T')
+        call results_writeDataset(current(ph)%T,group//'thermal','T', 'temperature','K')
 
     end select