anchor.diff:        real(rs) :: M(21,21)    ! Projector
         real(rs) :: C(6,6)
         real(rs) :: X(21),XH(21),Xin(21),CH(6,6)
  -      integer  :: i,j
  +      integer  :: i
         
   !  Input matrix and vector
         C = Cin
  
positive.code:       real(rs) :: C(6,6), CR(6,6) ! Voigt notation matrix
        real(rs) :: alp,bet,gam ! rotation (clockwise) about 1,2,3 axis respectively
        real(rs) :: a,b,g     
        real(rs) :: R(3,3), D(6,6), R1(3,3), R2(3,3), R3(3,3), R21(3,3)
        
        integer :: i,j,k,l,m,n,lp,lq,lt
              
  
negative.code: !-------------------------------------------------------------------------------
  
  !===============================================================================
  !   subroutine CIJ_isotropic_average(C,r,Ciso,Vp,Vs)
  !===============================================================================
  !      implicit none
  !      
  !      real(rs), intent(in) :: C(6,6), r
  !      real(rs), intent(out), optional :: Ciso(6,6), Vp, Vs
  !      
  !      write(0,'(a)') 'anisotropy_ajn: CIJ_isotropic_average is not working yet.'
  !      stop
  !      
  !   end subroutine CIJ_isotropic_average
  !-------------------------------------------------------------------------------
  
  !===============================================================================
  