TALK=T;RUN( 1, 1)
 
 ************************************************************
   Q1 created by VDI menu, Version 2014, Date 23/07/14
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************
  Echo DISPLAY / USE settings
  DISPLAY
 
   A box is filled with two immiscible fluids. The bottom fluid is
   denser than the upper. At time t=0, the interface between the
   two liquids is uniformly inclined, and is then let free to
   oscillate under the infulence of gravity.
 
   The vessel is 13m long, and 4m wide. The lower liquid is
   0.2m deep, the upper liquid is 0.2m deep. The original
   inclination of the interface is +/0.035m with respect to the
   average interface height.
 
   The top liquid contains two anode blocks. The upper fluid depth
   under the anodes is 0.04m
  ENDDIS
 ************************************************************
 IRUNN = 1 ;LIBREF = 498
 ************************************************************
  Group 1. Run Title
 TEXT(ESTER demo, Transient one-phase - Case 3)
 ************************************************************
  Group 2. Transience
 STEADY=F
    * Set overall time and no. of steps
 RSET(U,0.,400.,100)
    * Modify regions
 RSET(T,1,100,1.)
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,5,9,5)
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
    * Non-default variable names
 NAME(15)=HI ;NAME(147)=PRPS
 NAME(148)=NPOR ;NAME(149)=EPOR
 NAME(150) =CPOR
    * Solved variables list
 SOLVE(P1,U1,V1,W1)
    * Stored variables list
 STORE(CPOR,EPOR,NPOR,PRPS,HI,SURF)
    * Additional solver options
 SOLUTN(P1,Y,Y,Y,N,N,Y)
 
 ************************************************************
  Group 8. Terms & Devices
 U1AD =0. ;V1AD =0.
 W1AD = GRND
 ************************************************************
  Group 9. Properties
    * Domain material index is  51 signifying:
    * etal
 SETPRPS(1, 51,ESTR)
 ENUT =0.
 DRH1DP =1.0E-10
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(HI)=1.001E-10 ;FIINIT(PRPS)=-1.
 FIINIT(NPOR)=1. ;FIINIT(EPOR)=1.
 FIINIT(CPOR)=1.
   No PATCHes used for this Group
 
 
 INIADD = F
 ************************************************************
  Group 12. Convection and diffusion adjustments
   No PATCHes used for this Group
 ************************************************************
  Group 13. Boundary & Special Sources
 
   This patch is attached to object PREF1
 PATCH(PREF1, CELL, -1, 1, 0, 0, 0, 0, 1, 100)
 COVAL(PREF1, P1, 1., 0.)
 COVAL(PREF1, U1, 0., SAME)
 COVAL(PREF1, V1, 0., SAME)
 COVAL(PREF1, W1, 0., SAME)
 
   This patch is attached to object PREF2
 PATCH(PREF2, CELL, -1, 1, 0, 0, 0, 0, 1, 100)
 COVAL(PREF2, P1, 1., 0.)
 COVAL(PREF2, U1, 0., SAME)
 COVAL(PREF2, V1, 0., SAME)
 COVAL(PREF2, W1, 0., SAME)
 
 EGWF = T
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 60
 SARAH =1000.
 RESREF(P1)=1.0E-10 ;RESREF(U1)=1.0E-10
 RESREF(V1)=1.0E-10 ;RESREF(W1)=1.0E-10
 RESFAC =1.0E-04
 ************************************************************
  Group 16. Terminate Iterations
 ************************************************************
  Group 17. Relaxation
 RELAX(P1 ,LINRLX,0.9 )
 RELAX(U1 ,FALSDT,8. )
 RELAX(V1 ,FALSDT,8. )
 RELAX(W1 ,FALSDT,8. )
 ************************************************************
  Group 18. Limits
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 NAMGRD =ESTR
 PARSOL = F
 ISG62 = 1
 CSG8    ='ESTR'
 CSG10   ='ESTR'
 SPEDAT(SET,GXMONI,TRANSIENT,L,F)
 SPEDAT(SET,MATERIAL,51,L,T)
 SPEDAT(SET,MATERIAL,52,L,T)
 SPEDAT(SET,MATERIAL,100,L,T)
 SPEDAT(SET,MATERIAL,151,L,T)
 SPEDAT(SET,MATERIAL,198,L,T)
 SPEDAT(SET,MATERIAL,199,L,T)
 LG( 5) = T
 IG( 1) = 2
 IG( 2) = 4
 IG( 3) = 1
 IG( 4) = 1
 IG( 6) = 5
 RG( 1) =2300.
 RG( 2) =2100.
 RG( 3) =1.
 RG( 4) =4.1E+06
 RG( 5) =200.
 RG( 6) =2.33E+04
 RG( 7) =-9.81
 RG( 11) =1.
 RG( 12) =100.
 RG( 13) =1.
 RG( 14) =1.
 RG( 16) =-0.035
 RG( 17) =0.035
 RG( 18) =2.608696E-07
 RG( 19) =1.714286E-06
 RG( 31) =9.09E+04
 RG( 32) =999.869995
 RG( 33) =7800.
 RG( 34) =9.08E+05
 RG( 35) =2700.
 RG( 36) =1.98E+07
 RG( 37) =2700.
 RG( 38) =2.86E+07
 RG( 71) =7800.
 RG( 72) =2.17E+06
 RG( 73) =8950.
 RG( 74) =1.34E+07
 RG( 75) =7800.
 RG( 76) =6.45E+05
 RG( 77) =7800.
 RG( 78) =8.7E+05
 ************************************************************
  Group 20. Preliminary Printout
 ************************************************************
  Group 21. Print-out of Variables
 OUTPUT(HI  ,Y,N,Y,N,Y,Y)
 OUTPUT(PRPS,N,N,N,N,N,N)
 OUTPUT(NPOR,Y,N,N,N,N,N)
 OUTPUT(EPOR,Y,N,N,N,N,N)
 OUTPUT(CPOR,Y,N,N,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 2 ;IYMON = 2 ;IZMON = 1
 NPRMON = 100000
 NPRMNT = 1
 TSTSWP = -5
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000
 NTPRIN = 10 ;ISTPRF = 1 ;ISTPRL = 100000
 NPLT = 1
 ISWPRF = 1 ;ISWPRL = 100000
 IPROF = 3
 
 PATCH(PLOT1, PROFIL, 1, 1, 1, 1, 2, 2, 1, 100)
 PLOT(PLOT1, HI, 0., 0.)
 
 PATCH(PLOT2, PROFIL, 1, 1, 4, 4, 2, 2, 1, 100)
 PLOT(PLOT2, HI, 0., 0.)
 
 PATCH(PLOT3, PROFIL, 1, 1, 9, 9, 2, 2, 1, 100)
 PLOT(PLOT3, HI, 0., 0.)
 ************************************************************
  Group 24. Dumps For Restarts
 IDISPA = 10 ;IDISPB = 0 ;IDISPC = 0
 CSG1    ='c'
 
 GVIEW(P,0.57735,-0.57735,0.57735)
 GVIEW(UP,-0.408248,0.408248,0.816497)
 GVIEW(DEPTH,3674.429932)
 GVIEW(NEARPLANE,0.026)
 GVIEW(VDIS,6.531529)
 GVIEW(CENTRE,2.,6.5,0.2)
 
> DOM,    SIZE,        4.000000E+00, 1.300000E+01, 4.000000E-01
> DOM,    MONIT,       1.012500E+00, 1.057142E+00, 5.000001E-02
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    INCREMENT,   1.000000E-02, 1.000000E-02, 1.000000E-02
> GRID,   RSET_X_1,      1, 1.000000E+00
> GRID,   RSET_X_2,      1, 1.000000E+00
> GRID,   RSET_X_3,      1, 1.000000E+00
> GRID,   RSET_X_4,      1, 1.000000E+00
> GRID,   RSET_X_5,      1, 1.000000E+00
> GRID,   RSET_Y_1,      1, 1.000000E+00
> GRID,   RSET_Y_2,      7, 1.000000E+00
> GRID,   RSET_Y_3,      1, 1.000000E+00
> GRID,   RSET_Z_1,      2, 1.000000E+00
> GRID,   RSET_Z_2,      2, 1.000000E+00
> GRID,   RSET_Z_3,      1, 1.000000E+00
> DOM,    T_AMBIENT,   0.000000E+00
> DOM,    CARBON-BLOCKS,NO
> DOM,    ANODES,       YES
> DOM,    SIZE-X,        1.725000E+00, 2.500000E-01, 1.500001E-01
> DOM,    SIZE-Y,        1.270000E+01, 0.000000E+00, 1.499991E-01
> DOM,    SIZE-Z-1,      2.000000E-01, 1.600000E-01
> DOM,    SIZE-Z-2,      2.000000E-01, 4.000001E-02
> DOM,    NUM-X,             1,     1,     1
> DOM,    NUM-Y,             7,     0,     1
> DOM,    NUM-Z,             2,     2,     1,     1
> DOM,    NUM-ANODES,        2,     1
> DOM,    TAPPING-GAP,  NO
> DOM,    LORENTZ,      OFF
> DOM,    AUTOMESH,     ON
> DOM,    INTERFACE,    SLOPE_IN_Y
 
> OBJ,    NAME,        INIELEC
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 2.000000E-01
> OBJ,    SIZE,        4.000000E+00, 1.300000E+01, 2.000000E-01
> OBJ,    GEOMETRY,    default
> OBJ,    VISIBLE,     NO
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    COLOR-MODE,  USER
> OBJ,    OPAQUE,      50
> OBJ,    MATERIAL,    52,Bath
> OBJ,    TIME_LIMITS, ALWAYS_ACTIVE
 
> OBJ,    NAME,        PREF1
> OBJ,    POSITION,    1.500001E-01, 3.778571E+00, 0.000000E+00
> OBJ,    SIZE,        1.725000E+00, 1.814286E+00, 1.000000E-01
> OBJ,    GEOMETRY,    default
> OBJ,    TYPE,        CELLTYPE
> OBJ,    COLOR-MODE,  USER
> OBJ,    COLOR-VAL,   49
> OBJ,    OPAQUE,      50
> OBJ,    PATCHES,     PREF1
 
> OBJ,    NAME,        PREF2
> OBJ,    POSITION,    1.500001E-01, 3.778571E+00, 2.000000E-01
> OBJ,    SIZE,        1.725000E+00, 1.814286E+00, 2.000003E-02
> OBJ,    GEOMETRY,    default
> OBJ,    TYPE,        CELLTYPE
> OBJ,    COLOR-MODE,  USER
> OBJ,    COLOR-VAL,   49
> OBJ,    OPAQUE,      50
> OBJ,    PATCHES,     PREF2
 
> OBJ,    NAME,        ANOD1
> OBJ,    POSITION,    1.500001E-01, 1.499991E-01, 2.400000E-01
> OBJ,    SIZE,        1.725000E+00, 1.270000E+01, 1.600000E-01
> OBJ,    GEOMETRY,    default
> OBJ,    TYPE,        ANODE
> OBJ,    COLOR-MODE,  USER
> OBJ,    COLOR-VAL,   238
> OBJ,    OPAQUE,      50
 
> OBJ,    NAME,        ANOD2
> OBJ,    POSITION,    2.125000E+00, 1.499991E-01, 2.400000E-01
> OBJ,    SIZE,        1.725000E+00, 1.270000E+01, 1.600000E-01
> OBJ,    GEOMETRY,    default
> OBJ,    TYPE,        ANODE
> OBJ,    COLOR-MODE,  USER
> OBJ,    COLOR-VAL,   238
> OBJ,    OPAQUE,      50
 
> OBJ,    NAME,        INTERFAC
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 2.000000E-01
> OBJ,    SIZE,        4.000000E+00, 1.300000E+01, 0.000000E+00
> OBJ,    GEOMETRY,    cube
> OBJ,    TYPE,        NULL
> OBJ,    COLOR-MODE,  DEFAULT2
> OBJ,    OPAQUE,      50
 
> OBJ,    NAME,        FREESURF
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 4.000000E-01
> OBJ,    SIZE,        4.000000E+00, 1.300000E+01, 0.000000E+00
> OBJ,    GEOMETRY,    cube
> OBJ,    GRID,        N,N,N
> OBJ,    TYPE,        NULL
> OBJ,    COLOR-MODE,  DEFAULT2
> OBJ,    OPAQUE,      50
 ************************************************************
 MENSAV(S,MENU,ESTR,ESTR)
STOP