diff --git a/src/phase.f90 b/src/phase.f90 index 42e495483..c34455f2a 100644 --- a/src/phase.f90 +++ b/src/phase.f90 @@ -394,6 +394,7 @@ subroutine phase_init do ph = 1,phases%length phase => phases%get_dict(ph) + print'(a,i0,a)', ' phase ',ph,' '//material_references(phase) phase_lattice(ph) = phase%get_asString('lattice') if (all(phase_lattice(ph) /= ['cF','cI','hP','tI'])) & call IO_error(130,ext_msg='phase_init: '//phase%get_asString('lattice')) @@ -516,7 +517,7 @@ end subroutine phase_result !-------------------------------------------------------------------------------------------------- -!> @brief allocates and initialize per grain variables +!> @brief Allocate and initialize. !-------------------------------------------------------------------------------------------------- subroutine crystallite_init() diff --git a/src/phase_damage.f90 b/src/phase_damage.f90 index b21df97e9..7c5b31973 100644 --- a/src/phase_damage.f90 +++ b/src/phase_damage.f90 @@ -103,6 +103,7 @@ module subroutine damage_init() phase => phases%get_dict(ph) source => phase%get_dict('damage',defaultVal=emptyDict) if (source%length > 0) then + print'(a,i0,a)', ' phase ',ph,' '//material_references(source) damage_active = .true. param(ph)%mu = source%get_asFloat('mu') param(ph)%l_c = source%get_asFloat('l_c') diff --git a/src/phase_damage_anisobrittle.f90 b/src/phase_damage_anisobrittle.f90 index 6988ac39a..7fe11a71e 100644 --- a/src/phase_damage_anisobrittle.f90 +++ b/src/phase_damage_anisobrittle.f90 @@ -62,6 +62,8 @@ module function anisobrittle_init() result(mySources) associate(prm => param(ph)) + print'(a,i0,a)', ' phase ',ph,' '//material_references(src) + N_cl = src%get_as1dInt('N_cl',defaultVal=emptyIntArray) prm%sum_N_cl = sum(abs(N_cl)) diff --git a/src/phase_damage_isobrittle.f90 b/src/phase_damage_isobrittle.f90 index 562ffe6c2..0b1641d38 100644 --- a/src/phase_damage_isobrittle.f90 +++ b/src/phase_damage_isobrittle.f90 @@ -64,6 +64,8 @@ module function isobrittle_init() result(mySources) prm%W_crit = src%get_asFloat('G_crit')/src%get_asFloat('l_c') + print'(a,i0,a)', ' phase ',ph,' '//material_references(src) + #if defined (__GFORTRAN__) prm%output = output_as1dString(src) #else diff --git a/src/phase_mechanical_elastic.f90 b/src/phase_mechanical_elastic.f90 index d5ad9f916..91cb7319d 100644 --- a/src/phase_mechanical_elastic.f90 +++ b/src/phase_mechanical_elastic.f90 @@ -42,6 +42,7 @@ module subroutine elastic_init(phases) phase => phases%get_dict(ph) mech => phase%get_dict('mechanical') elastic => mech%get_dict('elastic') + print'(a,i0,a)', ' phase ',ph,' '//material_references(elastic) if (elastic%get_asString('type') /= 'Hooke') call IO_error(200,ext_msg=elastic%get_asString('type')) associate(prm => param(ph)) diff --git a/src/phase_mechanical_plastic_dislotungsten.f90 b/src/phase_mechanical_plastic_dislotungsten.f90 index 40c1743bb..dedd876e2 100644 --- a/src/phase_mechanical_plastic_dislotungsten.f90 +++ b/src/phase_mechanical_plastic_dislotungsten.f90 @@ -128,6 +128,8 @@ module function plastic_dislotungsten_init() result(myPlasticity) mech => phase%get_dict('mechanical') pl => mech%get_dict('plastic') + print'(a,i0,a)', ' phase ',ph,' '//material_references(pl) + #if defined (__GFORTRAN__) prm%output = output_as1dString(pl) #else diff --git a/src/phase_mechanical_plastic_dislotwin.f90 b/src/phase_mechanical_plastic_dislotwin.f90 index 52864ddf6..b1c97a9e1 100644 --- a/src/phase_mechanical_plastic_dislotwin.f90 +++ b/src/phase_mechanical_plastic_dislotwin.f90 @@ -181,6 +181,8 @@ module function plastic_dislotwin_init() result(myPlasticity) mech => phase%get_dict('mechanical') pl => mech%get_dict('plastic') + print'(a,i0,a)', ' phase ',ph,' '//material_references(pl) + #if defined (__GFORTRAN__) prm%output = output_as1dString(pl) #else diff --git a/src/phase_mechanical_plastic_isotropic.f90 b/src/phase_mechanical_plastic_isotropic.f90 index 64e3060ab..ae3fa2e80 100644 --- a/src/phase_mechanical_plastic_isotropic.f90 +++ b/src/phase_mechanical_plastic_isotropic.f90 @@ -86,6 +86,8 @@ module function plastic_isotropic_init() result(myPlasticity) mech => phase%get_dict('mechanical') pl => mech%get_dict('plastic') + print'(a,i0,a)', ' phase ',ph,' '//material_references(pl) + #if defined (__GFORTRAN__) prm%output = output_as1dString(pl) #else diff --git a/src/phase_mechanical_plastic_kinehardening.f90 b/src/phase_mechanical_plastic_kinehardening.f90 index 622de124a..8aa824a5f 100644 --- a/src/phase_mechanical_plastic_kinehardening.f90 +++ b/src/phase_mechanical_plastic_kinehardening.f90 @@ -84,6 +84,7 @@ module function plastic_kinehardening_init() result(myPlasticity) mech, & pl + myPlasticity = plastic_active('kinehardening') if (count(myPlasticity) == 0) return @@ -111,6 +112,8 @@ module function plastic_kinehardening_init() result(myPlasticity) mech => phase%get_dict('mechanical') pl => mech%get_dict('plastic') + print'(a,i0,a)', ' phase ',ph,' '//material_references(pl) + #if defined (__GFORTRAN__) prm%output = output_as1dString(pl) #else diff --git a/src/phase_mechanical_plastic_none.f90 b/src/phase_mechanical_plastic_none.f90 index 351d4244b..b79a61183 100644 --- a/src/phase_mechanical_plastic_none.f90 +++ b/src/phase_mechanical_plastic_none.f90 @@ -31,6 +31,7 @@ module function plastic_none_init() result(myPlasticity) phases => config_material%get_dict('phase') do ph = 1, phases%length if (.not. myPlasticity(ph)) cycle + print'(a,i0,a)', ' phase ',ph call phase_allocateState(plasticState(ph),count(material_ID_phase == ph),0,0,0) end do diff --git a/src/phase_mechanical_plastic_nonlocal.f90 b/src/phase_mechanical_plastic_nonlocal.f90 index baaecb7b9..f157b2339 100644 --- a/src/phase_mechanical_plastic_nonlocal.f90 +++ b/src/phase_mechanical_plastic_nonlocal.f90 @@ -234,13 +234,15 @@ module function plastic_nonlocal_init() result(myPlasticity) mech => phase%get_dict('mechanical') pl => mech%get_dict('plastic') - plasticState(ph)%nonlocal = pl%get_asBool('flux',defaultVal=.True.) + print'(a,i0,a)', ' phase ',ph,' '//material_references(pl) + #if defined (__GFORTRAN__) prm%output = output_as1dString(pl) #else prm%output = pl%get_as1dString('output',defaultVal=emptyStringArray) #endif + plasticState(ph)%nonlocal = pl%get_asBool('flux',defaultVal=.True.) prm%isotropic_bound = pl%get_asString('isotropic_bound',defaultVal='isostrain') prm%atol_rho = pl%get_asFloat('atol_rho',defaultVal=1.0_pReal) diff --git a/src/phase_mechanical_plastic_phenopowerlaw.f90 b/src/phase_mechanical_plastic_phenopowerlaw.f90 index c60e58411..0f8b9a299 100644 --- a/src/phase_mechanical_plastic_phenopowerlaw.f90 +++ b/src/phase_mechanical_plastic_phenopowerlaw.f90 @@ -122,6 +122,14 @@ module function plastic_phenopowerlaw_init() result(myPlasticity) mech => phase%get_dict('mechanical') pl => mech%get_dict('plastic') + print'(a,i0,a)', ' phase ',ph,' '//material_references(pl) + +#if defined (__GFORTRAN__) + prm%output = output_as1dString(pl) +#else + prm%output = pl%get_as1dString('output',defaultVal=emptyStringArray) +#endif + !-------------------------------------------------------------------------------------------------- ! slip related parameters N_sl = pl%get_as1dInt('N_sl',defaultVal=emptyIntArray) @@ -217,15 +225,6 @@ module function plastic_phenopowerlaw_init() result(myPlasticity) prm%h_0_tw_sl = 0.0_pReal end if slipAndTwinActive -!-------------------------------------------------------------------------------------------------- -! output pararameters - -#if defined (__GFORTRAN__) - prm%output = output_as1dString(pl) -#else - prm%output = pl%get_as1dString('output',defaultVal=emptyStringArray) -#endif - !-------------------------------------------------------------------------------------------------- ! allocate state arrays Nmembers = count(material_ID_phase == ph) diff --git a/src/phase_thermal.f90 b/src/phase_thermal.f90 index 1e44258ef..dc80aa423 100644 --- a/src/phase_thermal.f90 +++ b/src/phase_thermal.f90 @@ -107,6 +107,7 @@ module subroutine thermal_init(phases) ! ToDo: temperature dependency of K and C_p if (thermal%length > 0) then + print'(a,i0,a)', ' phase ',ph,' '//material_references(thermal) param(ph)%C_p = thermal%get_asFloat('C_p') param(ph)%K(1,1) = thermal%get_asFloat('K_11') if (any(phase_lattice(ph) == ['hP','tI'])) param(ph)%K(3,3) = thermal%get_asFloat('K_33') diff --git a/src/phase_thermal_dissipation.f90 b/src/phase_thermal_dissipation.f90 index ba2a9aafc..ffb48f80c 100644 --- a/src/phase_thermal_dissipation.f90 +++ b/src/phase_thermal_dissipation.f90 @@ -56,6 +56,7 @@ module function dissipation_init(source_length) result(mySources) if (mySources(so,ph)) then associate(prm => param(ph)) src => sources%get_dict(so) + print'(a,i0,a,i0,a)', ' phase ',ph,' source ',so,' '//material_references(src) prm%kappa = src%get_asFloat('kappa') Nmembers = count(material_ID_phase == ph) diff --git a/src/phase_thermal_externalheat.f90 b/src/phase_thermal_externalheat.f90 index 0e1d11b62..421117f9e 100644 --- a/src/phase_thermal_externalheat.f90 +++ b/src/phase_thermal_externalheat.f90 @@ -60,6 +60,7 @@ module function externalheat_init(source_length) result(mySources) source_thermal_externalheat_offset(ph) = so associate(prm => param(ph)) src => sources%get_dict(so) + print'(a,i0,a,i0,a)', ' phase ',ph,' source ',so,' '//material_references(src) prm%f = table(src,'t','f')