From 8fb8e9be6ea2a3fff9c2d037df9a11f831943d86 Mon Sep 17 00:00:00 2001
From: Martin Diehl <mail@martin-diehl.net>
Date: Tue, 10 Dec 2019 20:25:19 +0100
Subject: [PATCH] write temperature to DADF5

---
 src/homogenization.f90     |  7 +++++++
 src/thermal_adiabatic.f90  | 26 ++++++++++++++++++++++++++
 src/thermal_conduction.f90 | 27 +++++++++++++++++++++++++++
 3 files changed, 60 insertions(+)

diff --git a/src/homogenization.f90 b/src/homogenization.f90
index 5fa723eda..b014bebb5 100644
--- a/src/homogenization.f90
+++ b/src/homogenization.f90
@@ -779,6 +779,13 @@ subroutine homogenization_results
     
     group = trim(group_base)//'/thermal'
     call results_closeGroup(results_addGroup(group))
+    select case(thermal_type(p))
+      case(THERMAL_ADIABATIC_ID)
+        call thermal_adiabatic_results(p,group)
+      case(THERMAL_CONDUCTION_ID)
+        call thermal_conduction_results(p,group)
+    end select
+    
     
  enddo
 
diff --git a/src/thermal_adiabatic.f90 b/src/thermal_adiabatic.f90
index 2aa69bec5..5c30f280b 100644
--- a/src/thermal_adiabatic.f90
+++ b/src/thermal_adiabatic.f90
@@ -7,6 +7,7 @@ module thermal_adiabatic
   use config
   use numerics
   use material
+  use results
   use source_thermal_dissipation
   use source_thermal_externalheat
   use crystallite
@@ -37,6 +38,7 @@ module thermal_adiabatic
     thermal_adiabatic_getSourceAndItsTangent, &
     thermal_adiabatic_getSpecificHeat, &
     thermal_adiabatic_getMassDensity, &
+    thermal_adiabatic_results, &
     thermal_adiabatic_postResults
  
 contains
@@ -251,6 +253,30 @@ function thermal_adiabatic_getMassDensity(ip,el)
 end function thermal_adiabatic_getMassDensity
 
 
+!--------------------------------------------------------------------------------------------------
+!> @brief writes results to HDF5 output file
+!--------------------------------------------------------------------------------------------------
+subroutine thermal_adiabatic_results(homog,group)
+
+  integer,          intent(in) :: homog
+  character(len=*), intent(in) :: group
+#if defined(PETSc) || defined(DAMASK_HDF5)  
+  integer :: o, instance
+  
+  instance  = thermal_typeInstance(homog)
+
+  outputsLoop: do o = 1,thermal_adiabatic_Noutput(instance)
+     select case(thermal_adiabatic_outputID(o,instance))
+    
+      case (temperature_ID)
+        call results_writeDataset(group,temperature(homog)%p,'T',&
+                                  'temperature','K')
+    end select
+  enddo outputsLoop
+#endif
+
+end subroutine thermal_adiabatic_results
+
 !--------------------------------------------------------------------------------------------------
 !> @brief return array of thermal results
 !--------------------------------------------------------------------------------------------------
diff --git a/src/thermal_conduction.f90 b/src/thermal_conduction.f90
index e513d709f..5507d8a7d 100644
--- a/src/thermal_conduction.f90
+++ b/src/thermal_conduction.f90
@@ -7,6 +7,7 @@ module thermal_conduction
   use material
   use config
   use lattice
+  use results
   use crystallite
   use source_thermal_dissipation
   use source_thermal_externalheat
@@ -37,6 +38,7 @@ module thermal_conduction
     thermal_conduction_getSpecificHeat, &
     thermal_conduction_getMassDensity, &
     thermal_conduction_putTemperatureAndItsRate, &
+    thermal_conduction_results, &
     thermal_conduction_postResults
 
 contains
@@ -263,6 +265,31 @@ subroutine thermal_conduction_putTemperatureAndItsRate(T,Tdot,ip,el)
 
 end subroutine thermal_conduction_putTemperatureAndItsRate
  
+
+!--------------------------------------------------------------------------------------------------
+!> @brief writes results to HDF5 output file
+!--------------------------------------------------------------------------------------------------
+subroutine thermal_conduction_results(homog,group)
+
+  integer,          intent(in) :: homog
+  character(len=*), intent(in) :: group
+#if defined(PETSc) || defined(DAMASK_HDF5)  
+  integer :: o, instance
+  
+  instance  = thermal_typeInstance(homog)
+
+  outputsLoop: do o = 1,thermal_conduction_Noutput(instance)
+     select case(thermal_conduction_outputID(o,instance))
+    
+      case (temperature_ID)
+        call results_writeDataset(group,temperature(homog)%p,'T',&
+                                  'temperature','K')
+    end select
+  enddo outputsLoop
+#endif
+
+end subroutine thermal_conduction_results
+
  
 !--------------------------------------------------------------------------------------------------
 !> @brief return array of thermal results