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

YAML structure allows more than one mechanism

This commit is contained in:
Martin Diehl 2021-04-11 12:42:57 +02:00
parent 081d51f224
commit 9d90ed525e
1 changed files with 9 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
src/phase_damage.f90
View File

@ -105,7 +105,8 @@ module subroutine damage_init
class(tNode), pointer :: & class(tNode), pointer :: &
phases, & phases, &
phase, & phase, &
sources sources, &
source
logical:: damage_active logical:: damage_active
print'(/,a)', ' <<<+- phase:damage init -+>>>' print'(/,a)', ' <<<+- phase:damage init -+>>>'
@ -113,7 +114,6 @@ module subroutine damage_init
phases => config_material%get('phase') phases => config_material%get('phase')
allocate(current(phases%length)) allocate(current(phases%length))
allocate(damageState (phases%length)) allocate(damageState (phases%length))
allocate(param(phases%length)) allocate(param(phases%length))
@ -127,11 +127,14 @@ module subroutine damage_init
phase => phases%get(ph) phase => phases%get(ph)
sources => phase%get('damage',defaultVal=emptyList) sources => phase%get('damage',defaultVal=emptyList)
param(ph)%K(1,1) = sources%get_asFloat('K_11',defaultVal=0.0_pReal)
param(ph)%K(2,2) = sources%get_asFloat('K_22',defaultVal=0.0_pReal)
param(ph)%K(3,3) = sources%get_asFloat('K_33',defaultVal=0.0_pReal)
if (sources%length > 1) error stop if (sources%length > 1) error stop
if (sources%length == 1) damage_active = .true. if (sources%length == 1) then
damage_active = .true.
source => sources%get(1)
param(ph)%K(1,1) = source%get_asFloat('K_11',defaultVal=0.0_pReal)
param(ph)%K(2,2) = source%get_asFloat('K_22',defaultVal=0.0_pReal)
param(ph)%K(3,3) = source%get_asFloat('K_33',defaultVal=0.0_pReal)
endif
enddo enddo