2013_08_29
This model has MULTPV =100000 around zone that has Sw=1, and boundaries i,j
Also model has new cut-off that is defined by relation: sowcr + swl == 1
PCW is counted by initial permeability to take into account GIS distribution of reserves (if we will want to change permeability)

DATAfile
INCLUDE
-'INCL/grid_arith_0823_2.inc' / -- grid_arith
+'GEO_0822/multpv.inc' / -- multpv is non-ones only swat==1 cells
+
+INCLUDE
+'INCL/grid_arith_0829_1.inc' / -- grid_arith

-PCW(,,1:40)=21.72973*28*0.318316*((poro/permx)^0.5)
+PCW(,,1:40)=21.72973*28*0.318316*((poro/(1.00005E-04*((PORO*100)^5.1)))^0.5)
-PCW(,,41:76)=54.59103*28*0.318316*((poro/permx)^0.5)
+PCW(,,41:76)=54.59103*28*0.318316*((poro/(4.02456E-02*EXP(44.17*PORO)))^0.5)

 INCLUDE
-  'INCL\event_0823_2.inc' /
+  'INCL\event_0829_1.inc' /

grid_arithmetic
ARR001(,,1:40) = MAX((42.9+11.835*(PORO*100)-1.5613*(PORO*100)^2+0.06093*(PORO*100)^3-0.000782*(PORO*100)^4)*0.01,0.07)
ARR001(,,41:76) = 0.01*(31040*(PORO*100)^-2.97)
ARR002(,,1:40) = 0.01*(-1.966*LOG(PERMX/1000) + 28.084)
ARR002(,,41:76) = 0.01*(-1.724 * LOG(PERMX/1000)+33.227)
ARR003 = (ARR001 + ARR002)
ARR003 = 1 - ARR003
NTG = NTG* (ARR003 >= 0.01) + 0 * (ARR003 < 0.01)
PORO = PORO * (ARR003 >= 0.01) + 0* (ARR003 < 0.01)
PERMX = PERMX * (ARR003 >= 0.01) + 0* (ARR003 < 0.01)
MULTPV = 100000* (MULTPV!=1)+1* (MULTPV==1)