# Set mpirun, the full path to the flow executable, and simulator flags (except --output-dir)
flow = "flow --relaxed-max-pv-fraction=0 --enable-opm-rst-file=true --newton-min-iterations=1"

# Set the model parameters
regional_dims = [45000,15000,81] # Regional aquifer length, width, and depth [m]
regional_x_n = [3,5,5] # Variable array of x-refinement (Regional)
regional_y_n = [5,5,5] # Variable array of y-refinement (Regional)
regional_z_n = [1,1,1] # Variable array of z-refinement (Regional)
reference_x_n = [1,1,1,5,5,15,15,15,15,15,5,5,5,5,5] # Variable array of x-refinement (Reference)
reference_y_n = [25,25,25] # Variable array of y-refinement (Reference)
reference_z_n = [6,5,5,5,3,3,3,3,3] # Variable array of z-refinement (Reference)
site_location = [18000,5000,0,27000,10000,81] # Site xi, yi, zi, xf, yf, and zf box positions [m]
fault_regional = [10000,11000,0.01,10,22.5] # Regional fault x, and y positions [m], x and y multipliers for the trans, and height of the fault jump [m] 
fault_site = [[21583,5710],[24081,8233],[0,0]] # Site fault x, and y positions [m] (initial and final) and x and y multipliers for the trans
thickness = [9,9,9,9,9,9,9,9,9] # Thickness of the layers [m]
pressure = 2E2  # Pressure on the reservoir top [Bar]
temperature = [60,50] # Top and bottom temperatures [C]
rock_comp = 6.11423e-5 # Rock compressibility [1/Bar]
sensor_coords = [20000,8000,0] # Sensor position x, y, and z to assess the error over time w.r.t the reference solution [m]
z_xy = "(20-20*mt.sin((2*mt.pi*(x+y)/10000)))" # The function for the reservoir surface
hysteresis = ["Killough", "Both"] # Activate hysteresis: 1) Killough or Carlson and 2) Both, Pc, or Kr (HYSTOPT in EHYSTR keyword, see the OPM manual)
salinity = 2.92 # Add salinity (value [1e-3 kg-M/kg])

# Set the boundary conditions
# Use open/closed/porv for the Regional aquifer (if porv, enter the bottom, right, top, and left values (e.g, porv 1e8 1e7 1e6 1e5))
regional_bctype = ["open"]
# Use open/closed/porv/porvproj/flux/pres/pres2p/wells for the BC site (if porv; bottom, right, top, and left values (e.g, porv 1e4 1e3 1e2 1e1))
# For pres/flux BC's, add 'interp' to use linear interpolation in time
site_bctype = ["open"]

# Set the saturation functions
krw = "krw * ((sw - swi) / (1.0 - sni -swi)) ** nkrw"        # Wetting rel perm saturation function [-]
krn = "krn * ((1.0 - sw - sni) / (1.0 - sni - swi)) ** nkrn" # Non-wetting rel perm saturation function [-]
pcap = "pen * ((sw - swi) / (1.0 - swi)) ** (-(1.0 / npen))" # Capillary pressure saturation function [Bar]

# Properties sat functions: 1) swi [-], 2) sni [-], 3) krw [-], 4) krn [-], 5) pen [Bar], 6) nkrw [-], 7) nkrn [-],
# 8) npen [-], 9) threshold cP evaluation, and 10) npoints [-] (entry per layer, if hysteresis, additional entries per layer)
safu = [[0.12,0.10,1,1,3060e-5,2,2,2,1e-4,1000],
[0.12,0.10,1.0,1.0,3870.63e-5,2,2,2,1e-4,1000],
[0.12,0.10,1.0,1.0,3060.00e-5,2,2,2,1e-4,1000],
[0.12,0.10,1.0,1.0,3870.63e-5,2,2,2,1e-4,1000],
[0.32,0.10,1.0,1.0, 193531e-5,2,2,2,1e-4,1000],
[0.14,0.10,1.0,1.0,8654.99e-5,2,2,2,1e-4,1000],
[0.12,0.10,1.0,1.0,6120.00e-5,2,2,2,1e-4,1000],
[0.14,0.10,1.0,1.0,8654.99e-5,2,2,2,1e-4,1000],
[0.12,0.10,1.0,1.0,6120.00e-5,2,2,2,1e-4,1000],
[0.22,0.30,1.0,1.0,3060.00e-5,2,4,2,1e-4,1000],
[0.12,0.30,1.0,1.0,3870.63e-5,2,4,2,1e-4,1000],
[0.12,0.30,1.0,1.0,3060.00e-5,2,4,2,1e-4,1000],
[0.12,0.30,1.0,1.0,3870.63e-5,2,4,2,1e-4,1000],
[0.32,0.30,1.0,1.0, 193531e-5,2,4,2,1e-4,1000],
[0.14,0.30,1.0,1.0,8654.99e-5,2,4,2,1e-4,1000],
[0.12,0.30,1.0,1.0,6120.00e-5,2,4,2,1e-4,1000],
[0.14,0.30,1.0,1.0,8654.99e-5,2,4,2,1e-4,1000],
[0.12,0.30,1.0,1.0,6120.00e-5,2,4,2,1e-4,1000]]

# Properties rock: 1) Kxy [mD], 2) Kz [mD], and 3) phi [-] (entry per layer)
rock = [[1013.25,101.325,0.25],
[506.625,50.6625,0.20],
[1013.25,101.325,0.25],
[506.625,50.6625,0.20],
[0.10132,0.01013,0.10],
[101.324,10.1324,0.20],
[202.650,20.2650,0.20],
[101.324,10.1324,0.20],
[202.650,20.2650,0.20]]

# Wells position: 1) x, 2) y, 3) zi, and 4) zf positions [m] (entry per well)
well_coords = [[21180,7068,0,81],[24200,7800,15,65],[21718,7122,45,81],[14518,11377,0,50],[31679,8883,0,30],[28477,2732,0,81]]

# Define the injection values (entry per change in the schedule): 
# 1) injection time [d], 2) time step size to write results regional [d], 3) time step size to write results site/reference [d]
# 1) fluid (0 wetting, 1 non-wetting) well 0, 2) injection rates [kg/day] well 0, 3) fluid ... well n, injection, ...well n (as many as num of wells) 
# if 'wells' for site_bctype, then 1) bottom, 2) right, 3) top, and 4) left values (0(prod)/1(inj), pressure [Bar]))
# if --enable-tuning=1, then 1) for TUNING values as described in the OPM manual
inj = [[[365,73,73],[1,3e5,1,3e5,1,3e5,1,5e6,1,5e6,0,1e7]],[[365,73,73],[1,3e5,1,3e5,1,3e5,1,5e6,1,0,0,1e7]]]