From 9ff31d4f2bc9a3ffea3b4b17f739b5637c1b2312 Mon Sep 17 00:00:00 2001 From: Sharan Date: Tue, 22 Feb 2022 23:16:14 +0100 Subject: [PATCH] suggested changes from MR !530 --- .../homogenization/mechanical/RGC_8grains.yaml | 2 +- src/homogenization_mechanical.f90 | 2 +- src/homogenization_thermal_pass.f90 | 2 +- src/phase_thermal.f90 | 12 +++++------- 4 files changed, 8 insertions(+), 10 deletions(-) diff --git a/examples/config/homogenization/mechanical/RGC_8grains.yaml b/examples/config/homogenization/mechanical/RGC_8grains.yaml index 2506a9bcf..cb7fba07b 100644 --- a/examples/config/homogenization/mechanical/RGC_8grains.yaml +++ b/examples/config/homogenization/mechanical/RGC_8grains.yaml @@ -1,4 +1,4 @@ -# For Relaxed Grain Cluster homogenization, requires N_constituents = 8 +# For Relaxed Grain Cluster homogenization, requires N_constituents = prod(cluster_size) type: RGC D_alpha: [4.0e-06, 4.0e-06, 2.0e-06] a_g: [0.0, 0.0, 0.0] diff --git a/src/homogenization_mechanical.f90 b/src/homogenization_mechanical.f90 index c4b74a15d..f1ca24e9b 100644 --- a/src/homogenization_mechanical.f90 +++ b/src/homogenization_mechanical.f90 @@ -204,7 +204,7 @@ module subroutine mechanical_results(group_base,ho) 'deformation gradient','1') case('P') call results_writeDataset(reshape(homogenization_P,[3,3,discretization_nCells]),group,'P', & - 'deformation gradient','1') + 'first Piola-Kirchhoff stress','Pa') end select end do diff --git a/src/homogenization_thermal_pass.f90 b/src/homogenization_thermal_pass.f90 index ac12aa0b1..f06fe4cf8 100644 --- a/src/homogenization_thermal_pass.f90 +++ b/src/homogenization_thermal_pass.f90 @@ -11,7 +11,7 @@ module subroutine pass_init() print'(/,1x,a)', '<<<+- homogenization:thermal:pass init -+>>>' if (homogenization_Nconstituents(1) /= 1) & - call IO_error(211,ext_msg='N_constituents (pass)') + call IO_error(211,ext_msg='(pass) with N_constituents !=1') end subroutine pass_init diff --git a/src/phase_thermal.f90 b/src/phase_thermal.f90 index 924ae2f5d..2aab37529 100644 --- a/src/phase_thermal.f90 +++ b/src/phase_thermal.f90 @@ -388,7 +388,7 @@ end function thermal_active !---------------------------------------------------------------------------------------------- -!< @brief writes damage sources results to HDF5 output file +!< @brief writes thermal sources results to HDF5 output file !---------------------------------------------------------------------------------------------- module subroutine thermal_results(group,ph) @@ -398,18 +398,16 @@ module subroutine thermal_results(group,ph) integer :: ou - if (allocated(param(ph)%output)) then - call results_closeGroup(results_addGroup(group//'thermal')) - else - return - endif + if (.not. allocated(param(ph)%output)) return + + call results_closeGroup(results_addGroup(group//'thermal')) do ou = 1, size(param(ph)%output) select case(trim(param(ph)%output(ou))) case ('T') - call results_writeDataset(current(ph)%T,group//'thermal','T', 'temperature','T') + call results_writeDataset(current(ph)%T,group//'thermal','T', 'temperature','K') end select