TALK=F;RUN( 1, 1)
 
 ************************************************************
   Q1 created by VDI menu, Version 2007, Date 26/09/07
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************
  Echo DISPLAY / USE settings
  DISPLAY
  Problem: The quarter of the square channel.
           Inner pressure = 4.0e8, outer pressure = 0
   1. 2D 'PLANE-STRAIN' : z-surface fixed
   2. Comparison - a calculation with code ELCUT (www.tor.ru)
  ENDDIS
  PHOTON USE
  p;;;;
 
 
  set prop off
  msg
  msg V1 contours
  msg
  cont V1 z 1 x 1 20 y 1 20 fil;.0001
  pause
 
  cl
  msg
  msg U1 contours
  msg
  cont U1 z 1 x 1 20 y 2 21 fil;.0001
  pause
 
  cl
  msg
  msg EpsX contours
  msg
  cont EPSX z 1 x 1 20 y 1 20 fil;.0001
  pause
 
  cl
  msg
  msg EpsY contours
  msg
  cont EPSY z 1 x 1 20 y 1 20 fil;.0001
  pause
 
  cl
  msg
  msg EpsZ contours
  msg
  cont EPSZ z 1 x 1 20 y 1 20 fil;.0001
  pause
 
  cl
  msg
  msg STRX contours
  msg
  cont STRX z 1 x 1 20 y 2 20 fil;.0001
  pause
 
  cl
  msg
  msg STRY contours
  msg
  cont STRY z 1 x 1 20 y 2 20 fil;.0001
  pause
 
  cl
  msg
  msg STXY contours
  msg
  cont STXY z 1 x 1 20 y 2 20 fil;.0001
  pause
 
  cl
  msg Displacements fields
  gr ou z 1
  cont DELT z 1 x 1 20 y 2 20 fil;.0001
  vec z 1 x 1 20 y 1 20 col 0
  pause
 
 
  cl
  AUTOPLOT
  file
  phi 5
  EDZZ\S221_1.txt 2
  EDZZ\S221_2.txt 2
 
  data 1
  STRX X 14 Y 2 21
  col1 1
  data 2
  y SXX
  colA 2
  scale
  msg STRX(Y) for X=40 mm (WHITE-Phoenics, YELLOW - Elcut) ...
  plot 1 2
  pause
 
  clear
  data 1
  STRY X 14 Y 2 21
  col1 1
  data 2
  y SYY
  colA 2
  scale
  msg STRY(Y) for X=40 mm (WHITE-Phoenics, YELLOW - Elcut) ...
  plot 1 2
  pause
 
  clear
  data 1
  STXY X 14 Y 2 21
  col1 1
  data 2
  y TXY
  colA 2
  scale
  msg STXY(Y) for X=40 mm (WHITE-Phoenics, YELLOW - Elcut) ...
  plot 1 2
  pause
 
  clear
  data 1
  U1 X 13 Y 2 21
  col1 1
  data 3
  y dx
  colA 2
  scale
  msg U1(Y) for X=40 mm (WHITE-Phoenics, YELLOW - Elcut) ...
  plot 1 2
  pause
 
  clear
  data 1
  V1 X 14 Y 2 21
  col1 1
  data 3
  y dy
  colA 2
  scale
  msg V1(Y) for X=40 mm (WHITE-Phoenics, YELLOW - Elcut) ...
  plot 1 2
  pause
 
  clear
  data 1
  EPSX X 14 Y 2 21
  col1 1
  data 3
  y exx
  colA 2
  scale
  msg EpsX(Y) for X=40 mm (WHITE-Phoenics, YELLOW - Elcut) ...
  plot 1 2
  pause
 
  clear
  data 1
  EPSY X 14 Y 2 21
  col1 1
  data 3
  y eyy
  colA 2
  scale
  msg EpsY(Y) for X=40 mm (WHITE-Phoenics, YELLOW - Elcut) ...
  plot 1 2
  pause
 
 
  ENDUSE
 ************************************************************
 IRUNN   =         1 ;LIBREF =         221
 ************************************************************
  Group 1. Run Title
 TEXT(square channel;2D[xy]                   )
 ************************************************************
  Group 2. Transience
 STEADY  =    T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,21,21,1,1.000000E-05)
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Non-default variable names
 NAME(140) =DELY ; NAME(141) =DELX
 NAME(142) =DELT ; NAME(143) =STXY
 NAME(144) =EPSZ ; NAME(145) =EPSX
 NAME(146) =EPSY ; NAME(147) =STRZ
 NAME(148) =STRY ; NAME(149) =STRX
 NAME(150) =PRPS
    * Solved variables list
 SOLVE(P1  ,U1  ,V1  )
    * Stored variables list
 STORE(PRPS,STRX,STRY,STRZ,EPSY,EPSX,EPSZ,STXY)
 STORE(DELT,DELX,DELY)
 
 ************************************************************
  Echo InForm settings for Group  7
  INFORM7BEGIN

