anchor.diff:  
           subroutine MakeCircleCoord(coord, lat, lon, theta0, cinterval, cnum, &
                                      exitstatus)
  -            use ftypes
  +            integer, parameter :: dp = selected_real_kind(p=15)
               real(dp), intent(in) :: lat, lon, theta0
               real(dp), intent(out) :: coord(:,:)
               real(dp), intent(in), optional :: cinterval
  
positive.code:         end subroutine PlmBar_d1
  
          subroutine PlBar(p, lmax, z, exitstatus)
              use ftypes
              integer, intent(in) :: lmax
              real(dp), intent(out) :: p(:)
              real(dp), intent(in) :: z
  
negative.code: 
          subroutine SHCrossPowerSpectrumDensity(c1, c2, lmax, cspectra, &
                                                 exitstatus)
              integer, parameter :: dp = selected_real_kind(p=15)
              real(dp), intent(in) :: c1(:,:,:), c2(:,:,:)
              integer, intent(in) :: lmax
              real(dp), intent(out) :: cspectra(:)
  