Non-advection phase

The non-advaction phase 
is solved as

    

where  is F calculated from
the quantities at .

This is of first-order accuracy in time.