Currently I have 18 variables to output, and I did like this in the user element:
valp(1)=rho
valp(2)=rho_ele
valp(3)=rho_hol
.........
valp(17)=sigma(2)
valp(18)=sigma(3)
Now the output file works just fine, however when feap proceed to get vtk file for paraview, such an error occurs: *** Error in `feap': corrupted double-linked list: 0x0000000001113ce0 ***
when I change back to only 9 valps, it works again.
I did change the "stress" options in umacr04, like this:
if(foundnp) then
write(17,'(a,2i10,a)') 'stress',12,numnp,' float'
do n=1,numnp
write(17,'(12f30.11)') (hr(begnp-1+i*numnp+n),i=1,12)
enddo
endif
and all these efforts didn't work.
So, is there a way to give more than 9 valps?