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: 
                   if ( inertial ) then
  
                     vx = (-2.0d0 * M_PI / rot_period) * yy
                     vy = ( 2.0d0 * M_PI / rot_period) * xx
                     vz = ZERO
  
                   else
  