From 37a7364a3eb427c5d1258dc8c8c73df8bc75cde1 Mon Sep 17 00:00:00 2001 From: Pratheek Shanthraj Date: Wed, 25 Mar 2015 16:02:30 +0000 Subject: [PATCH] multi physics output now working for all solvers --- code/constitutive.f90 | 317 +++++++++++++++++++--------------------- code/crystallite.f90 | 22 ++- code/homogenization.f90 | 92 +++++------- 3 files changed, 203 insertions(+), 228 deletions(-) diff --git a/code/constitutive.f90 b/code/constitutive.f90 index 08cfe7902..518b704cb 100644 --- a/code/constitutive.f90 +++ b/code/constitutive.f90 @@ -242,164 +242,164 @@ subroutine constitutive_init(temperature_init) !-------------------------------------------------------------------------------------------------- ! write description file for constitutive phase output - call IO_write_jobFile(FILEUNIT,'outputConstitutive') - do phase = 1_pInt,material_Nphase - instance = phase_plasticityInstance(phase) ! which instance of a plasticity is present phase - knownPlasticity = .true. ! assume valid - select case(phase_plasticity(phase)) ! split per constititution - case (PLASTICITY_NONE_ID) - outputName = PLASTICITY_NONE_label - thisNoutput => null() - thisOutput => null() ! plastic_none_output - thisSize => null() ! plastic_none_sizePostResult - case (PLASTICITY_J2_ID) - outputName = PLASTICITY_J2_label - thisNoutput => plastic_j2_Noutput - thisOutput => plastic_j2_output - thisSize => plastic_j2_sizePostResult - case (PLASTICITY_PHENOPOWERLAW_ID) - outputName = PLASTICITY_PHENOPOWERLAW_label - thisNoutput => plastic_phenopowerlaw_Noutput - thisOutput => plastic_phenopowerlaw_output - thisSize => plastic_phenopowerlaw_sizePostResult - case (PLASTICITY_DISLOTWIN_ID) - outputName = PLASTICITY_DISLOTWIN_label - thisNoutput => plastic_dislotwin_Noutput - thisOutput => plastic_dislotwin_output - thisSize => plastic_dislotwin_sizePostResult - case (PLASTICITY_DISLOKMC_ID) - outputName = PLASTICITY_DISLOKMC_label - thisNoutput => plastic_dislokmc_Noutput - thisOutput => plastic_dislokmc_output - thisSize => plastic_dislokmc_sizePostResult - case (PLASTICITY_DISLOUCLA_ID) - outputName = PLASTICITY_DISLOUCLA_label - thisNoutput => plastic_disloucla_Noutput - thisOutput => plastic_disloucla_output - thisSize => plastic_disloucla_sizePostResult - case (PLASTICITY_TITANMOD_ID) - outputName = PLASTICITY_TITANMOD_label - thisNoutput => plastic_titanmod_Noutput - thisOutput => plastic_titanmod_output - thisSize => plastic_titanmod_sizePostResult - case (PLASTICITY_NONLOCAL_ID) - outputName = PLASTICITY_NONLOCAL_label - thisNoutput => plastic_nonlocal_Noutput - thisOutput => plastic_nonlocal_output - thisSize => plastic_nonlocal_sizePostResult - case default - knownPlasticity = .false. - end select - write(FILEUNIT,'(/,a,/)') '['//trim(phase_name(phase))//']' - if (knownPlasticity) then - write(FILEUNIT,'(a)') '(plasticity)'//char(9)//trim(outputName) - if (phase_plasticity(phase) /= PLASTICITY_NONE_ID) then - do e = 1_pInt,thisNoutput(instance) - write(FILEUNIT,'(a,i4)') trim(thisOutput(e,instance))//char(9),thisSize(e,instance) - enddo + if (worldrank == 0_pInt) then + call IO_write_jobFile(FILEUNIT,'outputConstitutive') + do phase = 1_pInt,material_Nphase + instance = phase_plasticityInstance(phase) ! which instance of a plasticity is present phase + knownPlasticity = .true. ! assume valid + select case(phase_plasticity(phase)) ! split per constititution + case (PLASTICITY_NONE_ID) + outputName = PLASTICITY_NONE_label + thisNoutput => null() + thisOutput => null() ! plastic_none_output + thisSize => null() ! plastic_none_sizePostResult + case (PLASTICITY_J2_ID) + outputName = PLASTICITY_J2_label + thisNoutput => plastic_j2_Noutput + thisOutput => plastic_j2_output + thisSize => plastic_j2_sizePostResult + case (PLASTICITY_PHENOPOWERLAW_ID) + outputName = PLASTICITY_PHENOPOWERLAW_label + thisNoutput => plastic_phenopowerlaw_Noutput + thisOutput => plastic_phenopowerlaw_output + thisSize => plastic_phenopowerlaw_sizePostResult + case (PLASTICITY_DISLOTWIN_ID) + outputName = PLASTICITY_DISLOTWIN_label + thisNoutput => plastic_dislotwin_Noutput + thisOutput => plastic_dislotwin_output + thisSize => plastic_dislotwin_sizePostResult + case (PLASTICITY_DISLOKMC_ID) + outputName = PLASTICITY_DISLOKMC_label + thisNoutput => plastic_dislokmc_Noutput + thisOutput => plastic_dislokmc_output + thisSize => plastic_dislokmc_sizePostResult + case (PLASTICITY_DISLOUCLA_ID) + outputName = PLASTICITY_DISLOUCLA_label + thisNoutput => plastic_disloucla_Noutput + thisOutput => plastic_disloucla_output + thisSize => plastic_disloucla_sizePostResult + case (PLASTICITY_TITANMOD_ID) + outputName = PLASTICITY_TITANMOD_label + thisNoutput => plastic_titanmod_Noutput + thisOutput => plastic_titanmod_output + thisSize => plastic_titanmod_sizePostResult + case (PLASTICITY_NONLOCAL_ID) + outputName = PLASTICITY_NONLOCAL_label + thisNoutput => plastic_nonlocal_Noutput + thisOutput => plastic_nonlocal_output + thisSize => plastic_nonlocal_sizePostResult + case default + knownPlasticity = .false. + end select + write(FILEUNIT,'(/,a,/)') '['//trim(phase_name(phase))//']' + if (knownPlasticity) then + write(FILEUNIT,'(a)') '(plasticity)'//char(9)//trim(outputName) + if (phase_plasticity(phase) /= PLASTICITY_NONE_ID) then + do e = 1_pInt,thisNoutput(instance) + write(FILEUNIT,'(a,i4)') trim(thisOutput(e,instance))//char(9),thisSize(e,instance) + enddo + endif endif - endif -#ifdef multiphysicsOut - instance = phase_damageInstance(phase) ! which instance of a plasticity is present phase - knownDamage = .true. - select case(phase_damage(phase)) ! split per constititution - case (LOCAL_DAMAGE_none_ID) - outputName = LOCAL_DAMAGE_NONE_label - thisNoutput => null() - thisOutput => null() - thisSize => null() - case (LOCAL_DAMAGE_isoBrittle_ID) - outputName = LOCAL_DAMAGE_isoBrittle_LABEL - thisNoutput => damage_isoBrittle_Noutput - thisOutput => damage_isoBrittle_output - thisSize => damage_isoBrittle_sizePostResult - case (LOCAL_DAMAGE_isoDuctile_ID) - outputName = LOCAL_DAMAGE_isoDuctile_LABEL - thisNoutput => damage_isoDuctile_Noutput - thisOutput => damage_isoDuctile_output - thisSize => damage_isoDuctile_sizePostResult - case (LOCAL_DAMAGE_anisoBrittle_ID) - outputName = LOCAL_DAMAGE_anisoBrittle_label - thisNoutput => damage_anisoBrittle_Noutput - thisOutput => damage_anisoBrittle_output - thisSize => damage_anisoBrittle_sizePostResult - case (LOCAL_DAMAGE_anisoDuctile_ID) - outputName = LOCAL_DAMAGE_anisoDuctile_LABEL - thisNoutput => damage_anisoDuctile_Noutput - thisOutput => damage_anisoDuctile_output - thisSize => damage_anisoDuctile_sizePostResult - case (LOCAL_DAMAGE_gurson_ID) - outputName = LOCAL_DAMAGE_gurson_label - thisNoutput => damage_gurson_Noutput - thisOutput => damage_gurson_output - thisSize => damage_gurson_sizePostResult - case (LOCAL_DAMAGE_phaseField_ID) - outputName = LOCAL_DAMAGE_phaseField_label - thisNoutput => damage_phaseField_Noutput - thisOutput => damage_phaseField_output - thisSize => damage_phaseField_sizePostResult - case default - knownDamage = .false. - end select - if (knownDamage) then - write(FILEUNIT,'(a)') '(damage)'//char(9)//trim(outputName) - if (phase_damage(phase) /= LOCAL_DAMAGE_none_ID) then - do e = 1_pInt,thisNoutput(instance) - write(FILEUNIT,'(a,i4)') trim(thisOutput(e,instance))//char(9),thisSize(e,instance) - enddo + instance = phase_damageInstance(phase) ! which instance of a plasticity is present phase + knownDamage = .true. + select case(phase_damage(phase)) ! split per constititution + case (LOCAL_DAMAGE_none_ID) + outputName = LOCAL_DAMAGE_NONE_label + thisNoutput => null() + thisOutput => null() + thisSize => null() + case (LOCAL_DAMAGE_isoBrittle_ID) + outputName = LOCAL_DAMAGE_isoBrittle_LABEL + thisNoutput => damage_isoBrittle_Noutput + thisOutput => damage_isoBrittle_output + thisSize => damage_isoBrittle_sizePostResult + case (LOCAL_DAMAGE_isoDuctile_ID) + outputName = LOCAL_DAMAGE_isoDuctile_LABEL + thisNoutput => damage_isoDuctile_Noutput + thisOutput => damage_isoDuctile_output + thisSize => damage_isoDuctile_sizePostResult + case (LOCAL_DAMAGE_anisoBrittle_ID) + outputName = LOCAL_DAMAGE_anisoBrittle_label + thisNoutput => damage_anisoBrittle_Noutput + thisOutput => damage_anisoBrittle_output + thisSize => damage_anisoBrittle_sizePostResult + case (LOCAL_DAMAGE_anisoDuctile_ID) + outputName = LOCAL_DAMAGE_anisoDuctile_LABEL + thisNoutput => damage_anisoDuctile_Noutput + thisOutput => damage_anisoDuctile_output + thisSize => damage_anisoDuctile_sizePostResult + case (LOCAL_DAMAGE_gurson_ID) + outputName = LOCAL_DAMAGE_gurson_label + thisNoutput => damage_gurson_Noutput + thisOutput => damage_gurson_output + thisSize => damage_gurson_sizePostResult + case (LOCAL_DAMAGE_phaseField_ID) + outputName = LOCAL_DAMAGE_phaseField_label + thisNoutput => damage_phaseField_Noutput + thisOutput => damage_phaseField_output + thisSize => damage_phaseField_sizePostResult + case default + knownDamage = .false. + end select + if (knownDamage) then + write(FILEUNIT,'(a)') '(damage)'//char(9)//trim(outputName) + if (phase_damage(phase) /= LOCAL_DAMAGE_none_ID) then + do e = 1_pInt,thisNoutput(instance) + write(FILEUNIT,'(a,i4)') trim(thisOutput(e,instance))//char(9),thisSize(e,instance) + enddo + endif endif - endif - instance = phase_thermalInstance(phase) ! which instance is present phase - knownThermal = .true. - select case(phase_thermal(phase)) ! split per constititution - case (LOCAL_THERMAL_isothermal_ID) - outputName = LOCAL_THERMAL_ISOTHERMAL_label - thisNoutput => null() - thisOutput => null() - thisSize => null() - case (LOCAL_THERMAL_adiabatic_ID) - outputName = LOCAL_THERMAL_ADIABATIC_label - thisNoutput => thermal_adiabatic_Noutput - thisOutput => thermal_adiabatic_output - thisSize => thermal_adiabatic_sizePostResult - case default - knownThermal = .false. - end select - if (knownThermal) then - write(FILEUNIT,'(a)') '(thermal)'//char(9)//trim(outputName) - if (phase_thermal(phase) /= LOCAL_THERMAL_isothermal_ID) then - do e = 1_pInt,thisNoutput(instance) - write(FILEUNIT,'(a,i4)') trim(thisOutput(e,instance))//char(9),thisSize(e,instance) - enddo + instance = phase_thermalInstance(phase) ! which instance is present phase + knownThermal = .true. + select case(phase_thermal(phase)) ! split per constititution + case (LOCAL_THERMAL_isothermal_ID) + outputName = LOCAL_THERMAL_ISOTHERMAL_label + thisNoutput => null() + thisOutput => null() + thisSize => null() + case (LOCAL_THERMAL_adiabatic_ID) + outputName = LOCAL_THERMAL_ADIABATIC_label + thisNoutput => thermal_adiabatic_Noutput + thisOutput => thermal_adiabatic_output + thisSize => thermal_adiabatic_sizePostResult + case default + knownThermal = .false. + end select + if (knownThermal) then + write(FILEUNIT,'(a)') '(thermal)'//char(9)//trim(outputName) + if (phase_thermal(phase) /= LOCAL_THERMAL_isothermal_ID) then + do e = 1_pInt,thisNoutput(instance) + write(FILEUNIT,'(a,i4)') trim(thisOutput(e,instance))//char(9),thisSize(e,instance) + enddo + endif endif - endif - instance = phase_vacancyInstance(phase) ! which instance is present phase - knownVacancy = .true. - select case(phase_vacancy(phase)) ! split per constititution - case (LOCAL_VACANCY_constant_ID) - outputName = LOCAL_VACANCY_constant_label - thisNoutput => null() - thisOutput => null() - thisSize => null() - case (LOCAL_VACANCY_generation_ID) - outputName = LOCAL_VACANCY_generation_label - thisNoutput => vacancy_generation_Noutput - thisOutput => vacancy_generation_output - thisSize => vacancy_generation_sizePostResult - case default - knownVacancy = .false. - end select - if (knownVacancy) then - write(FILEUNIT,'(a)') '(vacancy)'//char(9)//trim(outputName) - if (phase_vacancy(phase) /= LOCAL_VACANCY_constant_ID) then - do e = 1_pInt,thisNoutput(instance) - write(FILEUNIT,'(a,i4)') trim(thisOutput(e,instance))//char(9),thisSize(e,instance) - enddo + instance = phase_vacancyInstance(phase) ! which instance is present phase + knownVacancy = .true. + select case(phase_vacancy(phase)) ! split per constititution + case (LOCAL_VACANCY_constant_ID) + outputName = LOCAL_VACANCY_constant_label + thisNoutput => null() + thisOutput => null() + thisSize => null() + case (LOCAL_VACANCY_generation_ID) + outputName = LOCAL_VACANCY_generation_label + thisNoutput => vacancy_generation_Noutput + thisOutput => vacancy_generation_output + thisSize => vacancy_generation_sizePostResult + case default + knownVacancy = .false. + end select + if (knownVacancy) then + write(FILEUNIT,'(a)') '(vacancy)'//char(9)//trim(outputName) + if (phase_vacancy(phase) /= LOCAL_VACANCY_constant_ID) then + do e = 1_pInt,thisNoutput(instance) + write(FILEUNIT,'(a,i4)') trim(thisOutput(e,instance))//char(9),thisSize(e,instance) + enddo + endif endif - endif -#endif - enddo - close(FILEUNIT) + enddo + close(FILEUNIT) + endif constitutive_maxSizeDotState = 0_pInt constitutive_maxSizePostResults = 0_pInt @@ -1739,7 +1739,6 @@ function constitutive_postResults(Tstar_v, FeArray, ipc, ip, el) plastic_titanmod_postResults use plastic_nonlocal, only: & plastic_nonlocal_postResults -#ifdef multiphysicsOut use damage_isoBrittle, only: & damage_isoBrittle_postResults use damage_isoDuctile, only: & @@ -1756,23 +1755,17 @@ function constitutive_postResults(Tstar_v, FeArray, ipc, ip, el) thermal_adiabatic_postResults use vacancy_generation, only: & vacancy_generation_postResults -#endif implicit none integer(pInt), intent(in) :: & ipc, & !< grain number ip, & !< integration point number el !< element number -#ifdef multiphysicsOut real(pReal), dimension(plasticState(material_phase(ipc,ip,el))%sizePostResults + & damageState( material_phase(ipc,ip,el))%sizePostResults + & thermalState(material_phase(ipc,ip,el))%sizePostResults + & vacancyState(material_phase(ipc,ip,el))%sizePostResults) :: & constitutive_postResults -#else - real(pReal), dimension(plasticState(material_phase(ipc,ip,el))%sizePostResults) :: & - constitutive_postResults -#endif real(pReal), intent(in), dimension(3,3,homogenization_maxNgrains,mesh_maxNips,mesh_NcpElems) :: & FeArray !< elastic deformation gradient real(pReal), intent(in), dimension(6) :: & @@ -1806,7 +1799,6 @@ function constitutive_postResults(Tstar_v, FeArray, ipc, ip, el) plastic_nonlocal_postResults (Tstar_v,FeArray,ip,el) end select -#ifdef multiphysicsOut startPos = endPos + 1_pInt endPos = endPos + damageState(material_phase(ipc,ip,el))%sizePostResults select case (phase_damage(material_phase(ipc,ip,el))) @@ -1837,7 +1829,6 @@ function constitutive_postResults(Tstar_v, FeArray, ipc, ip, el) case (LOCAL_VACANCY_generation_ID) constitutive_postResults(startPos:endPos) = vacancy_generation_postResults(ipc, ip, el) end select -#endif end function constitutive_postResults diff --git a/code/crystallite.f90 b/code/crystallite.f90 index 3c340a511..f0d5c34c2 100644 --- a/code/crystallite.f90 +++ b/code/crystallite.f90 @@ -395,16 +395,18 @@ subroutine crystallite_init !-------------------------------------------------------------------------------------------------- ! write description file for crystallite output - call IO_write_jobFile(FILEUNIT,'outputCrystallite') + if (worldrank == 0_pInt) then + call IO_write_jobFile(FILEUNIT,'outputCrystallite') - do p = 1_pInt,material_Ncrystallite - write(FILEUNIT,'(/,a,/)') '['//trim(crystallite_name(p))//']' - do e = 1_pInt,crystallite_Noutput(p) - write(FILEUNIT,'(a,i4)') trim(crystallite_output(e,p))//char(9),crystallite_sizePostResult(e,p) + do p = 1_pInt,material_Ncrystallite + write(FILEUNIT,'(/,a,/)') '['//trim(crystallite_name(p))//']' + do e = 1_pInt,crystallite_Noutput(p) + write(FILEUNIT,'(a,i4)') trim(crystallite_output(e,p))//char(9),crystallite_sizePostResult(e,p) + enddo enddo - enddo - close(FILEUNIT) + close(FILEUNIT) + endif !-------------------------------------------------------------------------------------------------- ! initialize @@ -4229,18 +4231,12 @@ function crystallite_postResults(ipc, ip, el) ip, & !< integration point index ipc !< grain index -#ifdef multiphysicsOut real(pReal), dimension(1+crystallite_sizePostResults(microstructure_crystallite(mesh_element(4,el))) + & 1+plasticState(material_phase(ipc,ip,el))%sizePostResults + & damageState( material_phase(ipc,ip,el))%sizePostResults + & thermalState(material_phase(ipc,ip,el))%sizePostResults + & vacancyState(material_phase(ipc,ip,el))%sizePostResults) :: & crystallite_postResults -#else - real(pReal), dimension(1+crystallite_sizePostResults(microstructure_crystallite(mesh_element(4,el)))+ & - 1+plasticState(material_phase(ipc,ip,el))%sizePostResults) :: & - crystallite_postResults -#endif real(pReal), dimension(3,3) :: & Ee real(pReal), dimension(4) :: & diff --git a/code/homogenization.f90 b/code/homogenization.f90 index 0627d6171..1b910a240 100644 --- a/code/homogenization.f90 +++ b/code/homogenization.f90 @@ -228,47 +228,46 @@ subroutine homogenization_init() !-------------------------------------------------------------------------------------------------- ! write description file for homogenization output - call IO_write_jobFile(FILEUNIT,'outputHomogenization') - do p = 1,material_Nhomogenization - i = homogenization_typeInstance(p) ! which instance of this homogenization type - knownHomogenization = .true. ! assume valid - select case(homogenization_type(p)) ! split per homogenization type - case (HOMOGENIZATION_NONE_ID) - outputName = HOMOGENIZATION_NONE_label - thisNoutput => null() - thisOutput => null() - thisSize => null() - case (HOMOGENIZATION_ISOSTRAIN_ID) - outputName = HOMOGENIZATION_ISOSTRAIN_label - thisNoutput => homogenization_isostrain_Noutput - thisOutput => homogenization_isostrain_output - thisSize => homogenization_isostrain_sizePostResult - case (HOMOGENIZATION_RGC_ID) - outputName = HOMOGENIZATION_RGC_label - thisNoutput => homogenization_RGC_Noutput - thisOutput => homogenization_RGC_output - thisSize => homogenization_RGC_sizePostResult - case default - knownHomogenization = .false. - end select - write(FILEUNIT,'(/,a,/)') '['//trim(homogenization_name(p))//']' - if (knownHomogenization) then - write(FILEUNIT,'(a)') '(type)'//char(9)//trim(outputName) - write(FILEUNIT,'(a,i4)') '(ngrains)'//char(9),homogenization_Ngrains(p) - if (homogenization_type(p) /= HOMOGENIZATION_NONE_ID) then - do e = 1,thisNoutput(i) - write(FILEUNIT,'(a,i4)') trim(thisOutput(e,i))//char(9),thisSize(e,i) - enddo + if (worldrank == 0_pInt) then + call IO_write_jobFile(FILEUNIT,'outputHomogenization') + do p = 1,material_Nhomogenization + i = homogenization_typeInstance(p) ! which instance of this homogenization type + knownHomogenization = .true. ! assume valid + select case(homogenization_type(p)) ! split per homogenization type + case (HOMOGENIZATION_NONE_ID) + outputName = HOMOGENIZATION_NONE_label + thisNoutput => null() + thisOutput => null() + thisSize => null() + case (HOMOGENIZATION_ISOSTRAIN_ID) + outputName = HOMOGENIZATION_ISOSTRAIN_label + thisNoutput => homogenization_isostrain_Noutput + thisOutput => homogenization_isostrain_output + thisSize => homogenization_isostrain_sizePostResult + case (HOMOGENIZATION_RGC_ID) + outputName = HOMOGENIZATION_RGC_label + thisNoutput => homogenization_RGC_Noutput + thisOutput => homogenization_RGC_output + thisSize => homogenization_RGC_sizePostResult + case default + knownHomogenization = .false. + end select + write(FILEUNIT,'(/,a,/)') '['//trim(homogenization_name(p))//']' + if (knownHomogenization) then + write(FILEUNIT,'(a)') '(type)'//char(9)//trim(outputName) + write(FILEUNIT,'(a,i4)') '(ngrains)'//char(9),homogenization_Ngrains(p) + if (homogenization_type(p) /= HOMOGENIZATION_NONE_ID) then + do e = 1,thisNoutput(i) + write(FILEUNIT,'(a,i4)') trim(thisOutput(e,i))//char(9),thisSize(e,i) + enddo + endif endif - endif -#ifdef multiphysicsOut - write(FILEUNIT,'(a)') '(field)' - do e = 1_pInt,field_Noutput(p) - write(FILEUNIT,'(a,i4)') trim(field_output(e,p))//char(9),field_sizePostResult(e,p) + do e = 1_pInt,field_Noutput(p) + write(FILEUNIT,'(a,i4)') trim(field_output(e,p))//char(9),field_sizePostResult(e,p) + enddo enddo -#endif - enddo - close(FILEUNIT) + close(FILEUNIT) + endif !-------------------------------------------------------------------------------------------------- ! allocate and initialize global variables @@ -310,16 +309,12 @@ subroutine homogenization_init() enddo materialpoint_sizeResults = 1 & ! grain count + 1 + homogenization_maxSizePostResults & ! homogSize & homogResult -#ifdef multiphysicsOut + field_maxSizePostResults & ! field size & field result -#endif + homogenization_maxNgrains * (1 + crystallite_maxSizePostResults & ! crystallite size & crystallite results -#ifdef multiphysicsOut + + 1 + constitutive_maxSizePostResults & ! constitutive size & constitutive results + constitutive_damage_maxSizePostResults & + constitutive_thermal_maxSizePostResults & - + constitutive_vacancy_maxSizePostResults & -#endif - + 1 + constitutive_maxSizePostResults) ! constitutive size & constitutive results + + constitutive_vacancy_maxSizePostResults) allocate(materialpoint_results(materialpoint_sizeResults,mesh_maxNips,mesh_NcpElems)) mainProcess: if (worldrank == 0) then @@ -751,28 +746,21 @@ subroutine materialpoint_postResults thePos = thePos + theSize endif -#ifdef multiphysicsOut theSize = field_sizePostResults(mappingHomogenization(2,i,e)) if (theSize > 0_pInt) then ! any homogenization results to mention? materialpoint_results(thePos+1:thePos+theSize,i,e) = field_postResults(i,e) ! tell field results thePos = thePos + theSize endif -#endif materialpoint_results(thePos+1,i,e) = real(myNgrains,pReal) ! tell number of grains at materialpoint thePos = thePos + 1_pInt grainLooping :do g = 1,myNgrains -#ifdef multiphysicsOut theSize = 1 + crystallite_sizePostResults(myCrystallite) + & 1 + plasticState(material_phase(g,i,e))%sizePostResults + & !ToDo damageState(material_phase(g,i,e))%sizePostResults + & thermalState(material_phase(g,i,e))%sizePostResults + & vacancyState(material_phase(g,i,e))%sizePostResults -#else - theSize = (1 + crystallite_sizePostResults(myCrystallite)) + & - (1 + plasticState(material_phase(g,i,e))%sizePostResults) -#endif materialpoint_results(thePos+1:thePos+theSize,i,e) = crystallite_postResults(g,i,e) ! tell crystallite results thePos = thePos + theSize enddo grainLooping