anchor.diff:       &     infmat(1,imat),matgpt,elas_strs_1,                           ! materl
        &     gauss,sh,shj,infetype,netypes,                               ! eltype
        &     skew,nskdim,numrot,                                          ! skew
  -     &     getshape,bmatrix,idebug,idout,kto,kw,ierr)                   ! info
  +     &     getshape,bmatrix,                                            ! bbar
  +     &     ierr)                                                        ! errcode
           else if(matmodel.eq.2) then
             call elas_strs_cmp_ss(
        &     b,neq,                                                       ! force
  
positive.code:      &     infmat(1,imat),matgpt,elas_strs_6,                           ! materl
       &     gauss,sh,shj,infetype,netypes,                               ! eltype
       &     skew,nskdim,numrot,                                          ! skew
       &     getshape,bmatrix,idebug,idout,kto,kw,ierr)                   ! info
          else if(matmodel.eq.7) then
            call elas_strs_cmp_ss(
       &     b,neq,                                                       ! force
  
negative.code:      &     infmat(1,imat),matgpt,elas_strs_9,                           ! materl
       &     gauss,sh,shj,infetype,netypes,                               ! eltype
       &     skew,nskdim,numrot,                                          ! skew
       &     getshape,bmatrix,                                            ! bbar
       &     ierr)                                                        ! errcode
          else if(matmodel.eq.10) then
            call elas_strs_cmp_ss(
       &     b,neq,                                                       ! force
  