anchor.diff:                  endif
   
                ! If we're in the inertial reference frame, and we want to provide an
  -             ! initial orbital kick, use counter-clockwise rigid body rotation.
  +             ! initial orbital kick, use rigid body rotation with velocity omega x r.
   
                else if ((do_rotation .ne. 1) .and. (.not. no_orbital_kick) .and. (.not. single_star)) then
   
  -                ! x velocity is -omega * r * sin(theta) == -omega * y
  -                state(i,j,k,UMX+star_axis-1) = state(i,j,k,UMX+star_axis-1) &
  -                                             + state(i,j,k,URHO) * (-TWO * M_PI / rot_period) * loc(initial_motion_dir)
  -
  -                ! y velocity is +omega * r * cos(theta) == +omega * x
  -                state(i,j,k,UMX+initial_motion_dir-1) = state(i,j,k,UMX+initial_motion_dir-1) &
  -                                                      + state(i,j,k,URHO) * ( TWO * M_PI / rot_period) * loc(star_axis)
  +                state(i,j,k,UMX:UMZ) = state(i,j,k,UMX:UMZ) + state(i,j,k,URHO) * cross_product(omega, loc)
   
                endif
   
  
positive.code:      double precision :: xlo(3), xhi(3), time, delta(3)
       double precision :: state(state_l1:state_h1,state_l2:state_h2,state_l3:state_h3,NVAR)
  
       double precision :: loc(3)
       double precision :: dist_P(3), dist_S(3)
  
       type (eos_t) :: zone_state, ambient_state
  
negative.code:      ! or if we are in an ambient zone.
  
       call get_ambient(ambient_state)
  
       omega = get_omega()
  
       !$OMP PARALLEL DO PRIVATE(i, j, k, loc) &
       !$OMP PRIVATE(dist_P, dist_S, zone_state)
  