anchor.diff:  
           subroutine SHBiasKMask(tapers, lwin, k, incspectra, ldata, &
                              outcspectra, taper_wt, save_cg, exitstatus)
  -            use ftypes
  +            integer, parameter :: dp = selected_real_kind(p=15)
               real(dp), intent(in) :: tapers(:,:), incspectra(:)
               real(dp), intent(out) :: outcspectra(:)
               integer, intent(in) :: lwin, ldata, k
  
positive.code: 
          subroutine SHExpandGLQC(cilm, lmax, gridglq, w, plx, zero, norm, &
                                  csphase, lmax_calc, exitstatus)
              use ftypes
              real(dp), intent(in) :: w(:)
              complex(dp), intent(in) :: gridglq(:,:)
              real(dp), intent(in), optional :: plx(:,:), zero(:)
  
negative.code: 
          subroutine CilmPlusRhoH(cilm, gridin, lmax, nmax, mass, d, rho, &
                                  gridtype, w, zero, plx, n, dref, exitstatus)
              integer, parameter :: dp = selected_real_kind(p=15)
              real(dp), intent(in) :: gridin(:,:), mass, rho(:,:)
              real(dp), intent(in), optional :: w(:), zero(:), plx(:,:), dref
              real(dp), intent(out) :: cilm(:,:,:), d
  