Fixed a small bug: Dirichlet boundary cobndition was not initialized with .false.
Moved porosity from vertex to cell data. Inactive cells with zero porosity are now allowed. Boundary conditions still do not work properly - the outflow always gets oversaturated