real(PIN,POUT,ZWW)
PIN=4.0e8  ! H/m^2
POUT = 0.0

INTEGER(CASENO)       ! start of menu
     caseno 1 : default settings
     caseno 3 : Poisson's ratio = 0.0 ; should affect result
     caseno 5 : zwlast = 1.0 m ; unlikely to affect result
     caseno 7 : focus print-out on channel corner
     caseno 9 : outside fluid pressure = inside pressure
    
    *************** 
caseno=1
ZWW = 0.001
    ***************  

if(caseno.eq.9) then
 POUT=PIN
endif 

  
  cell-centre displacements, for comparison with ELCUT
(STORED VAR Delt IS SQRT(U1^2+V1^2) with imat>100)
 
  INFORM7END
 ************************************************************
  Group 8. Terms & Devices
 ************************************************************
  Group 9. Properties
    * Domain material index is 160 signifying:
    * NOTSET
 SETPRPS(1,160,Q1)
 ENUT    = 0.000000E+00
 DRH1DP  = 5.000000E-12
 DVO1DT  = 1.000000E-05
   *  List of user-defined materials to be read by EARTH
    MATFLG=T;IMAT=1
   *  Name
   *Ind. Dens.  Viscos.  Spec.heat  Conduct.  Expans.  Compr.
   *  
    160 7800.0 0.3 473.0 43.0 1.0E-5 0.5E-11

 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(P1  ) =  0.000000E+00 ;FIINIT(U1  ) =  0.000000E+00
 FIINIT(V1  ) =  0.000000E+00 ;FIINIT(DELY) =  1.001000E-10
 FIINIT(DELX) =  1.001000E-10 ;FIINIT(DELT) =  0.000000E+00
 FIINIT(STXY) =  1.001000E-10 ;FIINIT(EPSZ) =  1.001000E-10
 FIINIT(EPSX) =  1.001000E-10 ;FIINIT(EPSY) =  1.001000E-10
 FIINIT(STRZ) =  1.001000E-10 ;FIINIT(STRY) =  1.001000E-10
 FIINIT(STRX) =  1.001000E-10 ;FIINIT(PRPS) =  1.600000E+02
   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
 
 PATCH(UP,NWALL,4,0,0,0,0,0,1,1)
 COVAL(UP,V1,1,0.0)

 PATCH(LEFT,WWALL,5,0,0,0,0,0,1,1)
 COVAL(LEFT,U1,1,0.0)

 
 PATCH (IN      ,CELL,6,0,0,0,0,0,1,1)
 COVAL (IN      ,P1  , FIXVAL      , :PIN:)
 
 PATCH (UIN     ,CELL  ,7,0,0,0,0,0,1,1)
 COVAL (UIN     ,U1  , FIXVAL      , 0.000000E+00)
 
 PATCH (VIN     ,CELL,8,0,0,0,0,0,1,1)
 COVAL (VIN     ,V1  , FIXVAL      , 0.000000E+00)
 
 PATCH (OUT1    ,CELL,9,0,0,0,0,0,1,1)
 COVAL (OUT1    ,P1  , FIXVAL      , :POUT:)
 
 PATCH (OUT2    ,CELL  ,10,1,0,0,0,0,1,1)
 COVAL (OUT2    ,P1  , FIXVAL      , :POUT:)
 
 EGWF    =    T
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =       300
 RESREF(V1  ) =  0.000000E+00
 RESFAC  = 1.000000E-07
 ************************************************************
  Group 16. Terminate Iterations
 LITER (P1  ) =        2 ;LITER (U1  ) =       50
 LITER (V1  ) =       50
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 1.000000E+00)
 ************************************************************
  Group 18. Limits
 VARMAX(U1  ) = 1.000000E+06 ;VARMIN(U1  ) =-1.000000E+06
 VARMAX(V1  ) = 1.000000E+06 ;VARMIN(V1  ) =-1.000000E+06
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 USEGRD  =    T  ;USEGRX =    T
 STRA    =    T
 ASAP    =    T
 PARSOL  =    F
 ISG21   =       300
 ISG52   =         1
 CSG10   ='q1'
 SPEDAT(SET,BOUNDARY,ZCONST,R,1.00000E+20)
 ************************************************************
  Group 20. Preliminary Printout
 ECHO    =    T
 ************************************************************
  Group 21. Print-out of Variables
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =        19 ;IYMON  =         2 ;IZMON  =         1
 NPRMON  =    100000
 NPRMNT  =         1
 TSTSWP  =        -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT  =    100000
 NXPRIN  =         1
 IXPRF   =        17 ;IXPRL  =        21
 NYPRIN  =         1
 IYPRF   =         1 ;IYPRL  =         5
 ISWPRF  =         1 ;ISWPRL =    100000
   No PATCHes used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 
 GVIEW(P,0.000000E+00,0.000000E+00,1.000000E+00)
 GVIEW(UP,4.999979E-03,9.999875E-01,0.000000E+00)
 
