anchor.diff:      real(kind=kreal)   :: s11, s22, s33, s12, s23, s13, t11, t22, t33, t12, t23, t13, ps, smises, tmises
       integer(kind=kint), allocatable :: nnumber(:)
       type(fstr_solid_physic_val), pointer :: layer => null()
  +
  +    fstrSOLID%ESTRAIN = 0.0d0
  +    fstrSOLID%ESTRESS = 0.0d0
   
       n_totlyr = fstrSOLID%max_lyr
   
  
positive.code:     tnstrain => fstrSOLID%tnstrain
      testrain => fstrSOLID%testrain
  
      allocate( nnumber(hecMESH%n_node) )
      allocate( fstrSOLID%is_rot(hecMESH%n_node) )
      nnumber = 0
  
negative.code:     end if
    end subroutine make_principal
  
  end module m_fstr_NodalStress