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