diff --git a/src/DAMASK_marc.f90 b/src/DAMASK_marc.f90 index 7789fc58c..233e684ee 100644 --- a/src/DAMASK_marc.f90 +++ b/src/DAMASK_marc.f90 @@ -255,7 +255,7 @@ subroutine hypela2(d,g,e,de,s,t,dt,ngens,m,nn,kcus,matus,ndi,nshear,disp, & outdatedByNewInc = .false., & !< needs description CPFEM_init_done = .false. !< remember whether init has been done already class(tNode), pointer :: & - debug_Marc + debug_Marc ! pointer to Marc debug options defaultNumThreadsInt = omp_get_num_threads() ! remember number of threads set by Marc call omp_set_num_threads(1) ! no openMP diff --git a/src/constitutive.f90 b/src/constitutive.f90 index 41baab7ab..dc15644d5 100644 --- a/src/constitutive.f90 +++ b/src/constitutive.f90 @@ -37,37 +37,37 @@ module constitutive module subroutine plastic_none_init(debug_constitutive) class(tNode), pointer , intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options end subroutine plastic_none_init module subroutine plastic_isotropic_init(debug_constitutive) class(tNode), pointer , intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options end subroutine plastic_isotropic_init module subroutine plastic_phenopowerlaw_init(debug_constitutive) class(tNode), pointer , intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options end subroutine plastic_phenopowerlaw_init module subroutine plastic_kinehardening_init(debug_constitutive) class(tNode), pointer , intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options end subroutine plastic_kinehardening_init module subroutine plastic_dislotwin_init(debug_constitutive) class(tNode), pointer , intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options end subroutine plastic_dislotwin_init module subroutine plastic_disloUCLA_init(debug_constitutive) class(tNode), pointer , intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options end subroutine plastic_disloUCLA_init module subroutine plastic_nonlocal_init(debug_constitutive) class(tNode), pointer , intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options end subroutine plastic_nonlocal_init @@ -84,7 +84,7 @@ module constitutive instance, & of class(tNode), pointer , intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options end subroutine plastic_isotropic_LpAndItsTangent pure module subroutine plastic_phenopowerlaw_LpAndItsTangent(Lp,dLp_dMp,Mp,instance,of) @@ -175,7 +175,7 @@ module constitutive instance, & of class(tNode), pointer , intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options end subroutine plastic_isotropic_LiAndItsTangent @@ -239,7 +239,7 @@ module constitutive ip, & !< current integration point el !< current element number class(tNode), pointer , intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options end subroutine plastic_nonlocal_dotState @@ -268,7 +268,7 @@ module constitutive ip, & el class(tNode), pointer , intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options end subroutine plastic_nonlocal_dependentState @@ -279,7 +279,7 @@ module constitutive instance, & of class(tNode), pointer , intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options end subroutine plastic_kinehardening_deltaState module subroutine plastic_nonlocal_deltaState(Mp,instance,of,ip,el,debug_constitutive) @@ -291,7 +291,7 @@ module constitutive ip, & el class(tNode), pointer , intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options end subroutine plastic_nonlocal_deltaState diff --git a/src/constitutive_plastic_disloUCLA.f90 b/src/constitutive_plastic_disloUCLA.f90 index 3a2946435..b355ca551 100644 --- a/src/constitutive_plastic_disloUCLA.f90 +++ b/src/constitutive_plastic_disloUCLA.f90 @@ -77,7 +77,7 @@ contains module subroutine plastic_disloUCLA_init(debug_constitutive) class(tNode), pointer, intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options integer :: & Ninstance, & diff --git a/src/constitutive_plastic_dislotwin.f90 b/src/constitutive_plastic_dislotwin.f90 index c643376c7..4cedaca94 100644 --- a/src/constitutive_plastic_dislotwin.f90 +++ b/src/constitutive_plastic_dislotwin.f90 @@ -125,7 +125,7 @@ contains module subroutine plastic_dislotwin_init(debug_constitutive) class(tNode), pointer, intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options integer :: & Ninstance, & diff --git a/src/constitutive_plastic_isotropic.f90 b/src/constitutive_plastic_isotropic.f90 index 7e5736d31..adb293c64 100644 --- a/src/constitutive_plastic_isotropic.f90 +++ b/src/constitutive_plastic_isotropic.f90 @@ -52,7 +52,7 @@ contains module subroutine plastic_isotropic_init(debug_constitutive) class(tNode), pointer, intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options integer :: & Ninstance, & @@ -173,7 +173,7 @@ module subroutine plastic_isotropic_LpAndItsTangent(Lp,dLp_dMp,Mp,instance,of,de instance, & of class(tNode), pointer, intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options real(pReal), dimension(3,3) :: & Mp_dev !< deviatoric part of the Mandel stress @@ -236,7 +236,7 @@ module subroutine plastic_isotropic_LiAndItsTangent(Li,dLi_dMi,Mi,instance,of,de instance, & of class(tNode), pointer, intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options real(pReal) :: & tr !< trace of spherical part of Mandel stress (= 3 x pressure) diff --git a/src/constitutive_plastic_kinehardening.f90 b/src/constitutive_plastic_kinehardening.f90 index 6ff250057..645a3916b 100644 --- a/src/constitutive_plastic_kinehardening.f90 +++ b/src/constitutive_plastic_kinehardening.f90 @@ -61,7 +61,7 @@ contains module subroutine plastic_kinehardening_init(debug_constitutive) class(tNode), pointer, intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options integer :: & Ninstance, & @@ -324,7 +324,7 @@ module subroutine plastic_kinehardening_deltaState(Mp,instance,of,debug_constitu instance, & of class(tNode), pointer , intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options real(pReal), dimension(param(instance)%sum_N_sl) :: & gdot_pos,gdot_neg, & diff --git a/src/constitutive_plastic_none.f90 b/src/constitutive_plastic_none.f90 index d79b08d93..3b7a4ff28 100644 --- a/src/constitutive_plastic_none.f90 +++ b/src/constitutive_plastic_none.f90 @@ -15,7 +15,7 @@ contains module subroutine plastic_none_init(debug_constitutive) class(tNode), pointer, intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options integer :: & Ninstance, & diff --git a/src/constitutive_plastic_nonlocal.f90 b/src/constitutive_plastic_nonlocal.f90 index 596c62c86..a131f4517 100644 --- a/src/constitutive_plastic_nonlocal.f90 +++ b/src/constitutive_plastic_nonlocal.f90 @@ -166,7 +166,7 @@ contains module subroutine plastic_nonlocal_init(debug_constitutive) class(tNode), pointer, intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options integer :: & Ninstance, & @@ -536,7 +536,7 @@ module subroutine plastic_nonlocal_dependentState(F, Fp, instance, of, ip, el,de ip, & el class(tNode), pointer, intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options integer :: & no, & !< neighbor offset @@ -854,7 +854,7 @@ module subroutine plastic_nonlocal_deltaState(Mp,instance,of,ip,el,debug_constit ip, & el class(tNode), pointer, intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options integer :: & ph, & !< phase @@ -973,7 +973,7 @@ module subroutine plastic_nonlocal_dotState(Mp, F, Fp, Temperature,timestep, & ip, & !< current integration point el !< current element number class(tNode), pointer, intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options integer :: & ph, & @@ -1177,7 +1177,7 @@ function rhoDotFlux(F,Fp,timestep, instance,of,ip,el,debug_constitutive) ip, & !< current integration point el !< current element number class(tNode), pointer, intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options integer :: & ph, & diff --git a/src/constitutive_plastic_phenopowerlaw.f90 b/src/constitutive_plastic_phenopowerlaw.f90 index 76e2606f8..1d3bb9961 100644 --- a/src/constitutive_plastic_phenopowerlaw.f90 +++ b/src/constitutive_plastic_phenopowerlaw.f90 @@ -69,7 +69,7 @@ contains module subroutine plastic_phenopowerlaw_init(debug_constitutive) class(tNode), pointer, intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options integer :: & Ninstance, & diff --git a/src/crystallite.f90 b/src/crystallite.f90 index 23e44c70d..47450d7ee 100644 --- a/src/crystallite.f90 +++ b/src/crystallite.f90 @@ -131,7 +131,7 @@ subroutine crystallite_init class(tNode) , pointer :: & num_crystallite, & - debug_crystallite + debug_crystallite ! pointer to debug options for crystallite write(6,'(/,a)') ' <<<+- crystallite init -+>>>' diff --git a/src/debug.f90 b/src/debug.f90 index 32fc8433a..deb01ec01 100644 --- a/src/debug.f90 +++ b/src/debug.f90 @@ -14,8 +14,8 @@ module debug implicit none private - class(tNode), pointer, public :: & - debug_root + class(tNode), pointer, protected, public :: & + debug_root !< root pointer storing the debug YAML structure public :: debug_init diff --git a/src/grid/DAMASK_grid.f90 b/src/grid/DAMASK_grid.f90 index d10c4badc..ebc1b6643 100644 --- a/src/grid/DAMASK_grid.f90 +++ b/src/grid/DAMASK_grid.f90 @@ -94,7 +94,7 @@ program DAMASK_grid class (tNode), pointer :: & num_grid, & num_generic, & - debug_grid + debug_grid ! pointer to grid debug options !-------------------------------------------------------------------------------------------------- ! init DAMASK (all modules) diff --git a/src/grid/grid_mech_FEM.f90 b/src/grid/grid_mech_FEM.f90 index 8a8791151..a543f679c 100644 --- a/src/grid/grid_mech_FEM.f90 +++ b/src/grid/grid_mech_FEM.f90 @@ -506,7 +506,7 @@ subroutine formResidual(da_local,x_local, & itmax class(tNode), pointer :: & num_generic, & - debug_grid + debug_grid ! pointer to grid debug options !--------------------------------------------------------------------- ! debug pointer to grid diff --git a/src/grid/grid_mech_spectral_basic.f90 b/src/grid/grid_mech_spectral_basic.f90 index 0ae081e37..533448837 100644 --- a/src/grid/grid_mech_spectral_basic.f90 +++ b/src/grid/grid_mech_spectral_basic.f90 @@ -475,7 +475,7 @@ subroutine formResidual(in, F, & itmax class(tNode), pointer :: & num_generic, & - debug_grid + debug_grid ! pointer to constitutive debug options !--------------------------------------------------------------------- ! debug pointer to grid diff --git a/src/grid/grid_mech_spectral_polarisation.f90 b/src/grid/grid_mech_spectral_polarisation.f90 index 543a9da9c..0cc221e7c 100644 --- a/src/grid/grid_mech_spectral_polarisation.f90 +++ b/src/grid/grid_mech_spectral_polarisation.f90 @@ -535,7 +535,7 @@ subroutine formResidual(in, FandF_tau, & class(tNode), pointer :: & num_grid, & num_generic, & - debug_grid + debug_grid ! pointer to grid debug options real(pReal) :: & polarAlpha, & !< polarization scheme parameter 0.0 < alpha < 2.0. alpha = 1.0 ==> AL scheme, alpha = 2.0 ==> accelerated scheme polarBeta !< polarization scheme parameter 0.0 < beta < 2.0. beta = 1.0 ==> AL scheme, beta = 2.0 ==> accelerated scheme diff --git a/src/grid/spectral_utilities.f90 b/src/grid/spectral_utilities.f90 index 273f0dabf..e0d730343 100644 --- a/src/grid/spectral_utilities.f90 +++ b/src/grid/spectral_utilities.f90 @@ -195,7 +195,7 @@ subroutine spectral_utilities_init class (tNode) , pointer :: & num_grid, & num_generic, & - debug_grid + debug_grid ! pointer to grid debug options write(6,'(/,a)') ' <<<+- spectral_utilities init -+>>>' diff --git a/src/homogenization.f90 b/src/homogenization.f90 index 4729933fe..38ff4ffd7 100644 --- a/src/homogenization.f90 +++ b/src/homogenization.f90 @@ -56,18 +56,18 @@ module homogenization module subroutine mech_none_init(debug_homogenization) class(tNode), pointer, intent(in) :: & - debug_homogenization + debug_homogenization !< pointer to debug options for homogenization end subroutine mech_none_init module subroutine mech_isostrain_init(debug_homogenization) class(tNode), pointer, intent(in) :: & - debug_homogenization + debug_homogenization !< pointer to debug options for homogenization end subroutine mech_isostrain_init module subroutine mech_RGC_init(num_homogMech, debug_homogenization) class(tNode), pointer, intent(in) :: & num_homogMech, & - debug_homogenization + debug_homogenization !< pointer to debug options for homogenization end subroutine mech_RGC_init @@ -84,7 +84,7 @@ module homogenization instance, & of class(tNode), pointer, intent(in) :: & - debug_homogenization + debug_homogenization !< pointer to debug options for homogenization end subroutine mech_RGC_partitionDeformation @@ -112,14 +112,14 @@ module homogenization P,& !< partitioned stresses F,& !< partitioned deformation gradients F0 !< partitioned initial deformation gradients - real(pReal), dimension(:,:,:,:,:), intent(in) :: dPdF !< partitioned stiffnesses - real(pReal), dimension(3,3), intent(in) :: avgF !< average F - real(pReal), intent(in) :: dt !< time increment + real(pReal), dimension(:,:,:,:,:), intent(in) :: dPdF !< partitioned stiffnesses + real(pReal), dimension(3,3), intent(in) :: avgF !< average F + real(pReal), intent(in) :: dt !< time increment integer, intent(in) :: & - ip, & !< integration point number - el !< element number + ip, & !< integration point number + el !< element number class(tNode), pointer, intent(in) :: & - debug_homogenization + debug_homogenization !< pointer to debug options for homogenization end function mech_RGC_updateState @@ -152,8 +152,8 @@ subroutine homogenization_init debug_g, & debug_e - num_homog => numerics_root%get('homogenization',defaultVal=emptyDict) - num_homogMech => num_homog%get('mech',defaultVal=emptyDict) + num_homog => numerics_root%get('homogenization',defaultVal=emptyDict) + num_homogMech => num_homog%get('mech',defaultVal=emptyDict) num_homogGeneric => num_homog%get('generic',defaultVal=emptyDict) debug_homogenization => debug_root%get('homogenization',defaultVal=emptyList) diff --git a/src/homogenization_mech_RGC.f90 b/src/homogenization_mech_RGC.f90 index 2a0b1800c..0cada9240 100644 --- a/src/homogenization_mech_RGC.f90 +++ b/src/homogenization_mech_RGC.f90 @@ -79,8 +79,7 @@ module subroutine mech_RGC_init(num_homogMech,debug_homogenization) class(tNode), pointer, intent(in) :: & num_homogMech, & - debug_homogenization - + debug_homogenization !< pointer to debug options for homogenization integer :: & Ninstance, & h, & @@ -214,8 +213,7 @@ module subroutine mech_RGC_partitionDeformation(F,avgF,instance,of,debug_homogen instance, & of class(tNode), pointer, intent(in) :: & - debug_homogenization - + debug_homogenization !< pointer to debug options for homogenization real(pReal), dimension(3) :: aVect,nVect integer, dimension(4) :: intFace integer, dimension(3) :: iGrain3 @@ -668,7 +666,7 @@ module procedure mech_RGC_updateState real(pReal), dimension (:,:,:), intent(in) :: fDef !< deformation gradients real(pReal), dimension (3,3), intent(in) :: avgF !< initial effective stretch tensor integer, intent(in) :: ip,el,instance,of - class(tNode), pointer, intent(in) :: debug_homogenization + class(tNode), pointer, intent(in) :: debug_homogenization !< pointer to debug options for homogenization integer, dimension (4) :: intFace integer, dimension (3) :: iGrain3,iGNghb3,nGDim @@ -783,7 +781,7 @@ module procedure mech_RGC_updateState Ngrain, & instance, & of - class(tNode), pointer, intent(in) :: debug_homogenization + class(tNode), pointer, intent(in) :: debug_homogenization !< pointer to debug options for homogenization real(pReal), dimension(size(vPen,3)) :: gVol integer :: i diff --git a/src/homogenization_mech_isostrain.f90 b/src/homogenization_mech_isostrain.f90 index b41550cd6..46a4c1215 100644 --- a/src/homogenization_mech_isostrain.f90 +++ b/src/homogenization_mech_isostrain.f90 @@ -29,7 +29,7 @@ contains module subroutine mech_isostrain_init(debug_homogenization) class(tNode), pointer, intent(in) :: & - debug_homogenization + debug_homogenization !< pointer to debug options for homogenization integer :: & Ninstance, & diff --git a/src/homogenization_mech_none.f90 b/src/homogenization_mech_none.f90 index 247c9e1f7..01e2568b3 100644 --- a/src/homogenization_mech_none.f90 +++ b/src/homogenization_mech_none.f90 @@ -14,7 +14,7 @@ contains module subroutine mech_none_init(debug_homogenization) class(tNode), pointer, intent(in) :: & - debug_homogenization + debug_homogenization !< pointer to debug options for homogenization integer :: & Ninstance, & diff --git a/src/kinematics_cleavage_opening.f90 b/src/kinematics_cleavage_opening.f90 index 2ffb27f2e..c35c9b615 100644 --- a/src/kinematics_cleavage_opening.f90 +++ b/src/kinematics_cleavage_opening.f90 @@ -46,7 +46,7 @@ contains subroutine kinematics_cleavage_opening_init(debug_constitutive) class(tNode), pointer, intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options integer :: Ninstance,p integer, dimension(:), allocatable :: N_cl !< active number of cleavage systems per family diff --git a/src/kinematics_slipplane_opening.f90 b/src/kinematics_slipplane_opening.f90 index 278754b79..caeff8e56 100644 --- a/src/kinematics_slipplane_opening.f90 +++ b/src/kinematics_slipplane_opening.f90 @@ -48,7 +48,7 @@ contains subroutine kinematics_slipplane_opening_init(debug_constitutive) class(tNode), pointer , intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options integer :: Ninstance,p,i character(len=pStringLen) :: extmsg = '' diff --git a/src/kinematics_thermal_expansion.f90 b/src/kinematics_thermal_expansion.f90 index 1a84fb7b8..1d099451c 100644 --- a/src/kinematics_thermal_expansion.f90 +++ b/src/kinematics_thermal_expansion.f90 @@ -41,7 +41,7 @@ contains subroutine kinematics_thermal_expansion_init(debug_constitutive) class(tNode), pointer, intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options integer :: Ninstance,p,i real(pReal), dimension(:), allocatable :: temp diff --git a/src/material.f90 b/src/material.f90 index 3915d705b..7de8214cd 100644 --- a/src/material.f90 +++ b/src/material.f90 @@ -217,7 +217,7 @@ subroutine material_init(restart) CounterPhase, & CounterHomogenization class(tNode), pointer :: & - debug_material + debug_material ! pointer to material debug options write(6,'(/,a)') ' <<<+- material init -+>>>'; flush(6) diff --git a/src/mesh/FEM_utilities.f90 b/src/mesh/FEM_utilities.f90 index 94b69c073..d91236598 100644 --- a/src/mesh/FEM_utilities.f90 +++ b/src/mesh/FEM_utilities.f90 @@ -105,7 +105,7 @@ subroutine FEM_utilities_init class(tNode), pointer :: & num_mesh, & num_generic, & - debug_mesh + debug_mesh ! pointer to mesh debug options integer :: structOrder !< order of displacement shape functions character(len=pStringLen) :: & petsc_options, & diff --git a/src/source_damage_anisoBrittle.f90 b/src/source_damage_anisoBrittle.f90 index 5826e7160..7422fa113 100644 --- a/src/source_damage_anisoBrittle.f90 +++ b/src/source_damage_anisoBrittle.f90 @@ -56,7 +56,7 @@ contains subroutine source_damage_anisoBrittle_init(debug_constitutive) class(tNode), pointer , intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options integer :: Ninstance,sourceOffset,NipcMyPhase,p integer, dimension(:), allocatable :: N_cl diff --git a/src/source_damage_anisoDuctile.f90 b/src/source_damage_anisoDuctile.f90 index 73e68d021..d511e9dc3 100644 --- a/src/source_damage_anisoDuctile.f90 +++ b/src/source_damage_anisoDuctile.f90 @@ -49,7 +49,7 @@ contains subroutine source_damage_anisoDuctile_init(debug_constitutive) class(tNode), pointer, intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options integer :: Ninstance,sourceOffset,NipcMyPhase,p integer, dimension(:), allocatable :: N_sl diff --git a/src/source_damage_isoBrittle.f90 b/src/source_damage_isoBrittle.f90 index e1c0625e4..7abd1e494 100644 --- a/src/source_damage_isoBrittle.f90 +++ b/src/source_damage_isoBrittle.f90 @@ -48,7 +48,7 @@ contains subroutine source_damage_isoBrittle_init(debug_constitutive) class(tNode), pointer, intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options integer :: Ninstance,sourceOffset,NipcMyPhase,p character(len=pStringLen) :: extmsg = '' diff --git a/src/source_damage_isoDuctile.f90 b/src/source_damage_isoDuctile.f90 index 83be82300..4dd77bed5 100644 --- a/src/source_damage_isoDuctile.f90 +++ b/src/source_damage_isoDuctile.f90 @@ -48,7 +48,7 @@ contains subroutine source_damage_isoDuctile_init(debug_constitutive) class(tNode), pointer, intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options integer :: Ninstance,sourceOffset,NipcMyPhase,p character(len=pStringLen) :: extmsg = '' diff --git a/src/source_thermal_dissipation.f90 b/src/source_thermal_dissipation.f90 index e8dd2eddb..099dd1e80 100644 --- a/src/source_thermal_dissipation.f90 +++ b/src/source_thermal_dissipation.f90 @@ -41,7 +41,7 @@ contains subroutine source_thermal_dissipation_init(debug_constitutive) class(tNode), pointer, intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options integer :: Ninstance,sourceOffset,NipcMyPhase,p diff --git a/src/source_thermal_externalheat.f90 b/src/source_thermal_externalheat.f90 index 65b11a499..023b81d22 100644 --- a/src/source_thermal_externalheat.f90 +++ b/src/source_thermal_externalheat.f90 @@ -45,7 +45,7 @@ contains subroutine source_thermal_externalheat_init(debug_constitutive) class(tNode), pointer, intent(in) :: & - debug_constitutive + debug_constitutive !< pointer to constitutive debug options integer :: Ninstance,sourceOffset,NipcMyPhase,p