anchor.diff:       ! 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)
  
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:      use bl_constants_module
       use model_parser_module, only: idens_model, itemp_model, ipres_model, ispec_model
       use initial_model_module, only: interpolate_3d_from_1d
       use rot_sources_module, only: cross_product, get_omega
  
       implicit none
  
  