anchor.diff:        C = phi*A
         F = A - 2._rs*L
         C12 = A - 2._rs*N
  -      
  +
         CIJ_panning_VTI = reshape( &
               (/ A , C12, F, O, O, O, &
                 C12,  A , F, O, O, O, &
  
positive.code:             CR(n,m) = x
           enddo
        enddo
        
     end subroutine CIJ_rot3_old
  !-------------------------------------------------------------------------------
  
  
negative.code:          end if
         end do ! i
        end if
        ! FINISHED INNER REDUCTION
      end do ! k
      ! END OF MAIN REDUCTION LOOP
  
      !                              UNSCRAMBLE ROWS
      do j=1,n
        do i=1,n
  