Earlier ultrasonic measurements of C33 for Gd, which showed step changes at various temperatures and fields, are shown to be due to changes in the domain configuration. The transition fields calculated using domain phase theory are in good agreement with the experimental results.