223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

write all outputs to HDF5 

not 100% sure about the units
This commit is contained in:
Martin Diehl 2019-05-12 15:14:29 +02:00
parent 434813e6e5
commit 6e96a57781
1 changed files with 13 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
src/homogenization_mech_RGC.f90
View File

@ -1097,6 +1097,7 @@ end function homogenization_RGC_postResults
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief writes results to HDF5 output file !> @brief writes results to HDF5 output file
! ToDo: check wheter units are correct
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
subroutine mech_RGC_results(instance,group) subroutine mech_RGC_results(instance,group)
#if defined(PETSc) || defined(DAMASK_HDF5) #if defined(PETSc) || defined(DAMASK_HDF5)
@ -1116,12 +1117,20 @@ subroutine mech_RGC_results(instance,group)
call results_writeDataset(group,stt%work,'W',& call results_writeDataset(group,stt%work,'W',&
'work density','J/m³') 'work density','J/m³')
case (magnitudemismatch_ID) case (magnitudemismatch_ID)
call results_writeDataset(group,dst%mismatch,'M',& call results_writeDataset(group,dst%mismatch,'N',&
'average mismatch tensor','tbd') 'average mismatch tensor','m')
case (penaltyenergy_ID) case (penaltyenergy_ID)
call results_writeDataset(group,stt%penaltyEnergy,'R',& call results_writeDataset(group,stt%penaltyEnergy,'R',&
'mismatch penalty density','tbd') 'mismatch penalty density','J/m³')
case (volumediscrepancy_ID)
call results_writeDataset(group,dst%volumeDiscrepancy,'Delta_V',&
'volume discrepancy','m³')
case (maximumrelaxrate_ID)
call results_writeDataset(group,dst%relaxationrate_max,'max_alpha_dot',&
'maximum relaxation rate','m/s')
case (averagerelaxrate_ID)
call results_writeDataset(group,dst%relaxationrate_avg,'avg_alpha_dot',&
'average relaxation rate','m/s')
end select end select
enddo outputsLoop enddo outputsLoop
end associate end associate