> DOM,    SIZE,        6.059999E-02, 6.060000E-02, :ZWW:
> DOM,    MONIT,       5.550000E-02, 2.100000E-03, 5.000000E-04
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    INCREMENT,   1.000000E-02, 1.000000E-02, 1.000000E-02
> GRID,   BOUNDS,       F F F F F F
> GRID,   RSET_X_1,     10, 1.000000E+00
> GRID,   RSET_X_2,     10, 1.000000E+00
> GRID,   RSET_X_3,      1, 1.000000E+00
> GRID,   RSET_Y_1,      1, 1.000000E+00
> GRID,   RSET_Y_2,     10, 1.000000E+00
> GRID,   RSET_Y_3,     10, 1.000000E+00
  > GRID,   RSET_Z_1,      1, 1.000000E+00
 
> OBJ,    NAME,        INNER
> OBJ,    POSITION,    0.000000E+00, 3.060000E-02, 0.000000E+00
> OBJ,    SIZE,        3.000000E-02, 3.000000E-02, :ZWW:
> OBJ,    GEOMETRY,    cubet
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    0,  Air at 20 deg C, 1 atm, treated as i
> OBJ,    SCAL_FIXF,   Per-unit-volume
 
> OBJ,    NAME,        OUTER1
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        6.059999E-02, 5.999999E-04, :ZWW:
> OBJ,    GEOMETRY,    cubet
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    0,  Air at 20 deg C, 1 atm, treated as i
> OBJ,    SCAL_FIXF,   Per-unit-volume
 
> OBJ,    NAME,        OUTER2
> OBJ,    POSITION,    6.000000E-02, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        5.999953E-04, 6.059999E-02, :ZWW:
> OBJ,    GEOMETRY,    cubet
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    0,  Air at 20 deg C, 1 atm, treated as i
> OBJ,    SCAL_FIXF,   Per-unit-volume
 
> OBJ,    NAME,        UP
> OBJ,    POSITION,    3.000000E-02, 6.059999E-02, 0.000000E+00
> OBJ,    SIZE,        3.000000E-02, 0.000000E+00, :ZWW:
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 
> OBJ,    NAME,        LEFT
> OBJ,    POSITION,    0.000000E+00, 5.999999E-04, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 3.000000E-02, :ZWW:
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 
> OBJ,    NAME,        IN
> OBJ,    POSITION,    0.000000E+00, 3.060000E-02, 0.000000E+00
> OBJ,    SIZE,        3.000000E-02, 2.999999E-02, :ZWW:
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 
> OBJ,    NAME,        UIN
> OBJ,    POSITION,    0.000000E+00, 3.060000E-02, 0.000000E+00
> OBJ,    SIZE,        2.700000E-02, 2.999999E-02, :ZWW:
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    GRID,        NO
> OBJ,    TYPE,        USER_DEFINED
 
> OBJ,    NAME,        VIN
> OBJ,    POSITION,    0.000000E+00, 3.060000E-02, 0.000000E+00
> OBJ,    SIZE,        3.000000E-02, 2.999999E-02, :ZWW:
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    GRID,        NO
> OBJ,    TYPE,        USER_DEFINED
 
> OBJ,    NAME,        OUT1
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        6.059999E-02, 5.999999E-04, :ZWW:
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 
> OBJ,    NAME,        OUT2
> OBJ,    POSITION,    6.000000E-02, 5.999999E-04, 0.000000E+00
> OBJ,    SIZE,        5.999953E-04, 5.999999E-02, :ZWW:
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 ************************************************************
STOP