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:       ! Non-exponential display (all values on same scale)
        if (.not.iexpo) then
           ! If necessary, work out max for matrix
           if (autopower) ipower = log10(maxval(abs(C)))
           write(fmt,'(a,i0.0,a,i0.0,a)') '(6(f',indp+3,'.',indp,',1x))'
           write(iunit,'(a,i0.0)') 'x 10^',ipower
           write(iunit,fmt) C/(10.**ipower)
  
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
  !-------------------------------------------------------------------------------
  
  !===============================================================================
  