From 3b673d23dc1ecb8a191592bc272b25c0df38ea5d Mon Sep 17 00:00:00 2001
From: Martin Diehl <m.diehl@mpie.de>
Date: Tue, 3 May 2016 17:06:55 +0200
Subject: [PATCH] added code for solving heat equation

---
 code/DAMASK_marc.f90 | 27 +++++++++++++++++++++++++++
 1 file changed, 27 insertions(+)

diff --git a/code/DAMASK_marc.f90 b/code/DAMASK_marc.f90
index a4542f96a..bef54489e 100644
--- a/code/DAMASK_marc.f90
+++ b/code/DAMASK_marc.f90
@@ -384,6 +384,33 @@ subroutine hypela2(d,g,e,de,s,t,dt,ngens,m,nn,kcus,matus,ndi,nshear,disp, &
 end subroutine hypela2
 
 
+!--------------------------------------------------------------------------------------------------
+!> @brief thermomechanical coupling
+!--------------------------------------------------------------------------------------------------
+subroutine flux(f,ts,n,time)
+ use prec, only: &
+   pReal, &
+   pInt
+ use homogenization, only: &
+   materialpoint_heat
+ use mesh, only: &
+   mesh_FEasCP
+ 
+ implicit none
+ real(pReal),   dimension(6),           intent(in) :: &
+   ts
+ integer(pInt), dimension(10),          intent(in) :: &
+   n
+ real(pReal),                           intent(in) :: &
+   time
+ real(pReal),   dimension(2),           intent(out) :: &
+   f
+ 
+ f(1) = materialpoint_heat(n(3),mesh_FEasCP('elem',n(1)))
+ 
+end subroutine flux
+
+
 !--------------------------------------------------------------------------------------------------
 !> @brief sets user defined output variables for Marc
 !> @details select a variable contour plotting (user subroutine).