anchor.diff:                                          lo, hi, dx, time, Qtt) bind(C)
   
     use bl_constants_module, only: ZERO, THIRD, HALF, ONE, TWO
  -  use prob_params_module, only: center, problo
  +  use prob_params_module, only: center, problo, probhi, physbc_lo, physbc_hi, Symmetry
     use probdata_module, only: inertial_rotation, inertial_velocity
  -  
  +  use castro_util_module, only: position
  +
     implicit none
   
     integer         , intent(in   ) :: r_lo(3), r_hi(3)
  
positive.code:   double precision, intent(inout) :: h_plus_1, h_cross_1, h_plus_2, h_cross_2, h_plus_3, h_cross_3
    double precision, intent(in   ) :: Qtt(3,3)
    double precision, intent(in   ) :: time
    
    integer          :: i, j, k, l, dir
    double precision :: h(3,3), proj(3,3,3,3), delta(3,3), n(3), r
    double precision :: dist(3)
    
    ! Standard Kronecker delta.
  
    delta(:,:) = ZERO
  
negative.code:   use prob_params_module, only: problo
  
    implicit none
  
    integer          :: r_lo(3), r_hi(3)
    integer          :: v_lo(3), v_hi(3)
    integer          :: lo(3), hi(3)
  