Merge branch 'MarcThermoMech' into development

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 calculate internal heat generated due to inelastic energy dissipation
+!--------------------------------------------------------------------------------------------------
+subroutine flux(f,ts,n,time)
+ use prec, only: &
+   pReal, &
+   pInt
+ use thermal_conduction, only: &
+   thermal_conduction_getSourceAndItsTangent
+ 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
+
+ call thermal_conduction_getSourceAndItsTangent(f(1), f(2), ts(3), 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).