anchor.diff:       use interpolate_module
        use eos_module
        use meth_params_module, only : NVAR, URHO, UMX, UMY, UMZ, UTEMP, &
  -          UEDEN, UEINT, UFS
  +          UEDEN, UEINT, UFS, rot_period
        use network, only : nspec
        use bl_constants_module
        use model_parser_module, only: idens_model, itemp_model, ipres_model, ispec_model
  
positive.code:      use model_parser_module
       use bl_constants_module
       use fundamental_constants_module
       use meth_params_module, only: small_temp, small_pres, small_dens
       use eos_module
       use com, only: mass_p, mass_s, com_loc_p, com_loc_s
  
  
negative.code: 
               ! x velocity is -omega * r * sin(theta) == -omega * y
  
               state(i,j,k,UMX) = state(i,j,k,URHO) * (-2.0d0 * M_PI / rot_period) * yy
             
               ! y velocity is +omega * r * cos(theta) == +omega * x
  
               state(i,j,k,UMY) = state(i,j,k,URHO) * ( 2.0d0 * M_PI / rot_period) * xx
  
               ! Add corresponding kinetic energy
  
  