anchor.diff:      interface
   
           subroutine PlmBar(p, lmax, z, csphase, cnorm, exitstatus)
  -            use ftypes
  +            integer, parameter :: dp = selected_real_kind(p=15)
               integer, intent(in) :: lmax
               real(dp), intent(out) :: p(:)
               real(dp), intent(in) :: z
  
positive.code:         end subroutine SHMagPowerSpectrum
  
          function SHMagPowerL(c, a, r, l)
              use ftypes
              real(dp) :: SHMagPowerL
              real(dp), intent(in) :: c(:,:,:)
              real(dp), intent(in) :: a, r
  
negative.code: 
          subroutine SHCrossPowerSpectrumDensityC(c1, c2, lmax, cspectra,&
                                                  exitstatus)
              integer, parameter :: dp = selected_real_kind(p=15)
              complex(dp), intent(in) :: c1(:,:,:), c2(:,:,:)
              integer, intent(in) :: lmax
              complex(dp), intent(out) :: cspectra(:)
  