anchor.diff:        real(rs),intent(in) :: VF_in(:), C_in(:,:,:), rh_in(:)
         real(rs),intent(out) :: Cave(6,6), rhave
         integer :: i,n
  -!      real(rs),allocatable :: VF(:)
         
   !  Get size of arrays and check they're consistent
         n = size(VF_in)
  
positive.code:       real(rs),intent(in) :: VF_in(:), C_in(:,:,:), rh_in(:)
        real(rs),intent(out) :: Cave(6,6), rhave
        integer :: i,n
  !      real(rs),allocatable :: VF(:)
        real(rs) :: S(6,6), S_in(6,6), C_temp(6,6)
        
  !  Get size of arrays and check they're consistent
  
negative.code:       
  !  Find stiffness from compliance
        call inverse(6,6,S,Cave)
              
     end subroutine CIJ_Reuss_av
  !-------------------------------------------------------------------------------
  
  