version #3

New version including mass balance and top flux boundary condition
changes