anchor.diff:  !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) )
  
positive.code:     end if
    end subroutine make_principal
  
  end module m_fstr_NodalStress
  
negative.code:     tnstrain => fstrSOLID%tnstrain
      testrain => fstrSOLID%testrain
  
      fstrSOLID%STRAIN = 0.0d0
      fstrSOLID%STRESS = 0.0d0
      allocate( nnumber(hecMESH%n_node) )
      allocate( fstrSOLID%is_rot(hecMESH%n_node) )
      nnumber = 0
  