******************************************************************************
** STPYS001.DAT:       STARS + PyScript/PyControl Communication             **
**                                                                          **
******************************************************************************
**************************************************************************************
**                                                                                  **
** FILE:   STPYS001.DAT                                                             **
** MODEL:  11x10x10                                                                 **
**         STARS + PYSCRIPT/PYCONTROL COMMUNICATION                                 **
**                                                                                  **
**************************************************************************************
**                                                                                  **
** This dataset, together with the Python script stpyc001.py, can be loaded into    **
** the PyControl application to test the coupling between STARS, OUTBOARD, and      **
** PyControl. It is useful when a simulation run needs to be controlled externally  **
** via a Python script, allowing results and conditions obtained from an external   **
** application to be incorporated.                                                  **
**                                                                                  **
**************************************************************************************
** CONTACT:    CMG,   (403)531-1300;   282-6495 (fax);support@cmgl.ca    (Internet) **
**************************************************************************************

RESULTS SIMULATOR STARS

*****************************************************************************
** I/O control section
*****************************************************************************

*INTERRUPT *STOP

*TITLE1 'STARS PyControl Template No. 1'
*TITLE2 'SAGD Production - Annulus Only'

*INUNIT *SI
*DIM *MDICLU  20000
     *MDLU  400000

***WRST time
*WPRN *GRID *TIME

*WPRN *ITER 1
*OUTPRN *WELL *ALL
*OUTPRN *GRID *TEMP *SG *PRES *SO
*OUTPRN *ITER *NEWTON
*OUTPRN *RES *ALL

*PRNTORIEN 2  0

*WSRF *GRID 1
*OUTSRF *GRID *TEMP *PRES *SG *SO *SW *X *Y *VISO

*****************************************************************************
** Grid data
*****************************************************************************
GRID CART 11 10 10

KDIR UP

DI CON 2
DJ JVAR
300. 100. 100. 100. 100. 100. 100. 100. 100. 300.
DK KVAR
1. 1. 1. 1. 1. 2. 2. 2. 2. 2.

VAMOD 2 0.5 1 0.5 0.5  ** Trimmed side blocks
VAMOD 3 0.5 0.5 0.5 1  ** Trimmed DW blocks like *JDIR hygrid:  ai = radial, aj = axial
VATYPE IVAR
           2            1            1            1            1            1
           1            1            1            1            2
DEPTH TOP 1 1 10 280.0

POR CON 0.34

PERMI CON 2200
PERMJ EQUALSI
PERMK EQUALSI * 0.8

END-GRID

PRPOR    1000
ROCKTYPE 1
   CPOR 7.E-06
   ROCKCP 2.35E+06
   THCONR 6.6E+05
   THCONW 5.34E+04
   THCONO 1.15E+04
   THCONG 140
   HLOSSPROP OVERBUR  2.35E+06 1.728E+05
             UNDERBUR 2.35E+06 1.728E+05

THTYPE  CON 1

*****************************************************************************
** Fluid model, component properties
*****************************************************************************
*MODEL 3 3 3

*COMPNAME   'WATER' 'BITUMEN' 'CH4'
*KV1        0.0      0.00     3.191E+04
*KV4        0.0      0.00    -3.307E+01
*KV5        0.0      0.00    -2.771E+01
*CMM        0.0      0.48     0.01604
*PCRIT      0.0      1.36E+3  4.6E+3
*TCRIT      0.0      624.65  -82.55
*MOLDEN     0.0      2045.8    3.75E+4
*CP         0.0      2*5.5e-7
*CT1        0.0      8.000E-04 8.000E-04

*VISCTABLE
**  Temp    WATER BITUMEN        CH4
        10  0.0   5000000        200
        20  0.0   371535.2291    93.32543008
        30  0.0   85113.80382    46.08
        40  0.0   23442.28815    26.30267992
        50  0.0   8028           12.58925412
        60  0.0   3073           8.16
        70  0.0   1415.4         5.623413252
        80  0.0   680.6          4
        100 0.0   212.6          2.454708916
        125 0.0   70.94          1.659586907
        150 0.0   30.92          1.230268771
        175 0.0   16.35          0.912010839
        200 0.0   9.96           0.831763771
        225 0.0   6.787          0.758577575
        250 0.0   4.999          0.72443596
        275 0.0   3.73           0.68
        300 0.0   3              0.616595002
        350 0.0   2.570395783    0.588843655
        400 0.0   2.398832919    0.562341325

*****************************************************************************
** Rock-fluid data
*****************************************************************************
*ROCKFLUID

