Revert rtvd to calculating 2 steps inside relaxing_tvd

This is a controversial change. Nevertheless I'd like to consider it. 

Disadvantages:
 - problems with external gravitational potential as it will certainly screw *outd* and such

Advantages:
 - less communication
 - lesser memory footprint
 - make certain algorithms (FARGO I'm looking at you) easier