Ejemplo n.º 1
0
 def recount_property(cell_container, k, i, j):
     cell = cell_container.get_cell(k, i, j)
     cell_state_n_plus = cell.get_cell_state_n_plus()
     for component in Layer.components:
         component_index = Layer.components.index(component)
         ro = component.count_ro(
             cell_state_n_plus.get_components_pressure()[component_index])
         k_r = component.count_k_r(
             cell_state_n_plus.get_components_saturation()[component_index])
         cell_state_n_plus.set_ro(ro, component_index)
         cell_state_n_plus.set_k_r(k_r, component_index)
     fi = Layer.count_fi(cell_state_n_plus.get_pressure_oil())
     cell_state_n_plus.set_fi(fi)
     if cell.has_well:
         cell.well.count_well_index(cell, (k, i, j))