anchor.diff:          end subroutine PLegendreA_d1
   
           subroutine PLegendre(p,lmax,z, 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: 
          subroutine MakeGridDHC(griddh, n, cilm, lmax, norm, sampling, &
                                 csphase, lmax_calc, exitstatus)
              use ftypes
              complex(dp), intent(in) :: cilm(:,:,:)
              complex(dp), intent(out) :: griddh(:,:)
              integer, intent(in) :: lmax
  
negative.code:         end subroutine SHPowerSpectrumDensityC
  
          subroutine SHCrossPowerSpectrumC(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(:)
  