anchor.diff:  
        namelist /fortin/ &
             model_P_name, model_S_name, &
  -          period, &
             nsub, &
             inertial, &
             interp_temp, &
  
positive.code: 
                   if ( inertial ) then
  
                     vx = (-2.0d0 * M_PI / period) * yy
                     vy = ( 2.0d0 * M_PI / period) * xx
                     vz = ZERO
  
                   else
  
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
  
  