From a89efaa4a6b7db2a4a204cbfa662d780df77bb05 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Sat, 9 Feb 2013 08:23:47 +0000 Subject: [PATCH] added missing bracket and substituted one more forall by do loop --- code/mesh.f90 | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/code/mesh.f90 b/code/mesh.f90 index 3a579498e..633eb1fd7 100644 --- a/code/mesh.f90 +++ b/code/mesh.f90 @@ -3192,7 +3192,9 @@ use numerics, only: numerics_unitlength read (myUnit,610,END=670) line myPos = IO_stringPos(line,maxNchunks) m = mesh_FEasCP('node',IO_intValue(line,myPos,1_pInt)) - forall (j=1_pInt:3_pInt) mesh_node0(j,m) = numerics_unitlength * IO_floatValue(line,myPos,j+1_pInt) + do j=1_pInt, 3_pInt + mesh_node0(j,m) = numerics_unitlength * IO_floatValue(line,myPos,j+1_pInt) + enddo enddo endif enddo @@ -3321,7 +3323,7 @@ subroutine mesh_abaqus_build_elements(myUnit) if (e /= 0_pInt) then ! disregard non CP elems mesh_element(1,e) = IO_intValue(line,myPos,1_pInt) ! FE id mesh_element(2,e) = t ! elem type - do j=1_pInt,FE_Nnodes(FE_geomtype(t))) + do j=1_pInt,FE_Nnodes(FE_geomtype(t)) mesh_element(4_pInt+j,e) = IO_intValue(line,myPos,1_pInt+j) ! copy FE ids of nodes to position 5: enddo call IO_skipChunks(myUnit,FE_NoriginalNodes(t)-(myPos(1_pInt)-1_pInt)) ! read on (even multiple lines) if FE_NoriginalNodes exceeds required node count