coRotating changes Omega too
Setting coRotating referential changes Omega too. As a result theoreticalProfile's vtheta shouldn't be updated too ? since it is used in boundaries conditions, it seems that we should used the "corrected" values.
Which raises a question for 2D theoretical profiles. Right now, if modification of omega is detected, we recompute vtheta from scratch instead of just applying a purely geometrical transformation. Is that correct ?