A numerical model for {CO}2 gas migration in a fault zone