anchor.diff:      end if
     end subroutine make_principal
   
  -end module m_fstr_NodalStress
  +end module m_fstr_NodalStress
positive.code:     tnstrain => fstrSOLID%tnstrain
      testrain => fstrSOLID%testrain
  
      allocate( nnumber(hecMESH%n_node) )
      allocate( fstrSOLID%is_rot(hecMESH%n_node) )
      nnumber = 0
  
negative.code: !C** Shell33 variables
      integer(kind=kint) :: isect, ihead, ntot_lyr, nlyr, flag33, cid, truss
      real(kind=kreal)   :: thick, thick_lyr, dtot_lyr
  
      fstrSOLID%STRAIN  = 0.0d0
      fstrSOLID%STRESS  = 0.0d0
      fstrSOLID%MISES   = 0.0d0
      fstrSOLID%ESTRAIN = 0.0d0
      fstrSOLID%ESTRESS = 0.0d0
      fstrSOLID%EMISES  = 0.0d0
  
      allocate( nnumber(hecMESH%n_node) )
      allocate( fstrSOLID%is_rot(hecMESH%n_node) )
  