anchor.diff:       & infmat,numat,                                                    ! materl
        & gauss,sh,shj,infetype,netypes,                                   ! eltype
        & skew,nskdim,numrot,                                              ! skew
  -     & getshape,bmatrix,idebug,idout,kto,kw,ierr)                       ! info
  +     & getshape,bmatrix,                                                ! bbar
  +     & ierr)                                                            ! errcode
   c
   c...program to compute the total stress and strain for the current
   c   iteration (elastic case).  This driver is for the small strain case.
  
positive.code:      &     infmat(1,imat),matgpt,elas_strs_16,                          ! materl
       &     gauss,sh,shj,infetype,netypes,                               ! eltype
       &     skew,nskdim,numrot,                                          ! skew
       &     getshape,bmatrix,idebug,idout,kto,kw,ierr)                   ! info
          else if(matmodel.eq.17) then
            call elas_strs_cmp_ss(
       &     b,neq,                                                       ! force
  
negative.code:      &     infmat(1,imat),matgpt,elas_strs_6,                           ! materl
       &     gauss,sh,shj,infetype,netypes,                               ! eltype
       &     skew,nskdim,numrot,                                          ! skew
       &     getshape,bmatrix,                                            ! bbar
       &     ierr)                                                        ! errcode
          else if(matmodel.eq.7) then
            call elas_strs_cmp_ss(
       &     b,neq,                                                       ! force
  