import numpy as np

PUE = 2

I_E = R * PUE *eta_E * ( I_E_S + I_E_R )

I_M_TM_BU = R * n_M_R * c_ON * ( I_ON + n_M_OA  * I_OA)

I_M_TM = np.random.uniform(I_M_TM_BU,I_M_TM_EMP)

I_M_IP = PUE * eta_M * (R * n_M_R * I_M_R)
I_M_O = PUE * eta_M * I_M_TM
I_M = I_M_IP + I_M_O


I_C_TM_BU = R * n_C_R * ( I_ON + n_C_OA  * I_OA)

I_C_TM = np.random.uniform(I_C_TM_BU,I_C_TM_EMP)

I_C_IP = PUE * eta_C * ( R * n_C_R * I_C_R)
I_C_O = PUE * eta_C *  I_C_TM
I_C= I_C_IP + I_C_O

I_S= s_T/s_SA * I_OSA + R * PUE * 2 I_ST + s_T / C_U