anchor.diff:           stop
         endif
         
  -!  Allocate space for arrays
  -!      allocate(VF(n))
  -      
  -!  Normalise the volume fractions to sum to unity
  -!      VF = VF_in / sum(VF_in)
  -      
   !  Construct Voigt average
         Cave = 0.  ;  rhave = 0.
         do i=1,n
  
positive.code:       
  !  Find stiffness from compliance
        call inverse(6,6,S,Cave)
        
  !      deallocate(VF)
        
     end subroutine CIJ_Reuss_av
  !-------------------------------------------------------------------------------
  
  
negative.code:       real(rs),intent(in) :: VF_in(:), C_in(:,:,:), rh_in(:)
        real(rs),intent(out) :: Cave(6,6), rhave
        integer :: i,n
        real(rs) :: S(6,6), S_in(6,6), C_temp(6,6)
        
  !  Get size of arrays and check they're consistent
  