create_spin_rpa.for:  500 format ('re = ', f4.2, '  vel = ', f4.2, '  mlt = ', f5.2,