anchor.diff:        enddo
         do i=1,ntot_lyr
           call fstr_make_result_main( hecMESH, fstrSOLID, fstrRESULT, &
  -          & fstrSOLID%SHELL%LAYER(i)%PLUS,  nitem, iitem, ncomp, i+1, clyr(2*i-1) )
  +          & fstrSOLID%SHELL%LAYER(i)%PLUS,  nitem, iitem, ncomp, eitem, jitem, ecomp, i+1, clyr(2*i-1) )
           call fstr_make_result_main( hecMESH, fstrSOLID, fstrRESULT, &
  -          & fstrSOLID%SHELL%LAYER(i)%MINUS, nitem, iitem, ncomp, i+1, clyr(2*i  ) )
  +          & fstrSOLID%SHELL%LAYER(i)%MINUS, nitem, iitem, ncomp, eitem, jitem, ecomp, i+1, clyr(2*i  ) )
         enddo
         deallocate(clyr)
       endif
  
positive.code: 
    end subroutine fstr_make_result
  
    subroutine fstr_make_result_main( hecMESH, fstrSOLID, fstrRESULT, RES, nitem, iitem, ncomp, nlyr, clyr )
      use m_fstr
      use m_out
      use m_static_lib
  
negative.code:     character(len=HECMW_NAME_LEN)   :: s, label, nameID, addfname
      character(len=6)                :: clyr
      character(len=4)                :: cnum
      integer(kind=kint) :: i, j, k, ndof, mdof, id, nitem, eitem, nn, mm, ngauss, it
      integer(kind=kint) :: iitem, ncomp, jitem, ecomp, nlyr
  
      ndof = hecMESH%n_dof
  
  