anchor.diff:            uks%data(1)=fum%data(1)
           end function umks2uksbl3
   
  -        function umksh2uks(fum,r,x2,hin) result(uks)
  +        function umksh2uks(fum,r,x2) result(uks)
             ! Converts a Kerr-Schild spherical 4-velocity to a Boyer-Lindquist one.
             ! JAD 11/10/2008
             ! Do simple transformation first:
  
positive.code:           real, dimension(:), intent(in) :: r, x2
            type (four_vector), dimension(:), intent(in) :: fum
            type (four_vector), dimension(size(r)) :: uks
            real, intent(in), optional :: hin
            real :: hval
            real(kind=8), dimension(size(r)) :: ur, uth, uph
            real, dimension(size(r)) :: dthdx2
            write(6,*) 'read harm umksh2uks present h'
            if (present(hin)) then
               hval=hin
            else
               hval=0.3
            endif
            write(6,*) 'hval: ',present(hin), hval
            ur=dble(r)*fum%data(2)
            ! Compute partial derivatives:
  !          dthdx2=pi
  !          dthdx2=pi*(1.+(1.-hval)*cos(2.*pi*x2))
            dthdx2=pi/2.*(1.+(1.-hval)*cos(pi*(1.+x2)))
            uth=fum%data(3)*dble(dthdx2)
            uph=fum%data(4)
            uks=fum
  
negative.code: !        write(6,*) 'uniform phi'
  ! uniform in phi
          pd=(zphi-minph)/(uniqph(2)-uniqph(1))
          if(BL.eq.3) then
             td=abs(theta-calcthmksbl3(dble(uniqx2(lx2)),dble(zr)))/dth
          else
             td=abs(theta-uniqth(lx2))/dth
          end if
          rd=(zr-uniqr(lx1))/(uniqr(ux1)-uniqr(lx1))
  ! for now just do nearest neighbor
          rd(:)=1.; td(:)=1.; pd(:)=1.
  