Numerical modeling of the productivity of vertical to shallowly dipping fractured zones in crystalline rocks