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

added missing lp output to crystallite

This commit is contained in:
Martin Diehl 2013-12-20 11:13:12 +00:00
parent de3fc70a02
commit 831e905b16
2 changed files with 13 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
code/constitutive_phenopowerlaw.f90
View File

@ -273,7 +273,6 @@ subroutine constitutive_phenopowerlaw_init(fileUnit)
tag = IO_lc(IO_stringValue(line,positions,1_pInt)) ! extract key tag = IO_lc(IO_stringValue(line,positions,1_pInt)) ! extract key
select case(tag) select case(tag)
case ('plasticity','elasticity') case ('plasticity','elasticity')
cycle
case ('(output)') case ('(output)')
constitutive_phenopowerlaw_Noutput(i) = constitutive_phenopowerlaw_Noutput(i) + 1_pInt constitutive_phenopowerlaw_Noutput(i) = constitutive_phenopowerlaw_Noutput(i) + 1_pInt
constitutive_phenopowerlaw_output(constitutive_phenopowerlaw_Noutput(i),i) = & constitutive_phenopowerlaw_output(constitutive_phenopowerlaw_Noutput(i),i) = &
@ -305,16 +304,16 @@ subroutine constitutive_phenopowerlaw_init(fileUnit)
case ('lattice_structure') case ('lattice_structure')
structure = IO_lc(IO_stringValue(line,positions,2_pInt)) structure = IO_lc(IO_stringValue(line,positions,2_pInt))
select case(structure(1:3)) select case(structure(1:3))
case(LATTICE_iso_label) case(LATTICE_iso_label)
constitutive_phenopowerlaw_structureID(i) = LATTICE_iso_ID constitutive_phenopowerlaw_structureID(i) = LATTICE_iso_ID
case(LATTICE_fcc_label) case(LATTICE_fcc_label)
constitutive_phenopowerlaw_structureID(i) = LATTICE_fcc_ID constitutive_phenopowerlaw_structureID(i) = LATTICE_fcc_ID
case(LATTICE_bcc_label) case(LATTICE_bcc_label)
constitutive_phenopowerlaw_structureID(i) = LATTICE_bcc_ID constitutive_phenopowerlaw_structureID(i) = LATTICE_bcc_ID
case(LATTICE_hex_label) case(LATTICE_hex_label)
constitutive_phenopowerlaw_structureID(i) = LATTICE_hex_ID constitutive_phenopowerlaw_structureID(i) = LATTICE_hex_ID
case(LATTICE_ort_label) case(LATTICE_ort_label)
constitutive_phenopowerlaw_structureID(i) = LATTICE_ort_ID constitutive_phenopowerlaw_structureID(i) = LATTICE_ort_ID
end select end select
configNchunks = lattice_configNchunks(constitutive_phenopowerlaw_structureID(i)) configNchunks = lattice_configNchunks(constitutive_phenopowerlaw_structureID(i))
Nchunks_SlipFamilies = configNchunks(1) Nchunks_SlipFamilies = configNchunks(1)

3
code/crystallite.f90
View File

@ -308,6 +308,7 @@ subroutine crystallite_init(temperature)
if (IO_getTag(line,'[',']') /= '') then ! next section if (IO_getTag(line,'[',']') /= '') then ! next section
section = section + 1_pInt section = section + 1_pInt
output = 0_pInt ! reset output counter output = 0_pInt ! reset output counter
cycle ! skip to next line
endif endif
if (section > 0_pInt) then if (section > 0_pInt) then
positions = IO_stringPos(line,maxNchunks) positions = IO_stringPos(line,maxNchunks)
@ -343,6 +344,8 @@ subroutine crystallite_init(temperature)
crystallite_outputID(output,section) = fe_ID crystallite_outputID(output,section) = fe_ID
case ('fp') case ('fp')
crystallite_outputID(output,section) = fp_ID crystallite_outputID(output,section) = fp_ID
case ('lp')
crystallite_outputID(output,section) = lp_ID
case ('e') case ('e')
crystallite_outputID(output,section) = e_ID crystallite_outputID(output,section) = e_ID
case ('ee') case ('ee')