Upscaling For Reservoir Simulation Using Multi-Layer Directional 2 Phase Pseudofunctions