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

adapted to new default table style in 2007R1 

hypoelastic part of input file needs skipping an additional line if new style tables are active
This commit is contained in:
Philip Eisenlohr 2007-09-28 14:16:03 +00:00
parent 28fc9cc221
commit 892f33591c
1 changed files with 5 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
trunk/mesh.f90
View File

@ -135,6 +135,7 @@
implicit none implicit none
integer(pInt), parameter :: fileUnit = 222 integer(pInt), parameter :: fileUnit = 222
integer(pInt), parameter :: hypoelasticTableStyle = 0
mesh_Nelems = 0_pInt mesh_Nelems = 0_pInt
mesh_NcpElems = 0_pInt mesh_NcpElems = 0_pInt
@ -291,11 +292,13 @@ candidate: do i=1,minN ! iterate over lonelyNode's shared elements
pos = IO_stringPos(line,20) pos = IO_stringPos(line,20)
select case ( IO_lc(IO_Stringvalue(line,pos,1))) select case ( IO_lc(IO_Stringvalue(line,pos,1)))
case('table')
hypoelasticTableStyle = IO_IntValue (line,pos,5)
case('sizing') case('sizing')
mesh_Nelems = IO_IntValue (line,pos,3) mesh_Nelems = IO_IntValue (line,pos,3)
mesh_Nnodes = IO_IntValue (line,pos,4) mesh_Nnodes = IO_IntValue (line,pos,4)
case('hypoelastic') case('hypoelastic')
do i=1,4 do i=1,4+hypoTableStyle
read (unit,610,END=620) line read (unit,610,END=620) line
end do end do
pos = IO_stringPos(line,20) pos = IO_stringPos(line,20)
@ -447,7 +450,7 @@ candidate: do i=1,minN ! iterate over lonelyNode's shared elements
read (unit,610,END=620) line read (unit,610,END=620) line
pos = IO_stringPos(line,1) pos = IO_stringPos(line,1)
if( IO_lc(IO_stringValue(line,pos,1)) == 'hypoelastic' ) then if( IO_lc(IO_stringValue(line,pos,1)) == 'hypoelastic' ) then
do i=1,3 ! skip three lines do i=1,3+hypoelasticTableStyle ! skip three (or four if new table style!) lines
read (unit,610,END=620) line read (unit,610,END=620) line
end do end do
contInts = IO_continousIntValues(unit,mesh_NcpElems) contInts = IO_continousIntValues(unit,mesh_NcpElems)