! TODO: we repeat the center node N times because the corner results have not been calculated
WRITE(OP2) (EID_OUT_ARRAY(I,1)*10+DEVICE_CODE, CID, CEN_WORD, NNODES-1, &
! grid_id
! 21
(GID_OUT_ARRAY(I,J), &
! oxx txy s1 a1 a2 a3
REAL(OGEL(I,1),4), REAL(OGEL(I,4),4), REAL(OGEL(I,9), 4), 0., 0., 0., &
! p ovm
REAL(OGEL(I,12),4), REAL(OGEL(I,7),4), &
! syy tyz s2 b1 b2 b3
REAL(OGEL(I,2),4), REAL(OGEL(I,5),4), REAL(OGEL(I,10),4), 0., 0., 0., &
! szz txz s3 c1 c2 c3
REAL(OGEL(I,3),4), REAL(OGEL(I,6),4), REAL(OGEL(I,11),4), 0., 0., 0., &
J=1,NNODES), I=1,NUM)
J is not used to index OGEL, so it writes the same values for each iteration of the J loop. OGEL does contain corner stress which is used correctly for F06 in the code below that.
Problem should be in WRITE_ELEM_STRESSES.f90:
J is not used to index OGEL, so it writes the same values for each iteration of the J loop. OGEL does contain corner stress which is used correctly for F06 in the code below that.