*RPT 1
*SWT
**  Sw        Krw      Krow
** --------------------------------
0.190000000  0.000000000  1.000000000
0.200000000  0.000600000  0.960000000
0.250000000  0.003200000  0.850000000
0.300000000  0.008000000  0.720000000
0.350000000  0.013000000  0.600000000
0.400000000  0.025400000  0.470000000
0.450000000  0.044000000  0.350000000
0.500000000  0.069800000  0.240000000
0.550000000  0.104000000  0.165000000
0.600000000  0.148000000  0.100000000
0.650000000  0.204000000  0.060000000
0.700000000  0.271000000  0.030000000
0.750000000  0.352000000  0.015000000
0.800000000  0.447000000  0.005000000
0.850000000  0.559000000  0.000000000
0.900000000  0.687000000  0.000000000
0.950000000  0.834000000  0.000000000
1.000000000  1.000000000  0.000000000

*SLT
**  Sl        Krg      Krog
** --------------------------------
0.190000000  1.000000000  0.000000000
0.200000000  0.950000000  0.000600000
0.250000000  0.840000000  0.003200000
0.300000000  0.720000000  0.008000000
0.350000000  0.600000000  0.013000000
0.400000000  0.470000000  0.025400000
0.450000000  0.350000000  0.044000000
0.500000000  0.240000000  0.069800000
0.550000000  0.165000000  0.104000000
0.600000000  0.093000000  0.148000000
0.650000000  0.075000000  0.204000000
0.700000000  0.045000000  0.271000000
0.750000000  0.027000000  0.352000000
0.800000000  0.020000000  0.447000000
0.850000000  0.010000000  0.559000000
0.900000000  0.005000000  0.687000000
0.950000000  0.000000000  0.834000000
1.000000000  0.000000000  1.000000000

*****************************************************************************
** Initial conditions
*****************************************************************************
*INITIAL

*VERTICAL *DEPTH_AVE
*REFPRES 1340
*REFDEPTH 280

TEMP CON 10

MFRAC_OIL 'BITUMEN' CON 0.95
MFRAC_OIL 'CH4'     CON 0.05

*****************************************************************************
** Numerical control
*****************************************************************************
*NUMERICAL

*DTMAX 15

*NORM *PRESS 500  *TEMP 60 SATUR 0.2  Y 0.2  X  0.2  W  0.2

*PIVOT *ON
*SDEGREE 2
*ITERMAX 40
*AIM STAB

*****************************************************************************
** Recurrent data
*****************************************************************************
RUN

DATE  2024 1 1

   *DTWELL 0.01

   *WELL 1 'injtb'   *FRAC 0.5
   *WELL 2 'injan'   *FRAC 0.5
   *WELL 3 'prdtb'   *FRAC 0.5
   *WELL 4 'prdan'   *FRAC 0.5

*INJECTOR *MOBWEIGHT 'injtb'
   *INCOMP *WATER  1.0  0.0  0.0
   *TINJW 224.0
   *QUAL 0.9
   *OPERATE *MAX *BHP 5500.
   *OPERATE *MAX *STW 400.
   *GEOMETRY *J 0.14 0.229 1.0  0.0
   *PERF *GEO 'injtb'
          1 2:9 5    1.0

*PRODUCER 'injan'
   *OPERATE *MAX *STL 500.
   *OPERATE *MIN *BHP 2480.
   *PERF *GEO 'injan'
          1 2:9 5    1.0

*INJECTOR *MOBWEIGHT 'prdtb'
   *INCOMP *WATER  1.0  0.0  0.0
   *TINJW 224.0
   *QUAL 0.9
   *OPERATE *MAX *BHP 5500.
   *OPERATE *MAX *STW 400.
   *PERF *GEO 'prdtb'
          1 2:9  2    1.0

*PRODUCER 'prdan'
   *OPERATE *MAX *STL 500.
   *OPERATE *MIN *BHP 2480.
   *PERF *GEO 'prdan'
          1 2:9 2   1.0

*TIME 30.0
	*DTWELL  0.01

*TIME 40.0
*TIME 60.0
*TIME 75.0

*INJECTOR *MOBWEIGHT 'injtb'
   *INCOMP *WATER  1.0  0.0  0.0
   *TINJW 224.0
   *QUAL 0.9
   *OPERATE *MAX *BHP 5500.
   *OPERATE *MAX *STW 400.

*TIME 100.0
*TIME 110.0
*TIME 120.0
*TIME 130.0
*TIME 140.0
*TIME 150.0

*STOP
