anchor.diff:      double precision, intent (in), dimension (:, :) :: Atheta
       type (meshSphereEquator), intent (in) :: meshEquator
       type (vecSphere), intent (out), &
  -      dimension (size(vals, 1) - 1, size(vals, 2) - 1) :: res
  +      dimension (size(Atheta, 1) - 1, size(Atheta, 2) - 1 ) :: res
   
       integer :: iR, iPhi, nR, nPhi
       double precision :: r, dr, dPhi
  
positive.code:     nR = size(vals, 1)
      nPhi = size(vals, 2)
  
      checkGrid(meshEquator%nR, nR, 'r mesh for Ar')
      checkGrid(meshEquator%nPhi, nPhi, 'phi mesh for Ar')
  
      dr = calcDr(meshEquator)
      dPhi = calcDphi(meshEquator)
  
negative.code:                         ) / dr &
          !             + 1/r * (partial Aphi / partial phi)
                        + 1 / r * ( &
                            Aphi(iR, iPhi + 1) - Aphi(iR, iPhi) &
                          ) / dPhi
        end do
      end do
  