anchor.diff:               l2_norm_resid = l2_norm_resid + dV * (state(i,j,k,URHO) - old_rho)**2
                l2_norm_source = l2_norm_source + dV * old_rho**2
   
  -             kin_eng = kin_eng + HALF * omega**2 * state(i,j,k,URHO) * (r(1)**2 + r(2)**2) * dV
  +             kin_eng = kin_eng + HALF * sum(cross_product(omega, r)**2) * state(i,j,k,URHO) * dV
   
                pot_eng = pot_eng + HALF * state(i,j,k,URHO) * phi(i,j,k) * dV
   
  
positive.code:     use prob_params_module, only: center
      use probdata_module, only: star_axis, get_ambient, h_max_P, h_max_S, enthalpy_min
      use eos_module
  
      implicit none
      
  
negative.code:           do i = rloc(1,3), rloc(1,3) + 1
               if (i .ge. lo(1) .and. j .ge. lo(2) .and. k .ge. lo(3) .and. &
                   i .le. hi(1) .and. j .le. hi(2) .and. k .le. hi(3)) then
  
                  omegasq = omegasq + c(i-rloc(1,3),j-rloc(2,3),k-rloc(3,3),3) &
                          * TWO * phi(i,j,k) / (sin(theta3)**2 * sum(d_vector(:,3)**2) - sin(theta2)**2 * sum(d_vector(:,2)**2))
  
               endif
            enddo
         enddo
  