Function setupSystem.m
...discrete differential operators, div and grad, the gravity force contribution, dz, function s = setupSystem(G, rock, bc, param) cf = G.cells.faces(:,1); nf = G.faces.num; nc = G.cells.num; Compute pore volumes s.pv = poreVolume(G, rock); s.poro = rock.poro; Compute the half, and then the full, transmissibilities...