anchor.diff:       max_temperr_lev = -1
        max_tempgrad_lev = -1
   
  +     period = 1.0
        model_P_name = "P"
        model_S_name = "S"
   
  
positive.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
  
  
negative.code: 
                   if ( inertial ) then
  
                     vx = (-2.0d0 * M_PI / period) * yy
                     vy = ( 2.0d0 * M_PI / period) * xx
                     vz = ZERO
  
                   else
  