anchor.diff:  
       use bl_constants_module, only: ZERO, HALF, TWO
       use meth_params_module, only: NVAR, URHO, UTEMP
  -    use prob_params_module, only: center, probhi
  +    use prob_params_module, only: center, probhi, dim, Symmetry, physbc_lo, physbc_hi
       use probdata_module, only: max_tagging_radius, &
                                  max_stellar_tagging_level, &
                                  max_temperature_tagging_level, &
  
positive.code:                    ! radius to ensure there aren't any sharp gradients in regions of
                     ! greater than ambient density.
  
                     r_P = ( (x-com_P(1))**2 + (y-com_P(2))**2 + (z-com_P(3))**2 )**HALF
                     r_S = ( (x-com_S(1))**2 + (y-com_S(2))**2 + (z-com_S(3))**2 )**HALF
  
                     if (r_P <= roche_tagging_factor * roche_rad_P) then
                        tag(i,j,k) = set
  
negative.code:     double precision :: problo(3),dx(3),time
      integer          :: level,set,clear
  
      integer          :: i, j, k, n
      double precision :: loc(3), r, r_P, r_S
  
      logical          :: outer_boundary_test(3)
  
      do k = lo(3), hi(3)
         loc(3) = problo(3) + (dble(k) + HALF) * dx(3)
  
         do j = lo(2), hi(2)
            loc(2) = problo(2) + (dble(j) + HALF) * dx(2)
  
            do i = lo(1), hi(1)
               loc(1) = problo(1) + (dble(i) + HALF) * dx(1)
  
               if (level < max_stellar_tagging_level) then
  
  