Comment:
subroutine umatl1(eps,theta,td,d,ud,hn,h1,nh,ii,istrt, sig,dd,isw)
real*8 eps(*), theta(*),d(*),ud(*),hn(nh),h1(nh), sig(*),dd(6,*)
Uncomment:
c subroutine umatl1( f , detf,td,d,ud,hn,h1,nh,ii,istrt, sig,dd,isw)
c real*8 f(3,3,*),detf(*),d(*),ud(*),hn(nh),h1(nh), sig(*),dd(6,*)
Then access f as a two-dimensional array, so that
do k=1:3
c(i,j) = c(i,j) + f(k,i)*f(k,j)
end do