import numpy as np

M = 200 #[kg]
M_H = 1 #[g/mol]
f_H = 1.75 #[-]
K_H = 0.78e4 #[atm^0.5]
G_m = 4.5e-6 #[kmol/s]
p_in = 1 #[atm]
k_ts = 2.3e-3 #[m/s]
A_s = 0.17 #[m^2]
rho = 2350 #[kg/m^3]
k_t = 2.3e-4 #[m/s]
G_v = 22.4 * G_m #[m^3/s]
d_b = 0.004 #[m]
h = 0.5 #[m]

tau = h/( 0.32 * np.power(0.004/d_b, 1/8) )
A = (6 * G_v * tau)/d_b #[m^2]

C0 = -(100 * M_H / M)
C1 = C0 * ((2 * np.power(f_H * K_H, 2) * G_m) / p_in)
C2 = C0 * ((k_ts * A_s * rho) / (100 * M_H))

K = (rho * k_t * A * p_in) / (400 * M_H * G_m * np.power(f_H * K_H, 2) )

scipy_integrator_method = 'dop853'